annotate tools/sr_mapping/lastz_paired_reads_wrapper.xml @ 1:cdcb0ce84a1b

Uploaded
author xuebing
date Fri, 09 Mar 2012 19:45:15 -0500
parents 9071e359b9a3
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
1 <tool id="lastz_paired_reads_wrapper" name="Lastz paired reads" version="1.1.1">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
2 <description> map short paired reads against reference sequence</description>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
3 <command interpreter="python">lastz_paired_reads_wrapper.py
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
4 #if $seq_name.how_to_name=="yes":
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
5 --ref_name=$seq_name.ref_name
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
6 #end if
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
7 --ref_source=$source.ref_source
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
8 --input2=$input2
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
9 --input3=$input3
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
10 --input4=$input4
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
11 #if $source.ref_source=="history":
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
12 --input1=$source.input1
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
13 --ref_sequences=$input1.metadata.sequences
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
14 #else:
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
15 --input1="${ filter( lambda x: str( x[0] ) == str( $source.input1_2bit ), $__app__.tool_data_tables[ 'lastz_seqs' ].get_fields() )[0][-1] }"
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
16 #end if
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
17 --output=$output1
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
18 --lastz_seqs_file_dir=${GALAXY_DATA_INDEX_DIR}
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
19 </command>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
20 <inputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
21 <param name="input2" format="fasta" type="data" label="Align sequencing reads in" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
22 <conditional name="source">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
23 <param name="ref_source" type="select" label="Against reference sequences that are">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
24 <option value="cached">locally cached</option>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
25 <option value="history">in your history</option>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
26 </param>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
27 <when value="cached">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
28 <param name="input1_2bit" type="select" label="Using reference genome" help="If your genome of interest is not listed, contact the Galaxy team">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
29 <options from_data_table="lastz_seqs" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
30 </param>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
31 </when>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
32 <when value="history">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
33 <param name="input1" type="data" format="fasta" label="Select a reference dataset" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
34 </when>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
35 </conditional>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
36 <param name="input3" format="fasta" type="data" label="Linker file" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
37 <param name="input4" format="qual454" type="data" label="Select a base quality score 454 dataset" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
38 <conditional name="seq_name">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
39 <param name="how_to_name" type="select" label="Do you want to modify the reference name?">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
40 <option value="no">No</option>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
41 <option value="yes">Yes</option>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
42 </param>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
43 <when value="yes">
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
44 <param name="ref_name" type="text" size="25" value="Type sequence name here" label="Enter name for the Reference sequence"/>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
45 </when>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
46 <when value="no" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
47 </conditional>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
48 </inputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
49 <outputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
50 <data format="sam" name="output1" label="${tool.name} on ${on_string}: mapped reads" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
51 </outputs>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
52 <requirements>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
53 <requirement type="package">lastz</requirement>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
54 </requirements>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
55 <tests>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
56 <test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
57 <!--
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
58 input1: a reference genome ( 2bit or fasta )
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
59 input2: a collection of 454 paired end reads ( a fasta file )
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
60 input3: a linker sequence ( a very small fasta file )
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
61 input4: a base quality score 454 file ( qual454 )
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
62 -->
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
63 <param name="input2" value="lastz_paired_input2.fasta" ftype="fasta" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
64 <param name="ref_source" value="cached" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
65 <param name="input1_2bit" value="/galaxy/data/hg18/seq/chr21.2bit" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
66 <param name="input3" value="lastz_paired_input3.fasta" ftype="fasta" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
67 <param name="input4" value="lastz_paired_input4.qual454" ftype="qual454" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
68 <param name="how_to_name" value="no" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
69 <output name="output1" file="lastz_paired_out1.sam" />
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
70 </test>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
71 </tests>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
72 <help>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
73
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
74 **What it does**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
75
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
76 **LASTZ** is a high performance pairwise sequence aligner derived from BLASTZ. It is written by Bob Harris in Webb Miller's laboratory at Penn State University. Special scoring sets were derived to improve runtime performance and quality. This Galaxy version of LASTZ is geared towards aligning short (Illumina/Solexa, AB/SOLiD) and medium (Roche/454) paired reads against a reference sequence. There is excellent, extensive documentation on LASTZ available here_.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
77
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
78 .. _here: http://www.bx.psu.edu/miller_lab/dist/README.lastz-1.02.00/README.lastz-1.02.00.html
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
79
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
80 ------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
81
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
82 **Input formats**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
83
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
84 LASTZ accepts reference and reads in FASTA format. However, because Galaxy supports implicit format conversion the tool will recognize fastq and other method specific formats.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
85
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
86 ------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
87
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
88 **Outputs**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
89
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
90 This LASTZ tool produces a SAM file showing sequence alignments.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
91
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
92 **SAM output**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
93
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
94 SAM has 12 columns::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
95
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
96 1 2 3 4 5 6 7 8 9 10 11 12
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
97 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
98 HWI-EAS91_1_30788AAXX:1:2:1670:915 99 chr9 58119878 60 36M = 58120234 392 GACCCCTACCCCACCGTGCTCTGGATCTCAGTGTTT IIIIIIIIIIIIIIIIEIIIIIII7IIIIIIIIIII XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
99 HWI-EAS91_1_30788AAXX:1:2:1670:915 147 chr9 58120234 60 36M = 58119878 -392 ATGAGTCGAATTCTATTTTCCAAACTGTTAACAAAA IFIIDI;IIICIIIIIIIIIIIIIIIIIIIIIIIII XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
100
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
101
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
102 where::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
103
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
104 Column Description
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
105 --------- ---------------------------------------------------------------------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
106 1. QNAME Query (pair) NAME
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
107 2. FLAG bitwise FLAG
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
108 3. RNAME Reference sequence NAME
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
109 4. POS 1-based leftmost POSition/coordinate of clipped sequence
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
110 5. MAPQ MAPping Quality (Phred-scaled)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
111 6. CIGAR extended CIGAR string
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
112 7. MRNM Mate Reference sequence NaMe ('=' if same as RNAME)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
113 8. MPOS 1-based Mate POSition
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
114 9. ISIZE Inferred insert SIZE
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
115 10. SEQ query SEQuence on the same strand as the reference
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
116 11. QUAL query QUALity (ASCII-33 gives the Phred base quality)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
117 12. OPT variable OPTional fields in the format TAG:VTYPE:VALUE, tab-separated
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
118
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
119 The flags are as follows::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
120
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
121 Flag Description
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
122 ------ -------------------------------------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
123 0x0001 the read is paired in sequencing
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
124 0x0002 the read is mapped in a proper pair
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
125 0x0004 the query sequence itself is unmapped
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
126 0x0008 the mate is unmapped
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
127 0x0010 strand of the query (1 for reverse)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
128 0x0020 strand of the mate
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
129 0x0040 the read is the first read in a pair
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
130 0x0080 the read is the second read in a pair
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
131 0x0100 the alignment is not primary
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
132
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
133 ------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
134
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
135 **Do you want to modify the reference name?**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
136
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
137 This option allows you to set the name of the reference sequence manually. This is helpful when, for example, you would like to make the reference name compatible with the UCSC naming conventions to be able to display your lastz results as a custom track at the UCSC Genome Browser.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
138
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
139 ------
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
140
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
141 **LASTZ parameter list**
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
142
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
143 This is an exhaustive list of LASTZ options. Once again, please note that not all parameters are included in this interface. If you would like to make additional options available through Galaxy, e-mail us at galaxy-bugs@bx.psu.edu::
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
144
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
145 target[[s..e]][-] spec/file containing target sequence (fasta or nib)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
146 [s..e] defines a subrange of the file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
147 - indicates reverse-complement
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
148 (use --help=files for more details)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
149 query[[s..e]][-] spec/file containing query sequences (fasta or nib)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
150 if absent, queries come from stdin (unless they
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
151 aren't needed, as for --self or --tableonly)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
152 (use --help=files for more details)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
153 --self the target sequence is also the query
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
154 --quantum the query sequence contains quantum DNA
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
155 --seed=match&lt;length&gt; use a word with no gaps instead of a seed pattern
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
156 --seed=half&lt;length&gt; use space-free half-weight word instead of seed pattern
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
157 --match=&lt;reward&gt;[,&lt;penalty&gt;] set the score values for a match (+&lt;reward&gt;)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
158 and mismatch (-&lt;penalty&gt;)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
159 --[no]trans[ition][=2] allow one or two transitions in a seed hit
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
160 (by default a transition is allowed)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
161 --word=&lt;bits&gt; set max bits for word hash; use this to trade time for
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
162 memory, eliminating thrashing for heavy seeds
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
163 (default is 28 bits)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
164 --[no]filter=[&lt;T&gt;:]&lt;M&gt; filter half-weight seed hits, requiring at least M
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
165 matches and allowing no more than T transversions
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
166 (default is no filtering)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
167 --notwins require just one seed hit
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
168 --twins=[&lt;min&gt;:]&lt;maxgap&gt; require two nearby seed hits on the same diagonal
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
169 (default is twins aren't required)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
170 --notwins allow single, isolated seeds
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
171 --[no]recoverseeds avoid losing seeds in hash collisions. Cannot be used with --twins
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
172 --seedqueue=&lt;entries&gt; set number of entries in seed hit queue
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
173 (default is 262144)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
174 --anchors=&lt;file&gt; read anchors from a file, instead of discovering anchors
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
175 via seeding
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
176 --recoverhits recover hash-collision seed hits
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
177 (default is not to recover seed hits)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
178 --step=&lt;length&gt; set step length (default is 1)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
179 --maxwordcount=&lt;limit&gt; words occurring more often than &lt;limit&gt; in the target
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
180 are not eligible for seeds
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
181 --strand=both search both strands
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
182 --strand=plus search + strand only (matching strand of query spec)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
183 --strand=minus search - strand only (opposite strand of query spec)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
184 (by default both strands are searched)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
185 --ambiguousn treat N as an ambiguous nucleotide
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
186 (by default N is treated as a sequence splicing character)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
187 --[no]gfextend perform gap-free extension of seed hits to HSPs
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
188 (by default no extension is performed)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
189 --[no]chain perform chaining
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
190 --chain=&lt;diag,anti&gt; perform chaining with given penalties for diagonal and
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
191 anti-diagonal
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
192 (by default no chaining is performed)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
193 --[no]gapped perform gapped alignment (instead of gap-free)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
194 (by default gapped alignment is performed)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
195 --score[s]=&lt;file&gt; read substitution scores from a file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
196 (default is HOXD70)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
197 --unitscore[s] scores are +1/-1 for match/mismatch
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
198 --gap=&lt;[open,]extend&gt; set gap open and extend penalties (default is 400,30)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
199 --xdrop=&lt;score&gt; set x-drop threshold (default is 10*sub[A][A])
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
200 --ydrop=&lt;score&gt; set y-drop threshold (default is open+300extend)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
201 --infer[=&lt;control&gt;] infer scores from the sequences, then use them
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
202 --inferonly[=&lt;control&gt;] infer scores, but don't use them (requires --infscores)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
203 all inference options are read from the control file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
204 --infscores[=&lt;file&gt;] write inferred scores to a file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
205 --hspthresh=&lt;score&gt; set threshold for high scoring pairs (default is 3000)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
206 ungapped extensions scoring lower are discarded
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
207 &lt;score&gt; can also be a percentage or base count
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
208 --entropy adjust for entropy when qualifying HSPs in the x-drop extension
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
209 method
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
210 --noentropy don't adjust for entropy when qualifying HSPs
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
211 --exact=&lt;length&gt; set threshold for exact matches
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
212 if specified, exact matches are found rather than high
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
213 scoring pairs (replaces --hspthresh)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
214 --inner=&lt;score&gt; set threshold for HSPs during interpolation
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
215 (default is no interpolation)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
216 --gappedthresh=&lt;score&gt; set threshold for gapped alignments
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
217 gapped extensions scoring lower are discarded
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
218 &lt;score&gt; can also be a percentage or base count
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
219 (default is to use same value as --hspthresh)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
220 --ball=&lt;score&gt; set minimum score required of words 'in' a quantum ball
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
221 --[no]entropy involve entropy in filtering high scoring pairs
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
222 (default is "entropy")
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
223 --[no]mirror report/use mirror image of all gap-free alignments
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
224 (default is "mirror" for self-alignments only)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
225 --traceback=&lt;bytes&gt; space for trace-back information
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
226 (default is 80.0M)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
227 --masking=&lt;count&gt; mask any position in target hit this many times
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
228 zero indicates no masking
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
229 (default is no masking)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
230 --targetcapsule=&lt;capsule_file&gt; the target seed word position table and seed
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
231 (as well as the target sequence)are read from specified file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
232 --segments=&lt;segment_file&gt; read segments from a file, instead of discovering
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
233 them via seeding. Replaces other seeding or gap-free extension
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
234 options
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
235 --[no]census[=&lt;file&gt;] count/report how many times each target base aligns
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
236 (default is to not report census)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
237 --identity=&lt;min&gt;[..&lt;max&gt;] filter alignments by percent identity
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
238 0&lt;=min&lt;=max&lt;=100; blocks (or HSPs) outside min..max
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
239 are discarded
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
240 (default is no identity filtering)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
241 --coverage=&lt;min&gt;[..&lt;max&gt;] filter alignments by percentage pf query covered
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
242 0&lt;=min&lt;=max&lt;=100; blocks (or HSPs) outside min..max
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
243 are discarded
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
244 (default is no query coverage filtering)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
245 --notrivial do not output trivial self-alignment block if the target and query
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
246 sequences are identical. Using --self enables this option automatically
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
247 --output=&lt;output_file&gt; write the alignments to the specified file name instead of stdout
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
248 --code=&lt;file&gt; give quantum code for query sequence (only for display)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
249 --format=&lt;type&gt; specify output format; one of lav, axt, maf, maf+, maf-, text,
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
250 lav+text, cigar, text, rdplot, general, or general:&lt;fields&gt;
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
251 (by default output is LAV)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
252 --rdotplot=&lt;file&gt; create an additional output file suitable for plotting the alignments
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
253 with the R statistical package.
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
254 --markend Just before normal completion, write "# lastz end-of-file" to output file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
255 --census[=&lt;output_file&gt;] count and report how many times each target base aligns, up
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
256 to 255. Ns are included in the count
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
257 --census16[=&lt;output_file&gt;] count and report how many times each target base aligns, up
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
258 up 65 thousand
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
259 --census32[=&lt;output_file&gt;] count and report how many times each target bas aligns, up
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
260 to 4 billion
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
261 --writecapsule=&lt;capsule_file&gt; just write out a targegt capsule file and quit; don't
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
262 search for seeds or perform subsequent stages
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
263 --verbosity=&lt;level&gt; set info level (0 is minimum, 10 is everything)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
264 (default is 0)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
265 --[no]runtime report runtime in the output file
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
266 (default is to not report runtime)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
267 --tableonly[=count] just produce the target position table, don't
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
268 search for seeds
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
269 --[no]stats[=&lt;file&gt;] show search statistics (or don't)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
270 (not available in this build)
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
271 --version report the program version and quit
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
272 --help list all options
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
273 --help=files list information about file specifiers
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
274 --help=short[cuts] list blastz-compatible shortcuts
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
275 --help=yasra list yasra-specific shortcuts
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
276
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
277 </help>
9071e359b9a3 Uploaded
xuebing
parents:
diff changeset
278 </tool>