Hi, I have a question dealing with PE RNAseq data (100x2) from Illumina.
It is my understanding that Bowtie2 checks a read for possible alignments across the reference genome, and reports the best alignment, discarding the rest. However, I am interested in finding just the reads that align to exactly 2 places, and determining where both locations are on the reference. I want to be able to extract this information from the BAM file if possible. I first tried setting the -k parameter to 5, and then looking for reads that could only find two alignments, but am not sure where in the BAM to find this information.
Does anyone have experience looking for such reads in a SAM/BAM file?
Thanks!
It is my understanding that Bowtie2 checks a read for possible alignments across the reference genome, and reports the best alignment, discarding the rest. However, I am interested in finding just the reads that align to exactly 2 places, and determining where both locations are on the reference. I want to be able to extract this information from the BAM file if possible. I first tried setting the -k parameter to 5, and then looking for reads that could only find two alignments, but am not sure where in the BAM to find this information.
Does anyone have experience looking for such reads in a SAM/BAM file?
Thanks!
Comment