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
0
Entering edit mode
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

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}'


0
Entering edit mode

what do you mean by "base group"?

0
Entering edit mode

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