I understand that the "bamtools split" command retrieves those reads. However, I still don't know the FLAG setting when dividing into pair-end reads (LR:-> <-), mate pair (RL: <- ->) and others (RR or LL).
I would appreciate it if you could explain it in an easy-to-understand manner.
Extract reads based on the orientation of the mapped reads.
I want to extract reads from a BAM file based on the orientation of the mapped reads. I want to extract the reads of each category separately based on the orientation of the reads such as LR (Illumina pair-end reads), RL (Illumina mate-pair reads), RR or LL as depicted in IGV. Please tell me if there is any solution to that.
• 1,744 views
•
link
updated
by
colindaven
•
written
by
kaz0816 •
•
written
by
kaz0816 •
1 answer
You can use bamtools for this. Look at the tutorial PDF here:
https://github.com/pezmaster31/bamtools/wiki/Using-the-toolkit
You have to write a config file for it.
• 0 views
•
link
written
by
colindaven
Log in to answer this question.
More posts like this
-
In what case there are RR reads but no LL reads when detect inversions using IGV?
written by Jingfang •
-
Why this alignment is shown as "Tandem Duplication"/"Translocation on the Same Chromosome" when co…
written by agreatgoodfellow •
-
How to extend assembled PACBIO contigs for scaffolding
written by K •
-
ADR, ADF and dp4 counts do not count strands correct
written by g.t.los •
-
Extract reads based on pair orientation
written by Emilie •
-
extracting orphan reads from bam file
written by Assa Yeroslaviz
-
How to interpret reads mapped to genome from samtools flagstat vs TopHat align_summary.txt ?
written by Vasu
-
Interpreting Flags in SAM file
written by jayg •
-
How to align Trimmomatic unpaired reads with BWA?
written by mcff23 •
-
Mate pair Illumina reads - recommended mapper of choice, adapter removal and PE contamination
written by Biomonika (Noolean)
