fixing more problems with allele lengths < 1bp in the VCF

mgymrek · Jul 29, 2014 · 46dd0c8 · 46dd0c8
1 parent 5fd8768
commit 46dd0c8
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/Genotyper.cpp b/src/Genotyper.cpp
@@ -119,8 +119,9 @@ void Genotyper::GetAlleles(const list<AlignedRead>& aligned_reads,
   for (list<AlignedRead>::const_iterator it = aligned_reads.begin();
        it != aligned_reads.end(); it++) {
     if (it->mate) continue;
+    int refcopy = it->msEnd - it->msStart + 1;
     if (it->diffFromRef != 0 && std::find(alleles->begin(), alleles->end(), it->diffFromRef) == alleles->end() &&
-	it->diffFromRef+it->refCopyNum*it->period > 0) {
+	it->diffFromRef+refcopy >= 1) {
       alleles->push_back(it->diffFromRef);
     }
   }