Combine paired reads to contigs


Combines paired reads to contigs for each sample and then collects the sequences from all the samples in one fasta file.




This tool takes a Tar package of fastq files (you can create one using the tool Utilities / Make a Tar package). It extracts the sequence and quality data from the fastq files, creates a reverse complement of the reverse read, and joins the reads into contigs by performing a Needleman alignment. When it identifies positions where the two reads disagree, it proceeds as follows: Note that Mothur Make.contigs tool is quite "greedy", so it will assemble contigs even if the reads are not high quality. Therefore it is recommended to use this tool only for good quality reads which totally overlap. Dr Patrick Schloss, the developer of Mothur, says that even with the advance in Illumina MiSeq technology to paired 300 nt reads, they are sticking with the paired 250 nt reads to sequence the 16S rRNA V4 region.

The tool figures out by itself which fastq files belong to each sample. You can check this in the output file samples.fastqs.txt. If the file assignment to samples is not correct, you can make a samples-files.txt yourself and give it as input. This file has 3 columns: sample name, forward fastq file name (R1), and reverse fastq file name (R2). So it should look something like this:
sample1 sample1_L001_R1_001.fastq sample1_L001_R2_001.fastq
sample2 sample2_L001_R1_001.fastq sample2_L001_R2_001.fastq
sample3 sample3_L001_R1_001.fastq sample3_L001_R2_001.fastq

This tool is based on the make.contigs and make.file commands of the Mothur package.


The analysis output consists of the following:


Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol, 2009. 75(23):7537-41.