Mercurial > repos > lparsons > ea_utils

changeset 2:71eaf2e85ae7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author lparsons
date Tue, 20 Nov 2012 16:00:51 -0500
parents cf4b5125a835
children f0d19a935325
files .hgignore .hgtags README.txt fastq-join.xml sam-stats.xml test-data/test.sam test-data/test_read1.fastq test-data/test_read3.fastq test-data/testout.join.fastq test-data/testout.txt test-data/testout.un1.fastq test-data/testout.un2.fastq tool_dependencies.xml
diffstat 13 files changed, 212 insertions(+), 235 deletions(-) [+]
line wrap: on
line diff
--- a/.hgignore	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,3 +0,0 @@
-^package$
-syntax: glob
-*.pyc
--- a/.hgtags	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,1 +0,0 @@
-d7b1bb0ab083311487138fd8f0f7561091f7e74e 0.1
--- a/README.txt	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,19 +0,0 @@
-== ea-utils Galaxy Wrapper ==
-
-This is a Galaxy wrapper for some ea-utils tools, fastq-join and sam-stats.
-
-** Installation **
-
-Installation from a tool shed provides the necessary tool dependencies.
-
-Otherwise, make sure fastq-join and sam-stats are in the path.
-Move the test data files to your galaxy root test-data.
-Move the xml files to a subdirectory of your tools directory and add lines in tool_conf.xml to point to them.
-Restart the Galaxy server.
-
-** Attribution **
-
-The ea-utils package and associated documentation can be found at: http://code.google.com/p/ea-utils/
-
-The galaxy wrapper code was written by Lance Parsons (lparsons@princeton.edu), Lewis-Sigler Institute for Integrative Genomics, Princeton University.
-The code is housed on BitBucket at: https://bitbucket.org/lance_parsons/ea_utils_galaxy_wrapper
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/fastq-join.xml	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,65 @@
+<tool id="fastq_join" name="fastq-join" version="0.1.1">
+    <description> - Joins two paired-end reads on the overlapping ends</description>
+    <requirements>
+        <requirement type="package" version="1.1.2-484">ea-utils</requirement>
+    </requirements>
+    <command>
+        fastq-join 
+        -v '$splitChar'
+        -p $pctMaxDiff
+        -m $minOverlap
+        #if $stitchLengthReport:
+        -r $outputStitchLengthReport
+        #end if
+        $read1
+        $read2
+        -o $outputUnmatched1 -o $outputUnmatched2 -o $outputJoined 
+    </command>
+    <inputs>
+        <param format="fastq, fastqillumina, fastqsanger, fastqsolexa" name="read1" type="data" label="Read 1 Fastq" />
+        <param format="fastq, fastqillumina, fastqsanger, fastqsolexa" name="read2" type="data" label="Read 2 Fastq" />
+        <param name="splitChar" type="text" value=" " label="Split read ids on this character" help="Default is space ' ' for Illumina reads" />
+        <param name="pctMaxDiff" type="float" value="8" min="0" max="100" label="Maximum percentage difference between matching segments" />
+        <param name="minOverlap" type="integer" value="6" min="1" label="Minimum length of matching segements" />
+        <param name="stitchLengthReport" type="boolean" value="False" label="Output verbose stitch length report" />
+    </inputs>
+
+    <outputs>
+        <data format="input" format_source="read1" name="outputJoined" label="${tool.name} on ${on_string} (joined)"/>
+        <data format="input" format_source="read1" name="outputUnmatched1" label="${tool.name} on ${on_string} (unmatched1)"/>
+        <data format="input" format_source="read2" name="outputUnmatched2" label="${tool.name} on ${on_string} (unmatched2)"/>
+        <data format="tabular" name="outputStitchLengthReport" label="${tool.name} on ${on_string} (stitch length report)">
+            <filter>stitchLengthReport</filter>
+        </data>
+    </outputs>
+
+    <stdio>
+        <exit_code range="1:" level="fatal" description="Unknown error occurred" />
+    </stdio>
+
+    <tests>
+        <test>
+            <param name="read1" value="test_read1.fastq" />
+            <param name="read2" value="test_read3.fastq" />
+            <output name="outputJoined" file="testout.join.fastq" />
+            <output name="outputUnmatched1" file="testout.un1.fastq" />
+            <output name="outputUnmatched2" file="testout.un2.fastq" />
+        </test>
+    </tests>
+
+    <help>
+Overview
+--------
+fastq-join joins two paired-end reads on the overlapping ends.
+
+Split read ids character: Verifies that the 2 files probe id's match up to char C. Use ' ' for Illumina reads.
+
+Maximum difference is the maximum allowed percentage of bases that differ in the matching region.
+
+Minimum overlap is the minimum number of bases that must overlap (with no more than the maximum difference) for reads to be joined.
+
+Verbose stitch length report is a report for each joined paired of reads showing how large the overlapping section was.
+
+This tool uses sqr(distance)/len for anchored alignment quality algorithm. It's a good measure of anchored alignment quality, akin to squared-deviation for means.  This tool uses the fastq-join program that is part of the ea-utils suite. See http://code.google.com/p/ea-utils/wiki/FastqJoin for details.
+    </help>
+</tool>
--- a/sam-stats.xml	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,83 +0,0 @@
-<tool id="sam_stats" name="sam-stats" version="0.1">
-    <description> - Compute statistics from SAM or BAM files</description>
-    <requirements>
-        <requirement type="package" version="1.1.2-484">ea-utils</requirement>
-    </requirements>
-    <command>
-        sam-stats
-        $trackMultAlign
-        $reportAllChr
-        #if $rnaSeqStats:
-        -R $rnaSeqStatsFile
-        #end if
-        #if $input.extension == "bam":
-        -B
-        #end if
-        -S $histBinSize
-        $input
-        &gt; $samStats
-    </command>
-    <inputs>
-        <param format="sam, bam" name="input" type="data" label="SAM/BAM File" />
-        <param name="trackMultAlign" type="boolean" value="False" truevalue="-D" falsevalue="" label="Keep track of multiple alignments (slower!)" />
-        <param name="reportAllChr" type="boolean" value="False" truevalue="-A" falsevalue="" label="Report all chr sigs, even if there are more than 1000" />
-        <!-- <param name="numReadsSubsample" type="integer" value="1000000" min="1" max="1000000" label="Number of reads to sample for per-base statistics (max 1,000,000)" /> -->
-        <param name="histBinSize" type="integer" value="30" min="1" label="Number of bins per chromosome for reads by chromosome &quot;histogram&quot;" />
-        <param name="rnaSeqStats" type="boolean" value="False" label="Output RNA-Seq statistics (coverage, 3 prime bias, etc.)" />
-    </inputs>
-
-    <outputs>
-        <data format="tabular" name="samStats" label="${tool.name} on ${on_string}"/>
-        <data format="tabular" name="rnaSeqStatsFile" label="${tool.name} on ${on_string} (RNA-Seq Stats)"> <filter>rnaSeqStats</filter>
-        </data>
-    </outputs>
-
-    <stdio>
-        <exit_code range="1:" level="fatal" description="Unknown error occurred" />
-    </stdio>
-
-    <tests>
-        <test>
-            <param name="input" value="test.sam" />
-            <output name="samStats" file="testout.txt" />
-        </test>
-    </tests>
-
-    <help>
-Overview
---------
-sam-stats computes varius statics on SAM/BAM alignment files.
-
-Complete Stats::
-
-  &lt;STATS&gt;           : mean, max, stdev, median, Q1 (25 percentile), Q3
-  reads             : # of entries in the sam file, might not be # reads
-  phred             : phred scale used
-  bsize             : # reads used for qual stats
-  mapped reads      : number of aligned reads (unique probe id sequences)
-  mapped bases      : total of the lengths of the aligned reads
-  forward           : number of forward-aligned reads
-  reverse           : number of reverse-aligned reads
-  snp rate          : mismatched bases / total bases
-  ins rate          : insert bases / total bases
-  del rate          : deleted bases / total bases
-  pct mismatch      : percent of reads that have mismatches
-  len &lt;STATS&gt;       : read length stats, ignored if fixed-length
-  mapq &lt;STATS&gt;      : stats for mapping qualities
-  insert &lt;STATS&gt;    : stats for insert sizes
-  &lt;CHR&gt;           : percentage of mapped bases per chr, followed by a signature
-
-Subsampled stats (1M reads max)::
-
-  base qual &lt;STATS&gt; : stats for base qualities
-  A,T,C,G       : base percentages
-
-Meaning of the per-chromosome signature:
-
-  A ascii-histogram of mapped reads by chromosome position.  It is only output if the original SAM/BAM has a header. The values are the log2 of the # of mapped reads at each position + ascii '0'.
-
-See http://code.google.com/p/ea-utils/wiki/SamStatsDetails for more information on each stat, how it's calculated and what it means.
-
-This tool uses the sam-stats program that is part of the ea-utils suite. See http://code.google.com/p/ea-utils/wiki/SamStats for details.
-    </help>
-</tool>
--- a/test-data/test.sam	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,100 +0,0 @@
-GEN-SEQ-ANA_0001:1:1:17434:14109#0/1	0	PCRPRIMER1	5	255	33M	*	0	0	ATACGGCGACCACCGAGATCTACACTCTTTCCC	?<??>>BB>>A5AAA;;9:B<A?:AAA<B<AAA	XA:i:0	MD:Z:33	NM:i:0
-GEN-SEQ-ANA_0001:1:1:4183:5576#0/1	0	PCRPRIMER1+1	1	255	38M	*	0	0	AATGATACGGCGACCACCGAGATCTACACTCTTTCCCT	AA=A=7=>;<BB8;>5>>?:B>BAABBBBBBBBABBB;	XA:i:0	MD:Z:38	NM:i:0
-GEN-SEQ-ANA_0001:1:1:16837:9043#0/1	0	PCRPRIMER1+1	14	255	27M	*	0	0	CCACCGAGATCTACACTCTTTCCCTAC	?BB?ABB?BBB<7A5?7??A?<AA?A?	XA:i:0	MD:Z:27	NM:i:0
-GEN-SEQ-ANA_0001:1:1:6753:19078#0/1	0	PCRPRIMER1+1	58	255	25M	*	0	0	TAATGATACGGCGACCACCGAGATC	A:BABA:AA<?B<<?BBBBBB@BBB	XA:i:0	MD:Z:25	NM:i:0
-GEN-SEQ-ANA_0001:1:1:14641:6640#0/1	0	PCRPRIMER1+1	58	255	43M	*	0	0	TAATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACAC	AC@CB>:CC>CCCCCC?@CCBB<?@B@BAB@ABBBBBB5:?5?	XA:i:0	MD:Z:43	NM:i:0
-GEN-SEQ-ANA_0001:1:1:18498:19719#0/1	0	PCRPRIMER1+1	58	255	38M	*	0	0	TAATGATACGGCGACCACCGAGATCTACACTCTTTCCC	CCCCCCCCCBCCACCC@CCC@BCCCCCC@CCCCCCCCC	XA:i:0	MD:Z:38	NM:i:0
-GEN-SEQ-ANA_0001:1:1:5165:8271#0/1	0	PCRPRIMER1+2	1	255	37M	*	0	0	AATGATACGGCGACCACCGAGATCTACACTCTTTCCC	CC=CC>CC:@CCC@C>BAB<B>AAA88?>>ABBB:AA	XA:i:0	MD:Z:37	NM:i:0
-GEN-SEQ-ANA_0001:1:1:3890:20139#0/1	0	PCRPRIMER2+1	34	255	23M	*	0	0	TAATGATACGGCGACCACCGAGA	CCCCCC<CBB;;C@CC5<CABB7	XA:i:0	MD:Z:23	NM:i:0
-GEN-SEQ-ANA_0001:1:1:6236:11228#0/1	0	PCRPRIMER2+1	34	255	20M	*	0	0	TAATGATACGGCGACCACCG	CCCCCC8CCCCC=CC?AAAA	XA:i:0	MD:Z:20	NM:i:0
-GEN-SEQ-ANA_0001:1:1:13390:7511#0/1	0	PCRPRIMER2+1	34	255	48M	*	0	0	TAATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGC	CCACBCBBBCCCCCCC?CCC;C?CCC?CCCCCCCCCCCACCCC;=>8A	XA:i:0	MD:Z:48	NM:i:0
-GEN-SEQ-ANA_0001:1:1:16810:9063#0/1	0	PCRPRIMER2+1	34	255	45M	*	0	0	TAATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGA	CCCCCCCCBCCCBBCCCCCCCCCCCCCCBA=@CCCCBC5B@9@@@	XA:i:0	MD:Z:45	NM:i:0
-GEN-SEQ-ANA_0001:1:1:1127:14386#0/1	4	*	0	0	*	*	0	0	TGAATTCATGAATGGAAAAA	?ACC<@>69@6>@BB?@:==	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1210:17845#0/1	4	*	0	0	*	*	0	0	GAGATAATTCAAAATCTTGCCGAT	>===>>@@@=<89A<B:@:BAA@@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1260:17853#0/1	4	*	0	0	*	*	0	0	GTACAAATTCAGGACCACGCGTAAACGTCACC	CCCCCB;ACBBB6@BB@BBB=A?A>6>AB?>B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1227:4069#0/1	4	*	0	0	*	*	0	0	TCGAGAATTCTTTATCAACCTCGTT	=CCCCCBC6@BCCC=9BA7A>>BB>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1309:16895#0/1	4	*	0	0	*	*	0	0	GTACAAATTCCCTCATTCTTGGATTTGTA	CA>CC@BBC8C@C>=C5<CCC@?CC;?>;	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1290:1939#0/1	4	*	0	0	*	*	0	0	TACGTAATTCCGTGTCGGAGCCCGT	??B:6>B;>AB=>B>B:BB>?A<6?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1311:12611#0/1	4	*	0	0	*	*	0	0	AAATTCGTTTATTTTTGTATTT	BBBBB@@@@@@@B:====:BBA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1312:13847#0/1	4	*	0	0	*	*	0	0	CGATAAATTCCTTGCTTCAACGATGCA	?BB>=8CC@??CCCC@B?>CC:CCC:8	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1327:14127#0/1	4	*	0	0	*	*	0	0	ATAAATTCAGAATAGTTTCACT	=A:??B>=;B@A;97C>B@=6B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1328:16612#0/1	4	*	0	0	*	*	0	0	GAGATAATTCTTCTCATCGG	@@@8@;C@C@CC=BC=>=AA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1339:13574#0/1	4	*	0	0	*	*	0	0	TTCTGTGATAATTTCAATATTGCC	CB8B>=B=ABBBB:=BABBB@>=6	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1343:12008#0/1	4	*	0	0	*	*	0	0	TTCGGCACTTTTCAGCACAATTTTCA	CA=ACC@CC@BB=?AC??CC?CC?<C	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1353:17775#0/1	4	*	0	0	*	*	0	0	CGATAACTTCCCAACAGTGTTGTAAAG	>5CBC?8?C=?ACBC>CCBB@=CAA8:	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1358:3384#0/1	4	*	0	0	*	*	0	0	AACAAACCATGTTTTGCCCTCAAACAG	B8:BBBBB:BBB=BBBA@<BBB>BB:9	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1361:15054#0/1	4	*	0	0	*	*	0	0	ATGCTAATTCAGTAGCTATTGTTCTCCAAAT	BCC:?CA?C?=CC?ABC@<?B@C@CC@8C>C	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1362:13972#0/1	4	*	0	0	*	*	0	0	GTAATTCTTCGCGAGTCTTTG	@>C?CA?:@??@@@@CC?C?B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1363:2154#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCTAGAAATGAAGAATATTTTTATTTATTTTCAG	CB@C?B;CC=CAAC@@>:C?@:???CCBB>8>>>>?C>6>A?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1373:2799#0/1	4	*	0	0	*	*	0	0	GAGATAATTCAAAACTACTACTGTAGAGGGGTACTG	?BABACCCC@B=<=@@?8B><@B@>BAAB9<@<>@B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1374:3530#0/1	4	*	0	0	*	*	0	0	GAGATAATTCTTCCTCCTTAA	AAAC8A=CA>CC@CC:B7ACC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1379:7364#0/1	4	*	0	0	*	*	0	0	CGCCAACCAAACCGATTGGCC	=???A;CCC?CCABCB8>CCC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1382:18005#0/1	4	*	0	0	*	*	0	0	GTAGCGAGTGACTCCAGTATGGAT	?CBB>=BBB:@B;>><<><@>@?7	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1385:14548#0/1	4	*	0	0	*	*	0	0	TTCCTATTCAACTCCTCAATT	8=9=AA:A@=?AA@@@:@?AA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1387:6982#0/1	4	*	0	0	*	*	0	0	ATTCGCGCCGATTGGGATTTTTG	6B?>:A=:A?7:<A<BBB?;?6<	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1390:2086#0/1	4	*	0	0	*	*	0	0	AATAAAACTTAAAATCTCAG	BBA?B:::A6B>A>B=B=BB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1392:4030#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCGAGCCCATTGATGA	?@6C@C@<A??CCC=C?C?C:A?A	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1397:15413#0/1	4	*	0	0	*	*	0	0	TCTTCGAAAAAACAAGCATCTACGTCTTG	7?BA?<B?B@<:::;;<>>BB<<<;95BA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1403:14145#0/1	4	*	0	0	*	*	0	0	ATGCTAATTCTCGTGGAACTCG	CA>C=CC@CA?8???7C7@=??	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1404:13576#0/1	4	*	0	0	*	*	0	0	CGATAAATTCTCTGAAAATGCATA	>BB=:BBBB@B@=B@BB=BBAB;B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1408:13430#0/1	4	*	0	0	*	*	0	0	GAACAATAAAACCTCCTGTC	@8B5@AA?@A==AAA=AA6@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1410:6288#0/1	4	*	0	0	*	*	0	0	TCCGAAAACCAAATTTTACTG	9?>B5:CBC:@C@CCB?CCCB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1420:17657#0/1	4	*	0	0	*	*	0	0	TTCAATGAGAGAATCCGCGACGAAAATGG	CC?CCC?BB?C:CAAB5AB>BB>BA8<B8	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1421:18483#0/1	4	*	0	0	*	*	0	0	GTACAAATTCAGAAAATATTTTATT	BC?CC@CBCA@BC@AA:6B@=BAB8	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1422:15203#0/1	4	*	0	0	*	*	0	0	GTACAAATTCGTGTCGAAATTGAAAAAAAAAG	CCC?CCCBCACC@@@CCCACC@=C;8?;98>>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1423:18328#0/1	4	*	0	0	*	*	0	0	ATAAATTCATCGACAGGTACCTTCTAT	<@@<>?B>>AB=>BBB>B=BBBA7>>6	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1426:4162#0/1	4	*	0	0	*	*	0	0	GCTAATTCCTATGGTTATAACTCTAAAACGGGTC	7==8<>@9B@?=5ACACA>>BCC?>C=?B>9@8;	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1426:18217#0/1	4	*	0	0	*	*	0	0	TCCTTCAAAAATCATCACTGAAA	>7:>@>@AA=@A:A@6@7A;>55	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1428:17122#0/1	4	*	0	0	*	*	0	0	GTACAAATTCATCTAATGGTCTAACTTTATA	C?CC8CCCBA@@=C>AAABC8BC:CCBACCC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1434:13892#0/1	4	*	0	0	*	*	0	0	GTACAAATTCGCCAAGCGTTGGATTG	A???88@:;=>=?@@>BBB6A>AA5>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1441:16403#0/1	4	*	0	0	*	*	0	0	TGAATTCTGGCACTTCGGCGATT	B>B@?B<BB?A@>BB@@>BAAB:	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1445:8376#0/1	4	*	0	0	*	*	0	0	TGAATATGAAATTCGTAAATAGCGTGTT	@@CCC@>CBCC?CCC<CCBB?B?C@?@<	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1449:14365#0/1	4	*	0	0	*	*	0	0	TTCTCCAGGTACATCTGTCTGAAA	CC<CCBCCB@@BB>@BBBBBBBB>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1451:7415#0/1	4	*	0	0	*	*	0	0	CTCTTAATTCAAGTCCATATTCTCCTTGTAACCCCAAGTGATACAATATTTGATT	BBA@;CCAC>@@CC>CBCC?BCCBBC<B=CBC=?7=<BABCC;@<68AA;==:=?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1452:4973#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCTTGCTAAGCCGTGAATACAC	CA<C;>AAAABCCCBCCCBC@CACAC;CC<	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1455:18050#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCCAGACATTATAGAAAA	BCCCC@=B@>>?C>CCCC?:CC:C=>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1455:12928#0/1	4	*	0	0	*	*	0	0	TCGAGAATTCCACATCACTTCTAT	@?@@;CCCC@B:C:CBBB=6<B8>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1456:13493#0/1	4	*	0	0	*	*	0	0	TTCATATTGGGTTTCTACCGCGAC	AB>B=?AA>B@BB>=BAAB=BBB>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1456:4776#0/1	4	*	0	0	*	*	0	0	ATGCTAATTCACATTCATGATT	CCCA=A@C@<CCCCCCC?A@C=	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1464:19773#0/1	4	*	0	0	*	*	0	0	TAAATTCTACTTCAAGGCGTGAT	=@B>@@<BBBBB@BBBB@B=<<?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1469:15650#0/1	4	*	0	0	*	*	0	0	TACGTAATTCACTGCACTAATC	B=>C@C5C???BCCBCCC@?CC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1469:15421#0/1	4	*	0	0	*	*	0	0	CGATAAATTCTGCCGGCAGAACTACGCATTT	AC;>CACCB8?CCA;ACBCBBC?CCAC=CCC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1472:17572#0/1	4	*	0	0	*	*	0	0	TTCTACCCGGAACTCTTAAACAGAAGCGAGGC	CC>B@@:BB>>>B>B<BBB:B@B7@@=A<BBB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1473:13951#0/1	4	*	0	0	*	*	0	0	AAATTTTTAAAAAAAATGTTTG	BCCCCCCC>C=CBC@8>CCC=C	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1473:13352#0/1	4	*	0	0	*	*	0	0	TCTGTGGAATTAACAATCACAATAACTAAAG	C?@BC@B<C=CC;C=CCCBCB8CA7C?CCC?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1475:14575#0/1	4	*	0	0	*	*	0	0	TACGTAATTCGAAATTTGAG	?CB@C==AA?C=CCCBB:@B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1480:6616#0/1	4	*	0	0	*	*	0	0	GTACAAATTCGAGCGCAAAGG	C;CCCCCCC?C@C9CC?CCC:	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1481:5435#0/1	4	*	0	0	*	*	0	0	GTACAAATTATGAATAGAGAACTTTTA	B@B6BBB=B>@@=@=@@BBBB>B?>B=	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1481:5276#0/1	4	*	0	0	*	*	0	0	GTACAAATTCGAATGAAATCATA	CCCCBBCDB@C=CDACCC=C?C>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1482:10064#0/1	4	*	0	0	*	*	0	0	CGATAAATTCAATTTGAAAA	CC=CCC6C@;CCCBCAB?9B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1484:2325#0/1	4	*	0	0	*	*	0	0	AATTCCAACCACCATTGGAA	@A@=>@C??@C7CB@B6?:?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1487:6701#0/1	4	*	0	0	*	*	0	0	CTCATTCATCTAGATATTTCTGATTTATCGTTGTTC	@9<?C@?BCCBCB@6@@AAA<<<<AB:BB@?@B6>>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1487:3645#0/1	4	*	0	0	*	*	0	0	CTAATTCGTGTAATAATTGTGTT	>87?@@>><9C:@7<A<BB?;>B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1488:9645#0/1	4	*	0	0	*	*	0	0	ATGCTAATTCGAAGTGCGTATCTCTCGGT	ACCCCC8C@>>BBB=BB=A>BBBBB6B?>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1489:3252#0/1	4	*	0	0	*	*	0	0	TAAATTCTGAAAAATGTATTT	@<@@B@@=?C:C>>>B@A6A8	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1492:14525#0/1	4	*	0	0	*	*	0	0	AACTTCGAATATTAACAACA	@7:@>@@5C@>@@@C@@?B@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1493:11409#0/1	4	*	0	0	*	*	0	0	CGAAACGGGCGTTGAGGCGTTTT	<BAB?>>=B6BBBBAB@=BB>BB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1494:17779#0/1	4	*	0	0	*	*	0	0	TTTTCCTGAAAAATGTTCAACCTGATTATT	:;A;:?:A>A8A:A=?AAA@@:=55AA6@@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1495:12537#0/1	4	*	0	0	*	*	0	0	GTACAAATTCACAATAATTA	BCCCBCCBCA;CCCBB@>>8	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1495:14319#0/1	4	*	0	0	*	*	0	0	ATTCGCCAAATTTGGCGTTAG	<?A>CCCCACBBCB@CCCACA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1496:18187#0/1	4	*	0	0	*	*	0	0	AGGCGCGGCCAGCGGCTGGG	<:B<<6B6<>A?BBB:5A6A	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1496:10780#0/1	4	*	0	0	*	*	0	0	GAAATTCCCAAGAAACCCAAG	C:?7@??CACCBCCC=8AA:6	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1497:8211#0/1	4	*	0	0	*	*	0	0	CGATAAATTCCAAAGTTTGTGAATTT	BCCCCC?C@?<::A:?AA?AB?:B@<	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1498:13789#0/1	4	*	0	0	*	*	0	0	ATGCTAATTCTTTCGGAAACGAGACCATCTGA	>C8C>CC>B@==C@AC?C8@B?:6<AA?A?AB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1501:3386#0/1	4	*	0	0	*	*	0	0	GTACAAATTCAGTAGTTTTG	?<=8?76@>;<C>AA@BBA@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1502:11211#0/1	4	*	0	0	*	*	0	0	ATTCGTACTCGTGATTTTCGAACAATCCAGT	@AB=CCCCCCCBBA;CACCCB@A@C==C@CA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1506:3772#0/1	4	*	0	0	*	*	0	0	CCATAAAATTATTGTCAATTTTCAAT	??:@?:BC:C?CC>C<?=@>BCCC::	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1506:13304#0/1	4	*	0	0	*	*	0	0	AATTCGACATCGCAAGTGAG	AAA?>?CC>CBB?CBCBCA>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1506:17897#0/1	4	*	0	0	*	*	0	0	TGAATTCCGCCTTGGAGACAATGTTCT	B:6<BB8B?=B>>?:B=B=BB>=?BBA	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1507:13225#0/1	4	*	0	0	*	*	0	0	TGCTAATTCTAATTTTCATGTATACTAAACCCCCAC	@?C>B8@<9B?C=CB5B66A<<?6?A>AA8?A<<A6	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1508:6179#0/1	4	*	0	0	*	*	0	0	GTCAAAAGTAGAGAAAAAAAA	><B=9<@B?@@BAB@A@9=9=	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1508:8256#0/1	4	*	0	0	*	*	0	0	ATTCTCGCCAGCCTTTTTAAGTCGAGTTT	ABB<ABA>BAAA=?B???6BB=<BAB=B?	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1508:11605#0/1	4	*	0	0	*	*	0	0	TACGCAATTCTTCTCAAAAATTC	B==B9=:=B=CBCCCB:B5B5>>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1510:9250#0/1	4	*	0	0	*	*	0	0	AAAAACTGAAAAATTAAGATTAAAAGC	A=?A6AA@?AAA@@AA:A8A@:>8AA5	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1510:4740#0/1	4	*	0	0	*	*	0	0	CTCTTAATTCATTTTTTTTTCAGTTTTAAAGGTG	BC;BBCCCC@>CC;CCB>8>=BBB6BBBB=B?=B	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1511:8709#0/1	4	*	0	0	*	*	0	0	TTCACTTTCTATGGTTATTTT	CC@=ACCCCCB?CCBBCCC=>	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1513:14094#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCCAATTCTCTATCACGATT	BC8CCBCC@;CCC8CCCCCACCCAB@@@	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1515:16388#0/1	4	*	0	0	*	*	0	0	AGCTGAATTCGCCAAGCGTTGGAT	BC;CCC@C?;=CCCCCCC@C;C@C	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1516:14913#0/1	4	*	0	0	*	*	0	0	CTAATTCGTCCAGATTTCCA	C:BCCB?CCBBACCBBCCCC	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1517:17176#0/1	4	*	0	0	*	*	0	0	CGATAAATTCTAATTTTGCCA	>C<@CBA9B?BB?CBCABCCB	XM:i:0
-GEN-SEQ-ANA_0001:1:1:1520:6860#0/1	4	*	0	0	*	*	0	0	GTACAAATTCTTGAGAGACAGAATT	C??C;CCCC?CCCCCCC@BC@C?AC	XM:i:0
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_read1.fastq	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,40 @@
+@JLK5VL1:222:D1888ACXX:1:1101:1656:2143 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+?+:BDDDAD8BDD?@4+<?EECD?BBFIDFD<E:*?DI@?DDEDDCEECEIECDCCAD>5;6?;;@A;AABA>A@>AAAD;A>A>5=5>>>BBAA######
+@JLK5VL1:222:D1888ACXX:1:1101:1613:2167 1:N:0:
+GTGATAGAGATACTGAGCACAGAGCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@?BDFFFFHBHHHJJJGFEHHEIHEGHJJJJDI@DHIIHGGGGGEHIHGIIIIJJJGIB?AACBEDDFCDCCCCCEDC>@CD>CC@BA>CDEECD=BDDB>
+@JLK5VL1:222:D1888ACXX:1:1101:1927:2121 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@B@DFFFFHHHHHJJIHJIJJJJJJJJJJJJIJHJJIJIJIHHGIJIIJGIIJJJIIII?EEEDFFFFEEEEEDECDDDEDDDDCCBACDEDECC@BB9<A
+@JLK5VL1:222:D1888ACXX:1:1101:2000:2166 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACGTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+B?@DFFFFHHHHHGIIJIIJJJJJJIJJJJJIIHIJJIJJIIGGJJGHGHIGGIIJHIGAEEEDFFEEEEFEEDDDDDCC>CCCCCACDDEEDCA9>BBDC
+@JLK5VL1:222:D1888ACXX:1:1101:1763:2172 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAAGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@?@DDDBDHDFDHBFB3CCGFFHDHHIIIIIIG@<EFHA?DDFFDHHHHCHIIIIGGGH-;@5=CHHEADDEBEECCCCDCA@CCCB@CC@BC(;';BBBC
+@JLK5VL1:222:D1888ACXX:1:1101:1832:2198 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGCATATACCGCCAAT
++
+B=BDFFFFHGHHHJJIIJHIJIIJJJJJJJJIJIJJJIIJJHIIIJJJJIIIIJJIGJJ??CBBFFFFEDEEEEECDDCDDDCDCC@BDCDDEEC><>BD?
+@JLK5VL1:222:D1888ACXX:1:1101:2169:2161 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTCGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@?BFFFFFHHHHHJJJGIGIJIIJJJIGJJJJJJIIIIHGHIIGIJGIIIJJIJIIJJJEHDBDCDEEEEDDDDD@CCCDAA@CA@A:@CADD>C9<@B9C
+@JLK5VL1:222:D1888ACXX:1:1101:2236:2214 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@=@DFDFFHAHDHGFEE@FFHGGEDEGHGIJ@GFGHAGFIIGBDGHIIJIGJJJJJIIIACEBDEFDFDEEEEDEDDCCDCD;CDCB@CF@DD@C?B<@DC
+@JLK5VL1:222:D1888ACXX:1:1101:2090:2240 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@:?DDDDDF<FHHIIEEGBHGGDGGIIIGHIIIGGGHGEHFIHGIIIII<HHHIEGGIIHE=?BDECD@CDDCCCCCAAC@ACACCBCCDCDDDC######
+@JLK5VL1:222:D1888ACXX:1:1101:2952:2161 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@@@DFFFFHFHHGIJIJJJJJIJJJJIIIJJIIJJIJIEHIHIEHHIGJEHIIJIJDCGCE?BFFFDFEDDEEECACCCDDD@C>ABCCDCFFCA:9BBBC
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_read3.fastq	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,40 @@
+@JLK5VL1:222:D1888ACXX:1:1101:1656:2143 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+BB8=4ADDHHHHACG<FHIGGIIIIICAEEC;CDF0@DEG3?BBACCGGICFD@@;@7;AED<<=CC@>(;(;@C@C:<:>>@>>;B>B(44?>C@B9<<8
+@JLK5VL1:222:D1888ACXX:1:1101:1613:2167 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+@?=DDDFFDFHHHIGGHHIGGFEGIJJFAFEGGGH@AFFHGIIIGEFHBF>@FCGAE@DD=D'9@DCC@;C@C@CDDCD>CCD@>>:@@C@CDDACA<28?
+@JLK5VL1:222:D1888ACXX:1:1101:1927:2121 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+CCCFFFFFHHHHHJIIJJJJJJIJJJJIJJJJJJJIJJJJJJJJJJJJJJIHIJJIJEEHHF>BBDEFEDDDDDDDCDDDDCDDDEEEEDDDDDDDBDBDD
+@JLK5VL1:222:D1888ACXX:1:1101:2000:2166 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+B@@DDFFFGHHHDIFHIJIIIJJJJIJJJJJIJJIIJIJIGIIJGIJJJJDFGGGII:DHFF<BDECDEDD@CCCDDCDDDCEDDDDC@C@CCDCD?@8<<
+@JLK5VL1:222:D1888ACXX:1:1101:1763:2172 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+@;?DFDFFFDHHDIGACFHGGGGJJEH3@GDGCGHIEIIG?FDHCEGHEI9@C@C@@@9@=<',5;5;;;AC>3>C>CCCC(5(;>CDECCAC>@AB08AB
+@JLK5VL1:222:D1888ACXX:1:1101:1832:2198 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATGCTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+C@@FFFFFHHHHHJJJJJIJIJJJJIIJJIJIIFIGIJGHIDGIIIGIJJIHGHIIJGIIHF<<@DDDDEEDDCCCDDDCCCDDDEDDEDDCCDDCBB?AB
+@JLK5VL1:222:D1888ACXX:1:1101:2169:2161 3:N:0:
+TTTGCCCTATTTATAGGACCCTAGAAAGAAAAAACCTCAAATAAAATAAAAGTGATTGGCGGGATATACTTGTGAACAATGTGTATTTACTTACCACCACC
++
+#####################################################################################################
+@JLK5VL1:222:D1888ACXX:1:1101:2236:2214 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+@@@DDFFFHHHHFEEGIEHIJJGGFGIGIHCHJJJJJJIJIJGIGGIJIJGCDGEGG>EGHF=ABC@BDDCDDA@CCCCDCDE@CCD>CC@CDDDDBD?<<
+@JLK5VL1:222:D1888ACXX:1:1101:2090:2240 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGATAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+CC@FFFFFGFHDDHH>@FHHCHIIGIIEGGHIJIIBEEHIGGGGJIJJJJGBFGIGGCAHEA,9?B>B@@C>@@ACACCCC@C@@CDCDC@>CCCDD@(8?
+@JLK5VL1:222:D1888ACXX:1:1101:2952:2161 3:N:0:
+TGGATCCTATTAATAGGACCCTAGAAAGAAAAAACGCCAATCACAATAAAGTTGATTGGCGGTATATACTTGTGAACAATGTGTATTCACTAACAACGTCA
++
+CCCFFFFFHHHGHJJJJJIJJJJJIJJJHIJJIJJIIGIIIJEIIJIIJIG@FHIJJJIJHF<A?CBCC>CCDCCDDDDDDCDCCDDEACA@CDDDDDABB
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/testout.join.fastq	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,36 @@
+@JLK5VL1:222:D1888ACXX:1:1101:1656:2143
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+?+:BDDDAD8BDD?@4+<?EECD?BBFIDFD<E:*?DI@?DDEDDCEECEIECDCCADB@C>?;;@B>BABA@A@>AACDCA>A>5>@CC>BBDEA;7@;@@DFCIGGCCABB?3GED@0FDC;CEEACIIIIIGGIHF<GCAHHHHDDA4=8BB
+@JLK5VL1:222:D1888ACXX:1:1101:1613:2167
+GTGATAGAGATACTGAGCACAGAGCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@?BDFFFFHBHHHJJJGFEHHEIHEGHJJJJDI@DHIIHGGGGGEHIHGIIIIJJJGIBCADDCEDDFCDCCDCCEDCDDCDCCC@BCCDDEEDDDDDEBGCF@>FBHFEGIIIGHFFA@HGGGEFAFJJIGEFGGIHHGGIHHHFDFFDDD=?@
+@JLK5VL1:222:D1888ACXX:1:1101:1927:2121
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@B@DFFFFHHHHHJJIHJIJJJJJJJJJJJJIJHJJIJIJIHHGIJIIJGIIJJJIIIIDEEEDFFFFEEEEEDEDDDDEDDDDDDEFEDEDEFHHEEJIJJIHIJJJJJJJJJJJJJJIJJJJJJJIJJJJIJJJJJJIIJHHHHHFFFFFCCC
+@JLK5VL1:222:D1888ACXX:1:1101:2000:2166
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACGTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+B?@DFFFFHHHHHGIIJIIJJJJJJIJJJJJIIHIJJIJJIIGGJJGHGHIGGIIJHIGDEEEDFFEEEEFEEDDDDDDDCCCCDDEDDEEEDFFHDBIIGGGFDJJJJIGJIIGIJIJIIJJIJJJJJIJJJJIIIJIHFIDHHHGFFFDD@@B
+@JLK5VL1:222:D1888ACXX:1:1101:1763:2172
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAAGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@?@DDDBDHDFDHBFB3CCGFFHDHHIIIIIIG@<EFHA?DDFFDHHHHCHIIIIGGGHA@@CACHHECDDEBEECCCCDCA@CCCB@CC@BC<=@;BBBC@C@9IEHGECHDF?GIIEIHGCGDG@3HEJJGGGGHFCAGIDHHDFFFDFD?;@
+@JLK5VL1:222:D1888ACXX:1:1101:1832:2198
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGCATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+B=BDFFFFHGHHHJJIIJHIJIIJJJJJJJJIJIJJJIIJJHIIIJJJJIIIIJJIGJJCDDCCFFFFEEEEEEECDDDDDDDDEEDDDDDDEFHIIGJIIHGHIJJIGIIIGDIHGJIGIFIIJIJJIIJJJJIJIJJJJJHHHHHFFFFF@@C
+@JLK5VL1:222:D1888ACXX:1:1101:2236:2214
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@=@DFDFFHAHDHGFEE@FFHGGEDEGHGIJ@GFGHAGFIIGBDGHIIJIGJJJJJIIIDDEDDEFDFDEEEEDEDDCCDCDDDDDDBCFBDDFHGE>GGEGDCGJIJIGGIGJIJIJJJJJJHCHIGIGFGGJJIHEIGEEFHHHHFFFDD@@@
+@JLK5VL1:222:D1888ACXX:1:1101:2090:2240
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTATCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@:?DDDDDF<FHHIIEEGBHGGDGGIIIGHIIIGGGHGEHFIHGIIIII<HHHIEGGIIHECCBDEDDDCDDCCCCCCACAACACCBCCDCDDDEHACGGIGFBGJJJJIJGGGGIHEEBIIJIHGGEIIGIIHCHHF@>HHDDHFGFFFFF@CC
+@JLK5VL1:222:D1888ACXX:1:1101:2952:2161
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTTGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAATCAACTTTATTGTGATTGGCGTTTTTTCTTTCTAGGGTCCTATTAATAGGATCCA
++
+@@@DFFFFHFHHGIJIJJJJJIJJJJIIIJJIIJJIJIEHIHIEHHIGJEHIIJIJDDGDEDCFFFDFEDDEEEDDDDDDDDDCCACCCDCFFFHJIJJJIHF@GIJIIJIIEJIIIGIIJJIJJIHJJJIJJJJJIJJJJJHGHHHFFFFFCCC
--- a/test-data/testout.txt	Tue Nov 20 15:54:14 2012 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,29 +0,0 @@
-reads	100
-version	1.32
-mapped reads	11
-mapped bases	377
-phred	33
-forward	0
-reverse	11
-len max	48
-len mean	34.2727
-len stdev	9.4137
-mapq mean	255.0000
-mapq stdev	0.0000
-mapq Q1	255.00
-mapq median	255.00
-mapq Q3	255.00
-snp rate	0.000000
-pct mismatch	0.0000
-base qual mean	31.5650
-base qual stdev	3.2627
-%A	27.5862
-%C	34.4828
-%G	16.7109
-%T	21.2202
-%PCRPRIMER1	8.753316
-%PCRPRIMER1+1	45.358090
-%PCRPRIMER1+2	9.814324
-%PCRPRIMER2+1	36.074271
-num ref seqs	4
-num ref aligned	4
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/testout.un1.fastq	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,4 @@
+@JLK5VL1:222:D1888ACXX:1:1101:2169:2161 1:N:0:
+GTGATAGAGATACTGAGCACAGACCCTTATTAAGCCGAGGGTCACCTAGCCAACTGACGTCGTTAGTGAATACACATTGTTCACAAGTATATACCGCCAAT
++
+@?BFFFFFHHHHHJJJGIGIJIIJJJIGJJJJJJIIIIHGHIIGIJGIIIJJIJIIJJJEHDBDCDEEEEDDDDD@CCCDAA@CA@A:@CADD>C9<@B9C
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/testout.un2.fastq	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,4 @@
+@JLK5VL1:222:D1888ACXX:1:1101:2169:2161 3:N:0:
+TTTGCCCTATTTATAGGACCCTAGAAAGAAAAAACCTCAAATAAAATAAAAGTGATTGGCGGGATATACTTGTGAACAATGTGTATTTACTTACCACCACC
++
+#####################################################################################################
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml	Tue Nov 20 16:00:51 2012 -0500
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<tool_dependency>
+    <package name="ea-utils" version="1.1.2-484">
+        <install version="1.0">
+            <actions>
+                <action type="download_by_url">http://ea-utils.googlecode.com/files/ea-utils.1.1.2-484.tar.gz</action>
+                <!-- Unable to get varcall to compile properly we don't need it, so skip -->
+                <action type="shell_command">sed -i.bak 's/varcall$//' Makefile</action>
+                <action type="shell_command">PREFIX=$INSTALL_DIR make install</action>
+                <!-- For some reason, this is not installed by the Makefile, so move it ourselves -->
+                <action type="move_file">
+                    <source>sam-stats</source>
+                    <destination>$INSTALL_DIR/bin</destination>
+                </action>
+                <action type="set_environment">
+                    <environment_variable name="PATH" action="prepend_to">$INSTALL_DIR/bin</environment_variable>
+                </action>
+            </actions>
+        </install>
+        <readme>
+        </readme>
+    </package>
+</tool_dependency>