Hi there, I'm attempting to use the CalculateMixingFractions tool to estimate the mixture of genomes from a VCF in a pooled BAM file and getting all NaN as output.
The output looks like this:
SAMPLE MIXING_FRACTION GENO1 NaN GENO2 NaN GENO3 NaN etc...
And the headers of my VCF file is the following:
##fileformat=VCFv4.1 ##filedate=2017.8.8 ##source=Minimac3 ##FORMAT= ##FORMAT= ##FORMAT= ##INFO= ##INFO= ##INFO= #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT GENO1 GENO2 GENO3 chr10 67193 rs1111111:67193:C:T C T . . . GT:AS 0|0:0,0 0|0:0,0 0|0:4,0
And the header for the SAM file is:
@HD VN:1.4 SO:coordinate SQ SN:chr1 LN:249250621 SQ SN:chr2 LN:243199373 ...etc SL-HEL:C66HBACXX150512:C66HBACXX:6:1302:18557:85138 147 chr1 10000 0 25M = 10003 -22 ATAACCCTAACCCTAACCCTAACCC ##BA39A<@CBBCCBBBAAAA<==> MD:Z:25 PG:Z:MarkDuplicates NM:i:0 AS:i:25 XS:i:25 etc...