Dear all,
I'm trying to recover reads sequences from specific region in bam file.
Reads were align using bowtie2.
I have reads align to one specific sequence and I would like to get sequences of reads mapped to a specific region of 3 nucleotides, and count the occurence of each tri-nucleotide at this position.
How would I achieve this ?
Thanks for your help !
I'm trying to recover reads sequences from specific region in bam file.
Reads were align using bowtie2.
I have reads align to one specific sequence and I would like to get sequences of reads mapped to a specific region of 3 nucleotides, and count the occurence of each tri-nucleotide at this position.
How would I achieve this ?
Thanks for your help !
Comment