Skip to content
Snippets Groups Projects
Commit 3d9705fc authored by Brandi Cantarel's avatar Brandi Cantarel
Browse files

update filter cnvkit

parent 50154ca3
No related merge requests found
...@@ -56,7 +56,8 @@ while (my $line = <CNR>) { ...@@ -56,7 +56,8 @@ while (my $line = <CNR>) {
foreach my $gid (@ids) { foreach my $gid (@ids) {
next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/); next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/);
my ($gene,@other) = split(/:/,$gid); my ($gene,@other) = split(/:/,$gid);
$genes{$gene} = 1; my ($gname,@loc) = split(/_/,$gene);
$genes{$gname} = 1;
} }
} }
foreach $gene (keys %genes) { foreach $gene (keys %genes) {
...@@ -94,7 +95,8 @@ while (my $line = <IN>) { ...@@ -94,7 +95,8 @@ while (my $line = <IN>) {
foreach my $gid (@ids) { foreach my $gid (@ids) {
next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/); next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/);
my ($gene,@other) = split(/:/,$gid); my ($gene,@other) = split(/:/,$gid);
$genes{$gene} = 1; my ($gname,@loc) = split(/_/,$gene);
$genes{$gname} = 1;
} }
} }
my $len = sprintf("%.1f",($hash{end}-$hash{start})/1000); my $len = sprintf("%.1f",($hash{end}-$hash{start})/1000);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment