# HG changeset patch
# User iuc
# Date 1510656824 18000
# Node ID daad7189e48ea4a73f778a990d7b2e905bb9708b

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit ae3b9baaf7040ed8b165d17466b8b2fe016d3d48

diff -r 000000000000 -r daad7189e48e macros.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,13 @@
+<macros>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@VERSION@">bctools</requirement>
+        </requirements>
+    </xml>
+    <token name="@VERSION@">0.2.1</token>
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1016/j.molcel.2013.07.001</citation>
+        </citations>
+    </xml>
+</macros>
diff -r 000000000000 -r daad7189e48e merge_pcr_duplicates.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/merge_pcr_duplicates.xml	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,51 @@
+<tool id="bctools_merge_pcr_duplicates" name="Merge PCR duplicates" version="@VERSION@">
+    <description>according to UMIs</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+    <command detect_errors="exit_code"><![CDATA[
+        merge_pcr_duplicates.py
+        '$alignments_bed'
+        '$barcode_library'
+        --outfile '$events'
+    ]]></command>
+    <inputs>
+        <param name="alignments_bed" type="data" format="bed" label="BED6 file containing alignments" />
+        <param name="barcode_library" type="data" format="fastq" label="FASTQ barcode library" />
+    </inputs>
+    <outputs>
+        <data name="events" format="bed"/>
+    </outputs>
+    <tests>
+        <test>
+            <param name="alignments_bed" value="pcr_dupes_sorted_2.bed"/>
+            <param name="barcode_library" value="pcr_dupes_randomdict.fastq"/>
+            <output name="events" file="merged_pcr_dupes.bed"/>
+        </test>
+    </tests>
+    <help><![CDATA[
+
+bctools - Merge PCR duplicates according to UMIs
+================================================
+
+Merge PCR duplicates according to unique molecular identifier (UMI) library. All alignments with
+same outer coordinates and barcode will be merged into a single crosslinking
+event.
+
+Barcodes containing uncalled base 'N' are removed.
+
+Input:
+------
+
+* BED6 file containing alignments with FASTQ read-id in name field
+* FASTQ library of UMIs
+
+Output:
+-------
+
+* BED6 file with a read id from a representative alignment in the name field and number of PCR duplicates as score, sorted by fields chrom, start, stop, strand, name
+
+    ]]></help>
+    <expand macro="citations"/>
+</tool>
diff -r 000000000000 -r daad7189e48e test-data/converted_bcs.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/converted_bcs.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,400 @@
+@UID0
+NYYRRY
++
+#0<FFF
+@UID1
+RRYRYR
++
+BBBFFF
+@UID2
+RYYRRY
++
+B<BF<B
+@UID3
+YRYRYR
++
+BBBFFF
+@UID4
+NRRRRY
++
+#0<FFF
+@UID5
+RRRRRY
++
+BBBFFF
+@UID6
+YRRRYR
++
+BBBFFF
+@UID7
+NYYRRY
++
+#0<FFF
+@UID8
+RYRRRY
++
+BBBFFF
+@UID9
+RRRRRY
++
+BBBFFF
+@UID10
+NRRRRY
++
+#0<FFF
+@UID11
+NYRRYR
++
+#0<FFF
+@UID12
+NRYRRY
++
+#0<FFF
+@UID13
+NRRRRY
++
+#0<FFF
+@UID14
+YYRRYR
++
+BBBFFF
+@UID15
+YYRRRY
++
+BBBFFF
+@UID16
+YYYRYR
++
+BBBFFF
+@UID17
+YRRRYR
++
+BBBFFF
+@UID18
+RYYRRY
++
+BBBFFF
+@UID19
+RYYYYY
++
+BBBFFF
+@UID20
+RYRRYR
++
+BBBFFF
+@UID21
+RYRRYR
++
+BBBFFF
+@UID22
+YYRRYR
++
+BBBFFF
+@UID23
+YYYRYR
++
+BBBFFF
+@UID24
+YRRRYR
++
+BBBFFF
+@UID25
+RRYYRR
++
+BBBFFF
+@UID26
+YYRRRY
++
+BBBFFB
+@UID27
+RYYRRY
++
+BBBFFF
+@UID28
+RYYRYR
++
+BBBFFF
+@UID29
+NYYRYR
++
+#0<FFF
+@UID30
+NYRRYR
++
+#0<FFF
+@UID31
+YRYRRY
++
+BBBFFF
+@UID32
+YRYRRY
++
+BBBFFF
+@UID33
+RYYRRY
++
+BBBFFF
+@UID34
+RRRRYR
++
+BBBFFF
+@UID35
+NYYRRY
++
+#0<FFF
+@UID36
+YYRRYR
++
+BBBFFF
+@UID37
+RRRRRY
++
+BBBFFF
+@UID38
+RRRRRY
++
+BBBFFF
+@UID39
+YYRRRY
++
+BBBFFF
+@UID40
+YYRRRY
++
+BBBFFF
+@UID41
+RYYRYR
++
+BBBFFF
+@UID42
+YRYRYR
++
+BBBFFF
+@UID43
+RRRRYR
++
+BBBFFF
+@UID44
+RYRRYR
++
+BBBFFF
+@UID45
+NYYRRY
++
+#0<FBF
+@UID46
+YYYRYR
++
+BBBFFF
+@UID47
+RRRRRY
++
+BBBFFF
+@UID48
+RYRRRY
++
+BBBFFF
+@UID49
+RYRRYR
++
+BBBFFF
+@UID50
+RYYRRY
++
+BBBFFF
+@UID51
+YRYYRR
++
+BBBFFF
+@UID52
+RYRRRY
++
+BBBFFF
+@UID53
+RRRRRY
++
+BBBF<F
+@UID54
+YYRRRY
++
+BBBFFF
+@UID55
+RYYRRY
++
+BBBFFF
+@UID56
+YYYRRY
++
+BBBFFF
+@UID57
+NYYYRY
++
+#0BFFF
+@UID58
+NRRRRY
++
+#0<FFF
+@UID59
+NYRRRY
++
+#0<FFF
+@UID60
+RRYRRY
++
+BB<FFF
+@UID61
+RRRRYR
++
+BBBFFF
+@UID62
+YRRRRY
++
+BBBFFF
+@UID63
+RRRRRY
++
+BBBFFF
+@UID64
+YRRRYR
++
+BBBFFF
+@UID65
+RRYRYR
++
+BBBFFF
+@UID66
+RRRRRY
++
+BBBFFF
+@UID67
+RRYRYR
++
+BB<FFF
+@UID68
+NYYRRY
++
+#00FFF
+@UID69
+YYRRRY
++
+BBBFFF
+@UID70
+YYYRYR
++
+BBBFFF
+@UID71
+RRRRYR
++
+BBBFFF
+@UID72
+RRRRYR
++
+BBBFFF
+@UID73
+NRYRYR
++
+#0<BFF
+@UID74
+YYYRYY
++
+BBBFFF
+@UID75
+YYRRRY
++
+BBBFFF
+@UID76
+YRRRRY
++
+BBBFFF
+@UID77
+NRYRYR
++
+#0<FFF
+@UID78
+NRRRRY
++
+#0<FFF
+@UID79
+YYYRRY
++
+BBBFFF
+@UID80
+RYYRRY
++
+BBBFFF
+@UID81
+NYYRYR
++
+#0<FFF
+@UID82
+NRRRRY
++
+#0<FFF
+@UID83
+NYRRRY
++
+#0<FFF
+@UID84
+NRRRRY
++
+#0<FFF
+@UID85
+RYRRRY
++
+BBBFFF
+@UID86
+RRRRRY
++
+BBBFFF
+@UID87
+NRYRRY
++
+#0<FFF
+@UID88
+NRRRRY
++
+#0<FFF
+@UID89
+NYRRRY
++
+#0<FFF
+@UID90
+NRYRRY
++
+#0<FFF
+@UID91
+RRRRYR
++
+BBBFFF
+@UID92
+RRRRYR
++
+BBBFFF
+@UID93
+YRYRRY
++
+BBBFFF
+@UID94
+RRYRYR
++
+BBBFFF
+@UID95
+RYRRYR
++
+BBBFFF
+@UID96
+NYRRRY
++
+#0<FFF
+@UID97
+NYYRRY
++
+#0<FFF
+@UID98
+YYRRRY
++
+BBBFFF
+@UID99
+YRYRYR
++
+BBBFFF
diff -r 000000000000 -r daad7189e48e test-data/extracted_bcs.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/extracted_bcs.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,400 @@
+@UID0
+NCTAAT
++
+#0<FFF
+@UID1
+GGTACG
++
+BBBFFF
+@UID2
+GCTAAT
++
+B<BF<B
+@UID3
+CATACG
++
+BBBFFF
+@UID4
+NAAAAT
++
+#0<FFF
+@UID5
+AAGAAT
++
+BBBFFF
+@UID6
+CAAACG
++
+BBBFFF
+@UID7
+NCCGGT
++
+#0<FFF
+@UID8
+ACAAAT
++
+BBBFFF
+@UID9
+AAGAAT
++
+BBBFFF
+@UID10
+NAAAAT
++
+#0<FFF
+@UID11
+NCAACG
++
+#0<FFF
+@UID12
+NGCAAT
++
+#0<FFF
+@UID13
+NAGAAT
++
+#0<FFF
+@UID14
+TCAACG
++
+BBBFFF
+@UID15
+TCGAAT
++
+BBBFFF
+@UID16
+CCCGCG
++
+BBBFFF
+@UID17
+CAAGCG
++
+BBBFFF
+@UID18
+ACCAAT
++
+BBBFFF
+@UID19
+ACTCTT
++
+BBBFFF
+@UID20
+GCAACG
++
+BBBFFF
+@UID21
+ACGACG
++
+BBBFFF
+@UID22
+TCAGCG
++
+BBBFFF
+@UID23
+TTCACG
++
+BBBFFF
+@UID24
+TAAACG
++
+BBBFFF
+@UID25
+GATTGA
++
+BBBFFF
+@UID26
+TCAGGT
++
+BBBFFB
+@UID27
+ACTAAT
++
+BBBFFF
+@UID28
+ACTACG
++
+BBBFFF
+@UID29
+NTTGCG
++
+#0<FFF
+@UID30
+NCAACG
++
+#0<FFF
+@UID31
+CACAAT
++
+BBBFFF
+@UID32
+CACAAT
++
+BBBFFF
+@UID33
+GCCAAT
++
+BBBFFF
+@UID34
+AGAACG
++
+BBBFFF
+@UID35
+NTTAAT
++
+#0<FFF
+@UID36
+CTAGCG
++
+BBBFFF
+@UID37
+GAGAAT
++
+BBBFFF
+@UID38
+AGGAAT
++
+BBBFFF
+@UID39
+CCAAAT
++
+BBBFFF
+@UID40
+CCAAAT
++
+BBBFFF
+@UID41
+ACTGCG
++
+BBBFFF
+@UID42
+TACGCG
++
+BBBFFF
+@UID43
+GGGACG
++
+BBBFFF
+@UID44
+GCAACG
++
+BBBFFF
+@UID45
+NCTAAT
++
+#0<FBF
+@UID46
+CCCGCG
++
+BBBFFF
+@UID47
+AAGAAT
++
+BBBFFF
+@UID48
+ATAAAT
++
+BBBFFF
+@UID49
+ACAGCG
++
+BBBFFF
+@UID50
+ACTGGT
++
+BBBFFF
+@UID51
+CACTGA
++
+BBBFFF
+@UID52
+GCGAAT
++
+BBBFFF
+@UID53
+GAAAAT
++
+BBBF<F
+@UID54
+CCAAAT
++
+BBBFFF
+@UID55
+ACCAAT
++
+BBBFFF
+@UID56
+TCCAAT
++
+BBBFFF
+@UID57
+NCCCAC
++
+#0BFFF
+@UID58
+NAGAAT
++
+#0<FFF
+@UID59
+NTAAAT
++
+#0<FFF
+@UID60
+GACAAT
++
+BB<FFF
+@UID61
+GGGGCG
++
+BBBFFF
+@UID62
+CAAAAT
++
+BBBFFF
+@UID63
+AAAAAT
++
+BBBFFF
+@UID64
+TAGACG
++
+BBBFFF
+@UID65
+AGCACG
++
+BBBFFF
+@UID66
+AGGAAT
++
+BBBFFF
+@UID67
+AACACG
++
+BB<FFF
+@UID68
+NCCAAT
++
+#00FFF
+@UID69
+CTGGGT
++
+BBBFFF
+@UID70
+TCTACG
++
+BBBFFF
+@UID71
+AAAACG
++
+BBBFFF
+@UID72
+GGAATA
++
+BBBFFF
+@UID73
+NACGCG
++
+#0<BFF
+@UID74
+TCTACC
++
+BBBFFF
+@UID75
+CCAAAT
++
+BBBFFF
+@UID76
+TGAAAT
++
+BBBFFF
+@UID77
+NGTACG
++
+#0<FFF
+@UID78
+NAAAAT
++
+#0<FFF
+@UID79
+TTTAAT
++
+BBBFFF
+@UID80
+ACCAAT
++
+BBBFFF
+@UID81
+NCCACG
++
+#0<FFF
+@UID82
+NAGAAT
++
+#0<FFF
+@UID83
+NTAAAT
++
+#0<FFF
+@UID84
+NAAAAT
++
+#0<FFF
+@UID85
+ACAAAT
++
+BBBFFF
+@UID86
+AAGAAT
++
+BBBFFF
+@UID87
+NATAAT
++
+#0<FFF
+@UID88
+NAAAAT
++
+#0<FFF
+@UID89
+NCGAAT
++
+#0<FFF
+@UID90
+NGCAAT
++
+#0<FFF
+@UID91
+AAGGCG
++
+BBBFFF
+@UID92
+AAGACG
++
+BBBFFF
+@UID93
+CATAAT
++
+BBBFFF
+@UID94
+AATACG
++
+BBBFFF
+@UID95
+GCAACG
++
+BBBFFF
+@UID96
+NCAAAT
++
+#0<FFF
+@UID97
+NCCAAT
++
+#0<FFF
+@UID98
+CCAAAT
++
+BBBFFF
+@UID99
+TATGCG
++
+BBBFFF
diff -r 000000000000 -r daad7189e48e test-data/merged_pcr_dupes.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/merged_pcr_dupes.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,3 @@
+chr1	10	20	readid_11	10	+
+chr1	10	20	readid_0	5	-
+chrX	20	30	readid_10	7	+
diff -r 000000000000 -r daad7189e48e test-data/merged_pcr_dupes_clnts.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/merged_pcr_dupes_clnts.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,3 @@
+chr1	9	10	readid_11	10	+
+chr1	20	21	readid_0	5	-
+chrX	19	20	readid_10	7	+
diff -r 000000000000 -r daad7189e48e test-data/merged_pcr_dupes_spurious.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/merged_pcr_dupes_spurious.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,12 @@
+chr1	10	20	AAAAA	1	+
+chr1	10	20	AAAAA	2	+
+chr1	10	20	AAAAA	3	+
+chr1	10	20	AAAAA	4	+
+chr1	10	20	AAAAA	5	+
+chr1	10	20	AAAAA	10	+
+chr1	10	20	AAAAA	100	+
+chr1	10	20	AAAAA	3	-
+chr1	10	20	AAAAA	5	-
+chr1	10	20	AAAAA	4	-
+chr1	10	20	AAAAA	5	-
+chrX	20	30	TTTTT	7	+
diff -r 000000000000 -r daad7189e48e test-data/merged_pcr_dupes_spurious_filtered_thresh05.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/merged_pcr_dupes_spurious_filtered_thresh05.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,6 @@
+chr1	10	20	AAAAA	5	-
+chr1	10	20	AAAAA	5	-
+chr1	10	20	AAAAA	4	-
+chr1	10	20	AAAAA	3	-
+chr1	10	20	AAAAA	100	+
+chrX	20	30	TTTTT	7	+
diff -r 000000000000 -r daad7189e48e test-data/pcr_dupes_randomdict.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pcr_dupes_randomdict.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,88 @@
+@readid_0
+AAAAA
++
+BBBBB
+@readid_1
+AAAAA
++
+BBBBB
+@readid_2
+AAAAA
++
+BBBBB
+@readid_3
+TTTTT
++
+BBBBB
+@readid_4
+AAAAA
++
+BBBBB
+@readid_5
+AAAAA
++
+BBBBB
+@readid_6
+AAAAA
++
+BBBBB
+@readid_7
+TTTTT
++
+BBBBB
+@readid_8
+AAAAA
++
+BBBBB
+@readid_9
+AAAAA
++
+BBBBB
+@readid_10
+TTTTT
++
+BBBBB
+@readid_11
+AAAAA
++
+BBBBB
+@readid_12
+TTTTT
++
+BBBBB
+@readid_13
+TTTTT
++
+BBBBB
+@readid_14
+TTTTT
++
+BBBBB
+@readid_15
+AAAAA
++
+BBBBB
+@readid_16
+AAAAA
++
+BBBBB
+@readid_17
+AAAAA
++
+BBBBB
+@readid_18
+AAAAA
++
+BBBBB
+@readid_19
+AAAAA
++
+BBBBB
+@readid_20
+TTTTT
++
+BBBBB
+@readid_21
+AAAAA
++
+BBBBB
diff -r 000000000000 -r daad7189e48e test-data/pcr_dupes_sorted_2.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/pcr_dupes_sorted_2.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,22 @@
+chr1	10	20	readid_0	0	-
+chr1	10	20	readid_1	0	-
+chr1	10	20	readid_2	0	+
+chr1	10	20	readid_4	0	+
+chr1	10	20	readid_5	0	+
+chr1	10	20	readid_6	0	-
+chr1	10	20	readid_8	0	+
+chr1	10	20	readid_9	0	-
+chr1	10	20	readid_11	0	+
+chr1	10	20	readid_15	0	+
+chr1	10	20	readid_16	0	+
+chr1	10	20	readid_17	0	-
+chr1	10	20	readid_18	0	+
+chr1	10	20	readid_19	0	+
+chr1	10	20	readid_21	0	+
+chrX	20	30	readid_3	0	+
+chrX	20	30	readid_7	0	+
+chrX	20	30	readid_10	0	+
+chrX	20	30	readid_12	0	+
+chrX	20	30	readid_13	0	+
+chrX	20	30	readid_14	0	+
+chrX	20	30	readid_20	0	+
diff -r 000000000000 -r daad7189e48e test-data/reads.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/reads.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,400 @@
+@UID0
+NCTAGAATGAAAAATGTATGCTTGTGACCAATAAATATAAAAGGTTAAGATCGGAAGAGCACACGTCTG
++
+#0<FFFFFFFFFFFFIFIIIIIIIIIIIIIIFIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID1
+GGTTAACGCCAGTTGGTCCGAGTGTTGTGGGTTATTGTTAACTAAGATCGGAAGAGCACACGTCTGAAC
++
+BBBFFFFFFFFFFIIIIIIIFFFFFIIIIIIFIFIIIIIIIIIFFIFIIIIFFIIFFFFFFFFFFFBFB
+@UID2
+GCTAGAATCCTTTTCCTTTGCCCTGATCTTGTCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+B<B<BF<BFBFBFIFFFBFFFFFFFBBB0BBFFFBB<FB707'BFF0FFB'<B<BF#############
+@UID3
+CATTCACGTGGTTGGCGTTTGGTGTGTGGCTGTTTCTGCTCATTGATAGATCGGAAGAGCACACGTCTG
++
+BBBFFFFFFFFFFIIFFFFIIIFFFFFFIFIIFIIFFIFIIIFIIFFIIIIFIIFFFFFFFFFFFFBFF
+@UID4
+NAAAGAATACTGCTTTAATTCTAAATAAAAATTTATTGCCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+#0<BFFFFFFFFFIIIFFFIIIIIIFIIFIFFIFIIIIFFIIIIIIIIIIIBFFFIIIIIIIIFIIIFF
+@UID5
+AAGAGAATGAATTGCCCCCCCCCACAACCGCGCTTGACTAGCGTAGCAGATCGGAAGAGCACACGTCTG
++
+BBBFFFFFFFFFFI0<FFFFIIIIIIIIIIIFBFFBFFFFBBFBFFFBBFFFFFFFFBBBFBFFFBBBB
+@UID6
+CAATCACGAAGTGTCCTCAGCAAAAGCTGGGTGGTTACCAAGATCGGAAGAGCACACGTCTGAACTCCA
++
+BBBFFFFFFFFFFFFFFIIIIIIFIFFIIIIBFFIIBBF<BFBFFFFIII7FFBFFFFBBBBBBFF<<B
+@UID7
+NCCGAGGTCAATCCGAGTCACGGCACCATGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+#0<FFFFFFFFFFIIIIFIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFB
+@UID8
+ACAAGAATAATCTTCTTAGAGTTGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCA
++
+BBBFFFFFFFFFFIFFIFI<FBFFFFFFIFFFFFFIBFBBFBFFFFFFFFFBBFFBB<7BFFF<<BB<B
+@UID9
+AAGAGAATAAATTACTAAATGTGAAAAACATAATGTGAAAAGTATATGTTAAGATCGGAAGAGCACACG
++
+BBBFFFFFFFFFFIIIIFFIIFIFIFIIFIFFFFF<BFFFIFFFIIIIIIIFFFIFF<FIIIIIIFIFI
+@UID10
+NAAAGAATGACCTGTAAACAGAATGATAAAGTATGTGATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIFFIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID11
+NCATCACGAATCTGGCATGTTGGAACAATGTAGGTAAGGGAAGTGGTCAAGATCGGAAGAGCACACGTC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIFFIIIIFIIFIIFIIIIIIIIIIIIIIIIIIFFFF
+@UID12
+NGCAGAATGGTTAAAGTTGCTGGCCTAATGAGCTCAGTAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIFIIIIIIIF
+@UID13
+NAGAGAATAAATCTTATTAATTCCCTCAAAAGCAGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID14
+TCATCACGATGGTTGAGATCGGAAAGCACTATCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIFIIIIFIFIIIIIIIIIIIIIIIFFIFFFFFFFFFBFBFF
+@UID15
+TCGAGAATCCGCCCACCCAGGGACGCCAGGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFBFFFFFF
+@UID16
+CCCCGGCGCTATCCCACTTCTGACACCAGGCTGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFB
+@UID17
+CAACGGCGGAGTCTGTTATAAAATGGACTTAACTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIFIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFBFF
+@UID18
+ACCAGAATCCGTCTAGAATGTAGCATCTAGTGACTTTTTAAAGCCGGTTAAGATCGGAAGAGCACACGT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID19
+ACTTTCTTTCCTGGGCATCCCGGGGATCCCAGAGCCGGCCCAGGTAGGCCAAGATCGGAAGAGCACACG
++
+BBBFFFFFFFFFFIIFIIIIIIFFFIIIIIIIIIIIFIIBFIIIBBFFFFBFFFFFFFFFBFFBFBFFF
+@UID20
+GCATCACGGGCCTACTCTACTTTCTGTTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+BBBFFFFFFFFFFFIFIIIIIIIIIFFFFIIIIIIFFFIFF<FFFFFFF<<BFBFFFFFFFBBBBFFFB
+@UID21
+ACGTCACGCAATCATGTTCGCCTCACACGCGAAAGGTCTTGACAGATCGGAAGAGCACACGTCTGAACT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFF
+@UID22
+TCACGGCGCGACCTAATGTTAGTGAATATTACTAAGTTGCCAAGATCGGAAGAGCACACGTCTGAACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFF
+@UID23
+TTCTCACGGCGGTGGGATCCCGCGGCCACAGTAGATCGGGAGAGCACACGTCTGAACTCCAGTCACTCA
++
+BBBFFFFFFFFFBBFF0'BFF<B<B7B0<7<'<<7'7<<'0'7<7<<<7'0<<<<B'0<B'00<BB###
+@UID24
+TAATCACGCGACAGCCACAAGTTACATAGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATT
++
+BBBFFFFFFFFFFIIIIIIIIFIIIIIIIIIIIIIIIIIIIIFIIIIIIIFFFFFFFFFFFFFFFFFFF
+@UID25
+GATGCTGAACTCTCGTAGTTGGATCTTGGGAGCGGGCGGGCGGCTAAGATCGGAAGAGCACACGTCTGA
++
+BBBFFFFFFFBFFFFIIFFIFFFFFFFFFFFFFFFFBBF7BB7BFFBFBBBBBFB7BBFFFFBFB<BBB
+@UID26
+TCAGAGGTCATGCACGCGTGGGTTCGTATCCCGATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+BBBFFFFBFFFFFIIIIIFFIIBFFFFFIIIIIIIIFFIIFFFFFFFFFFFFFFBFBBBBBBBBBBBBB
+@UID27
+ACTAGAATCCATTTTTTCTTACTAGCTTGCAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFBFFFFFFFFIIIIFFIFFBFFFFFFFIBBFFIFFFBFFFIFF<<BBBFFFFFIFFFFFFF<
+@UID28
+ACTTCACGAATTCCCCGCGCAGGTTCGTATCCTGCCGACTACGCCAGTAACAGATCGGAAGAGCACACG
++
+BBBFFFFFFFFBFFFIIFIFIIIFFFIIFIFFIFFFF<BFFBFFFFB<BBFFBFB<B<7B<7BBBB<B<
+@UID29
+NTTCGGCGACGTTGCACCGTGGAGGCCACAGGAGCAGAAAGATTAAGATCGGAAGAGCACACGTCTGAA
++
+#0<FFFFFFFFFFIIIIIIFIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID30
+NCATCACGCATTTTAATCTAGGGTTATGAAACTCAATTGCTGAGATCGGAAGAGCACACGTCTGAACTC
++
+#0<FFFFFFFFFFIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID31
+CACAGAATCACAGTAAGTGTTAATAGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCAT
++
+BBBFFFFFFFFFFFFIIFFIFIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIB
+@UID32
+CACAGAATCCAGTTTTTAAAGTTAAACTGAAGACAATATCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+BBBFFFFFFFFFFIFFIIIIIIFIIIFIFIFFFFIIIIIIIIFIIIIIFFFFIIIIIIFFIIFIFIIII
+@UID33
+GCCAGAATAATTAAACAGGTATTTTGATATTGGAGTATCAAGATCGGAAGAGCACACGTCTGAACTCCA
++
+BBBFFFFFFFFFFIIIIIIIFFIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID34
+AGATCACGTAAGTATCCCCGCCTGTCACGCGTTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+BBBFFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFF
+@UID35
+NTTAGAATACATAATCTTAAACAGCGGATGATAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIFFFFFFFFFF
+@UID36
+CTACGGCGACATCCGGGTGCCCCCTCCAGGATTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBFFFFFFFFFFIIFFBFFIIIIIIIIIIIIIIIIIIIFFFFFFFFFBFBBFFFFFFFFFFFFFFFFF
+@UID37
+GAGAGAATACTACAAACATGAAAATAAAACGGTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIIIFIIIIIFIFIFFBFFFF
+@UID38
+AGGAGAATAATCCAGAAAACACTGTACTTACCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIFFFFFFF
+@UID39
+CCAAGAATCCCTGTAATTCTTTTATAAATCCCTTAGTGTACTTTACTAAGATCGGAAGAGCACACGTCT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFIIFFIIIFFIIIIIIIIIIIFIIIIIIIIFII
+@UID40
+CCAAGAATACAGTCAAGACGGAAGAGCACAACTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBFFFFFFFFFFIIIIFFIIIIIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIFFFFFFFFF
+@UID41
+ACTCGGCGACAAAAAAACCCATTGCTAATATAGCAATAAAACGAGATCGGAAGAGCACACGTCTGAACT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID42
+TACCGGCGATTGAGCCCTGCAGCACCCCAACATAATATGGTAGATCGGAAGAGCACACGTCTGAACTCC
++
+BBBFFFFFBFFFFIIIIIIIIIIIIIIIIFFFFFFFFFIIIIIIIIFFFFFFFFFFFFFFFFFBF<BBF
+@UID43
+GGGTCACGGAAGAGGACGACCATCCCCGATAGAGGAGGATCAACTAAGATCGGAAGAGCACACGTCTGA
++
+BBBFFFFFFFFFFIFIIIIIIIIIIIIIFFIFFIIIFFFFFFFFFFFFFFFFFBFFFFFFFFFFFFFFF
+@UID44
+GCATCACGCATACATTTGGACAGTCAAGGTGCACGGTCAAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+BBBFFFFFFFFFFIIIIFIIIFFFFIIIIFIIIIIIFFFIIIIIIIIIIIFIFIIIFFFFFFFFFFFFF
+@UID45
+NCTAGAATACGCCAGGGATTGTGGGTTCGTGTCCCACCTGGGGTACCACGGTAGATCGGAAGAGCACAC
++
+#0<BFFBFFFFBFFFFFFFFFBFFFBBFBBFFFFFFFFFFFII<FBFFFFFBB77<BBFBB7BBBBBBB
+@UID46
+CCCCGGCGCCTCTCCAGTCTCTTCGGGGGCGTGGGTTTGACAGATCGGAAGAGCACACGTCTGAACTCC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIFFBFFFBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID47
+AAGAGAATGCCCTACTATACTTTGGTTGGCAGTGCCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID48
+ATAAGAATAATGGTAGAATTCTAAAATGTGTAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBFFFFFFFFFFFFIFFIIFIFFFFFFFFIIIIIFBIIIIIFIIIFFIIFFFFFFBFFFIFIIIIFII
+@UID49
+ACACGGCGCGTGAGAACTCCTCGTCTGGAAGAGTAATAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIFIIIBFIIIFIIFFFIIIIIFFFFFFFFFFFBFFFBFBBFFFFFF
+@UID50
+ACTGAGGTTCTCGCCCGGGCGAATGGCAGACGATGATCCCTACGGCCGAGATCGGAAGAGCACACGTCT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID51
+CACGCTGAGGTGCAACCTCTGCCTCCCGGGTTCAAGCATGTCGAGATCGGAAAAGCACACGTCTGAACT
++
+BBBFFFFFFFFFFIIIFIIIIIIIIIIIIIBFIIIIIIIIIIIFIIIIIIIIIIFFFFFFFFFFFFFFF
+@UID52
+GCGAGAATAATTTCTTTCTAGCATGCTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIFIIIIFIIIIIIFIIIIIIIIIIIIIIIIIIIFFF
+@UID53
+GAAAGAATACACATTTGTCAATAATACGTCAAAAAAAGATTAAGATCGGAAGAGCACACGTCTGAACTC
++
+BBBFFF<FFFFFFIIIIFFIIIIIIIFIIIIIIIIFIIIIIFIIIIIIIFFIIIIIIIFFFFFFBFFFF
+@UID54
+CCAAGAATACCATATTTACTTCTTATACATTTGATAGGTTAAGATCGGAAGAGCACACGTCTGAACTCC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIIIIIFFIIIIIIII
+@UID55
+ACCAGAATAAGACCAAGATCGGAAGAGCTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFFIIIIFIIIIIFIIIIIFIIIIIFFIIIIIIIFIIIIFFFFFFFFFFFFFBFFFBF
+@UID56
+TCCAGAATACAAATAATAGATAAGAACCAATCAATAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+BBBFFFFFFBBFFFIIFIFFBFIBFIFIFBFFFIIIBBFFBFFFIFBFFIFFFFFFFFBFBBFFFBFBF
+@UID57
+NCCGTCACGCATTTGGATTCGAGATTTTGGGATTTGGGTTGCTCATTCAATAGATCGGAAGAGCACACG
++
+#0BFFFFFFFFFFIIIIIIIIFFFIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIFFFFFFFFF
+@UID58
+NAGAGAATAATTTTAAAACTAATAAATAGTGACTGTAATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIII
+@UID59
+NTAAGAATCCATCTCGGTGGAACCTCCAGGCTGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+#0<FFFFFFFFFFIIIIBFFFIIIIIIIIFIIIBFBFFFFIIIIIIIIFIIFFIIIIIIIFFFFFFFFF
+@UID60
+GACAGAATCGTAATAAATACCCAGCACTAACTAAGCCTATCGAGATCGGAAGAGCACACGTCTGAACTC
++
+BB<FFFFFFFFFFFFIIFIFFFFBFFFFFFFBB<FFF<FBFFFFFBFB<FFIIFFFIIIIFFFFFFFFF
+@UID61
+GGGCGGCGAGAAAAACCCAAAAAGAAATATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+BBBFFFFFFFFFFIFFBFFIIFFIIIIFFFFFFFFFFFFFBBBBBBBFFFBBBFBBFF<BBBBBBBFBB
+@UID62
+CAAAGAATAATAACAAAACATTATTTGCAAAACATTGATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIFBFFFFFFFFIIFFIFFFFFFFIFFFFFFFFFFFIIFBIIII
+@UID63
+AAAAGAATGATTCAGGCTCAGTCCCCTCCCGATGATTAAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+BBBFFFFFFFFFFIIFIIIFIIIIFIFFIIIIIFFIBIFIIIIIIIIIFIIIIIFIFFFIFFFFFFFFB
+@UID64
+TAGTCACGAGAAGTACACATTTAAAATCTAAGTTCCATTTTTTCGTCTGTCAAGATCGGAAGAGCACAC
++
+BBBFFFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFF
+@UID65
+AGCTCACGTATCCTTCCACCTCAGCCTCACGAGTAGCTGGTGGTTAAGATCGGAAGAGCACACGTCTGA
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIFIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID66
+AGGAGAATGATTTTGACAATGGACCAAAGACGATTATTGAGATCGGAAGAGCACACGTCGAACTCCAGT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFF
+@UID67
+AACTCACGGGCCCCCTGGAGGCTGGACCCCATTCCCAGTTGGCAGATCGGAAGAGCACACGTCTGAACT
++
+BB<FFFFFFFBFFFI<FIFFIFIFFFBFBFIIIFFF<BBB7BBFF<FBB<BBBFFBBFBBFBBFF7BBB
+@UID68
+NCCAGAATAAATCCCAGCGGTGCCTCCATATCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+#00<BFFFFF<F0<BB<7B<<BFFFBBFF<<<B'BBF<BB<FF7<BB<<'07<''00<B<<BB<00770
+@UID69
+CTGGAGGTATAAGAGCTTTAGCCGGGCCGTGATAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFB
+@UID70
+TCTTCACGTACTCCACTCCCAGCTGTTCAGCCTGGTTGAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID71
+AAATCACGCAAAAATCTGTTTCCATATTGCTTCAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFF
+@UID72
+GGAGAATAAATTTCCAACAGCTTATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCAT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFIIIF
+@UID73
+NACCGGCGCATCCCCGCGAGGGGCCCGGGGCTCAGCAGATCGTAATAGATCGGAAGAGCACACGTCTGA
++
+#0<FFBFFFFFFFIIIIIFFBFFFFFFFFFFFFFBFFFFFFFFFFFFFBBBFFBBBFBFFFBFFBBFBB
+@UID74
+TCTTCACCGCATCCAGATTCTGTAGAGCTGGACAGGCAGGTGGTAGATCGGAAGAGCACACGTCTGAAC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIBFFFFFIIIIIIIFIIIIIIFFFFFFFFF
+@UID75
+CCAAGAATCCTTAACATTTTATTATTTATTGTAAAATCAATTTGTGACAGATCGGAAGAGCACACGTCT
++
+BBBFFFFFFFFFFIIFFFIIIIIIFIIFFIIFFIIFBFF<FIBFBBBFB7BFFFB<FFIFIIFFFFIBF
+@UID76
+TGAAGAATGATTTTAACCTGCAATTTGCCTCAGCAACAGGTTAAGATCGGAAGAGCACACGTCTGAACT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID77
+NGTTCACGGGCTCCTCCCCTTCGTCTTCACAGAGCTGGTCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIFFFFFFFFFFFFF
+@UID78
+NAAAGAATGATTTATAACCACAACAAGAAGGGAATCAATGTTAAGATCGGAAGAGCACACGTCTGAACT
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID79
+TTTAGAATCAATGACCAATTACCAAGGAAGCTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFF
+@UID80
+ACCAGAATCCTCACCCTAACACCCTTCAATACAATGTTTACAGGGTTAAGATCGGAAGAGCACACGTCT
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIFFFF
+@UID81
+NCCTCACGACTTTGTTAAGTTGATTTAACATTGCCCCCCCCACAACTTGGCAGATCGGAAGAGCACACG
++
+#0<FFFFFFFFFFIIIIIIFIIIIIIIIIIIFFFIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID82
+NAGAGAATGATTTTAACTAAGGTGGTGACACTGGTAACAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+#0<FFFFFFFFFFIIIIIIIIIFFFFFFIIIIIIFFIIIIIIIIIIIFFFIIIIIFFFIIIIIIIIFFF
+@UID83
+NTAAGAATAAAATAAGTGTGATCAATCTTTTTGTCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+#0<FFFFFFFFFFIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFF
+@UID84
+NAAAGAATAATTTGACAACTTTAAATGAGTGTAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+#0<FFFFFFFFFFIIIIIIIFIIFIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID85
+ACAAGAATAAGTTATTCCTTTAAAAGAAATTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBBBBFFFFFFFBFFIFFFIFFFFFFBFFFFFIIFFFIFFIBFBBB<FFBBB<BFFFFB<BBFBBFFFB
+@UID86
+AAGAGAATAAGCTAGAGCAGGCCTGCTTGCGGCTCCCTGTGGTTAAGATCGGAAGAGCACACGTTGAAC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIFFFFFFFFFFFFFF
+@UID87
+NATAGAATAACTTAAATTGATTAAAACTTTGTAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIFIIFIIIIIIIIIFIIIIIIIIIIIIFIIIIIIIIIIF
+@UID88
+NAAAGAATCGATCTCGGTGGAACCTCCAGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFB
+@UID89
+NCGAGAATCCCTATTATTTTGTTGCTAAAAATTTAATAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIF
+@UID90
+NGCAGAATAATGTAATACCTCCAGATATGTAATTATTTGTAACAGATCGGAAGAGCACACGTCTGAACT
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID91
+AAGCGGCGAGATCCCACTCCTGACACCATCGGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIFFIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFB
+@UID92
+AAGTCACGGATCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGTGCTAAGATCGGAAGAGCACACG
++
+BBBFFFFFFFFFFIIIIIIFIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFBFFFFBFFFFFFFF
+@UID93
+CATAGAATAATTGATAGATCGGAAGAGCACACGTCTTGTAGCAGATCGGAAGAGCACACGTCTGAACTC
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIFFFFFFF
+@UID94
+AATTCACGCAAAATAGAAGTCCTGATAGCCCAGACAACTGAGATCGGAAGAGCACACGTCTGAACTCCA
++
+BBBFFFFFFFFFFFFFFFFBFFIFFFFFFBBFFFFFIBFFBFFFFIFFFBBBFBFBFFBFFF<B<FFFF
+@UID95
+GCATCACGACTCATCCTCTATGGACATTGGAGTTTTTCATTATGTGCATGGATCAAGATCGGAAGAGCA
++
+BBBFFFFFFFFFFIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID96
+NCAAGAATCACTGAATCCAGCGATCCGAGTTCTGTCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIIFIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFF
+@UID97
+NCCAGAATCCTTTCACATGTATAAGGTGGCTCCATAGCTTACTGAGATCGGAAGACACACGTCTGAACT
++
+#0<FFFFFFFFFFIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIII
+@UID98
+CCAAGAATCCGCGTGGTGGCATGTGCCTGTGGTCCCAGATAGCAGATCGGAAGAGCACACGTCTGAACT
++
+BBBFFFFFFFFFFFIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFF
+@UID99
+TATCGGCGAGAACGCATAAATCTTTCGCCTTTTACTAAAGATTTGTGGCAGATCGGAAGAGCACACGTC
++
+BBBFFFFFFFFFFIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIFFFFFFFFBFFFBFFBFFFFF
diff -r 000000000000 -r daad7189e48e test-data/readswithtail.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/readswithtail.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,8 @@
+@longenough
+ACTAGAAT
++
+F0<FFFFF
+@tooshort
+GGTTAAC
++
+AAAFFFF
diff -r 000000000000 -r daad7189e48e test-data/readswithtailremoved.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/readswithtailremoved.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,8 @@
+@longenough
+A
++
+F
+@tooshort
+N
++
+B
diff -r 000000000000 -r daad7189e48e test-data/result_original_head.fastq
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/result_original_head.fastq	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,400 @@
+@UID0
+AGGAAAAATGTATGCTTGTGACCAATAAATATAAAAGGTTAAGATCGGAAGAGCACACGTCTG
++
+FFFFFFFFFIFIIIIIIIIIIIIIIFIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID1
+TACCAGTTGGTCCGAGTGTTGTGGGTTATTGTTAACTAAGATCGGAAGAGCACACGTCTGAAC
++
+FFFFFFFIIIIIIIFFFFFIIIIIIFIFIIIIIIIIIFFIFIIIIFFIIFFFFFFFFFFFBFB
+@UID2
+AGCCTTTTCCTTTGCCCTGATCTTGTCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+<BFBFBFIFFFBFFFFFFFBBB0BBFFFBB<FB707'BFF0FFB'<B<BF#############
+@UID3
+TCTGGTTGGCGTTTGGTGTGTGGCTGTTTCTGCTCATTGATAGATCGGAAGAGCACACGTCTG
++
+FFFFFFFIIFFFFIIIFFFFFFIFIIFIIFFIFIIIFIIFFIIIIFIIFFFFFFFFFFFFBFF
+@UID4
+AGACTGCTTTAATTCTAAATAAAAATTTATTGCCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+BFFFFFFIIIFFFIIIIIIFIIFIFFIFIIIIFFIIIIIIIIIIIBFFFIIIIIIIIFIIIFF
+@UID5
+AGGAATTGCCCCCCCCCACAACCGCGCTTGACTAGCGTAGCAGATCGGAAGAGCACACGTCTG
++
+FFFFFFFI0<FFFFIIIIIIIIIIIFBFFBFFFFBBFBFFFBBFFFFFFFFBBBFBFFFBBBB
+@UID6
+TCAAGTGTCCTCAGCAAAAGCTGGGTGGTTACCAAGATCGGAAGAGCACACGTCTGAACTCCA
++
+FFFFFFFFFFFIIIIIIFIFFIIIIBFFIIBBF<BFBFFFFIII7FFBFFFFBBBBBBFF<<B
+@UID7
+GACAATCCGAGTCACGGCACCATGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIFIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFB
+@UID8
+AGAATCTTCTTAGAGTTGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCA
++
+FFFFFFFIFFIFI<FBFFFFFFIFFFFFFIBFBBFBFFFFFFFFFBBFFBB<7BFFF<<BB<B
+@UID9
+AGAAATTACTAAATGTGAAAAACATAATGTGAAAAGTATATGTTAAGATCGGAAGAGCACACG
++
+FFFFFFFIIIIFFIIFIFIFIIFIFFFFF<BFFFIFFFIIIIIIIFFFIFF<FIIIIIIFIFI
+@UID10
+AGGACCTGTAAACAGAATGATAAAGTATGTGATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+FFFFFFFIIIIIIIIIIIIIIIIIIFFIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID11
+TCAATCTGGCATGTTGGAACAATGTAGGTAAGGGAAGTGGTCAAGATCGGAAGAGCACACGTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIFFIIIIFIIFIIFIIIIIIIIIIIIIIIIIIFFFF
+@UID12
+AGGGTTAAAGTTGCTGGCCTAATGAGCTCAGTAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+FFFFFFFIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIFIIIIIIIF
+@UID13
+AGAAATCTTATTAATTCCCTCAAAAGCAGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID14
+TCATGGTTGAGATCGGAAAGCACTATCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+FFFFFFFIIIIIIIIIIIIIIIIFIIIIFIFIIIIIIIIIIIIIIIFFIFFFFFFFFFBFBFF
+@UID15
+AGCCGCCCACCCAGGGACGCCAGGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFBFFFFFF
+@UID16
+CGCTATCCCACTTCTGACACCAGGCTGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFB
+@UID17
+CGGAGTCTGTTATAAAATGGACTTAACTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+FFFFFFFIIIIIIIIIIIIIFIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFBFF
+@UID18
+AGCCGTCTAGAATGTAGCATCTAGTGACTTTTTAAAGCCGGTTAAGATCGGAAGAGCACACGT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID19
+TTTCCTGGGCATCCCGGGGATCCCAGAGCCGGCCCAGGTAGGCCAAGATCGGAAGAGCACACG
++
+FFFFFFFIIFIIIIIIFFFIIIIIIIIIIIFIIBFIIIBBFFFFBFFFFFFFFFBFFBFBFFF
+@UID20
+TCGGCCTACTCTACTTTCTGTTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+FFFFFFFFIFIIIIIIIIIFFFFIIIIIIFFFIFF<FFFFFFF<<BFBFFFFFFFBBBBFFFB
+@UID21
+TCCAATCATGTTCGCCTCACACGCGAAAGGTCTTGACAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFF
+@UID22
+CGCGACCTAATGTTAGTGAATATTACTAAGTTGCCAAGATCGGAAGAGCACACGTCTGAACTC
++
+FFFFFFFIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFF
+@UID23
+TCGCGGTGGGATCCCGCGGCCACAGTAGATCGGGAGAGCACACGTCTGAACTCCAGTCACTCA
++
+FFFFFFBBFF0'BFF<B<B7B0<7<'<<7'7<<'0'7<7<<<7'0<<<<B'0<B'00<BB###
+@UID24
+TCCGACAGCCACAAGTTACATAGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATT
++
+FFFFFFFIIIIIIIIFIIIIIIIIIIIIIIIIIIIIFIIIIIIIFFFFFFFFFFFFFFFFFFF
+@UID25
+GCACTCTCGTAGTTGGATCTTGGGAGCGGGCGGGCGGCTAAGATCGGAAGAGCACACGTCTGA
++
+FFFFBFFFFIIFFIFFFFFFFFFFFFFFFFBBF7BB7BFFBFBBBBBFB7BBFFFFBFB<BBB
+@UID26
+GACATGCACGCGTGGGTTCGTATCCCGATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+FFFFFFFIIIIIFFIIBFFFFFIIIIIIIIFFIIFFFFFFFFFFFFFFBFBBBBBBBBBBBBB
+@UID27
+AGCCATTTTTTCTTACTAGCTTGCAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFBFFFFFFFFIIIIFFIFFBFFFFFFFIBBFFIFFFBFFFIFF<<BBBFFFFFIFFFFFFF<
+@UID28
+TCAATTCCCCGCGCAGGTTCGTATCCTGCCGACTACGCCAGTAACAGATCGGAAGAGCACACG
++
+FFFFFBFFFIIFIFIIIFFFIIFIFFIFFFF<BFFBFFFFB<BBFFBFB<B<7B<7BBBB<B<
+@UID29
+CGACGTTGCACCGTGGAGGCCACAGGAGCAGAAAGATTAAGATCGGAAGAGCACACGTCTGAA
++
+FFFFFFFIIIIIIFIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID30
+TCCATTTTAATCTAGGGTTATGAAACTCAATTGCTGAGATCGGAAGAGCACACGTCTGAACTC
++
+FFFFFFFIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID31
+AGCACAGTAAGTGTTAATAGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCAT
++
+FFFFFFFFFIIFFIFIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIB
+@UID32
+AGCCAGTTTTTAAAGTTAAACTGAAGACAATATCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+FFFFFFFIFFIIIIIIFIIIFIFIFFFFIIIIIIIIFIIIIIFFFFIIIIIIFFIIFIFIIII
+@UID33
+AGAATTAAACAGGTATTTTGATATTGGAGTATCAAGATCGGAAGAGCACACGTCTGAACTCCA
++
+FFFFFFFIIIIIIIFFIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID34
+TCTAAGTATCCCCGCCTGTCACGCGTTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+FFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFF
+@UID35
+AGACATAATCTTAAACAGCGGATGATAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIFFFFFFFFFF
+@UID36
+CGACATCCGGGTGCCCCCTCCAGGATTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+FFFFFFFIIFFBFFIIIIIIIIIIIIIIIIIIIFFFFFFFFFBFBBFFFFFFFFFFFFFFFFF
+@UID37
+AGACTACAAACATGAAAATAAAACGGTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIIIFIIIIIFIFIFFBFFFF
+@UID38
+AGAATCCAGAAAACACTGTACTTACCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIFFFFFFF
+@UID39
+AGCCCTGTAATTCTTTTATAAATCCCTTAGTGTACTTTACTAAGATCGGAAGAGCACACGTCT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIFIIFFIIIFFIIIIIIIIIIIFIIIIIIIIFII
+@UID40
+AGACAGTCAAGACGGAAGAGCACAACTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+FFFFFFFIIIIFFIIIIIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIFFFFFFFFF
+@UID41
+CGACAAAAAAACCCATTGCTAATATAGCAATAAAACGAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID42
+CGATTGAGCCCTGCAGCACCCCAACATAATATGGTAGATCGGAAGAGCACACGTCTGAACTCC
++
+FFBFFFFIIIIIIIIIIIIIIIIFFFFFFFFFIIIIIIIIFFFFFFFFFFFFFFFFFBF<BBF
+@UID43
+TCGAAGAGGACGACCATCCCCGATAGAGGAGGATCAACTAAGATCGGAAGAGCACACGTCTGA
++
+FFFFFFFIFIIIIIIIIIIIIIFFIFFIIIFFFFFFFFFFFFFFFFFBFFFFFFFFFFFFFFF
+@UID44
+TCCATACATTTGGACAGTCAAGGTGCACGGTCAAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+FFFFFFFIIIIFIIIFFFFIIIIFIIIIIIFFFIIIIIIIIIIIFIFIIIFFFFFFFFFFFFF
+@UID45
+AGACGCCAGGGATTGTGGGTTCGTGTCCCACCTGGGGTACCACGGTAGATCGGAAGAGCACAC
++
+BFFFFBFFFFFFFFFBFFFBBFBBFFFFFFFFFFFII<FBFFFFFBB77<BBFBB7BBBBBBB
+@UID46
+CGCCTCTCCAGTCTCTTCGGGGGCGTGGGTTTGACAGATCGGAAGAGCACACGTCTGAACTCC
++
+FFFFFFFIIIIIIIIIIIIIIIIFFBFFFBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID47
+AGGCCCTACTATACTTTGGTTGGCAGTGCCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID48
+AGAATGGTAGAATTCTAAAATGTGTAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+FFFFFFFFFIFFIIFIFFFFFFFFIIIIIFBIIIIIFIIIFFIIFFFFFFBFFFIFIIIIFII
+@UID49
+CGCGTGAGAACTCCTCGTCTGGAAGAGTAATAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+FFFFFFFIIIIIIIIIIIFIIIBFIIIFIIFFFIIIIIFFFFFFFFFFFBFFFBFBBFFFFFF
+@UID50
+GATCTCGCCCGGGCGAATGGCAGACGATGATCCCTACGGCCGAGATCGGAAGAGCACACGTCT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID51
+GCGGTGCAACCTCTGCCTCCCGGGTTCAAGCATGTCGAGATCGGAAAAGCACACGTCTGAACT
++
+FFFFFFFIIIFIIIIIIIIIIIIIBFIIIIIIIIIIIFIIIIIIIIIIFFFFFFFFFFFFFFF
+@UID52
+AGAATTTCTTTCTAGCATGCTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIFIIIIFIIIIIIFIIIIIIIIIIIIIIIIIIIFFF
+@UID53
+AGACACATTTGTCAATAATACGTCAAAAAAAGATTAAGATCGGAAGAGCACACGTCTGAACTC
++
+FFFFFFFIIIIFFIIIIIIIFIIIIIIIIFIIIIIFIIIIIIIFFIIIIIIIFFFFFFBFFFF
+@UID54
+AGACCATATTTACTTCTTATACATTTGATAGGTTAAGATCGGAAGAGCACACGTCTGAACTCC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIFFFFIIIIIIIIIIIIIIIIIIFFIIIIIIII
+@UID55
+AGAAGACCAAGATCGGAAGAGCTGCCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFFIIIIFIIIIIFIIIIIFIIIIIFFIIIIIIIFIIIIFFFFFFFFFFFFFBFFFBF
+@UID56
+AGACAAATAATAGATAAGAACCAATCAATAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+FFFBBFFFIIFIFFBFIBFIFIFBFFFIIIBBFFBFFFIFBFFIFFFFFFFFBFBBFFFBFBF
+@UID57
+GTGCATTTGGATTCGAGATTTTGGGATTTGGGTTGCTCATTCAATAGATCGGAAGAGCACACG
++
+FFFFFFFIIIIIIIIFFFIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIFFFFFFFFF
+@UID58
+AGAATTTTAAAACTAATAAATAGTGACTGTAATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+FFFFFFFIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIII
+@UID59
+AGCCATCTCGGTGGAACCTCCAGGCTGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIBFFFIIIIIIIIFIIIBFBFFFFIIIIIIIIFIIFFIIIIIIIFFFFFFFFF
+@UID60
+AGCGTAATAAATACCCAGCACTAACTAAGCCTATCGAGATCGGAAGAGCACACGTCTGAACTC
++
+FFFFFFFFFIIFIFFFFBFFFFFFFBB<FFF<FBFFFFFBFB<FFIIFFFIIIIFFFFFFFFF
+@UID61
+CGAGAAAAACCCAAAAAGAAATATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCA
++
+FFFFFFFIFFBFFIIFFIIIIFFFFFFFFFFFFFBBBBBBBFFFBBBFBBFF<BBBBBBBFBB
+@UID62
+AGAATAACAAAACATTATTTGCAAAACATTGATAGATCGGAAGAGCACACGTCTGAACTCCAG
++
+FFFFFFFIIIIIIIIIIIIIIFBFFFFFFFFIIFFIFFFFFFFIFFFFFFFFFFFIIFBIIII
+@UID63
+AGGATTCAGGCTCAGTCCCCTCCCGATGATTAAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+FFFFFFFIIFIIIFIIIIFIFFIIIIIFFIBIFIIIIIIIIIFIIIIIFIFFFIFFFFFFFFB
+@UID64
+TCAGAAGTACACATTTAAAATCTAAGTTCCATTTTTTCGTCTGTCAAGATCGGAAGAGCACAC
++
+FFFFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFF
+@UID65
+TCTATCCTTCCACCTCAGCCTCACGAGTAGCTGGTGGTTAAGATCGGAAGAGCACACGTCTGA
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIFIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID66
+AGGATTTTGACAATGGACCAAAGACGATTATTGAGATCGGAAGAGCACACGTCGAACTCCAGT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFF
+@UID67
+TCGGCCCCCTGGAGGCTGGACCCCATTCCCAGTTGGCAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFBFFFI<FIFFIFIFFFBFBFIIIFFF<BBB7BBFF<FBB<BBBFFBBFBBFBBFF7BBB
+@UID68
+AGAAATCCCAGCGGTGCCTCCATATCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+<BFF<F0<BB<7B<<BFFFBBFF<<<B'BBF<BB<FF7<BB<<'07<''00<B<<BB<00770
+@UID69
+GAATAAGAGCTTTAGCCGGGCCGTGATAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFFFFFB
+@UID70
+TCTACTCCACTCCCAGCTGTTCAGCCTGGTTGAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFF
+@UID71
+TCCAAAAATCTGTTTCCATATTGCTTCAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFF
+@UID72
+GAAATTTCCAACAGCTTATCAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTCATTCAT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFIIIF
+@UID73
+CGCATCCCCGCGAGGGGCCCGGGGCTCAGCAGATCGTAATAGATCGGAAGAGCACACGTCTGA
++
+FFFFFFFIIIIIFFBFFFFFFFFFFFFFBFFFFFFFFFFFFFBBBFFBBBFBFFFBFFBBFBB
+@UID74
+TCGCATCCAGATTCTGTAGAGCTGGACAGGCAGGTGGTAGATCGGAAGAGCACACGTCTGAAC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIBFFFFFIIIIIIIFIIIIIIFFFFFFFFF
+@UID75
+AGCCTTAACATTTTATTATTTATTGTAAAATCAATTTGTGACAGATCGGAAGAGCACACGTCT
++
+FFFFFFFIIFFFIIIIIIFIIFFIIFFIIFBFF<FIBFBBBFB7BFFFB<FFIFIIFFFFIBF
+@UID76
+AGGATTTTAACCTGCAATTTGCCTCAGCAACAGGTTAAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID77
+TCGGCTCCTCCCCTTCGTCTTCACAGAGCTGGTCAAGATCGGAAGAGCACACGTCTGAACTCC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIFFFFFFFFFFFFF
+@UID78
+AGGATTTATAACCACAACAAGAAGGGAATCAATGTTAAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFF
+@UID79
+AGCAATGACCAATTACCAAGGAAGCTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFF
+@UID80
+AGCCTCACCCTAACACCCTTCAATACAATGTTTACAGGGTTAAGATCGGAAGAGCACACGTCT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIFFFF
+@UID81
+TCACTTTGTTAAGTTGATTTAACATTGCCCCCCCCACAACTTGGCAGATCGGAAGAGCACACG
++
+FFFFFFFIIIIIIFIIIIIIIIIIIFFFIIIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFF
+@UID82
+AGGATTTTAACTAAGGTGGTGACACTGGTAACAGATCGGAAGAGCACACGTCTGAACTCCAGT
++
+FFFFFFFIIIIIIIIIFFFFFFIIIIIIFFIIIIIIIIIIIFFFIIIIIFFFIIIIIIIIFFF
+@UID83
+AGAAAATAAGTGTGATCAATCTTTTTGTCGAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
++
+FFFFFFFIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFF
+@UID84
+AGAATTTGACAACTTTAAATGAGTGTAACAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+FFFFFFFIIIIIIIFIIFIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID85
+AGAAGTTATTCCTTTAAAAGAAATTGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACT
++
+BBFFFFBFFIFFFIFFFFFFBFFFFFIIFFFIFFIBFBBB<FFBBB<BFFFFB<BBFBBFFFB
+@UID86
+AGAAGCTAGAGCAGGCCTGCTTGCGGCTCCCTGTGGTTAAGATCGGAAGAGCACACGTTGAAC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIFFFFFFFFFFFFFF
+@UID87
+AGAACTTAAATTGATTAAAACTTTGTAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIFIIFIIIIIIIIIFIIIIIIIIIIIIFIIIIIIIIIIF
+@UID88
+AGCGATCTCGGTGGAACCTCCAGGTTAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFB
+@UID89
+AGCCCTATTATTTTGTTGCTAAAAATTTAATAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIF
+@UID90
+AGAATGTAATACCTCCAGATATGTAATTATTTGTAACAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIII
+@UID91
+CGAGATCCCACTCCTGACACCATCGGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIFFIIIIIIIFFFFFFFFFFFFFFFFFFFFFFFFFFFB
+@UID92
+TCGATCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGTGCTAAGATCGGAAGAGCACACG
++
+FFFFFFFIIIIIIFIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFFFBFFFFBFFFFFFFF
+@UID93
+AGAATTGATAGATCGGAAGAGCACACGTCTTGTAGCAGATCGGAAGAGCACACGTCTGAACTC
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIFFFFFFF
+@UID94
+TCCAAAATAGAAGTCCTGATAGCCCAGACAACTGAGATCGGAAGAGCACACGTCTGAACTCCA
++
+FFFFFFFFFFFFFBFFIFFFFFFBBFFFFFIBFFBFFFFIFFFBBBFBFBFFBFFF<B<FFFF
+@UID95
+TCACTCATCCTCTATGGACATTGGAGTTTTTCATTATGTGCATGGATCAAGATCGGAAGAGCA
++
+FFFFFFFIIIIIIIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIF
+@UID96
+AGCACTGAATCCAGCGATCCGAGTTCTGTCAAGATCGGAAGAGCACACGTCTGAACTCCAGTC
++
+FFFFFFFIIIIIIIIIIIIIIFIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFF
+@UID97
+AGCCTTTCACATGTATAAGGTGGCTCCATAGCTTACTGAGATCGGAAGACACACGTCTGAACT
++
+FFFFFFFIIIIIIIIIIIIIFIIIIIIIIIIIIIIIIIFFIIIIIIIIIIIIIIIIIIIIIII
+@UID98
+AGCCGCGTGGTGGCATGTGCCTGTGGTCCCAGATAGCAGATCGGAAGAGCACACGTCTGAACT
++
+FFFFFFFFIIFFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFFFFFFF
+@UID99
+CGAGAACGCATAAATCTTTCGCCTTTTACTAAAGATTTGTGGCAGATCGGAAGAGCACACGTC
++
+FFFFFFFIIFIIIIIIIIIIIIIIIIIIIIIIIIIIIIIFIIFFFFFFFFBFFFBFFBFFFFF
diff -r 000000000000 -r daad7189e48e test-data/twomates.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/twomates.sam	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,5 @@
+@HD	VN:1.0	SO:unsorted
+@SQ	SN:chr1	LN:249250621
+@PG	ID:bowtie2	PN:bowtie2	VN:2.2.0	CL:"/usr/local/user/bowtie2/bowtie2-2.2.0/bowtie2-align-s --wrapper basic-0 --very-sensitive --end-to-end --no-mixed --no-discordant --fr --omit-sec-seq --minins=0 --maxins=200 --threads 24 --time -x ../genome/hg19/hg19 --passthrough -1 /tmp/10410.inpipe1 -2 /tmp/10410.inpipe2"
+39V34V1:167:HAWHUADXX:1:1101:19513:10724	99	chr1	45244081	40	41M	=	45244081	-41	GTCCAGTTCTGCTACCGACAGTAAGTGAAGATAAAGTGTGT	FFFIIIFIIIIIIFIIIIIIIFIIIBFFIIIIIIIIFFFFF	AS:i:-6	XN:i:0	XM:i:1	XO:i:0	XG:i:0	NM:i:1	MD:Z:15T25	YS:i:-6	YT:Z:CP
+39V34V1:167:HAWHUADXX:1:1101:19513:10724	147	chr1	45244081	40	41M	=	45244081	-41	GTCCAGTTCTGCTACCGACAGTAAGTGAAGATAAAGTGTGT	IFIIFIIIFIIIIIIIIIIIIIIIIIIIIIIIIFFFFFFFF	AS:i:-6	XN:i:0	XM:i:1	XO:i:0	XG:i:0	NM:i:1	MD:Z:15T25	YS:i:-6	YT:Z:CP
diff -r 000000000000 -r daad7189e48e test-data/tworeads_aln_ends.bed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/tworeads_aln_ends.bed	Tue Nov 14 05:53:44 2017 -0500
@@ -0,0 +1,1 @@
+chr1	45244080	45244121	39V34V1:167:HAWHUADXX:1:1101:19513:10724	2	+