NAME

illuminapairedend - description of illuminapairedend

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.