Question: Call alternative base for each position in a bam file
1
gravatar for sacha
3 months ago by
sacha2.0k
France
sacha2.0k wrote:

Hi,

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:

enter image description here

#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
bam calling • 144 views
ADD COMMENTlink modified 3 months ago by dariober11k • written 3 months ago by sacha2.0k
2

Take a look at igvtools count included in IGV. Link to manual.

ADD REPLYlink modified 3 months ago • written 3 months ago by genomax91k
0
gravatar for dariober
3 months ago by
dariober11k
WCIP | Glasgow | UK
dariober11k wrote:

See also this thread Calculate The Frequency Of Nucleotides At Each Position In An Mpileup File I suggested there pysamstats and I'm still quite happy with it.

ADD COMMENTlink written 3 months ago by dariober11k

I tried pysamstats by the past . I was thinking it is only for read depths .

ADD REPLYlink modified 3 months ago • written 3 months ago by sacha2.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1664 users visited in the last hour