Mercurial > repos > iuc > pear

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/.shed.yml	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,3 @@
+# repository published to https://toolshed.g2.bx.psu.edu/repos/iuc/pear
+owner: iuc
+name: pear
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/pear.xml	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,204 @@
+<tool id="iuc_pear" name="Pear" version="0.9.6.0">
+    <description>Paired-End read merger</description>
+    <!--<version_command>bismark version</version_command>-->
+    <requirements>
+        <requirement type="package" version="0.9.6">pear</requirement>
+    </requirements>
+    <stdio>
+        <exit_code range="1:" />
+        <exit_code range=":-1" />
+        <regex match="Error:" />
+        <regex match="Exception:" />
+    </stdio>
+    <command>
+<![CDATA[
+    pear
+        -f $forward
+        -r $reverse
+        --output pear
+        --p-value $pvalue
+        --min-overlap $min_overlap
+        #if int($max_assembly_length) > 0:
+            --max-asm-length $max_assembly_length
+        #end if
+        --min-asm-length $min_assembly_length
+        --min-trim-length $min_trim_length
+        --quality-theshold $quality_threshold
+        --max-uncalled-base $max_uncalled_base
+        --test-method $test_method
+        --empirical-freqs $empirical_freqs
+        -j "\${GALAXY_SLOTS:-8}"
+        --score-method $score_method
+        --cap $cap
+        $nbase
+]]>
+    </command>
+    <inputs>
+        <param name="forward" type="data" format="fastqillumina, fastqsanger, fastq"
+            label="Specify the name of file that contains the forward paired-end reads" help="-f" />
+        <param name="reverse" type="data" format="fastqillumina, fastqsanger, fastq"
+            label="Specify the name of file that contains the reverse paired-end reads" help="-r" />
+
+        <!-- optional -->
+        <param name="pvalue" type="float" value="0.01" min="0" optional="True" max="1" label="Specify a p-value for the statistical test"
+            help="If the computed p-value of a possible assembly exceeds the specified p-value then the paired-end read will not be assembled. Setting 1.0 disables the test. (--p-value)" />
+
+        <param name="min_overlap" type="integer" value="10" optional="True" label="Specify the minimum overlap size"
+            help="The minimum overlap may be set to 1 when the statistical test is used. However, further restricting the minimum overlap size to a proper value may reduce false-positive assembles. (--min-overlap)" />
+
+        <param name="max_assembly_length" type="integer" value="0" optional="True" label="Specify the maximum possible length of the assembled sequences"
+            help="Setting this value to 0 disables the restriction and assembled sequences may be arbitrary long. (--max-assembly-length)" />
+
+        <param name="min_assembly_length" type="integer" value="50" optional="True" label="Specify the minimum possible length of the assembled sequences"
+            help="Setting this value to 0 disables the restriction and assembled sequences may be arbitrary short. (--min-assembly-length)" />
+
+        <param name="min_trim_length" type="integer" value="1" optional="True" label="Specify the minimum length of reads after trimming the low quality part"
+            help="See option -q. (--min-trim-length)" />
+
+        <param name="quality_threshold" type="integer" value="0" optional="True" label="Specify the quality score threshold for trimming the low quality part of a read"
+                help="If the quality scores of two consecutive bases are strictly less than the specified threshold, the rest of the read will be trimmed. (--quality-threshold)" />
+
+        <param name="max_uncalled_base" type="float" value="1.0" min="0" optional="True" max="1" label="Specify the maximal proportion of uncalled bases in a read"
+            help="Setting this value to 0 will cause PEAR to discard all reads containing uncalled bases. The other extreme setting is 1 which causes PEAR to process all reads independent on the number of uncalled bases. (--max-uncalled-base)" />
+
+        <param name="cap" type="integer" value="40" optional="True" label="Specify  the upper bound for the resulting quality score"
+            help="If set to zero, capping is disabled. (--cap)" />
+
+        <param name="test_method" type="select" label="Specify the type of statistical test" help="(--test-method)">
+            <option value="1" selected="True">Given the minimum allowed overlap, test using the highest OES (1)</option>
+            <option value="2">Use the acceptance probability (2)</option>
+        </param>
+
+        <param name="empirical_freqs" type="boolean" truevalue="-e" falsevalue="" checked="false"
+            label="Disable empirical base frequencies" help="(--empirical-freqs)" />
+        <param name="nbase" type="boolean" truevalue="--nbase" falsevalue="" checked="false"
+            label="Use N base if uncertain" help="When  merging a base-pair that consists of two non-equal bases out of which none is degenerate, set the merged base to N and use the highest quality score of the two bases. (--nbase)" />
+
+        <param name="score_method" type="select" label="Specify the scoring method" help="(--score-method)">
+            <option value="1">OES with +1 for match and -1 for mismatch</option>
+            <option value="2" selected="True">Assembly score (AS) use +1 for match and -1 for mismatch multiplied by base quality scores</option>
+            <option value="3">Ignore quality scores and use +1 for a match and -1 for a mismatch</option>
+        </param>
+
+        <param name="outputs" type="select" display="checkboxes" multiple="True" label="Output files">
+            <option value="assembled" selected="True">Assembled reads</option>
+            <option value="forward">Forward unassembled reads</option>
+            <option value="reverse">Reverse unassembled reads</option>
+            <option value="discarded">Discarded reads</option>
+
+            <validator type="no_options" message="Please select at least one output file." />
+        </param>
+    </inputs>
+    <outputs>
+        <data format="fastq" name="assembled_reads" from_work_dir="pear.assembled.fastq" label="${tool.name} on ${on_string}: Assembled reads">
+            <filter>'assembled' in outputs</filter>
+        </data>
+        <data format="fastq" name="unassembled_forward_reads" from_work_dir="pear.unassembled.forward.fastq" label="${tool.name} on ${on_string}: Unassembled forward reads">
+            <filter>'forward' in outputs</filter>
+        </data>
+        <data format="fastq" name="unassembled_reverse_reads" from_work_dir="pear.unassembled.reverse.fastq" label="${tool.name} on ${on_string}: Unassembled reverse reads">
+            <filter>'reverse' in outputs</filter>
+        </data>
+        <data format="fastq" name="discarded_reads" from_work_dir="pear.discarded.fastq" label="${tool.name} on ${on_string}: Discarded reads">
+            <filter>'discarded' in outputs</filter>
+        </data>
+    </outputs>
+    <tests>
+        <test>
+            <param name="forward" value="forward.fastq" ftype="fastq" />
+            <param name="reverse" value="reverse.fastq" ftype="fastq" />
+            <param name="min_overlap" value="10" />
+            <param name="min_assembly_length" value="50" />
+            <param name="cap" value="0" />
+            <param name="outputs" value="assembled,forward" />
+            <output name="assembled_reads" file="pear_assembled_results1.fastq" ftype="fastq"/>
+            <output name="unassembled_forward_reads" file="pear_unassembled_forward_results1.fastq" ftype="fastq"/>
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+
+**What it does**
+
+PEAR is an ultrafast, memory-efficient and highly accurate pair-end read merger.
+It is fully parallelized and can run with as low as just a few kilobytes of memory.
+
+
+**Parameters**
+
+
+Standard (mandatory):
+  -f, --forward-fastq         <str>     Forward paired-end FASTQ file.
+  -r, --reverse-fastq         <str>     Reverse paired-end FASTQ file.
+
+Optional:
+  -p, --p-value               <float>   Specify  a p-value for the statistical test. If the computed
+                                        p-value of a possible assembly exceeds the specified p-value
+                                        then  paired-end  read  will not be assembled. Valid options
+                                        are: 0.0001, 0.001, 0.01, 0.05 and 1.0. Setting 1.0 disables
+                                        the test. (default: 0.01)
+  -v, --min-overlap           <int>     Specify the minimum overlap size. The minimum overlap may be
+                                        set to 1 when the statistical test is used. However, further
+                                        restricting  the  minimum overlap size to a proper value may
+                                        reduce false-positive assembles. (default: 10)
+  -m, --max-assembly-length   <int>     Specify   the  maximum  possible  length  of  the  assembled
+                                        sequences.  Setting this value to 0 disables the restriction
+                                        and assembled sequences may be arbitrary long. (default: 0)
+  -n, --min-assembly-length   <int>     Specify   the  minimum  possible  length  of  the  assembled
+                                        sequences.  Setting this value to 0 disables the restriction
+                                        and  assembled  sequences  may be arbitrary short. (default:
+                                        50)
+  -t, --min-trim-length       <int>     Specify  the  minimum length of reads after trimming the low
+                                        quality part (see option -q). (default: 1)
+  -q, --quality-threshold     <int>     Specify  the  quality  score  threshold for trimming the low
+                                        quality  part  of  a  read.  If  the  quality  scores of two
+                                        consecutive  bases  are  strictly  less  than  the specified
+                                        threshold,  the  rest of the read will be trimmed. (default:
+                                        0)
+  -u, --max-uncalled-base     <float>   Specify  the maximal proportion of uncalled bases in a read.
+                                        Setting this value to 0 will cause PEAR to discard all reads
+                                        containing  uncalled  bases.  The other extreme setting is 1
+                                        which  causes  PEAR  to process all reads independent on the
+                                        number of uncalled bases. (default: 1)
+  -g, --test-method           <int>     Specify  the  type  of  statistical  test.  Two  options are
+                                        available. (default: 1)
+                                        1: Given the minimum allowed overlap, test using the highest
+                                        OES. Note that due to its discrete nature, this test usually
+                                        yields  a lower p-value for the assembled read than the cut-
+                                        off  (specified  by -p). For example, setting the cut-off to
+                                        0.05  using  this  test,  the  assembled reads might have an
+                                        actual p-value of 0.02.
+
+                                        2. Use the acceptance probability (m.a.p). This test methods
+                                        computes  the same probability as test method 1. However, it
+                                        assumes  that  the  minimal  overlap is the observed overlap
+                                        with  the  highest  OES, instead of the one specified by -v.
+                                        Therefore,  this  is  not  a  valid statistical test and the
+                                        'p-value'  is  in fact the maximal probability for accepting
+                                        the assembly. Nevertheless, we observed in practice that for
+                                        the case the actual overlap sizes are relatively small, test
+                                        2  can  correctly  assemble  more  reads  with only slightly
+                                        higher false-positive rate.
+  -e, --empirical-freqs                 Disable  empirical base frequencies. (default: use empirical
+                                        base frequencies)
+  -s, --score-method          <int>     Specify the scoring method. (default: 2)
+                                        1. OES with +1 for match and -1 for mismatch.
+                                        2: Assembly score (AS). Use +1 for match and -1 for mismatch
+                                        multiplied by base quality scores.
+                                        3: Ignore quality scores and use +1 for a match and -1 for a
+                                        mismatch.
+  -c, --cap                   <int>     Specify  the upper bound for the resulting quality score. If
+                                        set to zero, capping is disabled. (default: 40)
+  -z, --nbase                           When  merging  a  base-pair  that  consists of two non-equal
+                                        bases  out  of which none is degenerate, set the merged base
+                                        to N and use the highest quality score of the two bases
+
+
+
+
+]]>
+  </help>
+  <citations>
+      <citation type="doi">10.1093/bioinformatics/btt593</citation>
+  </citations>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/forward.fastq	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,56 @@
+@EU861894-140/1
+CAACCATTCCACTCGCCTCTGCCGCACTCGAGGCCCCCAGTATGAGAGGCAGTTACGAGGTTGAACCTCAGAATCTCACCCCACACTTAAGGGCGCGCCT
++
+@#CAAF4FHFG;DIJJI?JG#GIIIDECD#4H*1IBIIIH#?IDID*DFIIF?C#HGBJICJBC#JEGEA>ECH#=7#DDF@##D#D<@C#CC@#C?D@,
+@EU861894-138/1
+TCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAACGGCCCGACTCCGCGCTCCTTACGCCCAGTAAGTTCGGACAACGCTGGCCCCCTACGTACTAC
++
+@CCFFFBFFDHFHJJJ+<HJJGIGGBJJIDJ?JIJII#I5D*GB#EJ#FEJEIEI#GG#AIIDIEAD?@#AD3HG3;E;D7##FCDH0?#A,DD#A#ADA
+@EU861894-136/1
+GCTGATCGTCCTCTCAGACCACCTACCCGTCAAAGCCTTGGTGCGACGTAACTCCCCCAACTAGACGATAGGCCAAGGGCCTATCCCGAAGTGCCCGAAG
++
+@@C?F;FF#HGHGJC<JJJC<#CJG61GHIE#?J:JJHJJHFJ#?#1JJ#<J(DHFEIJJD=JF##DIECH)GD(#D<HC>#C@B9DEFDD@ED#EDBD>
+@EU861894-134/1
+TTGTGCAAGAGTCCTAATTGCTGCCCCCCGTAGGAGTCTGGGCCGTATCTCAGTCCCAATCCGGCTGTTCGTCCCCTCAGACCATGTACCCGTAGATGCC
++
+=CCFDFFFA?#BHIJJJCIJ<JIJIGHJ8JHCI9JJIDJIGIJEFI9@JIJ#FFI5J>JJJD;EEJB#@2GFIE#>A;C@FCA@##BD>AD=D#D?##A,
+@EU861894-132/1
+CGGAACACCCATGGCGAAGGCAGCTCGCTGGCCTGCAACTGACCCTGAGGCGCAAAAGCCGGGGGATCAAACAGGTTTATATACCCTGGTAGTCCCGAGC
++
+CCC4FFDFFFHH#AIJ<EJJJGGJJJJJIHEGE@GIIIBJ#JJ##GJ8IG#IJFHICFH#JJC?8I#HD;GD#ED#<B;##>DEDBC#H;#E8A(?>##D
+@EU861894-130/1
+TCGTCACTAACCGAAGAACCTTACCAGGGTTTGACATCCCCGGCCAGCACGAGAGATCGGGTGGTGCCCGCAAGGGGACCGGGTGACAGGTGTTGACTAG
++
+:B@FDFDFH<HAFIGFJJJ7JHJHJGFJHJJHHD?JBHIII1J#IAEI#7IIGFAEJGBCEAI#E#IHECC>CA@:CDFDEDAA#DDF:-<DA#D##CDD
+@EU861894-128/1
+CGGCTAACTACGTGCCAGCAGCCGCGGTTATACGTAGGGGGCCAGCGTTGTCCGTAGTTACCGGGAGTTAAGAGCGCGTAGCCAGGCCATTAAGTGGGGG
++
+@C?BFFFD#9HAHHHIJDJIJ<IJJJAB#GHFJJHGG*JG=IHI9#FJBDGBE8#B)@E7G#DFH##6#HD<D8IC;ACCD#D#;DA?#DACDE<CE>D,
+@EU861894-126/1
+GGAATGGTCGGTGTAGCGGTGAAATGCGTAGATATCGATCGGAACACCCATGGCTAAGGTAGCTTGCTGGCGTGCAACAGACGCTCAGGCGCAAAAGCCG
++
+BC?#FF=FBFHHHFGIGEIJ@IJJIFI@HJI9IJ;GHBGIIJG<CJIH8AEBI@#DGHF#G;)A#>HDC<>#9E7@CG#D(BHED#=;CD:;#>C#C?B,
+@EU861894-124/1
+CTGTGCCGCAGCTAACCCAAAAAGTATTCCGCCTGGGGAGTACGGCCGCAAGGTTGAAACTCAACGGAATTGACGGGGAACCGCACAAGCGGTGGACCAT
++
+CCCFFFFD,GB<CEI@#<I##JJG<HG<IBGFJ#DHI3IIHJII@GJAGGBGD#JEGBG'EJDB!HIEHE@AD#.DC###C<CF?C@DDFF?D@CD#BB,
+@EU861894-122/1
+TACACCGACCATTCCACTCGCCTCTCCCGCACTCAAGGCCCCCAGTATGATAGGCAGTTCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAAGTGTC
++
+B@@D=F2FHDFHFJJJIJ#JI@JHJGJJJEIHI0JIJJHG3JAGHJ@BIA!FDIJA#H2J#HDFH?@G)H;;EBD=@FC3#CCB>BD@8D>:BDD<###,
+@EU861894-120/1
+CTGTCACCCGGTCCGCTTGCGGGGAACACCCGATCTCTCGCGCTGTCCGGGGATGTCAAACCCTGGTAATGTTCTTCCGTTAGTGACGAATTAATGCAGA
++
+@C@FAFFFAHHHFA#JJ#J:IJ?H:HJJH9IJ#9I?IG:J#*IEBIC;I?AEJDIJJHHCICCHDJJ;C##C2EF?B#C>?DD6(B+#(D@E:DCDCC#,
+@EU861894-118/1
+CGAGCGCAACCCTGACCCAATGTTGGCAGCGTAAAGTCGGCAACTCTTTGGGAACTGCCCGTGGAAGCGGGAGGAAGGTGGGGATGACGTCAAATAGGCA
++
+@?;FF?EFFHHHH#HEAIGJ@J1JIJIJJGDGIF7GGH=J#JGIJJ#:GFJID#FJJ#JID4?#DBDJHEFID;EA5EDD;C?CCDCBDA>:CB>#DD#C
+@EU861894-116/1
+GAGGCAGTTCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAAGGGCCCGCCTACGCCCTATTTACGCCCATTAAATTCGGACAACGCTGGCCGCCTA
++
+C@@FFFD2GHHFFICFAGJJI3JJIIII8II#<*JJH#JJIID(EGEGJ;GE@D7IJHF#B;#C?H@GH;FCI#D<C#EA(>D@DE6CCDCDDC<#ACD4
+@EU861894-114/1
+TTACTAGGGATTCCGGCTTCATGCAGGCGAGTTGCAGCCTGCAATCCGAACTGAGACCATATTTGTGCGATTGGCTCCTCTACGCGGCATAGCGTCGCGC
++
+@C:FFDD#BHHAAIJ(GIFGIJGHBHDI9ICHGE7JHEJI4IJBBIDIAHJIBCBFIF#J=BDFFDD8HCHJ8E:@@C#D##BA@DD5=E#A5B#D#B#F
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pear_assembled_results1.fastq	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,56 @@
+@EU861894-140/1
+CAACCATTCCACTCGCCTCTGCCGCACTCGAGGCCCCCAGTATGAGAGGCAGTTACGAGGTTGAGCCTCAGAATTTCTCCCCCCACTTAAGGGCCCGCCTACTCGCTCTTTACGACCAGTAAAATCGGTCAAAGCTCACCCCCTACGTATTACCGC
++
+@#CAAF4FHFG;DIJJI?JG#GIIIDECD#4H*1IBIIIH#?IDID*DFIIF?C#HjZmKE^ed@LEiedWfekA=7#iFmcD=W%Qc\g7ekeDjf_bU)G#EE=?'>:EHDG#CIHHJIH!#J:BI#JJH#I<E##IHJJG)HFHHF#FDFCB;
+@EU861894-138/1
+TCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAACGGCCCGACTCCGCGCTCTTTACGCCCAGTAACTTCGGACAACGCCGGCCCCCTACGTATTACCGCGGCTGCTGGCACGTAGTTAGCCGGGGCTATTTAGAGGGGTACAGTCA
++
+@CCFFFBFFDHFHJJJ+<HJJGIGGBJJIDJ?JIJII#I5D*GB#EJ#FEldIEdCGjEYKIFhiAfA`;baUmhXQf;gYE#jEfHWhGCQXDKhG`jAIGJHH#IHJGGJBHGDGFHBA#IGGGJFG#@#!JJ!+B+#<HBF++F==@
+@EU861894-136/1
+GCTGATCGTCCTCTCAGACCACCTACCCGTCAAAGCCTTGGTGGGACGTTACCCCCCCAACTAGCTGATAGGCCGCGGGCCCATCCCGAAGTGCCCGAAGGCTTTGGATAAGCTCTGAGGGCAGAGCACATCCACATGCGGTA
++
+@@C?F;FF#HGHGJC<JJJC<#CJG61GHIE#?J:JJHJJHFJA_%T\^D^l8gHdVhLlD`ol(CFIgfl)meD;m<ngaIiif^dlgfj@gfLnmeicI>EGGGGCGJJH@*GJ;J#<J@?GEGGGHA#AHHH:FFDF@B1
+@EU861894-134/1
+TTGTGCAAGAGTCCTAATTGCTGCCCCCCGTAGGAGTCTGGGCCGTATCTCAGTCCCAATCCGGCTGATCGTCCCCTCAGACCAGCTACCCGTCGAAGCCTTGATGTGACAGTACCCTCCCAACTAGCTGAAAGGCCGCGGGCCCATCCCGAAGTGCCCGAAGGC
++
+=CCFDFFFA?#BHIJJJCIJ<JIJIGHJ8JHCI9JJIDJIGIJEFI9@JIJ#FFI5J>JJJD;EEiaAd4dFIG%ad^fciCdc6=ej_d^dlClbDAAT#7J#C>#J#F!#JDIBI(GGIIJCJHEH:JJ#IHJF=HJJHG+HIJI<FJICHCHADDF;+?+@@
+@EU861894-132/1
+CGGAACACCCATGGCGAAGGCAGCTCGCTGGCCTGCAACTGACCCTGAGGCGCAAAAGCCGGGGGAGCAAACAGGATTAGGTACCCTGGTAGTCCCGGCCGTAAACGATGAATACTAGGTGCTGGGGGTAGCGACCCCCCCTGTGCCGC
++
+CCC4FFDFFFHH#AIJ<EJJJGGJJJJJIHEGE@GIIIBJ#JJ##GJ8IMFiLHRjhZHDmmEaUlCjFIbg?eFC<dZD#f]gDkgElUDEY`*fd@?iJJFIDH8EFGBGJ?JHEB#CJ#+HHIE+FF!3FJGJ;H2HH)F=FD@1@
+@EU861894-130/1
+TCGTCACTAACCGAAGAACCTTACCAGGGTTTGACATCCCCGGCCAGCACGAGAGATCGGGTGGTGCCCGCAAGGGGACCGGGTGACAGGTGTTGCATAGTTGTCGTCAGCTACTATCGTGAGATGTTGGGTTAAGTCATGCAACAAGCGCAACCCTCACCCAATGTTG
++
+:B@FDFDFH<HAFIGFJJJ7JHJHJGFJHJJHHD?JBHIII1J#IAEI#7IIGFAEJGBCEAI#E#IHEec@eI]^egh^iiAAHh_j]P<aeEd3=CKfJ;HGEDABJG@J##H#)FHBJDJG?D#I#JI?IAC9FI##H::JJ+;IHFGE:GGGD+#HADFDFD:@1
+@EU861894-128/1
+CGGCTAACTACGTGCCAGCAGCCGCGGTTATACGTAGGGGGCCAGCGTTGTCCGAAGTTACTGGGCGTTAAGAGCGCGTAGGCGGGCCCTTAAGTGGGGGGTGAGATTCTGAGGCTCAACCTCGGAACTGCCTCCCATACTGGGGGCCTTG
++
+@C?BFFFD#9HAHHHIJDJIJ<IJJJAB#GHFJJHGG*JG=IHI9#FJBDGddW;e8WGUl&^ejBD>)kD_hQICPaCfaC`ERFfg@mAfjmeddfFPIIGI#HIG#HCI??JDIHIJBAG?F3?JJAJ*JJ!FG:?DD1HFDADF@@=
+@EU861894-126/1
+GGAATGGTCGGTGTAGCGGTGAAATGCGTAGATATCGATCGGAACACCCATGGCGAAGGCAGCTCGCTGGCCTGCAACTGACGCTGAGGCGCAAAAGCCGGGGGACCAACCAGGATTAGATACTCTGGTAGTCCCGGCCGTAAACGATGAATA
++
+BC?#FF=FBFHHHFGIGEIJ@IJJIFI@HJI9IJ;GHBGIIJG<CJIH8AEBIFDgeJHCG^+^A[HF`_>E;g7adIEkLfk_gFbbC][aDfeDlYkQHFGG'(HG?#@F2JGIIEDCJG>!J3EFGEJ3EGJGJHAJHDHF?:CADF8C@
+@EU861894-124/1
+CTGTGCCGCAGCTAACCCAAAAAGTATTCCGCCTGGGGAGTACGGCCGCAAGGTTGAAACTCAAAGGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATCTGCATTAATTCGTTACTAACCTAAGAACCTTAGCAGGGTTTGACATCCCCG
++
+CCCFFFFD,GB<CEI@#<I##JJG<HG<IBGFJ#DHI3IIHJII@GJAGGBGJ%JcfViIEiDBDJKhmebAaCPhEJHDc>dlclfkD`kAiiCD5hZS##HGBGJFFI6HFG#IJF:JHI#FIEJFG2#HI#HJ:<BHAHGHDDD#F@C=
+@EU861894-122/1
+TACACCGACCATTCCACTCGCCTCTCCCGCACTCAAGGCCCCCAGTATGAGAGGCAGTTCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAAGGGTCCGCGTACGCGCTCTTTACGCCCAGTATCTTCGGAC
++
+B@@D=F2FHDFHFJJJIJ#JI@JHJGJJJEIHI0JiJckeUJd`^`bcKADFipfc%mRn@kFekMejNn^ZmgkdMhj\=kkcffDY]lfcJll>GL#NJJG#E@#<:J2JJCAIJEFHHJA2CH!422FF@@+
+@EU861894-120/1
+CTGTCACCCGGTCCGCTTGCGGGGAACACCCGATCTCTCGCGCTGTCCGGGGATGTCAAACCCTGGTAATTTTCTTCGGTTAGTGACGAATTAATGCACATGCTCCTCCCCTTGTGCGGGCCGCCGTCAATTCCTTTGAGTTTCAACCTTGCG
++
+@C@FAFFFAHHHFA#JJ#J:IJ?H:HJJH9IJ#9I?IG:J#*IEBIC;I?AEJJIJljHCkbbf[L^[C##jUYide=i_dheYBcSI(id^BdCll`CE7H=JGA!9G*#>E#HGC@4IF:#DEJCI#GJBAJJEJJ<BHHHHHBDDFFB1B
+@EU861894-118/1
+CGAGCGCAACCCTGACCCAATGTTGGCAGCGTAAAGTCGGCAACTCTTTGGGAACTGCCCGTGTAAGCGGGAGGAAGGTGGGGATGACGTCAAATCGGCATGGCCCTTACGCCCTGGGCTTCACACATGCTACACTGG
++
+@?;FF?EFFHHHH#HEAIGJ@J1JIJIJJGDGIF7GGH\l%JIkcl%XdiglfFFclDmkD6^Hc^Dnidikk_ljIEFhRieeCkiBdeEXkDfIml6lGJIJEHJHJCJJGJI2D3J#I#G+AHFHBA?DD++CC@
+@EU861894-116/1
+GAGGCAGTTCCGAGGTTGAGCCTCAGAATTTCACCCCCCACTTAAGGGCCCGCCTACGCCCTTTTTACGCCCAGTAACTTCGGACAACGCTGGCCCCCTACGTATTATCGCGGCTGCTGGCACGTAGTTATCCAGGGCTTTTTCGTGTGGTACGGTCATT
++
+C@@FFFD2GHHFFICFAGJJI3JJIIII8II#<*JJH#JJIID(EGEGJ;GE@D7IJHF#eZ#b?fbiJXidlCD<fCid(^ibFeREiDfcfjcB^_lXE#EGIDB!EBJFCE?JH=FEGIG#FGGIIH#JC#@JJJ+J#F@EJ#G:HHHFB#DFFCBB
+@EU861894-114/1
+TTACTAGGGATTCCGGCTTCATGCAGGCGAGTTGCAGCCTGCAATCCGAACTGAGACCCTATTTGTGCGATTGGCTCCGCCTCGCGGCATCGCGTCGCTCTGATTAACGGCCATTGTAGCATGTATCAGGCCCAGGGCGTAAGGG
++
+@C:FFDD#BHHAAIJ(GIFGIJGHBHDI9ICHGE7JHEJI4IJBBOFhcjQleCfZni#LHBehdFjRmekmTP]_db7lI;je@D^]bnGg\bIDCcHa?D#JG:)IFJ#@;IB>GB4J;AJH#EF)#??BHHHHH;F=F=CC8
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/reverse.fastq	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,56 @@
+@EU861894-140/2
+GCGGTAATACGTAGGGGGTGAGCTTTGACCGATTTTACTGGTCGTAAAGAGCGAGTAGGCGGGCCCTTAAGTGGGGGGAATAATTCTGACGCTCAACCTC
++
+;BCFDF#FHHFH)GJJHI##E<I#HJJ#IB:J#!HIJHHIC#GDHE:>'?=EE#G)JC<HHDFIC5E=H.#4;DDH#F###ADCB:DAC##@BD5##D9D
+@EU861894-138/2
+TGACTGTACCCCTCTAAATAGCCCCGGCTAACTACGTGCCAGCAGCCGCGCTAATATGTAGGGCGCCGGCGCTGTCCGAAGTTAGTGGCCGTATAGCTCG
++
+@==F++FBH<#+B+!JJ!#@#GFJGGGI#ABHFGDGHBJGGJHI#HHJGI#G@GHI#5F#EJH#C#E#CCD!B7FBFC>B;A#C#E@###9CD#C<##@C
+@EU861894-136/2
+TACCGCATGTGGATGTGCTCTGCCCTCAGAGCTTATCCAAAGCCTTCGGGCGCTTCGGGATGGGCACGCGGTCTACCAGCTATTTGGGTGGGTAACGTCC
++
+1B@FDFF:HHHA#AHGGGEG?@J<#J;JG*@HJJGCGGGGE>IFFDJJJCC#GCBHAFEJGIDEG#J;DBG#EDC)#C(GFD#C#@2?#D8CCD53D#AA
+@EU861894-134/2
+GCCTTCGGGCACTTCGGGATGGGCCCGCGGCCTTTCAGCTAGTTGGGAGGGTACTGTCACATCAAGACTTCGACGGGTAGCTGCTCTGAGGGTTCGATCA
++
+@@+?+;FDDAHCHCIJF<IJIH+GHJJH=FJHI#JJ:HEHJCJIIGG(IBIDJ#!F#J#>C#J7#I#?DDICIH;DBGD=6DD#DDDDDD####>#EA@@
+@EU861894-132/2
+GCGGCACAGGGGGGGTCGCTACCCCCAGCACCTAGTATTCATCGTTTACGGCCGGGAATACCAGCGTACCTAGTCCTGTTTGCTCCCCCGTCTTTTGCGC
++
+@1@DF=F)HH2H;JGJF3!FF+EIHH+#JC#BEHJ?JGBGFE8HDIFJJF?@GH#@B#B;ECEJ#C:I#D@C#C#A=D</#CCD>C#DDB!5FB+##AD'
+@EU861894-130/2
+CAACATTGGGTGAGGGTTGCGCTTGTTGCATGACTTAACCCAACATCTCACGATAGTAGCTGACGACAACTGTGCAACCCCTGTCGACCGGTCCCCTTGC
++
+1@:DFDFDAH#+DGGG:EGFHI;+JJ::H##IF9CAI?IJ#I#D?GJDJBHF)#H##J@GJBADEGH;JC(!=3ACE>#DDE<EF##FE;CDCE>)C#AC
+@EU861894-128/2
+CAAGGCCCCCAGTATGGGAGGCAGTTCCGAGGTTGAGCCTCAGAATCTCACCCCCCACTTCAGGGCCCGCCTTCGGTCTCATAACGCCCAGTAACTTCGG
++
+=@@FDADFH1DD?:GF!JJ*JAJJ?3F?GABJIHIDJ??ICH#GIH#IGIIE#I@BJIGD#J@IF#8E=C>D#A6##:ED#D')BBC@;&F?#80D;@@C
+@EU861894-126/2
+TATTCATCGTTTACGGCCGGGACTACCAGAGTATCTAATCCTGGTTGGTCCCCCGGCTTTTGCGGCTCAGCGTCAGTTTCAGACCACCGAGCCGCCTTCG
++
+@C8FDAC:?FHDHJAHJGJGE3JEGFE3J!>GJCDEIIGJ2F@#?GH('GGFHFJ;JBCIBGB:#HFFD;DEEHE#BB#C#E#D>#!>A>#D#C##?DD'
+@EU861894-124/2
+CGGGGATGTCAAACCCTGCTAAGGTTCTTAGGTTAGTAACGAATTAATGCAGATGCGTCACCGTTTGTGCGGGCCCCCGTTAATTCCTCAGCGTTTCTAC
++
+=C@F#DDDHGHAHB<:JH#IH#2GFJEIF#IHJ:FJI#GFH6IFFJGBGH##H9G5##JF#F;!HGJEGB#ADHH#ECA>(CAFD##D##@6CC5@?##'
+@EU861894-122/2
+GTCCGAAGATACTGGGCGTAAAGAGCGCGTACGCGGCCCCTTAAGTGGAGGGTGAAATTCTGAGGCTCAACCTCGGAACTGCCACACATACTGTGGGCTT
++
++@@FF224!HC2AJHHFEJIACJJ2J:<#@E#GJJC!JG#II)JIIF:#EIBII;JHC.HHFI@DGFDF/D@#D>EAF#C=HF#D##BC77:D#C?D:#A
+@EU861894-120/2
+CGCAAGGTTGAAACTCAAAGGAATTGACGGCGGCCCGCACAAGGGGAGGAGCATGTGCGTTAATCCGTCACTAACCGAAGAAAGGTACCAGGGGATGGGA
++
+B1BFFDDBHHHHHB<JJEJJABJG#ICJED#:FI4@CGH#E>#*G9!AGJ=H7:C>JI#A):EF#GIB;DBEFBG=DFD5DH#!#A5#8?@@C##CC##'
+@EU861894-118/2
+CCAGTGTAGCATGTGTGAAGCCCAGGGCGTAAGGGCCATGCCGATTTGACCTCCTCCCCATCTTCCTCCCGTTTACATGGGCACTTCCCAAAGAGTCGCC
++
+@CC++DD?ABHFHA+G#I#J3D2IJGJJCJHJHEJIJGJ4IJII#I?(EA!GH#CGG8E##5JHEHCD@BE#=@H@##CDBC:&DCD>D>?#C:C###C@
+@EU861894-116/2
+AATGACCGTACCACACGAAAAAGCCCTGGATAACTACGTGCCAGCAGCCGCGATAATACGTAGGGGGCCAACGTTGTCCTAAGTGTCTGGGCGTATAAAG
++
+BBCFFD#BFHHH:G#JE@F#J+JJJ@#CJ#HIIGGF#GIGEF=HJ?ECFJBE!BDIGE#EEI=>BHHC@D#G#=A#CFA#DECD##CDBD>#CC?#@#@D
+@EU861894-114/2
+CCCTTACGCCCTGGGCCTGATACATGCTACAATGGCCGTTAATCAGAGCTACGCGATGCTACGAGGCGGAGCCAATCGCACAAGTAGGGTCACAGTTCGG
++
+8CC=F=F;HHHHHB??#)FE#HJA;J4BG>BI;@#JFI):GJ#D?<HBA#GAHGGJFI;##EI;II7@E@D,=DDCF;G#?CB#,##DF5E#DD(CC@#'
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml	Wed Feb 11 07:52:43 2015 -0500
@@ -0,0 +1,6 @@
+<?xml version="1.0"?>
+<tool_dependency>
+    <package name="pear" version="0.9.6">
+        <repository changeset_revision="e8448893d2be" name="package_pear_0_9_6" owner="iuc" toolshed="https://testtoolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>