Hi,
I've a question about intersectBed. Is it possible to extract only alignment like this :
chromosome ===============================================================
BED/BAM A ============== =================
BED FILE B ============
RESULT ==============
But no alignment like this (even if the read overlapp 100% of the feature, I don't want to extract these reads)
chromosome ===============================================================
BED/BAM A ========================= =================
BED FILE B =============
RESULT
So, only extracting reads that have 90-95% of its sequence overlapping 90-95% of the feature.
Is it clear ?
Thanks,
N.