Entering edit mode
9.0 years ago
andynkili
▴
10
Is it possible to get every reads that match several loci on the reference (using the read mappers bowtie2)?
I've looked for any options allowing to do so, but I'm only aware of the --al
option, saving all the reads that matched the reference sequence. The ultimate goal is to keep only the reads that have a unique position (match) in the reference sequence, i.e the reads that can be trusted regarding their position.
THANK YOU !