diff --git a/alignment/virusalign.sh b/alignment/virusalign.sh index 1158137cc20b6f28bf3f1a0f2bc4fb66c2402845..6343ef918fe68008cbfea7545041bd9c1632f474 100644 --- a/alignment/virusalign.sh +++ b/alignment/virusalign.sh @@ -49,5 +49,5 @@ samtools index ${pairid}.viral.bam samtools idxstats ${pairid}.viral.bam >${pairid}.viral.idxstats.txt if [[ $filter == 1 ]] then - perl $baseDir/filter_viral_idxstats.txt -p ${pairid} ${pairid}.viral.idxstats.txt + perl $baseDir/filter_viral_idxstats.pl -p ${pairid} ${pairid}.viral.idxstats.txt fi diff --git a/variants/filter_cnvkit.pl b/variants/filter_cnvkit.pl index 86545f5e029436d89e241c5b124a177aa950e35f..772feaf5159faa1cf9b276140ad0faa9e6268c8e 100755 --- a/variants/filter_cnvkit.pl +++ b/variants/filter_cnvkit.pl @@ -54,7 +54,7 @@ while (my $line = <CNR>) { } else { my @ids = split(/,/,$geneids); foreach my $gid (@ids) { - next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/); + next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion|MSI/); my ($gene,@other) = split(/:/,$gid); my ($gname,@loc) = split(/_/,$gene); $genes{$gname} = 1; @@ -93,7 +93,7 @@ while (my $line = <IN>) { } else { my @ids = split(/,/,$geneids); foreach my $gid (@ids) { - next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion/); + next if ($gid =~ /^rs\d+$|^SNP:rs\d+$|^-$|Fusion|MSI/); my ($gene,@other) = split(/:/,$gid); my ($gname,@loc) = split(/_/,$gene); $genes{$gname} = 1;