Question: 16S Paired-End Illumina MiSeq Data in Qiime
0
gravatar for Ruth M.
3.0 years ago by
Ruth M.10
USA
Ruth M.10 wrote:

I am learning to use Qiime to analyze 16S V3/V4 sequencing data from the Illumina Miseq (PE 2x300).

I'm wondering if there is any regular convention in deciding whether or not to pair ends before processing data. Because the second read from the Illumina 600v3 2x300 gives you poor quality after 200-250 bases, it seems that there would be many reads that don't get paired if using any kind of pairing program, right?

Is it possible to simply merge the two fastq files for read 1 and 2, or would that not work because of directionality?

I was going to use multiple_split_libraries_fastq.py, but if I have read 1 and 2 separate, that won't work either, correct?

Thanks in advance for any information!

paired-end miseq 16s illumina • 3.1k views
ADD COMMENTlink modified 3.0 years ago by Picasa440 • written 3.0 years ago by Ruth M.10
1
gravatar for Brian Bushnell
3.0 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

I suggest merging the reads with BBMerge (part of the BBMap package) prior to analysis. If the quality is so bad that you get a low merge rate, you can add the flag "loose" or "vloose" to increase it, at the expense of a slightly higher false-positive rate; or add "qtrim=r trimq=10" to do quality-trimming first. But overall, merged reads are easier to analyze and will have a lower error rate, so you should get more accurate clustering.

ADD COMMENTlink written 3.0 years ago by Brian Bushnell16k
1
gravatar for Charles Warden
3.0 years ago by
Charles Warden7.0k
Duarte, CA
Charles Warden7.0k wrote:

1) If the read merging options in QIIME's join_paired_ends.py don't work well, I've found that PEAR usually works pretty good, typically merging a greater percentage of reads that still pass the QIIME quality score filter at a similar rate.

2) You can try using RTAX (inside or outside QIIME), which is supposed to work without merging the reads.

3) If you are still having problems, you can probably just use the forward read and get similar classification results.

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by Charles Warden7.0k
1
gravatar for Picasa
3.0 years ago by
Picasa440
Picasa440 wrote:

What I usually do is using trimmomatic to perform QC processing then merging the 2 PE with FLASh.

https://ccb.jhu.edu/software/FLASH/

http://www.usadellab.org/cms/?page=trimmomatic

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by Picasa440
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: 829 users visited in the last hour