I've got a bam file of paired-end illumina sequences mapped to a reference sequence and I was wondering if when I call the samtools view command if the results represent only reads that completely span over the region of interest or if it also included reads that intersect it by one or two base pairs.
If not are there any existing tools out there to only read thats completely span over the region of interest?
with this visual example of what I mean I would want to only include sequences like case_2 and exclude case_1:
[===================Region-of-Interest=====================]
[---------------------------------------------------------------------------------------------------------] case_1
[-------------------------------------------------------------------------------------------------------------------] case_2
Thanks!
Hey h.mon,
bedtools is definitely the way to go, I ended up using
Thanks for the help!
Checkmark this as the accepted answer if they helped you :)