Illumina / GTCtoVCF

Script to convert GTC/BPM files to VCF
Apache License 2.0
41 stars 31 forks source link

Bug/baf mnv handling #57

Closed jjzieve closed 4 years ago

jjzieve commented 4 years ago

This PR is to handle a bug with MNV handling (via auxiliary loci file) and the BAlleleFrequencyFormat.

Previously, this class used the SNP string directly from a BPMRecord (e.g. [A/G]) when canonicalizing the BAF for the VCF. But the auxiliary loci logic only overrides the plus_strand_alleles which are essentially a canonical representation of the SNP based on the raw SNP and the ref_strand.

The change to the class, is to use the plus_strand_alleles, which is the approach used in the GenotypeFormat class.