Question: Extract reads based on pair orientation
0
gravatar for Emilie
10 months ago by
Emilie0
Emilie0 wrote:

Hi!

I have paired-end mapped reads (using bbmap) that I would like to separate based on the orientation of the read pairs. I see pairs as being able to have the following orientations:

LR: --> <--

LL or RR: --> --> or <-- <--

RL: <-- -->

I've managed to separate the reads in two, namely LR and one with the others ("discordant") using samtools:

samtools view -b infile.bam -F 1294 > outfile.bam

Which I picked up here , but am not too sure about... In IGV I can colour these discordant reads according to 'pair orientation'. I can see that they group together according to read orientation in clearly defined genomic positions.

But I have no idea how to split them based on orientation from here...

I tried several options, including a tool in the Jvarkit (unable to complete make of BioAlcidae so could not run) and other custom scripts that did not work due being dependent on older versions of samtools/htslib...

I hope someone can help me from here - thanks! ^@^

ADD COMMENTlink written 10 months ago by Emilie0

(unable to complete make of BioAlcidae so could not run

I need more information, or submit an issue: https://github.com/lindenb/jvarkit/issues

ADD REPLYlink written 10 months ago by Pierre Lindenbaum124k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1030 users visited in the last hour