Entering edit mode
2.9 years ago
sacha ★ 2.4k
From a bam file and a target bed file, I would like to extract each base count. I don't want to skip low frequency variant and sequencing errors. It seems I cannot do this using GATK or freebayes.
For instance, I would like to extract same data from IGV , even if I have a small allele frequency:
#CHROM POS DEPTH A C G T chr4 795890 40 2 0 38 0 .....
I tried with bcftools, but I didn't find this information:
chr4 795890 . G A,<*> 0 . DP=40;I16=18,15,1,0,1212,48324,66,4356,1980,118800,60,3600,574,12390,25,625;QS=0.952494,0.0475059,0;SGB=-0.379885;RPB=1;MQB=1;MQSB=1;BQB=1;MQ0F=0 PL 0,42,255,99,255,255
Take a look at
igvtools countincluded in IGV. Link to manual.