Hello, I used samtools to get the data for chromosome 22. The command is:
samtools view input.bam 'chr22' >raw_chr22.txt
Then I extracted the 3,4 and 10th columns from it. The partial data for demo purpose like
chr22 14430092 NNNNCNAGCNGAGCGNNTCTGGGNACCTCGAAGGCAGACATG chr22 14430092 NNNNNNNNNNNTNNNNNNNNANNGNNTNNNNNNNNNNNNNNN chr22 14430092 CCTCGCGGGACTGGTATGGGGACGGTCATGCAATCTGGACAA
The data has three columns, chromosome, position and sequence. My question is that they all have the same position but different sequence. So what is the real sequence in this specific position.