/tools/sr_mapping/PerM.xml

https://bitbucket.org/h_morita_dbcls/galaxy-central · XML · 368 lines · 311 code · 33 blank · 24 comment · 0 complexity · fc1ac66d5359f331235a69c8b2b79ba6 MD5 · raw file 

    

  1. <tool id="PerM" name="Map with PerM" version="1.0.0">
    2. 

  2.   <description>for SOLiD and Illumina</description>
    3. 

  3.   <!-- works with PerM version 0.2.6 -->
    4. 

  4.   <command>
    5. 

  5. PerM
    6. 

  6. #if $s.sourceOfRef.refSource == "history":
    7. 

  7.     $s.sourceOfRef.ref
    8. 

  8. #else:
    9. 

  9.     $s.sourceOfRef.index.value
    10. 

  10. #end if
    11. 

  11. #if $s.mate.singleOrPairs == "single":
    12. 

  12.     $s.mate.reads
    13. 

  13. #else:
    14. 

  14.     -1 $s.mate.reads1 -2 $s.mate.reads2
    15. 

  15.     -U $s.mate.upperbound
    16. 

  16.     -L $s.mate.lowerbound
    17. 

  17.     $s.mate.excludeAmbiguousPairs
    18. 

  18. #end if
    19. 

  19. #if $s.space == "color":
    20. 

  20.     --readFormat "csfastq"
    21. 

  21. #else:
    22. 

  22.     --readFormat "fastq"
    23. 

  23. #end if
    24. 

  24. #if $int($str($valAlign)) &gt;= 0:
    25. 

  25.     -v $valAlign
    26. 

  26. #end if
    27. 

  27. #if $align.options == "full":
    28. 

  28.     --seed $align.seed
    29. 

  29.     -$align.alignments
    30. 

  30.     #if $str($align.delimiter) != "None":
    31. 

  31.         --delimiter $align.delimiter
    32. 

  32.     #end if
    33. 

  33.     -T $align.sTrimL
    34. 

  34.     $align.includeReadsWN
    35. 

  35.     $align.statsOnly
    36. 

  36.     $align.ignoreQS
    37. 

  37. #end if
    38. 

  38. #if $str($bUnmappedRead) == "true" and $s.space == "color":
    39. 

  39.   -u $unmappedReadOutCS
    40. 

  40. #elif $str($bUnmappedRead) == "true" and $s.space == "base":
    41. 

  41.   -u $unmappedReadOut
    42. 

  42. #end if
    43. 

  43. -o $output --outputFormat sam --noSamHeader | tr '\r' '\n' | tr -cd "[:print:]\t\n " | grep "Reads\|Sub0\|Pairs\|single" | sed 's/.*Reads:,//' | sed 's/\/.*dat,_ Sub0/Sub0/'
    44. 

  44.   </command>
    45. 

  45.   <inputs>
    46. 

  46.     <conditional name="s">
    47. 

  47.       <param name="space" label="Is your data color space (SOLiD) or base space (Illumina)?" type="select">
    48. 

  48.         <option value="color">Color space</option>
    49. 

  49.         <option value="base">Base space</option>
    50. 

  50.       </param>
    51. 

  51.       <when value="color">
    52. 

  52.         <conditional name="sourceOfRef">
    53. 

  53.           <param name="refSource" label="Will you provide your own reference file from the history or use a built-in index?" type="select">
    54. 

  54.             <option value="indexed">Built-in index</option>
    55. 

  55.             <option value="history">Fasta file from history</option>
    56. 

  56.           </param>
    57. 

  57.           <when value="indexed">
    58. 

  58.             <param name="index" type="select" label="Select a reference genome (with seed and read length)" help="if your genome of interest is not listed - contact Galaxy team">
    59. 

  59.               <options from_file="perm_color_index.loc">
    60. 

  60.                 <column name="value" index="1" />
    61. 

  61.                 <column name="name" index="0" />
    62. 

  62.               </options>
    63. 

  63.             </param>
    64. 

  64.           </when>
    65. 

  65.           <when value="history">
    66. 

  66.             <param name="ref" format="fasta" type="data" label="Reference" />
    67. 

  67.           </when>
    68. 

  68.         </conditional>
    69. 

  69.         <conditional name="mate">
    70. 

  70.           <param name="singleOrPairs" label="Mate-paired?" type="select">
    71. 

  71.             <option value="single">Single-end</option>
    72. 

  72.             <option value="paired">Mate pairs</option>
    73. 

  73.           </param>
    74. 

  74.           <when value="single">
    75. 

  75.             <param format="fastqcssanger" name="reads" type="data" label="Reads" />
    76. 

  76.           </when>
    77. 

  77.           <when value="paired">
    78. 

  78.             <param name="reads1" format="fastqcssanger" label="Forward FASTQ file" type="data" />
    79. 

  79.             <param name="reads2" format="fastqcssanger" label="Reverse FASTQ file" type="data" />
    80. 

  80.             <param label="Upperbound of pairs separation (-U)" name="upperbound" type="integer" size="8" value="100000" />
    81. 

  81.             <param label="Lowerbound of pairs separation (-L)" name="lowerbound" type="integer" size="8" value="0" />
    82. 

  82.             <param label="Exclude ambiguous pairs (-e)" name="excludeAmbiguousPairs" type="boolean" checked="false" truevalue="-e" falsevalue="" />
    83. 

  83.           </when>
    84. 

  84.         </conditional>
    85. 

  85.       </when>
    86. 

  86.       <when value="base">
    87. 

  87.         <conditional name="sourceOfRef">
    88. 

  88.           <param name="refSource" label="Will you provide your own reference file from the history or use a built-in index?" type="select">
    89. 

  89.             <option value="indexed">Built-in index</option>
    90. 

  90.             <option value="history">Fasta file from history</option>
    91. 

  91.           </param>
    92. 

  92.           <when value="indexed">
    93. 

  93.             <param name="index" type="select" label="Select a reference genome with seed and read length" help="if your genome of interest is not listed - contact Galaxy team">
    94. 

  94.               <options from_file="perm_base_index.loc">
    95. 

  95.                 <column name="value" index="1" />
    96. 

  96.                 <column name="name" index="0" />
    97. 

  97.               </options>
    98. 

  98.             </param>
    99. 

  99.           </when>
    100. 

  100.           <when value="history">
    101. 

  101.             <param name="ref" format="fasta" type="data" label="Reference" />
    102. 

  102.           </when>
    103. 

  103.         </conditional>
    104. 

  104.         <conditional name="mate">
    105. 

  105.           <param name="singleOrPairs" label="Mate-paired?" type="select">
    106. 

  106.             <option value="single">Single-end</option>
    107. 

  107.             <option value="paired">Mate pairs</option>
    108. 

  108.           </param>
    109. 

  109.           <when value="single">
    110. 

  110.             <param format="fastqsanger" name="reads" type="data" label="Reads" />
    111. 

  111.           </when>
    112. 

  112.           <when value="paired">
    113. 

  113.             <param name="reads1" format="fastqsanger" label="Forward FASTQ file" type="data" />
    114. 

  114.             <param name="reads2" format="fastqsanger" label="Reverse FASTQ file" type="data" />
    115. 

  115.             <param label="Upperbound of pairs separation (-U)" name="upperbound" type="integer" size="8" value="100000" />
    116. 

  116.             <param label="Lowerbound of pairs separation (-L)" name="lowerbound" type="integer" size="8" value="0" />
    117. 

  117.             <param label="Exclude ambiguous pairs (-e)" name="excludeAmbiguousPairs" type="boolean" checked="false" truevalue="-e" falsevalue="" />
    118. 

  118.           </when>
    119. 

  119.         </conditional>
    120. 

  120.       </when>
    121. 

  121.     </conditional>
    122. 

  122.     <param label="Maximum number of mismatches permitted in one end of full read (-v)" name="valAlign" type="integer" size="5" value="2" />
    123. 

  123.     <conditional name="align">
    124. 

  124.       <param help="Use default setting or specify full parameters list" label="PerM settings to use" name="options" type="select">
    125. 

  125.         <option value="preSet">Commonly used</option>
    126. 

  126.         <option value="full">Full parameter list</option>
    127. 

  127.       </param>
    128. 

  128.       <when value="preSet"/>
    129. 

  129.       <when value="full">
    130. 

  130.         <param label="Whether or not to report all valid alignments per read (-A/-B/-E)" name="alignments" type="select">
    131. 

  131.           <option value="A">Report all valid alignments</option>
    132. 

  132.           <option value="B">Report the best alignments in terms of number of mismatches</option>
    133. 

  133.           <option value="E">Report only uniquely mapped reads</option>
    134. 

  134.         </param>
    135. 

  135.         <param label="Choose the seed full sensitive to different number of mismatches (--seed)" name="seed" type="select" >
    136. 

  136.           <option value="F2">2 mismatches</option>
    137. 

  137.           <option value="S11">1 SNP + 1 color error</option>
    138. 

  138.           <option value="F3">3 mismatches</option>
    139. 

  139.           <option value="F4">4 mismatches</option>
    140. 

  140.         </param>
    141. 

  141.         <param label="Choose the delimiter to identify read name (--delimiter)" name="delimiter" type="select">
    142. 

  142.           <option value="None">Tab/Space/Comma</option>
    143. 

  143.           <option value=":">Colon</option>
    144. 

  144.           <option value="_">Underscore</option>
    145. 

  145.         </param>
    146. 

  146.         <param label="Use the first n bases of each read for alignment (-T)" name="sTrimL" type="integer" size="5" value="50" />
    147. 

  147.         <param name="includeReadsWN" type="boolean" checked="true" truevalue="--includeReadsWN" falsevalue="" label="Include reads with 'N' or '.' by encoding '.' as 3, 'N' as 'A' (--includeReadsWN)" /> 
    148. 

  148.         <param name="statsOnly" type="boolean" checked="false" truevalue="--statsOnly" falsevalue="" label="output mapping stats only. Don't output alignments (--statsOnly)" />
    149. 

  149.         <param name="ignoreQS" type="boolean" checked="false" truevalue="--ignoreQS" falsevalue="" label="Ignore quality scores (--ignoreQS)" />
    150. 

  150.       </when>
    151. 

  151.     </conditional> <!-- options -->
    152. 

  152.     <param name="bUnmappedRead" type="select" label="Output the unmapped reads (-u)">
    153. 

  153.       <option value="true">Yes</option>
    154. 

  154.       <option value="false">No</option>
    155. 

  155.     </param>
    156. 

  156.   </inputs>
    157. 

  157.   <outputs>
    158. 

  158.     <data name="output" format="sam"/> 
    159. 

  159.     <data name="unmappedReadOut" format="fastqsanger">
    160. 

  160.       <filter>bUnmappedRead == "true" and s["space"] == "base"</filter>
    161. 

  161.     </data>
    162. 

  162.     <data name="unmappedReadOutCS" format="fastqcssanger">
    163. 

  163.       <filter>bUnmappedRead == "true" and s["space"] == "color"</filter>
    164. 

  164.     </data>
    165. 

  165.   </outputs>
    166. 

  166.   <tests>
    167. 

  167.     <test>
    168. 

  168.       <!--
    169. 

  169.       PerM command:
    170. 

  170.       PerM test-data/phiX.fasta 50 +seed F3 -m -s phiX_F3_50.index +readFormat .fastq
    171. 

  171.       PerM phiX_F3_50.index -1 test-data/perm_in1.fastqsanger -2 test-data/perm_in2.fastqsanger -U 100000 -L 0 -e +readFormat .fastq -v 0 +seed F2 -A -T 50 +includeReadsWN -o perm_out1.sam +outputFormat sam +noSamHeader | tr '\r' '\n' | tr -cd "[:print:]\t\n " | grep "Reads\|Sub0\|Pairs\|single" | sed 's/.*Reads:,//' | sed 's/\/.*dat,_ Sub0/Sub0/'
    172. 

  172.       You need to replace the + with 2 dashes.
    173. 

  173.       -->
    174. 

  174.       <param name="space" value="base" />
    175. 

  175.       <param name="refSource" value="indexed" />
    176. 

  176.       <param name="index" value="phiX_F3_50" />
    177. 

  177.       <param name="singleOrPairs" value="paired" />
    178. 

  178.       <param name="reads1" value="perm_in1.fastqsanger" ftype="fastqsanger" />
    179. 

  179.       <param name="reads2" value="perm_in2.fastqsanger" ftype="fastqsanger" />
    180. 

  180.       <param name="upperbound" value="100000" />
    181. 

  181.       <param name="lowerbound" value="0" />
    182. 

  182.       <param name="excludeAmbiguousPairs" value="true" />
    183. 

  183.       <param name="valAlign" value="0" />
    184. 

  184.       <param name="options" value="full" />
    185. 

  185.       <param name="alignments" value="A" />
    186. 

  186.       <param name="seed" value="F2" />
    187. 

  187.       <param name="delimiter" value="None" />
    188. 

  188.       <param name="sTrimL" value="50" />
    189. 

  189.       <param name="includeReadsWN" value="true" />
    190. 

  190.       <param name="statsOnly" value="false" />
    191. 

  191.       <param name="ignoreQS" value="false" />
    192. 

  192.       <param name="bUnmappedRead" value="false" />
    193. 

  193.       <output name="output" file="perm_out1.sam" ftype="sam" />
    194. 

  194.     </test>
    195. 

  195.     <test>
    196. 

  196.       <!--
    197. 

  197.       PerM command:
    198. 

  198.       PerM test-data/chr_m.fasta test-data/perm_in3.fastqsanger +readFormat .fastq -v 2 -u perm_out3.fastqsanger -o perm_out2.sam +outputFormat sam +noSamHeader | tr '\r' '\n' | tr -cd "[:print:]\t\n " | grep "Reads\|Sub0\|Pairs\|single" | sed 's/.*Reads:,//' | sed 's/\/.*dat,_ Sub0/Sub0/'
    199. 

  199.       You need to replace the + with 2 dashes.
    200. 

  200.       -->
    201. 

  201.       <param name="space" value="base" />
    202. 

  202.       <param name="refSource" value="history" />
    203. 

  203.       <param name="ref" value="chr_m.fasta" ftype="fasta" />
    204. 

  204.       <param name="singleOrPairs" value="single" />
    205. 

  205.       <param name="reads" value="perm_in3.fastqsanger" ftype="fastqsanger" />
    206. 

  206.       <param name="valAlign" value="2" />
    207. 

  207.       <param name="options" value="preSet" />
    208. 

  208.       <param name="bUnmappedRead" value="true" />
    209. 

  209.       <output name="output" file="perm_out2.sam" ftype="sam" />
    210. 

  210.       <output name="unmappedReadOut" file="perm_out3.fastqsanger" ftype="fastqsanger" />
    211. 

  211.     </test>
    212. 

  212.     <test>
    213. 

  213.       <!--
    214. 

  214.       PerM command:
    215. 

  215.       PerM test-data/phiX.fasta test-data/perm_in4.fastqcssanger +readFormat .csfastq -v 1 -o perm_out4.sam +outputFormat sam +noSamHeader | tr '\r' '\n' | tr -cd "[:print:]\t\n " | grep "Reads\|Sub0\|Pairs\|single" | sed 's/.*Reads:,//' | sed 's/\/.*dat,_ Sub0/Sub0/'
    216. 

  216.       You need to replace the + with 2 dashes.
    217. 

  217.       -->
    218. 

  218.       <param name="space" value="color" />
    219. 

  219.       <param name="refSource" value="history" />
    220. 

  220.       <param name="ref" value="phiX.fasta" ftype="fasta" />
    221. 

  221.       <param name="singleOrPairs" value="single" />
    222. 

  222.       <param name="reads" value="perm_in4.fastqcssanger" ftype="fastqcssanger" />
    223. 

  223.       <param name="valAlign" value="1" />
    224. 

  224.       <param name="options" value="preSet" />
    225. 

  225.       <param name="bUnmappedRead" value="false" />
    226. 

  226.       <output name="output" file="perm_out4.sam" ftype="sam" />
    227. 

  227.     </test>
    228. 

  228.     <test>
    229. 

  229.       <!--
    230. 

  230.       PerM command:
    231. 

  231.       PerM equCab2.fasta 50 +seed F4 -m -s equCab2_F3_50.index +readFormat .csfastq
    232. 

  232.       PerM equCab2_F3_50.index -1 test-data/perm_in5.fastqcssanger -2 test-data/perm_in6.fastqcssanger -U 90000 -L 10000 +readFormat .csfastq -v 3 -o perm_out5.sam +outputFormat sam +noSamHeader | tr '\r' '\n' | tr -cd "[:print:]\t\n " | grep "Reads\|Sub0\|Pairs\|single" | sed 's/.*Reads:,//' | sed 's/\/.*dat,_ Sub0/Sub0/'
    233. 

  233.       You need to replace the + with 2 dashes.
    234. 

  234.       hg19.fasta needs to be supplied.
    235. 

  235.       -->
    236. 

  236.       <param name="space" value="color" />
    237. 

  237.       <param name="refSource" value="indexed" />
    238. 

  238.       <param name="index" value="equCab2_chrM_F3_50" />
    239. 

  239.       <param name="singleOrPairs" value="paired" />
    240. 

  240.       <param name="reads1" value="perm_in5.fastqcssanger" ftype="fastqcssanger" />
    241. 

  241.       <param name="reads2" value="perm_in6.fastqcssanger" ftype="fastqcssanger" />
    242. 

  242.       <param name="upperbound" value="90000" />
    243. 

  243.       <param name="lowerbound" value="10000" />
    244. 

  244.       <param name="excludeAmbiguousPairs" value="false" />
    245. 

  245.       <param name="valAlign" value="3" />
    246. 

  246.       <param name="options" value="preSet" />
    247. 

  247.       <param name="bUnmappedRead" value="false" />
    248. 

  248.       <output name="output" file="perm_out5.sam" ftype="sam" />
    249. 

  249.     </test>
    250. 

  250.   </tests>
    251. 

  251.   <help>
    252. 

  252. **What it does**
    253. 

  253. 

  254. PerM is a short read aligner designed to be ultrafast with long SOLiD reads to the whole genome or transcriptions. PerM can be fully sensitive to alignments with up to four mismatches and highly sensitive to a higher number of mismatches.
    255. 

  255. 

  256. **Development team**
    257. 

  257. 

  258. PerM is developed by Ting Chen's group, Center of Excellence in Genomic Sciences at the University of Southern California. If you have any questions, please email yanghoch at usc.edu or check the `project page`__.
    259. 

    260. 

  260.  .. __: http://code.google.com/p/perm/
    261. 

  261. 

  262. **Citation**
    263. 

  263. 

  264. PerM: Efficient mapping of short sequencing reads with periodic full sensitive spaced seeds. Bioinformatics, 2009, 25 (19): 2514-2521.
    265. 

  265. 

  266. **Input**
    267. 

  267. 

  268. The input files are read files and a reference. Users can use the pre-indexed reference in Galaxy or upload their own reference.
    269. 

  269. 

  270. The uploaded reference file should be in the fasta format. Multiple sequences like transcriptions should be concatenated together separated by a header line that starts with the ">" character.
    271. 

  271. 

  272. Reads files must be in either fastqsanger or fastqcssanger format to use in PerM. However, there are several possible starting formats that can be converted to one of those two: fastq (any type), color-space fastq, fasta, csfasta, or csfasta+qualsolid. 
    273. 

  273. 

  274. An uploaded base-space fastq file MUST be checked/transformed with FASTQGroomer tools in Galaxy to be converted to the fastqsanger format (this is true even if the original file is in Sanger format).
    275. 

  275. 

  276. Uploaded fasta and csfasta without quality score files can be transformed to fastqsanger by the FASTQGroomer, with pseudo quality scores added.
    277. 

  277. 

  278. An uploaded csfasta + qual pair can also be transformed into fastqcssanger by solid2fastq.
    279. 

  279. 

  280. **Outputs**
    281. 

  281. 

  282. The output mapping result is in SAM format, and has the following columns::
    283. 

  283. 

  284.     Column  Description
    285. 

  285.   --------  --------------------------------------------------------
    286. 

  286.    1 QNAME  Query (pair) NAME
    287. 

  287.    2 FLAG   bitwise FLAG
    288. 

  288.    3 RNAME  Reference sequence NAME
    289. 

  289.    4 POS    1-based leftmost POSition/coordinate of clipped sequence
    290. 

  290.    5 MAPQ   MAPping Quality (Phred-scaled)
    291. 

  291.    6 CIGAR  extended CIGAR string
    292. 

  292.    7 MRNM   Mate Reference sequence NaMe ('=' if same as RNAME)
    293. 

  293.    8 MPOS   1-based Mate POSition
    294. 

  294.    9 ISIZE  Inferred insert SIZE
    295. 

  295.   10 SEQ    query SEQuence on the same strand as the reference
    296. 

  296.   11 QUAL   query QUALity (ASCII-33 gives the Phred base quality)
    297. 

  297.   12 OPT    variable OPTional fields in the format TAG:VTYPE:VALUE
    298. 

  298.   12.1 NM   Number of mismatches (SOLiD-specific)
    299. 

  299.   12.2 CS   Reads in color space (SOLiD-specific)
    300. 

  300.   12.3 CQ   Bases quality in color spacehidden="true" (SOLiD-specific)
    301. 

  301. 

  302. The flags are as follows::
    303. 

  303. 

  304.     Flag  Description
    305. 

  305.   ------  -------------------------------------
    306. 

  306.   0x0001  the read is paired in sequencing
    307. 

  307.   0x0002  the read is mapped in a proper pair
    308. 

  308.   0x0004  the query sequence itself is unmapped
    309. 

  309.   0x0008  the mate is unmapped
    310. 

  310.   0x0010  strand of the query (1 for reverse)
    311. 

  311.   0x0020  strand of the mate
    312. 

  312.   0x0040  the read is the first read in a pair
    313. 

  313.   0x0080  the read is the second read in a pair
    314. 

  314.   0x0100  the alignment is not primary
    315. 

  315. 

  316. Here is some sample output::
    317. 

  317. 

  318.   Qname	FLAG	Rname	POS	MAPQ	CIAGR	MRNM	MPOS	ISIZE	SEQ	QUAL	NM	CS	CQ
    319. 

  319.   491_28_332_F3   16      ref-1   282734  255     35M     *       0       0       AGTCAAACTCCGAATGCCAATGACTTATCCTTAGG    #%%%%%%%!!%%%!!%%%%%%%%!!%%%%%%%%%%      NM:i:3  CS:Z:C0230202330012130103100230121001212        CQ:Z:###################################
    320. 

  320.   491_28_332_F3   16      ref-1   269436  255     35M     *       0       0       AGTCAAACTCCGAATGCCAATGACTTATCCTTAGG    #%%%%%%%!!%%%!!%%%%%%%%!!%%%%%%%%%%      NM:i:3  CS:Z:C0230202330012130103100230121001212        CQ:Z:###################################
    321. 

  321. 

  322. The user can check a checkbox for optional output containing the unmmaped reads in fastqsanger or fastqcssanger. The default is to produce it.
    323. 

  323. 

  324. **PerM parameter list**
    325. 

  325. 

  326. Below is a list of PerM command line options for PerM. Not all of these are relevant to Galaxy's implementation, but are included for completeness.
    327. 

    328. 

  328. The command for single-end::
    329. 

  329. 

  330.   PerM [ref_or_index] [read] [options]
    331. 

  331. 

  332. The command for paired-end::
    333. 

  333. 

  334.   PerM [ref_or_index] -1 [read1] -2 [read1] [options]
    335. 

  335. 

  336. The command-line options::
    337. 

  337. 

  338.   -A                Output all alignments within the given mismatch threshold, end-to-end.
    339. 

  339.   -B                Output best alignments in terms of mismatches in the given mismatch threshold. [Default]
    340. 

  340.   -E                Output only the uniquely mapped reads in the given mismatch threshold.
    341. 

  341.   -m                Create the reference index, without reusing the saved index.
    342. 

  342.   -s PATH           Save the reference index to accelerate the mapping in the future. If PATH is not specified, the default path will be used.
    343. 

  343.   -v INT            Where INT is the number of mismatches allowed in one end. [Default=2]
    344. 

  344.   -T INT            Where INT is the length to truncate read length to, so 30 means use only first 30 bases (signals). Leave blank if the full read is meant to be used.
    345. 

  345.   -o PATH           Where PATH is for output the mapping of one read set. PerM's output are in .mapping or .sam format, determined by the ext name of PATH. Ex: -o out.sam will output in SAM format; -o out.mapping will output in .mapping format.
    346. 

  346.   -d PATH           Where PATH is the directory for multiple read sets.
    347. 

  347.   -u PATH           Print the fastq file of those unmapped reads to the file in PATH.
    348. 

  348.   --noSamHeader     Print no SAM header so it is convenient to concatenate multiple SAM output files.
    349. 

  349.   --includeReadsWN  Encodes N or "." with A or 3, respectively.
    350. 

  350.   --statsOnly       Output the mapping statistics in stdout only, without saving alignments to files.
    351. 

  351.   --ignoreQS        Ignore the quality scores in fastq or QUAL files.
    352. 

  352.   --seed {F2 | S11 | F3 | F4}    Specify the seed pattern, which has a specific full sensitivity. Check the algorithm page (link below) for seed patterns to balance the sensitivity and running time.
    353. 

  353.   --readFormat {fasta | fastq | csfasta | csfastq}    Read in reads in the specified format, instead of guessing according to the extension name.
    354. 

  354.   --delimiter CHAR  Which is a character used as the delimiter to separate the the read id, and the additional info in the line with ">" in fasta or csfasta.
    355. 

  355. 

  356. Paired reads options::
    357. 

  357. 

  358.   -e        Exclude ambiguous paired.
    359. 

  359.   -L INT    Mate-paired separate lower bound.
    360. 

  360.   -U INT    Mate-paired separate upper bound.
    361. 

  361.   -1 PATH   The forward reads file path.
    362. 

  362.   -2 PATH   The reversed reads file path.
    363. 

  363. 

  364. See the PerM `algorithm page`__ for information on algorithms and seeds.
    365. 

  365. 

  366.  .. __: http://code.google.com/p/perm/wiki/Algorithms
    367. 

  367.   </help>
    368. 

  368. </tool>
    369.