Entering edit mode
9.7 years ago
mchone69
•
0
My mpileup output looks like this:
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT HS
Chr1 1954 . G A 155 . DP=21;VDB=0.0401;AF1=1;AC1=2;DP4=0,0,13,7;MQ=42;FQ=-87 GT:PL:DP:SP:GQ 1/1:188,60,0:20:0:99
Chr1 23056 . C A 169 . DP=22;VDB=0.0388;AF1=1;AC1=2;DP4=0,0,12,10;MQ=42;FQ=-93 GT:PL:DP:SP:GQ 1/1:202,66,0:22:0:99
Chr1 31071 . A G 182 . DP=31;VDB=0.0333;AF1=1;AC1=2;DP4=0,0,15,15;MQ=42;FQ=-117 GT:PL:DP:SP:GQ 1/1:215,90,0:30:0:99
Chr1 31478 . C T 169 . DP=26;VDB=0.0394;AF1=1;AC1=2;DP4=0,0,15,11;MQ=42;FQ=-105 GT:PL:DP:SP:GQ 1/1:202,78,0:26:0:99
Chr1 33667 . A G 172 . DP=24;VDB=0.0399;AF1=1;AC1=2;DP4=0,0,10,14;MQ=42;FQ=-99 GT:PL:DP:SP:GQ 1/1:205,72,0:24:0:99
Could someone tell me what's the relationship between QUAL and FQ, or how can I get the FQ value if I only have QUAL value?
THANKS A LOT.