.\" Man page generated from reStructuredText. . .TH "NGSFILTER" "1" "Jul 27, 2019" " 1.02 13" "OBITools" .SH NAME ngsfilter \- description of ngsfilter . .nr rst2man-indent-level 0 . .de1 rstReportMargin \\$1 \\n[an-margin] level \\n[rst2man-indent-level] level margin: \\n[rst2man-indent\\n[rst2man-indent-level]] - \\n[rst2man-indent0] \\n[rst2man-indent1] \\n[rst2man-indent2] .. .de1 INDENT .\" .rstReportMargin pre: . RS \\$1 . nr rst2man-indent\\n[rst2man-indent-level] \\n[an-margin] . nr rst2man-indent-level +1 .\" .rstReportMargin post: .. .de UNINDENT . RE .\" indent \\n[an-margin] .\" old: \\n[rst2man-indent\\n[rst2man-indent-level]] .nr rst2man-indent-level -1 .\" new: \\n[rst2man-indent\\n[rst2man-indent-level]] .in \\n[rst2man-indent\\n[rst2man-indent-level]]u .. .sp To distinguish between sequences from different PCR products pooled in the same sequencing library, pairs of small DNA sequences (call tags, see the \fBoligoTag\fP command and its associated paper for more informations on the design of such tags) can be concatenated to the PCR primers. .sp \fI\%ngsfilter\fP takes as input sequence record files and a file describing the DNA tags and primers sequences used for each PCR sample. \fI\%ngsfilter\fP allows to demultiplex sequence records file by identifying these DNA tags and the primers. .sp \fI\%ngsfilter\fP requires a sample description file containing the description of the primers and tags associated to each sample (specified by option \fB\-t\fP). The sample description file is a text file where each line describes one sample. Columns are separated by space or tab characters. Lines beginning with the ‘#’ character will be considered as commentary lines and will simply be ignored by \fI\%ngsfilter\fP\&. .sp Here is an example of a sample description file: .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C #exp sample tags forward_primer reverse_primer extra_information gh 01_11a cacgcagtc:cacgcatcg GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=1; extraction=1; gh 01_12a cacgcatcg:cacgcagtc GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=1; extraction=2; gh 01_21a cacgcgcat:cacgctact GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=2; extraction=1; gh 01_22a cacgctact:cacgcgcat GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=2; extraction=2; gh 02_11a cacgctgag:cacgtacga GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=1; extraction=1; gh 02_12a cacgtacga:cacgctgag GGGCAATCCTGAGCCAA CCATTGAGTCTCTGCACCTATC F @ community=Festuca; bucket=1; extraction=2; .ft P .fi .UNINDENT .UNINDENT .sp The results consist of sequence records, printed on the standard output, with their sequence trimmed of the primers and tags and annotated with the corresponding experiment and sample (and possibly some extra informations). Sequences for which the tags and primers have not been well identified, and which are thus unassigned to any sample, are stored in a file if option \fB\-u\fP is specified and tagged as erroneous sequences (\fBerror\fP attribute) by \fI\%ngsfilter\fP\&. .SH NGSFILTER SPECIFIC OPTIONS .INDENT 0.0 .TP .B \-t, \-\-tag\-list Used to specify the file containing the samples description (with tags, primers, sample names,…) .UNINDENT .INDENT 0.0 .TP .B \-u, \-\-unidentified Filename used to store the sequences unassigned to any sample .UNINDENT .INDENT 0.0 .TP .B \-e, \-\-error Used to specify the number of errors allowed for matching primers [default = 2] .UNINDENT .SH OPTIONS TO SPECIFY INPUT FORMAT .SS Restrict the analysis to a sub\-part of the input file .INDENT 0.0 .TP .B \-\-skip The N first sequence records of the file are discarded from the analysis and not reported to the output file .UNINDENT .INDENT 0.0 .TP .B \-\-only Only the N next sequence records of the file are analyzed. The following sequences in the file are neither analyzed, neither reported to the output file. This option can be used conjointly with the \fI–skip\fP option. .UNINDENT .SS Sequence annotated format .INDENT 0.0 .TP .B \-\-genbank Input file is in genbank format. .UNINDENT .INDENT 0.0 .TP .B \-\-embl Input file is in embl format. .UNINDENT .SS fasta related format .INDENT 0.0 .TP .B \-\-fasta Input file is in fasta format (including OBITools fasta extensions). .UNINDENT .SS fastq related format .INDENT 0.0 .TP .B \-\-sanger Input file is in Sanger fastq format (standard fastq used by HiSeq/MiSeq sequencers). .UNINDENT .INDENT 0.0 .TP .B \-\-solexa Input file is in fastq format produced by Solexa (Ga IIx) sequencers. .UNINDENT .SS ecoPCR related format .INDENT 0.0 .TP .B \-\-ecopcr Input file is in ecoPCR format. .UNINDENT .INDENT 0.0 .TP .B \-\-ecopcrdb Input is an ecoPCR database. .UNINDENT .SS Specifying the sequence type .INDENT 0.0 .TP .B \-\-nuc Input file contains nucleic sequences. .UNINDENT .INDENT 0.0 .TP .B \-\-prot Input file contains protein sequences. .UNINDENT .SH OPTIONS TO SPECIFY OUTPUT FORMAT .SS Standard output format .INDENT 0.0 .TP .B \-\-fasta\-output Output sequences in \fBOBITools\fP fasta format .UNINDENT .INDENT 0.0 .TP .B \-\-fastq\-output Output sequences in Sanger fastq format .UNINDENT .SS Generating an ecoPCR database .INDENT 0.0 .TP .B \-\-ecopcrdb\-output= Creates an ecoPCR database from sequence records results .UNINDENT .SS Miscellaneous option .INDENT 0.0 .TP .B \-\-uppercase Print sequences in upper case (default is lower case) .UNINDENT .SH COMMON OPTIONS .INDENT 0.0 .TP .B \-h, \-\-help Shows this help message and exits. .UNINDENT .INDENT 0.0 .TP .B \-\-DEBUG Sets logging in debug mode. .UNINDENT .SH NGSFILTER ADDED SEQUENCE ATTRIBUTES .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 2.0 .IP \(bu 2 avg_quality .IP \(bu 2 complemented .IP \(bu 2 cut .IP \(bu 2 direction .IP \(bu 2 error .IP \(bu 2 experiment .IP \(bu 2 forward_match .IP \(bu 2 forward_primer .UNINDENT .INDENT 2.0 .IP \(bu 2 forward_score .IP \(bu 2 forward_tag .IP \(bu 2 head_quality .IP \(bu 2 mid_quality .IP \(bu 2 partial .IP \(bu 2 reverse_match .IP \(bu 2 reverse_primer .UNINDENT .INDENT 2.0 .IP \(bu 2 reverse_score .IP \(bu 2 reverse_tag .IP \(bu 2 sample .IP \(bu 2 seq_length .IP \(bu 2 seq_length_ori .IP \(bu 2 status .IP \(bu 2 tail_quality .UNINDENT .UNINDENT .UNINDENT .UNINDENT .SH AUTHOR The OBITools Development Team - LECA .SH COPYRIGHT 2019 - 2015, OBITool Development Team .\" Generated by docutils manpage writer. .