SNPs coordinates for each read
2
0
Entering edit mode
2 days ago
mante93 • 0

Hi all. I need to find the coordinates for every point mutation for EACH READ in a bam file. The coordinates should be expressed as both (1) the position of the mutation in relation to the read length and (2) the position of the mutation in relation to the chromosome coordinates.

Also, for each point mutation I need to output the reference base pair and the mutated base pair.

Samtools mpileup outputs interesting info, but I would need them for EACH READ rather than piling up reads at a single genomic position.

Would you suggest to write my own script or are there already functions I can use?

Thank you.

SNPs mutation Coordinate BAM read • 117 views
ADD COMMENT
1
Entering edit mode
ADD COMMENT
0
Entering edit mode

Thank you for for quick reply! Very nice software you made. Could you please explain me why you used the CIGAR the rater than the MD tag? Thanks!

ADD REPLY
0
Entering edit mode
2 days ago
heskett ▴ 60

I don't know why you would want to do this, but it should be easy once you have the BAM and the list of variant positions

1) call variants to get a list of a variant positions 2) Filter the bam to keep all reads that overlap any variant 3) write a little script to parse each read to find the mutation position relative to that read, given that you already know the genomic position and the alternate allele

ADD COMMENT

Login before adding your answer.

Traffic: 1611 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6