Provided by: mothur_1.42.1-1build1_amd64 bug

NAME

       uchime - reads a fasta file and reference file and outputs potentially chimeric sequences

DESCRIPTION

       The  chimera.uchime  command  reads  a  fasta  file  and  reference file and outputs potentially chimeric
       sequences. The original uchime program was written by Robert C. Edgar.

SYNOPSIS

       uchime --input query.fasta [--db db.fasta] [--uchimeout results.uchime]

              [--uchimealns results.alns]

OPTIONS

       --input filename

              Query sequences in FASTA format.  If the --db option  is  not  specificed,  uchime  uses  de  novo
              detection. In de novo mode, relative abundance must be given by a string /ab=xxx/ somewhere in the
              label, where xxx is a floating-point number, e.g. >F00QGH67HG/ab=1.2/.

       --db filename

              Reference database in FASTA format.  Optional, if not specified uchime uses de novo mode.

              ***WARNING***  The  database  is  searched  ONLY  on  the   plus   strand.    You   MUST   include
              reverse-complemented sequences in the database if you want both strands to be searched.

       --abskew x

              Minimum abundance skew. Default 1.9. De novo mode only.  Abundance skew is:

              min [ abund(parent1), abund(parent2) ] / abund(query).

       --uchimeout filename

              Output  in  tabbed  format  with  one record per query sequence.  First field is score (h), second
              field is query label.  For details, see manual.

       --uchimealns filename

              Multiple alignments of query sequences to parents in humanreadable format. Alignments show columns
              with differences that support or contradict a chimeric model.

       --minh h

              Minimum  score  to  report  chimera.  Default 0.3. Values from 0.1 to 5 might be reasonable. Lower
              values increase sensitivity but may report more false positives. If you  decrease  --xn,  you  may
              need to increase --minh, and vice versa.

       --mindiv div

              Minimum  divergence  ratio,  default 0.5. Div ratio is 100% - %identity between query sequence and
              the closest candidate for being a parent. If you don't care about very close  chimeras,  then  you
              could  increase  --mindiv  to, say, 1.0 or 2.0, and also decrease --min h, say to 0.1, to increase
              sensitivity.  How well this works will depend on your data. Best is to tune parameters on  a  good
              benchmark.

       --xn beta

              Weight  of  a  no  vote,  also  called the beta parameter. Default 8.0.  Decreasing this weight to
              around 3 or 4 may give better performance on denoised data.

       --dn n

              Pseudo-count prior on number of no votes. Default 1.4. Probably no  good  reason  to  change  this
              unless  you  can  retune  to a good benchmark for your data. Reasonable values are probably in the
              range from 0.2 to 2.

       --xa w

              Weight of an abstain vote. Default 1. So far, results do not seem to be  very  sensitive  to  this
              parameter,  but  if  you  have  a good training set might be worth trying. Reasonable values might
              range from 0.1 to 2.

       --chunks n

              Number of chunks to extract from the query sequence when searching for parents. Default 4.

       --[no]ovchunks

              [Do not] use overlapping chunks. Default do not.

       --minchunk n

              Minimum length of a chunk. Default 64.

       --idsmoothwindow w

              Length of id smoothing window. Default 32.

       --minsmoothid f

              Minimum factional identity over smoothed window of candidate parent.  Default 0.95.

       --maxp n

              Maximum number of candidate parents to consider. Default 2. In tests  so  far,  increasing  --maxp
              gives only a very small improvement in sensivity but tends to increase the error rate quite a bit.

       --[no]skipgaps --[no]skipgaps2

              These  options  control  how gapped columns affect counting of diffs.  If --skipgaps is specified,
              columns containing gaps do not found  as  diffs.   If  --skipgaps2  is  specified,  if  column  is
              immediately  adjacent  to  a  column  containing  a  gap, it is not counted as a diff.  Default is
              --skipgaps --skipgaps2.

       --minlen L --maxlen L

              Minimum and maximum sequence length. Defaults 10, 10000.  Applies  to  both  query  and  reference
              sequences.

       --ucl

              Use  local-X  alignments.  Default  is  global-X. On tests so far, global-X is always better; this
              option is retained because it just might work well on some future type of data.

       --queryfract f

              Minimum fraction of the query sequence that must be covered by a local-X alignment.  Default  0.5.
              Applies only when --ucl is specified.

       --quiet

              Do not display progress messages on stderr.

       --log filename

              Write  miscellaneous  information  to  the  log  file.  Mostly  of  interest  to me (the algorithm
              developer). Use --verbose to get more info.

       --self

              In reference database mode, exclude a reference sequence if it has the same label  as  the  query.
              This is useful for benchmarking by using the ref db as a query to test for false positives.

       --abskew <float>
              help

       --absort <str>
              help

       --abx <float>
              help

       --allpairs <str>
              help

       --alpha <str>
              help

       --band <uint>
              help

       --blast6out <str>
              help

       --[no]blast_termgaps
              help

       --blastout <str>
              help

       --bump <uint>
              help

       --[no]cartoon_orfs
              help

       --cc <str>
              help

       --chain_evalue <float>
              help

       --chain_targetfract <float>
              help

       --chainhits <str>
              help

       --chainout <str>
              help

       --chunks <uint>
              help

       --clstr2uc <str>
              help

       --clump <str>
              help

       --clump2fasta <str>
              help

       --clumpfasta <str>
              help

       --clumpout <str>
              help

       --cluster <str>
              help

       --compilerinfo
              Write info about compiler types and #defines to stdout.

       --computekl <str>
              help

       --db <str>
              help

       --dbstep <uint>
              help

       --[no]denovo
              help

       --derep
              help

       --diffchar <str>
              help

       --dn <float>
              help

       --doug <str>
              help

       --droppct <uint>
              help

       --evalue <float>
              help

       --evalue_g <float>
              help

       --exact
              help

       --[no]fastalign
              help

       --fastapairs <str>
              help

       --fastq2fasta <str>
              help

       --findorfs <str>
              help

       --[no]flushuc
              help

       --frame <int>
              help

       --fspenalty <float>
              help

       --gapext <str>
              help

       --gapopen <str>
              help

       --getseqs <str>
              help

       --global
              help

       --hash help

       --hashsize <uint>
              help

       --help Display command-line options.

       --hireout <str>
              help

       --hspalpha <str>
              help

       --id <float>
              help

       --idchar <str>
              help

       --iddef <uint>
              help

       --idprefix <uint>
              help

       --ids <str>
              help

       --idsmoothwindow <uint>
              help

       --idsuffix <uint>
              help

       --indexstats <str>
              help

       --input <str>
              help

       --[no]isort
              help

       --k <uint>
              help

       --ka_dbsize <float>
              help

       --ka_gapped_k <float>
              help

       --ka_gapped_lambda <float>
              help

       --ka_ungapped_k <float>
              help

       --ka_ungapped_lambda <float>
              help

       --[no]label_ab
              help

       --labels <str>
              help

       --[no]leftjust
              help

       --lext <float>
              help

       --local
              help

       --log <str>
              Log file name.

       --[no]log_hothits
              help

       --[no]log_query
              help

       --[no]logmemgrows
              help

       --logopts
              Log options.

       --[no]logwordstats
              help

       --lopen <float>
              help

       --makeindex <str>
              help

       --match <float>
              help

       --matrix <str>
              help

       --max2 <uint>
              help

       --maxaccepts <uint>
              help

       --maxclump <uint>
              help

       --maxlen <uint>
              help

       --maxovd <uint>
              help

       --maxp <uint>
              help

       --maxpoly <uint>
              help

       --maxqgap <uint>
              help

       --maxrejects <uint>
              help

       --maxspan1 <uint>
              help

       --maxspan2 <uint>
              help

       --maxtargets <uint>
              help

       --maxtgap <uint>
              help

       --mcc <str>
              help

       --mergeclumps <str>
              help

       --mergesort <str>
              help

       --minchunk <uint>
              help

       --mincodons <uint>
              help

       --mindiffs <uint>
              help

       --mindiv <float>
              help

       --minh <float>
              help

       --minhsp <uint>
              help

       --minlen <uint>
              help

       --minorfcov <uint>
              help

       --minspanratio1 <float>
              help

       --minspanratio2 <float>
              help

       --[no]minus_frames
              help

       --mismatch <float>
              help

       --mkctest <str>
              help

       --[no]nb
              help

       --optimal
              help

       --orfstyle <uint>
              help

       --otusort <str>
              help

       --output <str>
              help

       --[no]output_rejects
              help

       --probmx <str>
              help

       --query <str>
              help

       --queryfract <float>
              help

       --querylen <uint>
              help

       --quiet
              Turn off progress messages.

       --randseed <uint>
              help

       --realign
              help

       --[no]rev
              help

       --[no]rightjust
              help

       --rowlen <uint>
              help

       --secs <uint>
              help

       --seeds <str>
              help

       --seedsout <str>
              help

       --seedt1 <float>
              help

       --seedt2 <float>
              help

       --self help

       --[no]selfid
              help

       --simcl <str>
              help

       --[no]skipgaps
              help

       --[no]skipgaps2
              help

       --sort <str>
              help

       --sortuc <str>
              help

       --sparsedist <str>
              help

       --sparsedistparams <str>
              help

       --split <float>
              help

       --[no]ssort
              help

       --sspenalty <float>
              help

       --[no]stable_sort
              help

       --staralign <str>
              help

       --stepwords <uint>
              help

       --strand <str>
              help

       --targetfract <float>
              help

       --targetlen <uint>
              help

       --tmpdir <str>
              help

       --[no]trace
              help

       --tracestate <str>
              help

       --[no]trunclabels
              help

       --[no]twohit
              help

       --uc <str>
              help

       --uc2clstr <str>
              help

       --uc2fasta <str>
              help

       --uc2fastax <str>
              help

       --uchime <str>
              help

       --uchimealns <str>
              help

       --uchimeout <str>
              help

       --[no]ucl
              help

       --uhire <str>
              help

       --ungapped
              help

       --userfields <str>
              help

       --userout <str>
              help

       --usersort
              help

       --uslink <str>
              help

       --[no]usort
              help

       --utax <str>
              help

       --[no]verbose
              help

       --version
              Show version and exit.

       --w <uint>
              help

       --weak_evalue <float>
              help

       --weak_id <float>
              help

       --[no]wordcountreject
              help

       --[no]wordweight
              help

       --xa <float>
              help

       --xdrop_g <float>
              help

       --xdrop_nw <float>
              help

       --xdrop_u <float>
              help

       --xdrop_ug <float>
              help

       --xframe <str>
              help

       --xlat help

       --xn <float>
              help

AUTHOR

       Robert C. Edgar

SEE ALSO

       http://www.drive5.com/uchime