File: denoise_wrapper.rst

package info (click to toggle)
qiime 1.4.0-2
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 29,704 kB
  • sloc: python: 77,837; haskell: 379; sh: 113; makefile: 103
file content (64 lines) | stat: -rw-r--r-- 1,966 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
.. _denoise_wrapper:

.. index:: denoise_wrapper.py

*denoise_wrapper.py* -- Denoise a flowgram file
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

**Description:**

This script will denoise a flowgram file in .sff.txt format, which is the output of sffinfo.


**Usage:** :file:`denoise_wrapper.py [options]`

**Input Arguments:**

.. note::

	
	**[REQUIRED]**
		
	-i, `-`-input_file
		Path to flowgram files (.sff.txt), comma separated
	-f, `-`-fasta_file
		Path to fasta file from `split_libraries.py <./split_libraries.html>`_
	
	**[OPTIONAL]**
		
	-o, `-`-output_dir
		Path to output directory [default: denoised_seqs/]
	-n, `-`-num_cpus
		Number of CPUs [default: 1]
	`-`-force_overwrite
		Overwrite files in output directory [default: False]
	-m, `-`-map_fname
		Name of mapping file, Has to contain field LinkerPrimerSequence. [REQUIRED unless --primer specified]
	-p, `-`-primer
		Primer sequence [REQUIRED unless --map_fname specified]
	`-`-titanium
		Select Titanium defaults for denoiser, otherwise use FLX defaults [default: False]


**Output:**

This script results in a OTU like mapping file along with a sequence file of denoised (FASTA-format). Note that the sequences coming from denoising are no real OTUs, and have to be sent to `pick_otus.py <./pick_otus.html>`_ if the users wishes to have a defined similarity threshold.


**Example:**

Denoise flowgrams in file 454Reads.sff.txt, discard flowgrams not in seqs.fna, and extract primer from map.txt:

::

	denoise_wrapper.py -i 454Reads.sff.txt -f seqs.fna -m map.txt

**Multi-core Example:**

Denoise flowgrams in file 454Reads.sff.txt using 2 cores on your machine in parallel:

::

	denoise_wrapper.py -n 2 -i 454Reads.sff.txt -f seqs.fna -m map.txt