.\" Man page generated from reStructuredText. . .TH "OBISELECT" "1" "Jul 27, 2019" " 1.02 13" "OBITools" .SH NAME obiselect \- description of obiselect . .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 \fI\%obiselect\fP command allows to select a subset of sequences records from a sequence file by describing sequence record groups and defining how many and which sequence records from each group must be retrieved. .sp In each group as definied by a set of \fI\-c\fP options, sequence records are ordered according to a score function. The \fIN\fP first sequences (\fIN\(gais selected using the \(ga\-n\fP option) are kept in the result subset of sequence records. .sp By default the score function is a random function and one sequence record is retrieved per group. This leads to select randomly one sequence per group. .SH OBISELECT SPECIFIC OPTIONS .INDENT 0.0 .TP .B \-c , \-\-category\-attribute= .INDENT 7.0 .INDENT 3.5 Attribute used to categorize the sequence records. Several \fB\-c\fP options can be combined. .sp \fBTIP:\fP .INDENT 0.0 .INDENT 3.5 The \fB\fP can be simply the key of an attribute, or a \fIPython\fP expression similarly to the \fB\-p\fP option of \fBobigrep\fP\&. .UNINDENT .UNINDENT .UNINDENT .UNINDENT .sp \fIExample:\fP .INDENT 7.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C > obiselect \-c sample \-c seq_length seq.fasta .ft P .fi .UNINDENT .UNINDENT .UNINDENT .UNINDENT .sp This command select randomly one sequence record per sample and sequence length from the sequence records included in the \fIseq.fasta\fP file. The selected sequence records are printed on the screen. .UNINDENT .UNINDENT .UNINDENT .INDENT 0.0 .TP .B \-n , \-\-number= .INDENT 7.0 .INDENT 3.5 Indicates how many sequence records per group have to be retrieved. If the size of the group is lesser than this \fINUMBER\fP, the whole group is retrieved. .UNINDENT .UNINDENT .sp \fIExample:\fP .INDENT 7.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C > obiselect \-n 2 \-c sample \-c seq_length seq.fasta .ft P .fi .UNINDENT .UNINDENT .UNINDENT .UNINDENT .sp This command has the same effect than the previous example except that two sequences are retrieved by class of sample/length. .UNINDENT .UNINDENT .UNINDENT .INDENT 0.0 .TP .B \-\-merge= Attribute to merge. .sp \fIExample:\fP .INDENT 7.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C > obiselect \-c seq_length \-n 2 \-m sample seq1.fasta > seq2.fasta .ft P .fi .UNINDENT .UNINDENT .sp This command keeps two sequences per sequence length, and records how many times they were observed for each sample in the new attribute \fBmerged_sample\fP\&. .UNINDENT .UNINDENT .UNINDENT .INDENT 0.0 .TP .B \-\-merge\-ids Adds a \fBmerged\fP attribute containing the list of sequence record ids merged within this group. .UNINDENT .INDENT 0.0 .TP .B \-m, \-\-min Sets the function used for scoring sequence records into a group to the minimum function. The minimum function is applied to the values used to define categories (see option \fI\-c\fP). Sequences will be ordered according to the distance of their values to the minimum value. .UNINDENT .INDENT 0.0 .TP .B \-M, \-\-max Sets the function used for scoring sequence records into a group to the maximum function. The maximum function is applied to the values used to define categories (see option \fI\-c\fP). Sequences will be ordered according to the distance of their values to the maximum value. .UNINDENT .INDENT 0.0 .TP .B \-a, \-\-mean Sets the function used for scoring sequence records into a group to the mean function. The mean function is applied to the values used to define categories (see option \fI\-c\fP). Sequences will be ordered according to the distance of their values to the mean value. .UNINDENT .INDENT 0.0 .TP .B \-\-median Sets the function used for scoring sequence records into a group to the median function. The median function is applied to the values used to define categories (see option \fI\-c\fP). Sequences will be ordered according to the distance of their values to the median value. .UNINDENT .INDENT 0.0 .TP .B \-f FUNCTION, \-\-function=FUNCTION Sets the function used for scoring sequence records into a group to a user define function. The user define function is declared using \fIPython\fP syntax. Attribute keys can be used as variables. An extra \fIsequence\fP variable representing the full sequence record is available. If option for loading a taxonomy database is provided, a \fItaxonomy\fP variable is also available. The function is estimated for each sequence record and the minimum value of this function in each group. Sequences will be ordered in each group according to the distance of their function estimation to the minimum value of their group. .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 TAXONOMY RELATED OPTIONS .INDENT 0.0 .TP .B \-d , \-\-database= ecoPCR taxonomy Database name .UNINDENT .INDENT 0.0 .TP .B \-t , \-\-taxonomy\-dump= NCBI Taxonomy dump repository name .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 OBISELECT ADDED SEQUENCE ATTRIBUTES .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .IP \(bu 2 class .IP \(bu 2 distance .IP \(bu 2 merged .IP \(bu 2 class .IP \(bu 2 merged_* .IP \(bu 2 select .UNINDENT .UNINDENT .UNINDENT .SH OBISELECT USED SEQUENCE ATTRIBUTE .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .IP \(bu 2 taxid .UNINDENT .UNINDENT .UNINDENT .SH AUTHOR The OBITools Development Team - LECA .SH COPYRIGHT 2019 - 2015, OBITool Development Team .\" Generated by docutils manpage writer. .