fix VCF indel position issue

OpenGene · Oct 22, 2017 · 7e2e3eb · 7e2e3eb
1 parent ade5922
commit 7e2e3eb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/mutation.cpp b/src/mutation.cpp
@@ -252,7 +252,7 @@ vector<Mutation> Mutation::parseVcf(string vcfFile, string refFile) {
             ss<<cds<<"_";
         ss<<chrom<<"_"<<v.pos<<"_"<<v.ref<<">"<<v.alt;
         string name = ss.str();
-        string left = ref[chrom].substr(v.pos-25-v.ref.length(), 25);
+        string left = ref[chrom].substr(v.pos-25-1, 25);
         str2upper(left);
         string center = v.alt;
         string right = ref[chrom].substr(v.pos+v.ref.length()-1, 25);