How to fix mates in PE Illumina data after alignment and discard reads with no mates?
How to fix mates in PE Illumina data after alignment and discard reads with no mates?
• 2,030 views
•
link
updated
by
Ram
•
written
by
lakhujanivijay
•
written
by
lakhujanivijay
1 answer
I'm sure there are bigger-name pieces of software that will do it, but these scripts by the MacDonald lab are pretty quick and easy and have good instructions. You need Picard installed to use them but they're written in Perl
http://www.mcdonaldlab.biology.gatech.edu/bioinformatics.htm
edit: no affiliation, I just found them very useful when I was getting started.
• 0 views
•
link
updated
by
Ram
•
written
by
memory_donk
•
written
by
memory_donk
Log in to answer this question.
More posts like this
-
masurca invalid forward file for PE library
written by m.t.lorenc •
-
Extracting Read and mates having variant positions Using BAM file vcf file in Python using pysam/ b…
written by Hassan •
-
Masurca hybrid assembly throwing Error correction of PE reads failed error
written by Pac314 •
-
It is possible to improve my assembly?
written by emile.cavalet_giorsa •
-
Does Last aligner only outputs the mapped reads to a bam file ?
written by pinn
-
Minimum length in fastq quality trimmer with FASTX toolkit
written by Lucila •
-
PBJelly error at mapping level
written by SJ Basu •
-
The comparison between HISAT2 and Tophat2
written by Peng Huang
-
Masurca assemble questions
written by tigerhui587 •
-
SAM format paired-end alignment handling
written by darxsys
It's better to do that before alignment. Indeed, it is best to use data that is not corrupted. So the real question is, how did you manage to destroy the pairing of the reads during preprocessing? You need to redo the preprocessing with a pair-aware tool.