Repository 'kat_sect'
hg clone https://toolshed.g2.bx.psu.edu/repos/nml/kat_sect

Changeset 1:ab2cf85b9ae8 (2017-03-13)
Previous changeset 0:699eef7398c1 (2017-02-06) Next changeset 2:7f59ee641f70 (2017-03-13)
Commit message:
planemo upload commit f86fb54b98f0dbed260735675219e64882f08efe
modified:
kat_sect.xml
test-data/db.fasta
test-data/input.fastq
test-data/kat-sect-counts.cvg
test-data/kat-sect-stats.tsv
added:
macros.xml
removed:
test-data/input-counts.cvg
test-data/input-stats.tsv
test-data/matched.fastq
b
diff -r 699eef7398c1 -r ab2cf85b9ae8 kat_sect.xml
--- a/kat_sect.xml Mon Feb 06 12:36:09 2017 -0500
+++ b/kat_sect.xml Mon Mar 13 15:36:43 2017 -0400
[
@@ -1,14 +1,16 @@
-<tool id="kat_sect" name="KAT Sect" version="1.3.0">
-  <description> SEquence Coverage estimator Tool. Estimates the coverage of each sequence in a file using K-mers from another sequence file.</description>
-  <requirements>
-    <requirement type="package" version="2.3.1">kat</requirement>
-  </requirements>
-  <stdio>
-    <exit_code range="1:" />
-  </stdio>
+<?xml version='1.0' encoding='utf-8'?>
+<tool id="kat_@EXECUTABLE@" name="KAT @EXECUTABLE@" version="@VERSION@.1">
+  <description> Estimates the coverage of each sequence in a file using K-mers from another sequence file(s).</description>
+  <macros>
+    <token name="@EXECUTABLE@">sect</token>
+    <import>macros.xml</import>
+  </macros>
+  <expand macro="requirements" />
+  <expand macro="stdio" />
+  <expand macro="version_command" />
   <command><![CDATA[
 
-  ln "$db" db.fasta &&
+  cp "$db" db.fasta &&
   #if $single_or_paired.type == "single"
     ln "$input_se" "$single_or_paired.input_se.element_identifier".fastq && 
     kat sect 'db.fasta'
@@ -22,45 +24,46 @@
     kat sect 'db.fasta'
     "$single_or_paired.fastq_collection.forward.name"_1.fastq "$single_or_paired.fastq_collection.forward.name"_2.fastq
   #end if  
-  --threads \${GALAXY_SLOTS:-4}
-  -o reads
+
+  @THREADS@
+  
+  $keep_counts
+  
+  -o results
   
   ]]></command>
   <inputs>
-    <conditional name="single_or_paired">
-      <param name="type" type="select" label="Read type">
-        <option value="single">Single-end</option>
-        <option value="paired">Paired-end</option>
-        <option value="collection">Collection Paired-end</option>
-      </param>
-      <when value="single">
-        <param name="input_se" type="data" format="fastqsanger" label="Single end read file(s)"/>
-      </when>
-      <when value="paired">
-        <param name="forward_pe" type="data" format="fastqsanger" label="Forward paired-end read file"/>
-        <param name="reverse_pe" type="data" format="fastqsanger" label="Reverse paired-end read file"/>
-      </when>
-      <when value="collection">
-        <param name="fastq_collection" type="data_collection" label="Paired-end reads collection" optional="false" format="fastqsanger" collection_type="paired" />
-      </when>
-    </conditional>
+    <expand macro="macro_input" />
+
     <param name="db" type="data" format="fasta" label="Fasta file of targeted regions"/>
+    <expand macro="macro_kmers" />
+    <param name="keep_counts" type="boolean" falsevalue='--no_count_stats' truevalue='' checked='False' label="Fasta Count File" help="Produce a fasta style representation of the input sequence file containing K-mer coverage counts mapped across each sequence"/>
   </inputs>
   <outputs>
-    <data format="tabular" name="reads" label='Reads' from_work_dir="reads-stats.tsv" />
+    <data format="tabular" name="stats" label='Stats' from_work_dir="results-stats.tsv" />
+    <data format="txt" name="counts" label='Counts' from_work_dir="results-counts.cvg">
+      <filter>keep_counts</filter>
+    </data>
   </outputs>
   <tests>
     <test>
-      <param name="input" value="input.fastq" />
+      <param name="input_se" value="input.fastq" />
       <param name="db" value="db.fasta" />
-      <output file="matched.fastq" ftype="fastqsanger" name="reads" />
+      <param name="keep_counts" value="False" />
+      <output file="kat-sect-stats.tsv" name="stats" />
     </test>
+    <test>
+      <param name="input_se" value="input.fastq" />
+      <param name="db" value="db.fasta" />
+      <param name="keep_counts" value="True" />
+      <output file="kat-sect-stats.tsv" name="stats" />
+      <output file="kat-sect-counts.cvg" name="counts" />
+    </test>    
   </tests>
   <help><![CDATA[
-http://kat.readthedocs.io/en/latest/using.html#filtering-tools
+  https://kat.readthedocs.io/en/latest/using.html#sect
 ]]>
   
   </help>
-  <citations>
-  </citations>
+  <expand macro="citations" />
 </tool>
b
diff -r 699eef7398c1 -r ab2cf85b9ae8 macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Mon Mar 13 15:36:43 2017 -0400
b
@@ -0,0 +1,50 @@
+<macros>
+  <token name="@VERSION@">2.3</token>
+  <xml name="requirements">
+    <requirements>
+      <requirement type="package" version="2.3.1">kat</requirement>
+    </requirements>
+  </xml>
+  <xml name="stdio">
+    <stdio>
+      <exit_code range="1:" />
+      <exit_code range=":-1" />
+      <regex match="Error:" />
+      <regex match="Exception:" />
+    </stdio>
+  </xml>
+  <xml name="version_command">
+    <version_command>kat --version</version_command>
+  </xml>
+  <xml name="citations">
+    <citations>
+      <citation type="doi">10.1101/064733</citation>
+      <yield />
+    </citations>
+  </xml>
+  <token name="@THREADS@">
+    --threads \${GALAXY_SLOTS:-4}
+  </token>
+  <xml name="macro_input">
+    <conditional name="single_or_paired">
+      <param name="type" type="select" label="Read type">
+        <option value="single">Single-end</option>
+        <option value="paired">Paired-end</option>
+        <option value="collection">Collection Paired-end</option>
+      </param>
+      <when value="single">
+        <param name="input_se" type="data" format="fastqsanger" label="Single end read file(s)"/>
+      </when>
+      <when value="paired">
+        <param name="forward_pe" type="data" format="fastqsanger" label="Forward paired-end read file"/>
+        <param name="reverse_pe" type="data" format="fastqsanger" label="Reverse paired-end read file"/>
+      </when>
+      <when value="collection">
+        <param name="fastq_collection" type="data_collection" label="Paired-end reads collection" optional="false" format="fastqsanger" collection_type="paired" />
+      </when>
+    </conditional>
+  </xml>
+  <xml name="macro_kmers"> 
+    <param name="kmer" type="integer" value='27' label="Kmer size" help="The kmer length to use in the kmer hashes.  Larger values will provide more discriminating power between kmers but at the expense of additional memory and lower coverage."/>
+  </xml>
+</macros>
b
diff -r 699eef7398c1 -r ab2cf85b9ae8 test-data/input-counts.cvg
--- a/test-data/input-counts.cvg Mon Feb 06 12:36:09 2017 -0500
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
b'@@ -1,3032 +0,0 @@\n->abcZ_1\n-0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 161 159 159 157 162 161 161 161 160 170 170 168 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 143 144 126 129 129 130 131 131 137 143 143 143 138 137 139 139 138 140 140 155 155 155 155 156 165 161 162 159 158 161 160 160 161 160 158 147 146 145 148 148 148 149 147 145 142 143 147 149 153 154 153 154 163 165 161 160 158 158 158 157 159 164 160 161 162 161 161 147 144 144 140 141 140 138 139 135 134 128 128 131 129 128 121 114 114 112 111 106 107 107 109 112 117 117 111 110 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 140 142 142 142 136 135 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 117 119 115 104 100 100 98 99 101 99 100 103 104 106 104 106 107 107 107 105 109 111 111 124 121 122 123 119 124 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 112 110 106 105 107 105 100 103 99 101 98 95 94 97 109 108 106 105 92 96 95 94 102 97 104 109 108 111 110 111 115 115 120 120 123 126 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 142 143 145 140 136 133 133 133 136 136 139 139 141 134 134 135 135 138 135 131 118 117 118 118 119 119 126 126 120 120 114 112 114 110 110 112 109 116 110 111 109 106 108 109 108 112 110 116 116 116 119 116 116 116 115 115 114 114 112 112 112\n->abcZ_2\n-0 0 0 0 0 0 124 124 125 127 128 125 125 112 112 113 112 111 114 124 127 128 129 128 129 132 132 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 161 159 159 157 162 161 161 161 160 170 170 168 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 143 144 126 129 129 130 131 131 137 143 143 143 138 137 139 139 138 140 140 155 155 155 155 156 165 161 162 159 158 161 160 160 161 160 158 147 146 145 148 148 148 149 147 145 142 143 147 149 153 154 153 154 163 165 161 160 158 158 158 157 159 164 160 161 162 161 161 147 144 144 140 141 140 138 139 135 134 128 128 131 129 128 121 114 114 112 111 106 107 107 109 112 117 117 111 110 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 140 142 142 142 136 135 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 117 119 115 104 100 100 98 99 101 99 100 103 104 106 104 106 107 107 107 105 109 111 111 124 121 122 123 119 124 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 112 110 106 105 107 105 100 103 99 101 98 95 94 97 109 108 106 105 92 96 95 94 102 97 104 109 108 111 110 111 115 115 120 120 123 126 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 142 143 145 140 136 133 133 133 136 136 139 139 141 134 134 135 135 138 135 131 118 117 118 118 119 119 126 126 120 120 114 112 114 110 110 112 109 116 110 111 109 106 108 109 108 112 110 116 116 116 119 116 116 116 115 115 114 114 112 112 112\n->abcZ_3\n-0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 143 144 126 129 129 130 131 131 137 143 143 143 138 137 139 139 138 140 140 155 155 155 155 156 165 161 162 159 158 161 160 160 161 160 158 147 146 145 148 148 148 149 147 145 142 143 147 149 153 154 153 154 163 165 161 160 158 158 158 157 159 164 160 161 162 161 161 147 144 144 140 141 140 138 139 135 134 128 128 131 129 128 121 114 114 112 111 106 107 107 109 112 117 117 111 110 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 140 142 142 142 136 135 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 117 119 115 104 100 100 98 99 101 99 100 103 104 106 104 106 107 107 107 105 109 111 111 124 121 122 123 119 124 0 0 0 0 0 0'..b'0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 99 99 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 81 82 84 85 85 70 70 67 66 66 62 61 69 70 70 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 103 102 107 107 106 105 105 110 110 112 117 117 119 119 125 125 133 132 132 124 124 129 129 129 127 127 141 140 146 143 143 135 130 127 134 134 144 141 140 141 141 145 156 155 150 151 148 148 158 156 159 159 157 156 151 152 153 153 149 142 144 143 142 140 139 141 141 141 141 140 140 140 140 140 139 139 139 135 135 139 142 142 138 135 128 124 125 124 125 120 120 116 116 116 120 121 117 114 114 113 115 105 105 100 101 101 101 100 99 94 95 87 87 89 91 91 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 81 79 79 78 78 74 71 74 71 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 78 80 80 79 76 79 80 81 81 82 82 73 70 70 91 90 92 91 91 93 94 93 92 87 91 90 88 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0\n->lhkA_163\n-47 61 61 61 62 62 64 64 64 66 66 67 67 70 72 72 71 71 72 71 81 77 77 77 77 79 81 82 82 83 83 83 83 83 92 92 90 96 101 101 101 99 99 99 114 113 116 117 116 120 121 122 122 124 124 124 120 120 119 117 117 117 116 116 116 122 122 120 120 120 117 117 114 114 118 115 102 102 102 102 102 100 102 101 99 99 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 81 82 84 85 85 70 70 67 66 66 62 61 69 70 70 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 103 102 107 107 106 105 105 110 110 112 117 117 119 119 125 125 133 132 132 124 124 129 129 129 127 127 141 140 146 143 143 135 130 127 134 134 144 141 140 141 141 145 156 155 150 151 148 148 158 156 159 159 157 156 151 152 153 153 149 142 144 143 142 140 139 141 141 141 141 140 140 140 140 140 139 139 139 135 135 139 142 142 138 135 128 124 125 124 125 120 120 116 116 116 120 121 117 114 114 113 115 105 105 100 101 101 101 100 99 94 95 87 87 89 91 91 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 76 76 77 79 80 81 81 79 79 78 78 74 71 74 71 71 71 70 70 70 70 70 70 69 69 70 66 65 66 67 74 74 72 71 72 73 73 74 75 77 77 79 78 80 80 79 76 79 80 81 81 82 82 73 70 70 91 90 92 91 91 93 94 93 92 87 91 90 88 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0\n->lhkA_164\n-47 61 61 61 62 62 64 64 64 66 66 67 67 70 72 72 71 71 72 71 81 77 77 77 77 79 81 82 82 83 83 83 83 83 92 92 90 96 101 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 122 120 120 120 117 117 114 114 118 115 102 102 102 102 102 100 102 101 99 99 100 102 102 101 100 100 100 100 101 92 94 94 94 96 93 93 92 92 92 93 94 94 94 91 90 90 85 81 81 82 84 85 85 70 70 67 66 66 62 61 69 70 70 75 75 75 82 82 88 92 98 99 101 102 97 105 105 105 106 105 107 107 107 105 106 110 110 110 109 109 109 106 106 107 107 114 112 108 109 110 114 115 114 112 113 114 114 122 125 127 125 124 124 124 123 122 121 121 111 113 119 117 120 119 118 120 118 118 118 118 127 127 127 126 127 128 128 128 126 127 128 122 123 116 116 113 110 108 108 107 100 103 102 107 107 106 105 105 110 110 112 117 117 119 119 125 125 133 132 132 124 124 129 129 129 127 127 141 140 146 143 143 135 130 127 134 134 144 141 140 141 141 145 156 155 150 151 148 148 158 156 159 159 157 156 151 152 153 153 149 142 144 143 142 140 139 141 141 141 141 140 140 140 140 140 139 139 139 135 135 139 142 142 138 135 128 124 125 124 125 120 120 116 116 116 120 121 117 114 114 113 115 105 105 100 101 101 101 100 99 94 95 87 87 89 91 91 87 77 75 75 76 85 88 77 77 77 78 79 79 76 75 77 79 85 85 80 83 83 83 83 77 77 76 76 76 77 79 80 81 81 79 79 78 78 74 71 74 71 71 71 70 70 70 70 70 70 69 69 70 66 65 66 67 74 74 72 71 72 73 73 74 75 77 77 79 78 80 80 79 76 79 80 81 81 82 82 73 70 70 91 90 92 91 91 93 94 93 92 87 91 90 88 89 104 105 107 109 111 114 116 115 113 113 112 115 116 116 118 118 117 117 114 114 114 114 114 113\n'
b
diff -r 699eef7398c1 -r ab2cf85b9ae8 test-data/input-stats.tsv
--- a/test-data/input-stats.tsv Mon Feb 06 12:36:09 2017 -0500
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
b'@@ -1,1517 +0,0 @@\n-seq_name\tmedian\tmean\tgc%\tseq_length\tinvalid_bases\t%_invalid\tnon_zero_bases\t%_non_zero\t%_non_zero_corrected\n-abcZ_1\t1\t60.48532\t0.37616\t537\t0\t0.00000\t280\t54.79452\t54.79452\n-abcZ_2\t96\t65.53620\t0.37803\t537\t0\t0.00000\t301\t58.90411\t58.90411\n-abcZ_3\t0\t49.75147\t0.37058\t537\t0\t0.00000\t237\t46.37965\t46.37965\n-abcZ_4\t100\t68.67906\t0.38547\t537\t0\t0.00000\t304\t59.49119\t59.49119\n-abcZ_5\t129\t130.83366\t0.37616\t537\t0\t0.00000\t511\t100.00000\t100.00000\n-abcZ_6\t116\t88.17613\t0.37058\t537\t0\t0.00000\t354\t69.27593\t69.27593\n-abcZ_7\t127\t121.54207\t0.37616\t537\t0\t0.00000\t478\t93.54207\t93.54207\n-abcZ_8\t124\t114.06458\t0.37430\t537\t0\t0.00000\t451\t88.25832\t88.25832\n-abcZ_9\t129\t125.37573\t0.37430\t537\t0\t0.00000\t484\t94.71624\t94.71624\n-abcZ_10\t0\t57.11742\t0.36872\t537\t0\t0.00000\t247\t48.33659\t48.33659\n-abcZ_11\t105\t70.80822\t0.37616\t537\t0\t0.00000\t301\t58.90411\t58.90411\n-abcZ_12\t0\t58.73973\t0.37058\t537\t0\t0.00000\t252\t49.31507\t49.31507\n-abcZ_13\t100\t68.67906\t0.38361\t537\t0\t0.00000\t304\t59.49119\t59.49119\n-abcZ_14\t129\t124.47554\t0.37803\t537\t0\t0.00000\t484\t94.71624\t94.71624\n-abcZ_15\t129\t124.71820\t0.37430\t537\t0\t0.00000\t484\t94.71624\t94.71624\n-abcZ_16\t126\t115.42661\t0.37430\t537\t0\t0.00000\t451\t88.25832\t88.25832\n-abcZ_17\t124\t114.28180\t0.37803\t537\t0\t0.00000\t451\t88.25832\t88.25832\n-abcZ_18\t0\t56.73777\t0.38175\t537\t0\t0.00000\t248\t48.53229\t48.53229\n-abcZ_19\t106\t74.37573\t0.37244\t537\t0\t0.00000\t300\t58.70841\t58.70841\n-abcZ_20\t99\t67.11546\t0.38920\t537\t0\t0.00000\t283\t55.38160\t55.38160\n-abcZ_21\t124\t114.21918\t0.37430\t537\t0\t0.00000\t451\t88.25832\t88.25832\n-abcZ_22\t1\t59.81018\t0.38361\t537\t0\t0.00000\t258\t50.48924\t50.48924\n-abcZ_23\t124\t96.90802\t0.37989\t537\t0\t0.00000\t371\t72.60274\t72.60274\n-abcZ_24\t101\t68.25049\t0.36499\t537\t0\t0.00000\t266\t52.05479\t52.05479\n-abcZ_25\t0\t20.09980\t0.35196\t537\t0\t0.00000\t82\t16.04697\t16.04697\n-abcZ_26\t0\t20.09980\t0.35568\t537\t0\t0.00000\t82\t16.04697\t16.04697\n-abcZ_27\t0\t60.90215\t0.37803\t537\t0\t0.00000\t249\t48.72798\t48.72798\n-abcZ_28\t0\t20.09980\t0.35382\t537\t0\t0.00000\t82\t16.04697\t16.04697\n-abcZ_29\t116\t79.26027\t0.37430\t537\t0\t0.00000\t304\t59.49119\t59.49119\n-abcZ_30\t127\t122.29354\t0.37616\t537\t0\t0.00000\t484\t94.71624\t94.71624\n-abcZ_31\t1\t60.60665\t0.38361\t537\t0\t0.00000\t277\t54.20744\t54.20744\n-abcZ_32\t106\t73.50881\t0.37616\t537\t0\t0.00000\t318\t62.23092\t62.23092\n-abcZ_33\t115\t84.36204\t0.36872\t537\t0\t0.00000\t342\t66.92759\t66.92759\n-abcZ_34\t101\t69.33659\t0.38361\t537\t0\t0.00000\t307\t60.07828\t60.07828\n-abcZ_35\t0\t52.48532\t0.37430\t537\t0\t0.00000\t253\t49.51076\t49.51076\n-abcZ_36\t0\t24.93738\t0.35754\t537\t0\t0.00000\t100\t19.56947\t19.56947\n-abcZ_37\t1\t61.95499\t0.38920\t537\t0\t0.00000\t258\t50.48924\t50.48924\n-abcZ_38\t0\t60.28767\t0.38547\t537\t0\t0.00000\t248\t48.53229\t48.53229\n-abcZ_39\t1\t65.35421\t0.37430\t537\t0\t0.00000\t272\t53.22896\t53.22896\n-abcZ_40\t0\t65.14481\t0.38175\t537\t0\t0.00000\t253\t49.51076\t49.51076\n-abcZ_41\t105\t72.06067\t0.37244\t537\t0\t0.00000\t300\t58.70841\t58.70841\n-abcZ_42\t114\t85.59883\t0.37803\t537\t0\t0.00000\t368\t72.01566\t72.01566\n-abcZ_43\t100\t69.15264\t0.38547\t537\t0\t0.00000\t285\t55.77299\t55.77299\n-abcZ_44\t111\t77.55773\t0.37244\t537\t0\t0.00000\t336\t65.75342\t65.75342\n-abcZ_45\t0\t48.97451\t0.36940\t536\t0\t0.00000\t233\t45.68627\t45.68627\n-abcZ_46\t1\t57.47554\t0.37989\t537\t0\t0.00000\t275\t53.81605\t53.81605\n-abcZ_47\t0\t58.05871\t0.35940\t537\t0\t0.00000\t227\t44.42270\t44.42270\n-abcZ_48\t100\t67.31311\t0.36127\t537\t0\t0.00000\t263\t51.46771\t51.46771\n-abcZ_49\t0\t57.30333\t0.38734\t537\t0\t0.00000\t229\t44.81409\t44.81409\n-abcZ_50\t115\t85.36008\t0.37616\t537\t0\t0.00000\t334\t65.36204\t65.36204\n-abcZ_51\t1\t65.09393\t0.37803\t537\t0\t0.00000\t281\t54.99022\t54.99022\n-abcZ_52\t0\t57.37182\t0.37616\t537\t0\t0.00000\t239\t46.77104\t46.77104\n-abcZ_53\t1\t65.42270\t0.38175\t537\t0\t0.00000\t280\t54.79452\t54.79452\n-abcZ_54\t112\t80.39530\t0.36872\t537\t0\t0.00000\t329\t64.38356\t64.38356\n-abcZ_55\t116\t87.63992\t0.36872\t537\t0\t0.00000\t352\t68.88454\t68.88454\n-abcZ_56\t102\t68.86497\t0.37430\t537\t0\t0.00000\t288\t56.36008\t56.36008\n-abcZ_57\t128\t123.44814\t0.37803\t537\t0\t0.00000\t484\t94.71624\t94.71624\n-abcZ_58\t110\t76.94912\t0.37430\t537\t0\t0.00000\t326\t63.79648\t63.79648\n-abcZ_59\t98\t66.45205\t0.38361\t537\t0\t0.00000\t270\t52.83757\t52.83757\n-abcZ_60\t1\t56.99804\t0.37616\t'..b'kA_101\t0\t0.78194\t0.35417\t480\t0\t0.00000\t3\t0.66079\t0.66079\n-lhkA_102\t0\t3.23568\t0.37083\t480\t0\t0.00000\t22\t4.84581\t4.84581\n-lhkA_103\t99\t86.44714\t0.37292\t480\t0\t0.00000\t373\t82.15859\t82.15859\n-lhkA_104\t102\t93.24229\t0.37083\t480\t0\t0.00000\t400\t88.10573\t88.10573\n-lhkA_105\t0\t0.00000\t0.37500\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_106\t0\t3.01542\t0.37708\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_107\t0\t3.01542\t0.37083\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_108\t0\t0.86344\t0.37292\t480\t0\t0.00000\t21\t4.62555\t4.62555\n-lhkA_109\t0\t3.01542\t0.37083\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_110\t0\t3.01322\t0.37708\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_111\t0\t3.01322\t0.36875\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_112\t0\t3.01322\t0.36875\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_113\t0\t3.01322\t0.37083\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_114\t0\t3.01322\t0.37292\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_115\t0\t43.10132\t0.36458\t480\t0\t0.00000\t173\t38.10573\t38.10573\n-lhkA_116\t0\t43.10132\t0.36458\t480\t0\t0.00000\t173\t38.10573\t38.10573\n-lhkA_117\t100\t96.24890\t0.37292\t480\t0\t0.00000\t427\t94.05286\t94.05286\n-lhkA_118\t0\t0.00000\t0.38542\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_119\t0\t44.53084\t0.36667\t480\t0\t0.00000\t179\t39.42731\t39.42731\n-lhkA_120\t0\t0.00000\t0.38750\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_121\t0\t0.00000\t0.41042\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_122\t0\t0.00000\t0.44167\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_123\t0\t0.00000\t0.38958\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_124\t0\t0.00000\t0.35208\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_125\t0\t0.00000\t0.35625\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_126\t0\t0.00000\t0.34792\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_127\t0\t40.73568\t0.36250\t480\t0\t0.00000\t162\t35.68282\t35.68282\n-lhkA_128\t0\t0.00000\t0.34583\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_129\t105\t98.56167\t0.37708\t480\t0\t0.00000\t427\t94.05286\t94.05286\n-lhkA_130\t71\t58.76211\t0.36875\t480\t0\t0.00000\t256\t56.38767\t56.38767\n-lhkA_131\t0\t43.10132\t0.36875\t480\t0\t0.00000\t173\t38.10573\t38.10573\n-lhkA_132\t0\t36.10793\t0.36875\t480\t0\t0.00000\t146\t32.15859\t32.15859\n-lhkA_133\t0\t42.57048\t0.36667\t480\t0\t0.00000\t171\t37.66520\t37.66520\n-lhkA_134\t0\t2.22467\t0.34375\t480\t0\t0.00000\t9\t1.98238\t1.98238\n-lhkA_135\t105\t96.63656\t0.37292\t480\t0\t0.00000\t413\t90.96916\t90.96916\n-lhkA_136\t0\t38.44053\t0.36458\t480\t0\t0.00000\t148\t32.59912\t32.59912\n-lhkA_137\t0\t0.00000\t0.34583\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_138\t0\t2.22467\t0.34792\t480\t0\t0.00000\t9\t1.98238\t1.98238\n-lhkA_139\t0\t48.22247\t0.36667\t480\t0\t0.00000\t189\t41.62996\t41.62996\n-lhkA_140\t100\t96.41189\t0.37292\t480\t0\t0.00000\t427\t94.05286\t94.05286\n-lhkA_141\t0\t36.07048\t0.36250\t480\t0\t0.00000\t144\t31.71806\t31.71806\n-lhkA_142\t94\t87.83260\t0.37500\t480\t0\t0.00000\t400\t88.10573\t88.10573\n-lhkA_143\t0\t0.00000\t0.35208\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_144\t0\t3.01542\t0.37083\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_145\t0\t2.18943\t0.36667\t480\t0\t0.00000\t26\t5.72687\t5.72687\n-lhkA_146\t0\t40.17401\t0.36250\t480\t0\t0.00000\t156\t34.36123\t34.36123\n-lhkA_147\t0\t3.01542\t0.36875\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_148\t0\t0.72026\t0.35000\t480\t0\t0.00000\t3\t0.66079\t0.66079\n-lhkA_149\t100\t96.39868\t0.37500\t480\t0\t0.00000\t427\t94.05286\t94.05286\n-lhkA_150\t0\t3.01322\t0.37500\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_151\t0\t50.26211\t0.37083\t480\t0\t0.00000\t200\t44.05286\t44.05286\n-lhkA_152\t0\t2.22467\t0.34792\t480\t0\t0.00000\t9\t1.98238\t1.98238\n-lhkA_153\t0\t3.01542\t0.37292\t480\t0\t0.00000\t30\t6.60793\t6.60793\n-lhkA_154\t0\t43.10132\t0.36667\t480\t0\t0.00000\t173\t38.10573\t38.10573\n-lhkA_155\t0\t2.15639\t0.37083\t480\t0\t0.00000\t11\t2.42291\t2.42291\n-lhkA_156\t0\t0.00000\t0.35625\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_157\t0\t0.00000\t0.35417\t480\t0\t0.00000\t0\t0.00000\t0.00000\n-lhkA_158\t0\t3.68943\t0.36875\t480\t0\t0.00000\t33\t7.26872\t7.26872\n-lhkA_159\t0\t2.22467\t0.35000\t480\t0\t0.00000\t9\t1.98238\t1.98238\n-lhkA_160\t0\t3.01322\t0.37292\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_161\t0\t3.01322\t0.37500\t480\t0\t0.00000\t29\t6.38767\t6.38767\n-lhkA_162\t0\t42.83700\t0.36458\t480\t0\t0.00000\t172\t37.88546\t37.88546\n-lhkA_163\t76\t65.49119\t0.37083\t480\t0\t0.00000\t288\t63.43612\t63.43612\n-lhkA_164\t101\t96.63656\t0.37500\t480\t0\t0.00000\t427\t94.05286\t94.05286\n'
b
diff -r 699eef7398c1 -r ab2cf85b9ae8 test-data/matched.fastq
--- a/test-data/matched.fastq Mon Feb 06 12:36:09 2017 -0500
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
b'@@ -1,24288 +0,0 @@\n-@HWI-ST1348:99:H0D4GADXX:1:1101:12325:8854/1\n-CTTATAATGAACACATTGATCGCTTATATGCTAGTGCAGCAAAAATTGACTTAGTTATTCCTTATTCCAAAGAAGAGCTACGTGAATTACTTGAAAAATTA\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJIIJJJJJJJJJJJJJJJJJJIJJJJJJJJJJJJJJJJJJJJJJHHHHHFFFFEEEEEEEDDDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:16009:11807/1\n-TTTCCGTACTTCCATCGCTTGGTCCCGTATCTTCCCAAATGGTGACGAAACTGAGCCAAACGAAGCAGGACTACAATTTTACGATGATTTATTCGATGAAC\n-+\n-CCCFFFFFHHHHHJJJJJJJJJIIJJJHIJJJJJJJJJJJJJFHIIJJJJJJJJJJJJJIJHHFFFEEDDDDDDDDDDDDFDDBDDDFEEDEDDDEEEDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:13122:12268/1\n-GAAAAGGACTTTCCGTTCAAGATGTTACTCCAAAAGGCGGATTCGGTCACATTACTGACGGTCCAACACCAGATAACTTAAAATTAGAAGGAATCGACTTC\n-+\n-CCCFFFFFHHHHHJJIJJJJJJJJIJJJJJJJJJJJJJJJJJJJJJHIIJJJJHHHHHHFFCDDEDDDDDDDDDDDDDEDDDDEDDDDDDDDDDDDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:3847:12540/1\n-CTTGGTTTAGTGTATTGGAAATATTATCTAAGTCGCTAGTGAAGCGGCTAAGCATTTCGCCATCATTGCGGCTATCGAAGAAACGAATCGATAGTTTTTCC\n-+\n-CCCFFDDFHHHHHJJJIIJJJJJJJJJIJJJJIHJIIIJDHGJJJJJJJJJJJJJJJJIJJIGHHHHFFFDDDDDDDDDDDDDDDDDDDDDDDD:BDDDCD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:15695:12654/1\n-GGCAGTTACAACAACTAGGTCCGCATCGTGGCAGTCGCTATAATTTGCTGAGTAGATTTTCTTCGGTGTAGAAAATGGAACGGCATGGCTTAAATCCATTG\n-+\n-CCCFFFFFHHHHHJJJJJJHHIJJJJJJJIJJJJIJIIIJJJJJJJJJJJGIIJJGJJJJJHHHHF?DDEEEEEEDDDCDDDDDDDDDDDDDDDECDDDED\n-@HWI-ST1348:99:H0D4GADXX:1:1101:11189:13900/1\n-GTTTAGTGTATTGGAAATATTATCTAAGTCGCTAGTGAAGCGGCTAAGCATTTCGCCATCATTGCGGCTATCGAAGAAACGAATCGATAGTTTTTCCATCT\n-+\n-CCCFFFEFHHHHHJGJIJJJJJJJJJJIHJJJJJJGGHIIJJJJJJJJJJJJJIIIIIJJJJIHHHFFDDDDDDDDDDDDDBDBDDDDDDCEEDDDDEEDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:10723:17926/1\n-TCTCTATTTCTGCGGATGATAAGAAAAAGCTTGAAAGTGGTGAAACGGTTAGTAAGAAAATAGATAATCGCTTTAACAAAGAAATGACAATTGTGTACGTC\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJIJJJIJJJIJJJJIFHIIJJIJJJJIIGIFHIJJJJJJJJJHGHHHHFFFDDEDEDEDDDDDDDDDDDDDDDDEEDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:20826:19324/1\n-GTTCACCGATGATTAAACCATCTAAATCATCTGCGTAACCTTTTTGTGTTAGTTGTTCTGCTCCAAGTTCACCGATTTCTTCCCCAACTGTTGCTAATAAT\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJIJJJJJJJJJJHIJJIJIJJJJJJJJIJJJJHHHHHHFFEDDDDEDDDDDDDDDDDDDDDEDDDDE\n-@HWI-ST1348:99:H0D4GADXX:1:1101:11646:21325/1\n-GAGTAACGAAGGCTAGTTCCACGTTTTGTTGGAACATCATGATGATAACACCAATCATTAGCGCGACGTTGGATAGTACTTGGATCAATGCTTGGTTTAGT\n-+\n-@@@?DEFFHHHGHIJIIJIGIIJFGIJJJJJIIGIIJJJJIEGGIJJJJJIJJJJJJIIIIIHIGDDDDBDDCCEDDEDEEDDDDDDDDDDDECD@BDCDC\n-@HWI-ST1348:99:H0D4GADXX:1:1101:13141:21911/1\n-AGCCTTGATACAGCACGTTTCCGTATGTCAATTGCTGACTATCTAAAAGTAGATGCTCGTAACGTCCATGGTTACATCCTTGGCGAACACGGCGATACAGA\n-+\n-CCCFFFFFHGHGHJJJJJJJJJJIIGIJI:::?CGHIJJJJJJJJJIJJDHHIJGHIIJBGHIGHIIJJJHEHHHHFFFFFFEEDDDBDDDDDDDDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:10800:28631/1\n-CTAGTGAAGCGGCTAAGCATTTCGCCATCATTGCGGCTATCGAAGAAACGAATCGATAGTTTTTCCATCTTGCGGAAAAGCCCTATACGCATTCTGTTCGT\n-+\n-@CCFDFFFHHHHHJJJIJJJJJJJJJJJJJJJJJJJJJIJJJIJJJJJJHHHHFFFFDDEEEDDDDEEDDDDDDDDDDDDDDDDDDDDDDDDDDEEEDDDB\n-@HWI-ST1348:99:H0D4GADXX:1:1101:20491:28620/1\n-ATTCAAGGCAATTACAACGCTCGCTTGAAGGAAAATAATTTTGATGAAATTGGTGCACTCGCCATTGATTTCAATAAAATGACACAAACCCTTGAAACATC\n-+\n-CCCFFFFFHHHHHJJIJIJJIJJJJJJJJJJJJIJJJJIJJJJJJJJJJIJJJFHIIJJJIHHHFFFFFFEEEEEEDDDEDCDDDDDDBDDDDDDDDDDDC\n-@HWI-ST1348:99:H0D4GADXX:1:1101:7329:32108/1\n-ATAGTATCCATTGCACCTTGTTCATCTTCGCTAATCCATTCAGTAATTGGAAGGCCGCCGACAGTTGTGTGCCTCCATGCTGGGAATTCTGTATCGCCGTG\n-+\n-BCCFDFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJIIIJJJJJIJJJJJJJJHHFFFEEEDEDDDDDDDDDDDDDDDDDDDEEEEEEDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:6537:33082/1\n-AGTTCAATCATTGCAATAACCATCGCCGCTAGACCTGACTTCATATCCGTTGCGCCGCGTCCGTATAGTTTCCCTTCATGCTCTGTCGCTTCAAAAGGTGG\n-+\n-BCBFFFFFHHHHHJJJJJJJJJJJIIJIJJJIJJJJJIJJJJJJJJJJJIJIJJIJHFDDDDDBDDDEDFEEDDDDDDEDDDDDDDDDDDDDDDDDDD@CB\n-@HWI-ST1348:99:H0D4GADXX:1:1101:8396:34357/1\n-CCGCATGCGTCCATAATACACCATCTTTAATAATAGAAACGTTTGAAGCAGAACATTCTGTTACTTGTTCCCCGCGATGTAAAATAGCTTCCAAAGCATTT\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJIJJJJJJJJJJJJJJJJIJJJJJJJJJJJJJJJJJHHHFDDDDDE@C@;ACDDDDDDDDDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:1:1101:4455:40892/1\n-GGTAAATATTGATTCGACTAATGGGCATGAAGAACAAGTTGCGAACTATTTGCAAAAGTTGTTAGCTGAACACGGTATTGAGTCCGAAAAGGTACAA'..b'CTTGAAACATCTCAAGAAGAAATTGAACGACAAGAGAAACGGAGACGCCAGTTTATTGCTGA\n-+\n-BB1BDFFFFHHHFHIJJJIJJJJJJIJJJIGHIIIJJJJJEIJJJJGGFEGIGIIIJJIJGHIGGIIGHGFFFED@CECDD?@???BDDDCDEECDCDCC@\n-@HWI-ST1348:99:H0D4GADXX:2:2216:7148:70790/1\n-CTTGTCCAAACTACAAGTAAATCGAATCAACAAACTACGAGCAGCGACAAAAGACGTTATTCAAGGCAATTACAACGCTCGCTTGAAGGAAAATAATTTTG\n-+\n-?@@DFFFEDFDBHEHGGGGIIIGHGFGIJGIJGGGACHH@GGGHGIG6FEF=A1CEAEEFFFFEFC?AACCDCCDC;;??;@?BBDDDD8>><CCDDDDD:\n-@HWI-ST1348:99:H0D4GADXX:2:2216:3356:73770/1\n-AACGATAGCTGTATTTAATAGGGAAATACCTTGCATCATTGGGAAAAGTAAACCGGAGTAAACTTGCCCTTTATAAGTGGCGTTTTTAACGATATTGTTTT\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJIJJJJJIJJJJJJJJJIJJJJJJJIFHIJJJJJJJGHHIJIJJJHHHHHFFFFFADEDDDDDDDDDDDDDDDEEDDDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:5940:75462/1\n-GCTGAAATGGGCTTCAAGGTTTTCCGTACTTCCATCGCTTGGTCCCGTATCTTCCCAAATGGTGACGAAACTGAGCCAAACGAAGCAGGACTACAATCTGT\n-+\n-@@@FDDDDHGHHGJIIIJJHIGIIJJFFHIEHJIJIIJGGIJDDHGIGIJJIIJJJIJJIHI7@CHHEEFFFEACCCCDDDD=BB@ACDBBCACCCADCD:\n-@HWI-ST1348:99:H0D4GADXX:2:2216:16564:76685/1\n-GTATTGGAAATATTATCTAAGTCGCTAGTGAAGCGGCTAAGCATTTCGCCATCATTGCGGCTATCGAAGAAACGAATCGATAGTTTTTCCATCTTGCGGAA\n-+\n-C@@FFFFFHFBFDIIGGIJDHIIGIIIIHIG9GHIIGGEFHEEIIIGGGEGGJIJEGGHEEBBCCDDCBBDCCBBD??BBB25:@CCD?A@ACCCD@BB>>\n-@HWI-ST1348:99:H0D4GADXX:2:2216:11696:77729/1\n-TTCGATGAACTTCTAGCACATAATATCGAACCACTAATCACTTTATCTCACTATGAAACACCACTTCACTTATCGAAAACTTACGACGGCTGGGTAAATAG\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJJJJJIJJJJJJJJJJIJJJJJIJJIJJJJIIJIJJJJJJJIJJJJJIJJJJHHHHFFFFDEDDDDDDDBDDA@9::@>\n-@HWI-ST1348:99:H0D4GADXX:2:2216:1712:78154/1\n-CAATGATGGCGAAATGCTTAGCCGCTTCACTAGCGACTTAGATAATATTTCCAATACACTAAACCAAGCATTGATCCAAGTACTATCCAACGTCGCGCTAA\n-+\n-CCCFFEFFHHHHHIJJIJIJJJIJJJJIJJJGIIJGIIIIIJJJIIIIJJIIIIJJIFHGHHHFFFFECEEEDDEDDDDDDDDEEDDDDCCBBBDDDDDDB\n-@HWI-ST1348:99:H0D4GADXX:2:2216:8423:78377/1\n-GTTGATGCGGGTGATGTATCTAAGTGGAAGTTCCCACCTTTTGAAGCGACAGAGCATGAAGGGAAACTATACGGACGCGGCGCAACGGATATGAAGTCAGG\n-+\n-CBCFFFFFHHH@FGIJHJJJJJJJIHIIJJHJJJJIJJJJJIJJJJIJIJJJIIHHHHFFFFFDEDEEEDDDDDDDDDDDBDDBDDDDDDDDDEEDDDEDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:4026:81369/1\n-CTGTATCGCCGTGTTCGCCAAGGATGTAACCATGGACGTTACGAGCATCTACTTTTAGATAGTCAGCAATTGACATACGGAAACGTGCTGTATCAAGGCTT\n-+\n-CCCFFFFFHHHHHJJJJJJJJJJIJJJJJJJJJJJIJIJJJJJJIJJJIIJHHHHHHDEEFFFFEEEEEEDDDDDDDDDDDDDDDADDDDDDEEEDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:9367:83208/1\n-CATCTAAATCATCTGCGTAACCTTTTTGTGTTAGTTGTTCTGCCCCAAGTTCACCGATTTCTTCCCCAACTGTTGCTAATAATCTGATCTTGCCGTTTAGT\n-+\n-CCCFFFFFHHHGHIJJJGGJJJJJJJIIFHGHIIIJJJJIJIJJJJJJJJJIIJJJHHIGIJJJHHHHFFFFDEEEDEDDCDEDDECCDDDDDDDDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:9137:83619/1\n-ATTAAATACGGTGCGGACATAGTTTTCATAGAAGTCGATCATTTTTCTATTTACCCAGCCGTCGTAAGTTTTCGATAAGTGAAGTGGTGTTTCATAGTGAG\n-+\n-@@CDFFDFHFHFDHHHGJIJJGGGIIGIJG?BFDFEDGIIIJJJIIIGIJIJJIIIIGHHIBEEFDECBCECCBBFBDDDCDEDCDDC@ABDCBCDDEEEC\n-@HWI-ST1348:99:H0D4GADXX:2:2216:4540:85153/1\n-GTGTTTCATAGTGAGATAAAGTGATTAGTGGTTCGATATTATGTGCTAGAAGTTCATCGAATAAATCATCGTAAAATTGTAGTCCTGCTTCGTTTGGCTCA\n-+\n-BB@FFFFFHHHHHGIJIGIJJ?FHIIJJEGHIJIJIIIJJHHIFGIJJJGGIGGIJJIJJJIJJJJIJJIJJJJJJJHHEHHHHFFFFFEEDEEDDDDCDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:14434:93168/1\n-GTTTTTCCATCTTGCGGAAAAGCCCTATACGCATTCTGTTCGTTGATTTACCAGCTACCCCTGTCATAATGATACTTTGAATAAAAGTGGCAGCAGCGAGC\n-+\n-CCCFFFFFHHHHHJIJJJJJJIJJJJJIIJJJJJJJJJJJJJIIJJJJJJJJIJIIJJJJHHGEFFFFFFFEDDEEEEDADDEDDDDACDD?BBDDDDBDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:8036:96051/1\n-CATCTGGAAACTTAATCGCATCACGAATGAAGAATACTGGCAAATTATTTCCGACAAAGTCATAATTTCCCTCTTCCGTATAAAACTTAACGGAGAAACCG\n-+\n-@CCFFFFFHHHHHJJJJJJJJJJJIJJJJJJJJIJJJJIIHGIJIIIIJJJJGIIIJJHC@EGHHHHHHHFFFFFFEDEEDDEEDDDCDDDBBDBDDDDDD\n-@HWI-ST1348:99:H0D4GADXX:2:2216:6952:98103/1\n-ATTTGTGATAATGATTTTTTGTCCAGAGATTTTTTCGTCAATGTAGCCGTTAAGTACGCCTAGTTCATCTTGTTGAACATCAACGAATTTACGTGCTTTTC\n-+\n-+18BD,2AC>DDDDAEEICE??49C@D?)41?C4?BDDDD6BC<9*/8:7;;;;77@DDD25?;>?>@A;;>AA?AAAAAA3:,3==59A###########\n-@HWI-ST1348:99:H0D4GADXX:2:2216:11200:100390/1\n-AACGTGCTGTATCAAGGCTTGTTCCAGAACCGATAACACGTTCTTTTGGAAGACCTGAGAATTTCCATGTAGCGTAAGTTAAGATGTCTACTGGGTTAGAA\n-+\n-:?14ADDFHHHHHJJJJJJJJJJJJJJJJJJJHJJJJEGIHHIJJJJJJJJIIJJIJIGHHHHHFFFFFFFEEDDDDDDDEEDDDDCCCDDDDDDCBDDDC\n'