How can I extract base and read name based on position from bam file?
0
0
Entering edit mode
8 months ago
gh • 0

Hello.

Sadly, I can't not use English well....

I want to make primer for many cultivar of my subject. In my primer position, there are two bialleic SNP.

ex) .....A/C..........C/G........

I think reason of this problem is that my subject is polyploid. I already checked that there are two bases on that posistion using bam file, IGV and bam-readcount.

What I want to know is base group(?).

ex)

A..........G (A,G)

C..........C (C,C)

or

A..........C (A,C)

C..........G (C,G)

If base group of reference is (A,G) and base groups that checked are (A,G) , (C,C), I can use primer that I maked using reference sequence.

I think it is similar work of haplotype phasing. Of course I'm not sure. I can't check information like this using bam-readcount. I can check it using IGV, but it takes too long. because I want check many position.

Is there any work that I check base on my target position and read name at the same time? Or is there a better way?

Thank you.

bam sam genotype • 682 views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Thank you! Unfortunately, I searched information but I don't know how can I see specific region that I want using sam2tsv. T-T

ADD REPLY
0
Entering edit mode

I'm trying to do this using this command line.

Of course It was unsuccessful yet...

samtools view /my_bam_file.bam scaffold158:794058-794078 | awk '{print $1"\t"$3"\t"$4"\t"$6"\t"$10}'

result : read_name contig_name read_start_position CIAGAR sequence

ADD REPLY
0
Entering edit mode

what do you mean by "base group"?

ADD REPLY
0
Entering edit mode

Ummm.. pair of base on the position that I want to see in per read. I try to explain it.

enter image description here

ADD REPLY

Login before adding your answer.

Traffic: 1799 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