]> git.donarmstrong.com Git - samtools.git/commitdiff
Speedup with -i, no need to query the reference all the time
authorPetr Danecek <pd3@sanger.ac.uk>
Wed, 14 Apr 2010 15:44:50 +0000 (15:44 +0000)
committerPetr Danecek <pd3@sanger.ac.uk>
Wed, 14 Apr 2010 15:44:50 +0000 (15:44 +0000)
misc/sam2vcf.pl

index 38ab64dc827d80757b193317c5d0302f7fb8efdc..874ad5c19462b0e874b19cbdf6ab575271e53672 100755 (executable)
@@ -131,7 +131,13 @@ sub do_pileup_to_vcf
         if ( $ref eq '*' )
         {
             # An indel is involved.
-            if ( $ignore_indels ) { next; }
+            if ( $ignore_indels )
+            { 
+                $prev_ref = $ref;
+                $prev_pos = $pos;
+                $prev_chr = $chr;
+                next; 
+            }
 
             if (!defined $prev_chr || $chr ne $prev_chr || $pos ne $prev_pos) 
             {
@@ -171,7 +177,13 @@ sub do_pileup_to_vcf
         }
         else
         {
-            if ( $ignore_snps ) { next; }
+            if ( $ignore_snps ) 
+            { 
+                $prev_ref = $ref;
+                $prev_pos = $pos;
+                $prev_chr = $chr;
+                next; 
+            }
 
             # SNP
             ($alt,$gt) = iupac_to_gtype($ref,$cons);