Hi
I would like to filter a bam file, keeping only reads overlapping with genomic intervals from a bed file. I used samtools for this:
samtools view -b -h -L bedfile.bed bamfile.bam However the -L option does not seem to take into account the strand information.
Do you know if there is another option or way to do it that would keep strand information?