illuminapairedend: aligns paired-end Illumina reads

Important

illuminapairedend replaces solexapairend.

illuminapairedend aims at aligning the two reads of a pair-end library sequenced using an Illumina platform.

  • If the two reads overlap, it returns the consensus sequence together with its quality
  • Otherwise, it concatenates sequence merging the forward read and the reversed-complemented reverse read.

The program uses as input one or two fastq sequences reads files.

  • If two files are used one of them must be specified using the -r option. Sequence records corresponding to the same read pair must be in the same order in the two files.
  • If just one file is provided, sequence records are supposed to be all of the same length. The first half of the sequence is used as forward read, the second half is used as the reverse read.

illuminapairedend align the forward sequence record with the reverse complement of the reverse sequence record. The alignment algorithm takes into account the base qualities.

Example:

> illuminapairedend -r seq3P.fastq seq5P.fastq > seq.fastq

The seq5P.fastq sequence file contains the forward sequence records. The seq3P.fastq sequence file contains the reverse sequence records. Pairs of reads are aligned together and the consensus sequence is stored in the `` seq.fastq`` file.

illuminapairedend specific options

-r <FILENAME>, --reverse-reads=<FILENAME>

Filename points to the file containing the reverse reads.

--index-file=<FILENAME>
Filename points to the file containing the illumina index reads
--score-min=<FLOAT>

minimum score for keeping alignment. If the alignment score is below this threshold both the sequences are just concatenated. The mode attribute is set to the value joined.

Options to specify input format

Options to specify output format

Standard output format

--fasta-output

Output sequences in OBITools fasta format

--fastq-output

Output sequences in Sanger fastq format

Generating an ecoPCR database

--ecopcrdb-output=<PREFIX_FILENAME>

Creates an ecoPCR database from sequence records results

Miscellaneous option

--uppercase

Print sequences in upper case (default is lower case)

Common options

-h, --help

Shows this help message and exits.

--DEBUG

Sets logging in debug mode.