I have a bam file that aligns a number of reads against a single reference. And I want to extract the fastq of each read in the exact location that those reads overlap a specific position in the reference.
The reference is ~3kb in length, and I want to extract the sequence from each read where it overlaps the positions 946-1047 in the reference. What is the best way to do that?
Thanks,
John
The reference is ~3kb in length, and I want to extract the sequence from each read where it overlaps the positions 946-1047 in the reference. What is the best way to do that?
Thanks,
John
Comment