How can I extract base and read name based on position from bam file?

Question

[gh]

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.

[GenoMax]

See if this helps: How can I extract base on specific loci from all reads overlapping that?

[gh]

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

[gh]

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

[sina.majidian]

what do you mean by "base group"?

[gh]

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