This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get All Alleles Of 1 Column In Bam File

Hi all, I have bam files from which I would like to extract 1 specific column (ie chr:position). I won't to have all alleles, even if they come from low quality reads. And even if the position is not a SNP (ie all alleles are the same or almost).

I tried to use "samtools view" giving the position but then the output is all reads that overlap the position so it's not exactly what I am looking for.

Does anyone has any idea how I could solve my problem ?

Thanks a lot. M.D

bam position extraction read
• 6,401 views • link
updated by guns.guru • guns.guru • written by Manon Delahaye • Manon Delahaye

4 answers

I managed to do what I was hoping for.

samtools mpileup -Bf myReference.fasta -r chr9:myPos-myPos [?]

-B allows for all location (not only SNPs)

Then with some combination of grep and others I could count my alleles.

Thank you for your help ! M.D

• 0 views • link
written by Manon Delahaye • Manon Delahaye

I think you are looking for:

samtools mpileup ....

There are a number of flags that you may want to tweak to get exactly the output you want.

• 1 views • link
written by Sean Davis Sean Davis

Hi, Thanks I will have a look. However, I thought that this command is for SNP calling. And I want the alleles at the position even if there is no variant at that position.

Do you think that is included in one of the option ?

I will have a better and closer look tomorrow.

Thanks. MD

• 0 views • link
written by Manon Delahaye • Manon Delahaye

samtools pileup -cf reference.fasta Input.bam

• 0 views • link
written by Rm Rm

Just a note that pileup is deprecated (or not even available) in newer samtools versions. Mpileup is the replacement.

• 0 views • link
written by Sean Davis Sean Davis

Hi M.D

It's been a long time, since you posted this but I have the exact same problem i.e., to extract the allele at particular positions, whether they are SNPs or not.

How did you achieve this?

Thanks

• 0 views • link
written by guns.guru • guns.guru

Please start a new thread. Also, the answer is samtools mpileup, as the others mentioned above.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this