Hi folks,
I am new to NGS, but I am about to work on amplicon sequencing with Miseq (PE). I am now studying an article from Balint et al. (2014) An Illumina metabarcoding pepline for fungi, which mentions 'It is important to preserve the order of the reads in both forward and reverse read files: the paired-end read assembler needs corresponding read orders in both files'.
This might be an innocent question, but I would like to ask how I obtain the files with the right order of the reads preserved in both forward and reverse read files? Do I get the files from a sequencing facility with specific requests or do I need a program or a script to work on afterwards?
Thank you in advance!
Johan
•
written
by
User 4014 •
2 answers
Your FASTQ data will come off the MiSeq like this:
[Sample Name]_S[Sample number]_L001_R[read number]_001.fastq.gz
# ^^^^^^^^^^^^^^
Read 1 data will have R1 for the highlighted portion of the file name.
Read 2 data will have R2 for the highlighted portion.
What the paper is saying is not to combine those files before feeding them to the assembler. There are some (mostly old) programs (trimmers, groomers, etc) that combine paired-end data into a single file (this is sometimes termed an interleaved dataset) in the course of processing.
If you have paired reads, they should come from the facility in two files, both in the same order (or in one file with read 1 and read 2 alternating). To keep them in the same order during filtering or trimming, you just need to use a pair-aware program such as BBDuk that processes both reads at the same time.
•
written
by
Brian Bushnell
Log in to answer this question.
-
Understanding vg filter nodes output in terms of forward and reverse compliment strand origin
written by cassiwatt •
-
vg convert does not preserve the order of the reads
written by Yihang •
-
When to merge sequencing data from multiple lanes (FastQToSam, SamToFastq, BWA, MergeBamFiles, or a…
written by Jordi •
-
Dealing with non overlapping reads
written by Adrien.andre •
-
Parsing out very similar sequences from paired amplicon sequencing
written by geo •
-
Forward and Reverse Reads Out of Order
written by dw845316 •
-
Low quality of reads in reverse files toward the last bases
written by deepti1rao •
-
Double peaks on FASTQC
written by User 4014 •
-
Paired-end reads - what is what?
written by darxsys
-
Sequencing And Strands
written by christiantheilhave •
Thank you very much for super quick replies! :)