File: filter_tree.rst

package info (click to toggle)
qiime 1.8.0%2Bdfsg-4
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 130,508 kB
  • ctags: 10,145
  • sloc: python: 110,826; haskell: 379; sh: 169; makefile: 125
file content (60 lines) | stat: -rw-r--r-- 1,845 bytes parent folder | download
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
.. _filter_tree:

.. index:: filter_tree.py

*filter_tree.py* -- This script prunes a tree based on a set of tip names
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

**Description:**

This script takes a tree and a list of OTU IDs (in one of several supported formats) and outputs a subtree retaining only the tips on the tree which are found in the inputted list of OTUs (or not found, if the --negate option is provided).


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

**Input Arguments:**

.. note::

	
	**[REQUIRED]**
		
	-i, `-`-input_tree_filepath
		Input tree filepath
	-o, `-`-output_tree_filepath
		Output tree filepath
	
	**[OPTIONAL]**
		
	-n, `-`-negate
		If negate is True will remove input tips/seqs, if    negate is False, will retain input tips/seqs [default: False]
	-t, `-`-tips_fp
		A list of tips (one tip per line) or sequence identifiers   (tab-delimited lines with a seq identifier in the first field)   which should be retained   [default: None]
	-f, `-`-fasta_fp
		A fasta file where the seq ids should be retained [default: None]


**Output:**

Output is a pruned tree in newick format.


**Prune a tree to include only the tips in tips_to_keep.txt:**

::

	filter_tree.py -i rep_seqs.tre -t tips_to_keep.txt -o pruned.tre

**Prune a tree to remove the tips in tips_to_remove.txt. Note that the -n/--negate option must be passed for this functionality:**

::

	filter_tree.py -i rep_seqs.tre -t tips_to_keep.txt -o negated.tre -n

**Prune a tree to include only the tips found in the fasta file provided:**

::

	filter_tree.py -i rep_seqs.tre -f fast_f.fna -o pruned_fast.tre