changeset 0:55107d4a728c draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sam_to_fastq commit 0651eb8c86d890e4b223fec82ab3980932710030
author drosofff
date Mon, 21 Mar 2016 17:33:26 -0400
parents
children 8acb72264319
files sam_to_fastq.py sam_to_fastq.xml test-data/input.sam test-data/output.fastq
diffstat 4 files changed, 282 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/sam_to_fastq.py	Mon Mar 21 17:33:26 2016 -0400
@@ -0,0 +1,33 @@
+#!/usr/bin/python
+#
+import sys
+import argparse
+
+def Parser():
+    the_parser = argparse.ArgumentParser()
+    the_parser.add_argument(
+        '--input', action="store", type=str, help="input SAM file")
+    the_parser.add_argument(
+        '--output', action="store", type=str, help="output FASTQ file")
+    args = the_parser.parse_args()
+    return args
+    
+        
+def print_fastq_sequence(samline, file):
+  samfields = samline[:-1].split("\t")
+  file.write ( '@%s\n%s\n+\n%s' % (samfields[0], samfields[9], samfields[10]) )
+
+def main(input, output):
+    infile = open (input, "r")
+    outfile = open (output, "w")
+    with open (input, "r") as infile:
+        with open (output, "w") as outfile:
+            for line in infile:
+                if line[0] == "@":
+                    continue
+                if line.split("\t")[1] != "4":
+                    print_fastq_sequence (line, outfile)
+
+if __name__ == "__main__":
+    args = Parser()
+    main (args.input, args.output)
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/sam_to_fastq.xml	Mon Mar 21 17:33:26 2016 -0400
@@ -0,0 +1,28 @@
+<tool id="sam_to_fastq" name="SAM to FASTQ" version="0.1">
+  <description></description>
+  <command interpreter="python">sam_to_fastq.py --input $input --output $output</command>
+  <inputs>
+       <param name="input" type="data" format="sam" label="SAM file for extraction of aligned reads in FASTQ format"/>
+   </inputs>
+
+ <outputs>
+   <data format="fastqsanger" name="output" label="FASTQ extraction" />
+</outputs>
+
+    <tests>
+        <test>
+            <param ftype="sam" name="input" value="input.sam" />
+            <output file="output.fastq" name="output" />
+        </test>
+    </tests>
+
+
+<help>
+
+**What it does**
+
+Extract sequence and sequence quality of aligned reads in a SAM alignment file and return a FASTQ file containing those reads
+
+</help>
+
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/input.sam	Mon Mar 21 17:33:26 2016 -0400
@@ -0,0 +1,201 @@
+@HD	VN:1.0	SO:unsorted
+@SQ	SN:chr1	LN:249250621
+@SQ	SN:chr2	LN:243199373
+@SQ	SN:chr3	LN:198022430
+@SQ	SN:chr4	LN:191154276
+@SQ	SN:chr5	LN:180915260
+@SQ	SN:chr6	LN:171115067
+@SQ	SN:chr7	LN:159138663
+@SQ	SN:chr8	LN:146364022
+@SQ	SN:chr9	LN:141213431
+@SQ	SN:chr10	LN:135534747
+@SQ	SN:chr11	LN:135006516
+@SQ	SN:chr12	LN:133851895
+@SQ	SN:chr13	LN:115169878
+@SQ	SN:chr14	LN:107349540
+@SQ	SN:chr15	LN:102531392
+@SQ	SN:chr16	LN:90354753
+@SQ	SN:chr17	LN:81195210
+@SQ	SN:chr18	LN:78077248
+@SQ	SN:chr19	LN:59128983
+@SQ	SN:chr20	LN:63025520
+@SQ	SN:chr21	LN:48129895
+@SQ	SN:chr22	LN:51304566
+@SQ	SN:chrX	LN:155270560
+@SQ	SN:chrY	LN:59373566
+@SQ	SN:chrM	LN:16571
+@SQ	SN:chr1_gl000191_random	LN:106433
+@SQ	SN:chr1_gl000192_random	LN:547496
+@SQ	SN:chr4_gl000193_random	LN:189789
+@SQ	SN:chr4_gl000194_random	LN:191469
+@SQ	SN:chr7_gl000195_random	LN:182896
+@SQ	SN:chr8_gl000196_random	LN:38914
+@SQ	SN:chr8_gl000197_random	LN:37175
+@SQ	SN:chr9_gl000198_random	LN:90085
+@SQ	SN:chr9_gl000199_random	LN:169874
+@SQ	SN:chr9_gl000200_random	LN:187035
+@SQ	SN:chr9_gl000201_random	LN:36148
+@SQ	SN:chr11_gl000202_random	LN:40103
+@SQ	SN:chr17_gl000203_random	LN:37498
+@SQ	SN:chr17_gl000204_random	LN:81310
+@SQ	SN:chr17_gl000205_random	LN:174588
+@SQ	SN:chr17_gl000206_random	LN:41001
+@SQ	SN:chr18_gl000207_random	LN:4262
+@SQ	SN:chr19_gl000208_random	LN:92689
+@SQ	SN:chr19_gl000209_random	LN:159169
+@SQ	SN:chr21_gl000210_random	LN:27682
+@SQ	SN:chrUn_gl000211	LN:166566
+@SQ	SN:chrUn_gl000212	LN:186858
+@SQ	SN:chrUn_gl000213	LN:164239
+@SQ	SN:chrUn_gl000214	LN:137718
+@SQ	SN:chrUn_gl000215	LN:172545
+@SQ	SN:chrUn_gl000216	LN:172294
+@SQ	SN:chrUn_gl000217	LN:172149
+@SQ	SN:chrUn_gl000218	LN:161147
+@SQ	SN:chrUn_gl000219	LN:179198
+@SQ	SN:chrUn_gl000220	LN:161802
+@SQ	SN:chrUn_gl000221	LN:155397
+@SQ	SN:chrUn_gl000222	LN:186861
+@SQ	SN:chrUn_gl000223	LN:180455
+@SQ	SN:chrUn_gl000224	LN:179693
+@SQ	SN:chrUn_gl000225	LN:211173
+@SQ	SN:chrUn_gl000226	LN:15008
+@SQ	SN:chrUn_gl000227	LN:128374
+@SQ	SN:chrUn_gl000228	LN:129120
+@SQ	SN:chrUn_gl000229	LN:19913
+@SQ	SN:chrUn_gl000230	LN:43691
+@SQ	SN:chrUn_gl000231	LN:27386
+@SQ	SN:chrUn_gl000232	LN:40652
+@SQ	SN:chrUn_gl000233	LN:45941
+@SQ	SN:chrUn_gl000234	LN:40531
+@SQ	SN:chrUn_gl000235	LN:34474
+@SQ	SN:chrUn_gl000236	LN:41934
+@SQ	SN:chrUn_gl000237	LN:45867
+@SQ	SN:chrUn_gl000238	LN:39939
+@SQ	SN:chrUn_gl000239	LN:33824
+@SQ	SN:chrUn_gl000240	LN:41933
+@SQ	SN:chrUn_gl000241	LN:42152
+@SQ	SN:chrUn_gl000242	LN:43523
+@SQ	SN:chrUn_gl000243	LN:43341
+@SQ	SN:chrUn_gl000244	LN:39929
+@SQ	SN:chrUn_gl000245	LN:36651
+@SQ	SN:chrUn_gl000246	LN:38154
+@SQ	SN:chrUn_gl000247	LN:36422
+@SQ	SN:chrUn_gl000248	LN:39786
+@SQ	SN:chrUn_gl000249	LN:38502
+@PG	ID:bowtie2	PN:bowtie2	VN:2.2.6	CL:"/home/galaxy/tool_dependency/bowtie2/2.2.6/iuc/package_bowtie_2_2_6/0d9cd7487cc9/bin/bowtie2-align-s --wrapper basic-0 -p 8 -x /home/galaxy/galaxy-dist/tool-data/hg19/bowtie2_index/hg19.fa -k 1 --passthrough -U /home/galaxy/galaxy-dist/database/files/000/162/dataset_162397.dat"
+1/1	4	*	0	0	*	*	0	0	GGGTTCTTGGCATCAGTTGCTGAATGAAGCCTCTCTGATGTAGCTTGGTCTTCTTCAGGGACTCCTAACAGTGGAGTGGAGGCTGTCTATGACTCTTTTAC	BBCDFFFFHHHHHJJJJJJJJJJIJIJJJJJIIIIJIHIIIIJJJJJJIIIJIJJJIIJJJJJJJJJGIGHGHJHDEHHHFFDDCAAC>@@;ACCDDDDCC	YT:Z:UU
+4/2	4	*	0	0	*	*	0	0	CACATCAACATGATTCATTACCATCATTTCATCAATAATCACTTTATTTATCTTGTTTCAACTAAAAATTTCTTTACAATCATTTCCTCTATCACCCTCAT	CCCFFFFDHHHHDIJJJJJJJJJJIJJJJJJJJJIJJJJJJJJJJJJJJJGIJJIJJJJJJJJJJJJHIIJJJJJJIIIJJJJJJHHHHHHHFFFFFDDED	YT:Z:UU
+1/2	4	*	0	0	*	*	0	0	CTGCTGCTTCTATGATTGCCATCTCCTGGAACTCAGACCCCTTCTAGAGTCTTTCTCCATTTAAAAAAGTTGTCTCATATGTATATTACATCCCAACCACA	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJGJJJJJJJJJJJJJJJIHIIJJJJJJJJJJJIJJJJHIIHHHHHHFFFFEDEFFEEEEEEDDDDDDDD	YT:Z:UU
+5/1	4	*	0	0	*	*	0	0	GTTTCTATCTTACTAATTTGTCTTGAGTTTTATTATTTCTTGCCATGTAAACATTTATGTGTTAGTTTTTTTGTGAGTGTGCATGTATACATGCATGTGGG	CCCFFFFFHHHHHJJIJJJJJIJJJHJHJJJHIJIJJJIJJJJJJJJJJJIJEIJJIIJIIJJJIIJJJJJJJHGHEEHHFFFFFFFFFEEEEECEDACDD	YT:Z:UU
+5/2	4	*	0	0	*	*	0	0	TGCGCAGCTCAAAAGGTACAAAAATTTGAAATAATACACTGTATCCTTTCTGACCACAATTTAAACATTAGTTCAAAATAAATATCAAAAACACCAGGAAG	CCCFFFFFHHHHHJJJHIJJJJJJJJJIJJJJJJJJJJJJJHIJIIJJJJGIJJJJJJJJJJJJGHIJJJJIIJHHGHHHFFFFFFFEEEDDBD@DDDDD<	YT:Z:UU
+2/2	4	*	0	0	*	*	0	0	GAAATATTCCCCATCTTCAGAGCCATGCAAAAAGGCCTCATTGACCAAGACACAGGCCTTGTGCTCCTGGAATCCCAGATTATCACGTCTGGCCTCATTGC	@@CFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJGIJJGIJJIJIJJJIJJJJJJJHHHHHFFFFFEEECCEEDDDDDDDDDDDCEDC	YT:Z:UU
+6/1	4	*	0	0	*	*	0	0	CCCTCCTATCTGATTCATCCACATGATGTTGACTCCAGTGGATTAACCTTGGCAGTTTATGCAGGTCACCCATTCCGGCCGCTCGTTGTAGTGGACCTCTT	CCCFFFFFHGHHHJJJJJJJJJJJJJJJJJJJJJJJIJIJJJJJJIJJJJJIJJJJJJJJJJJIJJJJJJJJJJJJJHHFDDDDDDDDDDDEDDDDDDDDD	YT:Z:UU
+6/2	4	*	0	0	*	*	0	0	CCGTGAGGCGCCGGCGGCCGGGCGCCCGGAGCGGCGGTAGCCGGCGGCCTACGGTTCTCCTCCCCCGCCGCCCTCCCACCCGAGCTGCGGCGGGGCTTGGC	CCCFFFFFHHHGHJJJJBHADFDBBBDBD6@@BDBBD-70(80757B-75>?5@+88>>C:3?BDD;;B59BD9ABB9<BD0&0509?9B>59@D######	YT:Z:UU
+7/1	4	*	0	0	*	*	0	0	CAATTATACAAATAAGGAAATAGGTCAACATAGAATCTCACATAAATCTAATTACTGGTCAATATAAAAAGGTAACAGTTCAAGGTCATATGTTGCCTAGC	CCCFFFFFHHHHHJJJJJJJJJJJHIJJJJJJJJJJJJIJJIJJJJJJJJIJJJJJJJJJIIIJJIIJJJJJ@GHHIIHIJGEHHHHGFFFFFFFECECCE	YT:Z:UU
+7/2	4	*	0	0	*	*	0	0	ATTTTAGTGTATAGGTAATGAATACTATATGTCAGTATGCCAGCATGTTCTGGGCTAAAGAAAGCAACATTTGACAGAAGAGAAATTTTGGTTATTGTCGG	?BBFFFFDHHHHHJJEGIJJJJJJJJJJJJIIJJJIIJJJJJJJJJJJIJIIJJJJJIJJJJJJJJIJJJJJJIJGIGIIIHHHHHFFFFDCDEEEACDDB	YT:Z:UU
+8/1	4	*	0	0	*	*	0	0	CACCTTCCTCCCGCCTTTTTCTACCGGATCAAACTGATGTTGACCCGTGAAAGGCAAATACACCGCCAGAGACGGCTGAATACCTTTCACACGACTGACTT	CCCFFFFFHHHHHJJIJJJJJJJJIIJJJJJJJJJJIIJJJGIJJJJJJJJJJJJJHHHGHFECDCD@DDDDDDDDDDDDDEDDDCDDDDDCDBDDDDDDD	YT:Z:UU
+8/2	4	*	0	0	*	*	0	0	CGAAGAGGGTGACAATCCCGTACTTGACGCGACAACCGGTGCTCTGTGATACGTTCTCAATGAGTCGAGTTGTTTGGGAATGCAGCTCTAAATGGGTGGTA	@@;DDDDDFADDDHGHIEFHGIIJJGGGIIGBGGHFGBAHABEEEC?CFDDACECAB>@CCC;C;:?B?ABCDBDCDBBDDCCCDDDDDDDDDEDD?BD9C	YT:Z:UU
+3/2	4	*	0	0	*	*	0	0	GCAAGCAGCCAGCAAAAGGAAATAGCCCAGCGACAGCCACGAAGAGGAAAATGGCAAACGTGTCGGAAAAGAAGAGAAAAAAGAAAATGTAACTCACATGG	CCCFFFFFHHHHHJJJJJJJIJJJIJJJJJJJJJJJJJJJJIIJJJJJJIIHHHHHFFFFEDEEDDDBBDDBDDDDDDDDDDDDDDDDDDCDDDDDDDCDC	YT:Z:UU
+9/2	4	*	0	0	*	*	0	0	GTTGAAGTTGCCTTTGGAAACTTTGGACAAACAGGAAGAGAAACAAACCACCCCAAGGATAATAGCCATCACTGTTGGCTTAAATCAGATTTGGTGAATTT	CCCFFFFFHHHHHJJJJJGJJJJJJJJJJJJJJJJIIJJJIJJJJIJJJJIJJJJJJJIJJJJJJHHHHHFFFFDEEDEEEDDDEDDDDDDDDDDDDDFED	YT:Z:UU
+9/1	4	*	0	0	*	*	0	0	GTGTAGAGTATTTATTCTACCTTAAAGGATACAGAAGTATAAATACAGCTTAGAACTTCAACTGCTCACTCAGGAAGAAACAAAGGTGACCACAACTGGGA	BBCFFFFFFHHHHIJJJJJJJJJHJJJJJJJJJJJJJFHIJJJJJJJJIJJJJIJJJJJJJIIIJJJJJJJJJJJIIJJJJHHFHGBEFFFEDDDDDDDDD	YT:Z:UU
+10/1	4	*	0	0	*	*	0	0	CCTGACCACTTCCACAAGCAAGGCCATGCAGACTGTCTTCTCCTGAATTGTCACTGAATCTATACACTAACTCAGCAACAATGGCATTTAATAATAACCAG	CCCFFFFDHHHHHJJJJJJJJJJJJJIJJIJJFIIIIJJGGIIIJGIIJBGIGIHIJJJJIJJJJJJJIIJJJFJEHHHFFFFFDCEEEEEDEDDDDDCC@	YT:Z:UU
+10/2	4	*	0	0	*	*	0	0	ATTGGGTTTAGCCCTCTGTGGTAGGTATCTGCACTCCAAGAGACAACTGGGTTAGTAAATAATATCCGTTCTATAGGAAGAAAGAAGTCAATCATCACCGG	@CCFFFDFHHHHHJJJJJHIJGHIIHHIGIJIIGIIJJJJJIGIJJJJJJHGGGIIJJJJJIJJJJIJJJHHHHHFFDEFDEEDEEDDDCADDDDDDDDD6	YT:Z:UU
+11/1	4	*	0	0	*	*	0	0	GGCGCGGCGCAGCATGGCGGGGGTGTGTGGTAAGCCCCACATGTCACCGGGATCGCTGCCCGGGAAACGACGTTTGGAACCTGATCAGGAGCTGCTGATAC	CCCFFFFFHHGHHJJJJIJJHDD)75?8AD?CCC@CDD-7<CDCCDDDDDD5:?BBB@<BDDDD>BDCB@;;9??<BCCDDDD3::ADC>@?@C>CC:A>4	YT:Z:UU
+11/2	4	*	0	0	*	*	0	0	GCCTCTCATCTTCAACCTGGCAATGACACACATGTGCTGCTTCAGTCATCTATCTCAGACCCCCAGATTGTTCTTCCAGATGGGGGGGAGAGCCTAGGCCA	CCCFFFFFHGHHHJJJIJJGHIJFGJGHFHIIIIGIJJIIJJIGHIGIIIHIIJJJJHIJJJJJIIHHIJHEHHFGHFF@CFFDDDDBBDDDBDD>CCCDB	YT:Z:UU
+12/1	4	*	0	0	*	*	0	0	CTTCTGTGGGCTTCACAGGCCACAGAGTGGATCTCAAGCCAGGCACCAGGACAAGCCACAGCCTTGCTCACTGTTGGCATGAAGCATCTACCGAGGTATTC	@BCFFFFFHHHHHJJJJJJJJJJJJHJGHIJJJJJJJIJJJJIJIJJJJJIJHGHIJJJJJHHHHHFFFFFFEEEEDEDDDDDDDDDDDDDDDDDD>CEED	YT:Z:UU
+13/1	4	*	0	0	*	*	0	0	CCGGCCGTTCATGGGCTGGCAGCGGCAAGATCATGGGCCGAGCGTGCAGGGCGCGATCGAGCGGGCGGCGCACGACATCACCGGCGAGACGATAACGGTGC	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJJJJJHHFDACDDDDDBDDDBDDDDDDDDDCDDDDDDDDDBDDDDDDDDDDDDDDDDDDDDDD@BB	YT:Z:UU
+14/1	4	*	0	0	*	*	0	0	GCCGTCAGGTAGATGGTGCCGCCGCGCACGGGGTTGCCGGCCTTGAAGTCCTGCGCACGGCGCGGATCGATCAGCCTGGCGCGCGAGGCCAGGTAGGCGTC	@CCFFFFFHFHHHJJJHIJJJJJJJJIIJJJJJAHHFFFDDDDDDDDDDDDDDDDDDDDDDDDDDBDDDDDDDDDCDDDDDDDDBBDDDDDBD>@CD>>>@	YT:Z:UU
+12/2	4	*	0	0	*	*	0	0	TGGGATTTGAACTCAGGACCTTCGGGAGAGCAGTCGGTGCTCTTAACTGCTGAGCCATCTTACCAGCCTGGGTCTCCCATAAGTTAATAACTGTTTACAGT	BCCFFFFFHHHHHJJJJIIJJJJJJJJJIJJJIHIJJIIJJJJJJJJJJJIJJJJJJJJJHHHHHHFFFFDDCEDDDDDDCDDDDCDEDCDDDDDDCDCCA	YT:Z:UU
+13/2	4	*	0	0	*	*	0	0	CGACGGCGTCGAGTGGCTGCGGCACGCGCCAAGCTAACCCCGCCTCCCACGTCAACGGCGCACGGCGGGTGACGATGTGGTATTCGTAATGCCGCGCCAGA	CCCFFFFFHHHHHJJJJJJIJJFGHIGEEFFDCCDDC?ADDDBDDDDDDDBBDBBDBDDDDDDDDDDDD0;<@<8<@8:C3>ADCCCCD?>@C<>59B###	YT:Z:UU
+14/2	4	*	0	0	*	*	0	0	GCGGCCACACGCTGCACGAGATCCCGCCGGGCGGGCAAGGCATTGCCGCGCAGATTGCGCTGGGCATTCTGTCCAATTTCGACCTGGCCAGCCTGCCGGTG	CC@FFFFFHHHHHJJJJJJJJJJJJJJJIHIJHDDDDDDDDDDDDDDDDDBDDDDDDDDDDDDDDDDDDEEEDDDDDDDDCDDDDDDDDDDDDDDDDDDBD	YT:Z:UU
+15/1	4	*	0	0	*	*	0	0	TGCCTCTCTCACTTCTCGCTCCCCCTTCCACGCAAGCCGGCGGGCACGCGGCGCCGGCGGCCGCAGCGCCCCGCGGGCGGGCCCCCCCCTCCCCCCCCGGC	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJIHJJJIIJFDDDDDDBB57<B@BB@B5;7557@505@D-0)5<559B05@BD&0909@DDD######	YT:Z:UU
+16/1	4	*	0	0	*	*	0	0	GGCCGATCTACTTGTTTGATTCTCTATCAAAGAGGCTGTTGGTCAATGAGAATAATAATAACGAAACTATGCAAGTGCCTTTTCTGTGTATTTCAATGATT	CCCFFFFFHHHHHJIJJJIJJJJIIIJIJJJJJJJJJJJJJIIJJJJJIJJIIJJIGJJJJJJJJJJHHHHHHFFFCEFEEEEEEEEDDEFFFECCDDDDE	YT:Z:UU
+15/2	4	*	0	0	*	*	0	0	CGCTTGACGATCAATCGCCCGCCCCTCCGCGGTGTGGCGCGGCTGGGAGTCTTGCTCGCAGGGCCCCTCCCCGGGGCCCTCCGTCTCCCGAAGTACAGACG	CCCFFFFFGHHHHJJJJJJJJJJJJJJJJJJJ@GGIJJHHFDDDDDDDD??CDCDDDDDDDDD?BDDDDDDDDBDDBDDDDDDBDDDDDD99BDDD@CCDD	YT:Z:UU
+4/1	4	*	0	0	*	*	0	0	GTTGAAATGAGTGAAGACGGTTACTGATTAAGCGTTCGGATGATGGAAATAGGATTGATGGGAATATAATAATGGTAATGACAATTGGTAGACCTATTACT	CCCFFFFFHHHHHJIJJJJJHIIJJJIJJJIJJJJJJJJGIIIIJJJJIJIIGHIJJIJJJJIIHHHHHHHFFFFFFFEEEEECDDDDCDDDDDDCDDEDC	YT:Z:UU
+17/1	4	*	0	0	*	*	0	0	GGGGTGTAGCCTTGTTAAAGTAGGAGTGGCTTTCTTGGAGGAAGTGTATCACTAGGGGTGGGCCTTGAGGTTTCAAATGCCTGCACCCGTGCTCCCTCAAC	CCCFDFFFHHHHHJJJJIJIHIJJIJFHIJJJJIJJJJFHJIIJIIIJJJJJJJJJJJGIGFFFFEEEEE?BCDDCDDDDDDDDDDDDDDDDDDDDDDDDD	YT:Z:UU
+17/2	4	*	0	0	*	*	0	0	TTGAACTCTAACTTTCTCTGTCGTTTCTCCCTGTAGCTTTCTTCTTCTGCACTTCCTTAGTTATGTGAGTCAGGTCCCTTTCTCACTCTGCATGAATTGGT	CC@FFFFFHHHHHJJJJIJJJJJJJJGIIJIJJIIJJJJJGJJJIJIJJIIHIJJJJIJJJJJJJJJIIIJJJJIIJIJJHHHGCCHFFFEDEECEEEEC>	YT:Z:UU
+19/1	4	*	0	0	*	*	0	0	GAAAGTAAGGCACACCCAAGAATACGGGCGTGACTCTGCCTCCAGGGCCAACCACACTAGCTAAGCCTTCTATGGAAGTTAGAATTCACTACCTTGAGACT	BCCFFDFFHHHHHIJIJJJJJJJJJJJJJJGHJJJJJIJJJJJJJJJJJJHHHHFFFEEEEEEDDDDDDDDDEEDDDCCDACCDDDCDDDDDDDDDCDBDD	YT:Z:UU
+19/2	4	*	0	0	*	*	0	0	CTGGGCTGTGTGTCTGAACGTGAGTGTGACTGGTTAGCCTGCTGTGCTTAGCCCGCTCTTATCCAAATTGTCAAGGAAGAGTGACTCCCATGTGGCAAGAT	BBCFFFFFHHHHHIIJJJJJIIJJHGIJBHIIJGHIIJJIJJJJJJJJJJJJJJIJJJJJJJJJHHHHHHFFFFFFEEEDEDDCDDDDDDDDDDCDDDBDD	YT:Z:UU
+16/2	4	*	0	0	*	*	0	0	CAGTCCAAGGTGGGCACGATGATTGCAGCTCATTGTTAGAAATTTGATGAATACAATACTTAATTCTACTCTGAGTTGAAAACCCTGCAAAACAGCATGGC	?BBFFFFFHHHHHJJJJJJJJJJJJJJJJJIIIJJIJGIJJIJJJJIJJJJIJJGIJJJJJJJJJJJJJJJJJJHHHHHHHFFFFEEDEEDDDDDDDDDDD	YT:Z:UU
+18/2	4	*	0	0	*	*	0	0	GCCTGTCAGCATACAAGTACCACAGTGAGCAGATGAAGGTCAGAGGACAACTTATGGGTGTTGATTCCTTCCATCTACTATGTAGGCTCTAGTTTTGCAGC	CCCFFFFFHHHHHJJJJIJJJJJJJHIJJIJJIJJIIJJEIIJJJJJJJJIIJGIJJJFHIJJJJJJJJJJHHHHHHFFFFFFFFEEDEEDDEEEDDDDDC	YT:Z:UU
+20/2	4	*	0	0	*	*	0	0	CGGGTGGGTAGCCGACATCGCCGCCGACCCCGTGCGCTCGGCTTCGTCCCGCCCCCCGGGAGGAGGGAGGGACGCGTGACCGACGGCCCCCCGGGCCCGAC	BCCFDFFFHHHGHJJJJJJJJJJJIJJJJJHHEFFDDDDDDDDDDDDDDDDDDDDDDDDDBDD?DDDBDDDBBDBDDDDDDDBDDDBBDDDDDDDDDDD9<	YT:Z:UU
+21/1	4	*	0	0	*	*	0	0	TTCGTCGCGCACCCGGGCCAGGGCCCGGACCCCGCGCAGGGGGTCACCCTCGCTGTGCAGATCCACATGGGTGCGCACCGCCGTGGTGCCGCTGGCCACCA	BCCFFFFFHHHHHJJJJJJIJJJIJJJIGHHFFCDDDBDDDDDBDBDDDDDDDDDDDDEDDDDDDDCDDDDCBDDDBDDDDDDDDDBBCDDDDDDDDDDDB	YT:Z:UU
+21/2	4	*	0	0	*	*	0	0	GGGCGAAACACCGACCCGTGACACTGATCTGTATTGACGAGGTCAATATCCCGCGTCCTGGTCGAAAGCGCAAGGCACCACACCGACCCGGCCGGGGACGG	CCCFFFFFHHHHHJJJJJJIJJ>GIIIJIJJIIIIIJJJJJJHJJJHHHHHHFFDBDDDCDCCDDDDDDDDDDDDBBDDDDDDDBBDBDDDDDDDDB>BBB	YT:Z:UU
+22/1	4	*	0	0	*	*	0	0	GCGAGGCCTACGGGAGAAGGCAGCAAGTGGTGGGCGCTCACCCATCGAGGAAGAGCCTTCTTGCCACTCGGGCAAGATCCGCTTTCTTCTCGACCCCCCTT	CCCFFFFFHHHHHIIIIHIJJGIIJJJHIIHIIJDFGIIGIIHHGHFEFDDEEDDDDDDDDDCDDDDDDDDDDDDDDDDDDDDDDDDDCCDDCDDDDBBDD	YT:Z:UU
+22/2	4	*	0	0	*	*	0	0	TGGACGCGTCAGGCTCTCGATTGTCCTTACGCCGTCCGTGTCGCTCTCGTGTCCCCTTCCCCGAGCCGGGGAGGGGTGGCTGCTCTGTTGCGGTGCATACC	CCCFFFFFHHHHHJJIJIJIJIHGIIIJEIIJJIGHIIHEGIHHHFFFFDDDEDDDDDDDDDDBDDDBDDD@BDDD9@DBDDBDCCDCDCCDD@B>CDCAC	YT:Z:UU
+23/2	4	*	0	0	*	*	0	0	GGAAGGATTAACTTGGAATAAAGTTTACATCCTGTTTAGAAGACAGCACTAGCACTCTCTGAAATACCAGTTATGACAGACAGATGCTGCCTGGGAGCCCT	BCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJIJJJJGIJJJJJJIJJJJJJJJJJJJJJJJJJIIJJHIJJJJIJHHHHFDFFFFDEEEEDDDDDDDC	YT:Z:UU
+24/1	4	*	0	0	*	*	0	0	GCGATGGAGATCCTGAAGACTCAAGTGGTGCCCGCGCTGGGGGCTCGCGAGGGGGTGGGCGACGCGAAAATTGCCAAGGTGTCTATCGTTGGCATCGGCAT	@@@DDDFDFDAFHHIGJIEHGGGIGEFHAFCHIGGEGEFIJJJEE?<9A8=B'5<07BD?B@7;5>;BDD:>:@CCCBD4>4:A@CC@<A?B:>@ABB###	YT:Z:UU
+23/1	4	*	0	0	*	*	0	0	CGGAAACTCTTCCCCCAAAGTTCTACTGTTTAAATGTGTGGAAAAATTAAAACCGGGATCAGAGTCTAGAAGTCTTGATCCAGTCTGACCAAGTTTTCATT	CCCFFFFFHHHHHJJJJJJJJJIJGIJGJJJJJJJJIGHIJJJJJIJJJJJJJHGIJJJJHHHHHHFFFFFFCDDCECEDDDDDDDDDDCDDDDDDDDDEE	YT:Z:UU
+24/2	4	*	0	0	*	*	0	0	CTCTCGGTCACGTCTCCAGCAACAGCAAGCATTCTAGCTCAGAAATTCAAGCCCTCTGCCTGTTGTGCTCAAAACAACAGATCAGGCCTTGGCCTCATCCA	@CCFFFFFHHHHHJJJJJJJIJJJIJJJJIJJJJJJJJJJJJJJGIJJJIJJJJJJJJJJJJJJJJHHHHHHFFFFDDDEDDDDDDDDDDDDDDDDCDDDC	YT:Z:UU
+25/1	4	*	0	0	*	*	0	0	ATGGCAAGCTTCTCTGAAGTCATCCAGAGTTTCTTGGATCATGTTCTGAATGAAGCGGATCTGGACTGAGGTAAGTGGCGCACTTGGCCGACTGTTCCCAA	BCCFFFFFHHHHHJJJJJJJJJJJJJIIIHIIIJJIJIJJJJJJJJJJJIJJJIJJJJHIJJJJJJJJHHHHFEFFFFDDDDDDDDDDDDDDDDDDDDDDC	YT:Z:UU
+26/1	4	*	0	0	*	*	0	0	AGCTTTGGGTGGCCCTCTCTCTAAAACAAAACCTACTGGCCTTGTAATGTGCTGTGAGCCCGTTTTTTTCTGTTATAAAGAACTCTCCTCCACATGCTCAT	@CCFFFFFHHHHHJJJJJJJJJJJJJIJJIJJJJJJJJJJJJJJIJJJJIJJJJJJJJJIJIHHHHFDDDDDEDDCDEEDDDDDDDDDDDCDDDDDDDDDD	YT:Z:UU
+27/1	4	*	0	0	*	*	0	0	GAGGGAATGCAATTAATTCATGTTGTTCTGGTTCCTGTCTCCATTTAAGCTTCCTCAGTCCCCACTAATTCCCATCAAGTGCCTTAATTAAACAAAGCTTA	CBCFFFFFHHHHHJJJJJJJJJIIJJJJJJJIJJJJJJIJJJJJJJJJJJJJJJJJFIIJIJJJJJGIJJJIIIIHIIHHHHHHFFFFFFEEEDDCDCDDD	YT:Z:UU
+2/1	4	*	0	0	*	*	0	0	GTCCAGATTAATACAGTTCTCACTAGGGGGGAAACTGAAACCTCCAGTTGCCAGGTGAATTTCTAAGACTTTTAATCCAAGTCTCTCACTGATTATTTTCT	CCCFFFFFHHHHHJJJJHHJIIJJIJJJJJDBDDDDDDDDDDDDDDDDDDDDDDDCCCCDEDDEDEDDCCDDDDDCDCDCCCDEDDDDDDDDDDDDEEEEC	YT:Z:UU
+27/2	4	*	0	0	*	*	0	0	TTTTTTTCTTTTTTTGATTGGGTCTTCCTGTTTTGGAGTCTAGGCTGTCCCGGATTCTGTAGATCAGGCTAGCCCCCAAATTCAGAGAGATCCTCTTACCT	BBBFFFFFHHHHHJJJIJJJJJGIIJIIIIJIJJJJGJIIJJJJJJIJJHHHFEDDEEEEEDEEDDDDDDDDDDDDDDD@CCDCDDDDDDDDDDDDDACDC	YT:Z:UU
+25/2	4	*	0	0	*	*	0	0	GGTGCTGCACTCCAGCCCTCTTAACGTCCTCATGGGGTCTCCAGGGAAAGAAGAAAATGAAAGCCATGAGCTGTCAGCCGAGACTAAGAGAGCATACCTAG	CC@FFFFFHHHHHIJJJJIJJJJJJJJJIIJJJJJJJDGIJJJJJJGIIIGGGIIJJIIJJJJJHHHHHFFFFFEEEEDBDDBDDDDC@CDDDDDCDDCCD	YT:Z:UU
+26/2	4	*	0	0	*	*	0	0	CTGGATTTTAGTTGAAAGTTGTCAAGGAATGTATTTTGGAACTGTTTTGACCATGAGACCATCTTGTATAAAGGTTGTTTTAGTCACCTACTGACTAAGCT	CCCFFFFFHHHHHJJJJJIJJJJJJJJJJJJHIIJJJJJIJJJJIJJJIIDHIJJJJJJJJJJJJJIJJJIJJJCDHAHHHHFFFFFFEEEEDEDDCDDDD	YT:Z:UU
+28/2	4	*	0	0	*	*	0	0	GTTCTATACATAACACTGGCCACAATCTGAAAGGACTTTCAAATTTGCATTTTACTCTTGTACCTTTGGGTAAGATGATCTTTTGAGTCCTGTTCATTTTC	CBBFFFFFHHHHFJJJJIJJJGJJJJIJJIJIJJJJJJJJJJJJJJJIIJIJJJJJJJJJJJJJJJJJJJJJJJJJJJHHHHHHHFFDFEEEEEEEEEEED	YT:Z:UU
+3/1	4	*	0	0	*	*	0	0	GCCATCTGCCAAAATAATTTTATCATTCAACAAAACAAAGCACCACCTAAAAGTGGTTCTTTTTTTAAGCTGAACATTTTCTCCAGAAGGAGAAATTTTTT	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJIJJJIIGIJIGIJIJJJJJJJJJDGIHIJJJJJJJHHFFFFFFEECEEFEEDDDDDDDDDDDBDCDDDDDD	YT:Z:UU
+29/1	4	*	0	0	*	*	0	0	GTGGGGCCAGATGTAGAACTTTTCCAGTACCATGCCTGCCTGTGTACCACCATGTTCTCAGCCATGATAATTGACTAAATCTCTAAAACTTTAAGCAATTT	B@CFFFFFHHHHHIIJJJJIJJJJJJJJIJGGGFJIIIJJFIHIHGIJGIGGIIJJGGEHIIIIHI@GIGHIGHHEHHEDDFFDFFEEDEEECCA@CCCDA	YT:Z:UU
+29/2	4	*	0	0	*	*	0	0	CCTTGGTCAAGTTACCTAGGTTTGTTGTTCCTGTATTAAATTTTCACACAAAAATAATTGGGTAAAAGCAGGGAAGAAAGAGACAAGACCATCACAGCACC	CCCFFFFFHHHHHJJJJJJJHIIJJJJIJJJJIFHIIJIJJJJIJJJJJIIHIJJJJJJJJJGGIEHIJGIJJEHHFFFFFEEEEDDDDCCDCDCCDDDDB	YT:Z:UU
+28/1	4	*	0	0	*	*	0	0	CTCAGATATTAACAGATGTTTACACAGGGTAAGGCAATCCATTCAAAGGAACAGGGTATTTGGGGGCTGGATAAATTTTGCAAAGCTGGTTTACAAGGTTT	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJFFHIJJJJJJJJJJJJJJJIJJJIJJJFGHIJJJIJJHFFDEEEEE@CDDDDEDDDDDDDDDDCDDDDDDDD	YT:Z:UU
+30/1	4	*	0	0	*	*	0	0	GCCAAGTGTGGATGTAAACTGAGAAGGGAGATAGGAACACATATAGTCCATTACAGTTCCCTAAGAGCTTTGATCGCATTCTCCTAAGTTCATGCATGTGT	CCCFFFFFHHHHHJIJIJJJJHIJJJJJFHHIJJIJJJJJJJJJIJIIIGIIJJJJJJJAHHJJJJJJJJIIJJJHHHFFFFFEEEEEEDEEEDDDEDDDE	YT:Z:UU
+31/2	4	*	0	0	*	*	0	0	TGTGCCTTGTGAACTGCTTGCTTGTGGCACGGAAATTTCCTTCCCGGCACCTGTAGAAACTGAGCAAGCAAGATGCTTGTTCCAATCAAAGGAAGTGGAAG	BCCFFFFFHHHHHJJJJJJJJJJJHIJJJJJJJJJJJJJJJJJJJJJJJJHIIHGIJJJJJFHHHGFFFF@CDEEEDDDDDDDEDDDDDDDDDDDCCDDD>	YT:Z:UU
+32/1	4	*	0	0	*	*	0	0	GTTGTGTTTAAGTTCTAACAATGAATTTTTACATATTAGTTGAATTTTAATTCTATTAATGTATGTTCAATTAACACGTTTTACGCCGAGAATAATTAATT	CCCFFFFFHHHHHJJJIJJJJJJIJJJJJJEIIJIIJJIJIGIIIJJJJHIJIJIJJJIJJIJJJJJJJJJJJJJFHIIJJJJJGEHFDDDDDDDEFECED	YT:Z:UU
+32/2	4	*	0	0	*	*	0	0	CCTGGCCTTATAATTAATTGGAGGTAAGATTACACATGCAAACATCCATAAACCGGTGTAGAATCCCTTAAAGAACTATTAAAATTTAAGGAGAGGGTATC	@CCDDFFFHHHHHJJJJJJJJJJJFFHIIIJIJJJJJJJJIJJJJJJJJJIJJJJJIIIHIJHJJIIJJHHHHFFFFFFFEFEEEEEDEDDDDDDDD?AC@	YT:Z:UU
+33/2	4	*	0	0	*	*	0	0	CGGCATTCACAGCATGCCAGCAGTTGTCTCTGAAACTGTAATTTCTGTTTGCCAGTGTTTGTTATTTATCCTATAGTTCTCTGTCTGGACTGCAGTTTAAT	@BCFFFFFHHHHGJJJJJGJJJJJJJJJJJJJIGIJJJIGHJJJJJJJJJIJJJJHIIJIJHIGIJJJJJJJIJJJJJJJJICHHHGHFFFFFEEEEEEDD	YT:Z:UU
+30/2	4	*	0	0	*	*	0	0	CCCAGATCAACCACAAGCTTGAGGTGAGCCTCGGGGCTGGCCGCAATGGCATTCAAACCCCGGGCTGTCCGAGGGCGTTAGATCCCATTTAAACGCCCTAG	CCCFFFFFHGHHHJJJJJJJJIJJHIIJJJJJJJJJJIIJIJJJGHHHHCDEDFFFEEDDDDDDDDDDDDDBBBDDDBB>BDDDDDDDDCCDACB>BDBCD	YT:Z:UU
+33/1	4	*	0	0	*	*	0	0	AGTGGCTTCACGTGGGTCTCCGTTGGTATTTGAGTGCATAATTCTTGATTCTTCCATCTTCCAAACTGTAAGTTTTAAGATATATTCTCCTGGCTTTCAGA	@CCFFFFFHHHHHJJJFHIJJJIJJJGEIIIDFDBGHIJIIJJFIJIGIJDHHHDHHGIJGIJGGIIJEHEEHHHEFHFFFFFFFFEEEEEEDDDDDDDDD	YT:Z:UU
+31/1	4	*	0	0	*	*	0	0	TTCAAGTTTAGGATTATTGTCTCCCGAACAGTTTACTAGTGAAATGCTAAGTTCTTCACCCAAAGGCTCCAGATGGCTCTCACTGGCTGCAGCTGCTTTTG	CCCFFFFFHHHHHJJJJJIJJIJJJJJJIJJJJJIIIJIHIJJJJJIJJIJJJJJJJJJIJJJJJJJJJJJJHHHHHFFFFFEEEEDDDDADDDDDDDDDC	YT:Z:UU
+34/2	4	*	0	0	*	*	0	0	GTCAGGGGCTCGTCCGAAAGCCGCCGTGGCGCAATGAAGGTGAAGGGCCCCGTCCGGGGGCCTGAGGTGGGATCCCGAGGCCTCTCCAGTCCGCCGAGGGC	BBCFFFFFHHHHHIJIIJJGGIEGIJGIJGHGHHHHBDFF;BCDDEDDDDBD@DDBDDDDBDDDD(<9?BDBDDDDDBDD@BCDDDDDDCCDBDDBBBBD9	YT:Z:UU
+35/1	4	*	0	0	*	*	0	0	TGGCTGTCTAATAAGAAAATGAAACATGAAGTCTCCAAGTACCACCTGCACAGAAATGCTCTAGCAATTTTGCTCAAGTATTGTTGGCTTACTTCTTCCAT	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJJJJHIJHIJJJIJJJJJJJJJJJJJJJJIJJIIIIJJIIHIJFHHHHHHFFFFEEEEEEEDDDDE	YT:Z:UU
+36/1	4	*	0	0	*	*	0	0	GCCCTCCTACTCGTCGCGGCGTAGCGTCCGCGGGGCCCGACGCCGCGGGGGCGAGACCCGGCGCGCGGGGGAGGGCGGCCGGAACCGACACCCCCCCCGAC	CCCFFFFFHHHHHIJJJJJJJJJJJIJIJJGFDDDDDDD6BBBDDBBDDDDDDBD<BD?7@<BBBDDDDDDBDDDBB@9@>B@BB?B3905@BDDDDD###	YT:Z:UU
+37/1	4	*	0	0	*	*	0	0	GCGGCGACGGGTATCAGGCTTCCTCGGCCCCGGGATTCGGCGAACGCTGCTGCCGGGGGGCTGTAACACTCGGGGCGGGGTGGGTCGGCCCCACGGGGGAC	@C@FFFFFHHF??DGIGFCGHIJIGGGHGIIGHFCDCC?ABDD19BDDB8CCCD<888<<B@7AC@CCCDCDD>55<<09&5@##################	YT:Z:UU
+38/1	4	*	0	0	*	*	0	0	CATTCACTCACCAGTACATACTGGGAACTGGCTAGGTGCCGGGCCTAAAGATACAATGTTGAATAGACATATTTGTGGAACTTATGAAACCAGACTGACTA	CCCFFFFFHHHHHJIJJJJJJJJJJJJJJJJJJJJJDHIJJJJJJJJIJJJJJJJJJJJJJJHHHHHHFFFFFFFFEDEEDDDDDDDEDDDDD?BDDDDDD	YT:Z:UU
+35/2	4	*	0	0	*	*	0	0	CTCAAAGTAAACTTCATTTTCTGCCAATGATAAGAAAAAAAAGAAAAACATGATAATCAGTTGTATGCATCAGTGATGAAATTTGGCAGTATGAGTTCAAG	CCCFFFFEHHHHHJJJJJJJJJJJJIJJJJJJJJJJJIIIJJJJJJJJJJHIJJJJJJJJHHHHHHHFFFFFFEEEFEECDDDDDDDDDDDEEDDDEDEED	YT:Z:UU
+38/2	4	*	0	0	*	*	0	0	ATGCGGTTCACCCTAAGGCCAGGCCCAACCCATCTCTTACCTATCCCAGCAGGTCACCCAGGTCACTCCCAAGCACTCTGGCCACTACACAGGTAATCTCA	@BCFFFFFHHHHHJJJJJJJJJJGHIJJJJIJJJJJJJJJJJJJJJJJJJHIJIJJJJJHHHHFFFFFFDEEDDDDDDDDDDDDDDDDDDDDDCDDEECDD	YT:Z:UU
+39/1	4	*	0	0	*	*	0	0	AGCTTAAAGAACCACCACAGATGCCTTCTTTAGGAAACACTTAATGGAAGATAACAAGTTTTCCAGGGGAGGAAGGTAGAAGTGACTTCTATTTTTTCAAT	C@CFFFFFHHHHFIJIGHGIIGGIIJJIJGGGIIGEIJGGIGHFHGFEEGGHIGHCIJCFHHJFFHHIJBHFDFD@>>@ACC;AAACCCCCCCDDDDDCDC	YT:Z:UU
+20/1	4	*	0	0	*	*	0	0	GGCGCCGAACCACCCCGCCCCGAGTGTTACAGCCCCCCGGCAGCAGCGTTCGCCGAATCCCGGGGCCGAGGAAGCCTGATACCCGTCGCCGCGCTCTCCCC	CCCFFFFFHHHHHJJJJJJJJJJJFHIIJIJIJJJHHDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDBDDDCDDDDDDDDDDBDDDDDDDD	YT:Z:UU
+39/2	4	*	0	0	*	*	0	0	GGAAACCCCCACAGGCGCGCTCGCGGGGCGGCCTCCGCGGGGCGAGTCCACCGCGCGGCCGGCCTGGTCCTCAACTCCCGCCTGCCGAGCGCCCGCGGCGC	B@?DFFFFHHHHHGGHIJJJJJJIJGIDDB@BD1AD@BDBDDDD8B@CCCCBBBDDDBBBDD@B@BDCCCAACCDCACBBBDBAAB<5<BDDBD@<>BDDD	YT:Z:UU
+40/2	4	*	0	0	*	*	0	0	GTCATCCAAGAAGAGAAACCACTCTTGACATGATGCCACCAGAAGATTGGCCTATCGGAAGCCGGAGTGGCATTTCCTTAATTAGTGATTGATGTGTGATG	CCCFFFFFHHHHHJJJJJJJJJJJJJGJJJJIJJJJJJJJJJIJJHHIIJHIJJJJJJJJJJJJHFF;ACDEDDDEEDDDDEEED>ACDECCDDDDDBCDC	YT:Z:UU
+34/1	16	chrUn_gl000220	158989	255	32M3D11M4D58M	*	0	0	GTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCCGAAACGATCTCAAC	DBDBDDDDDEDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD?FFHHHJJJJJJIHJJJJIJJJJIJJJJJJIIJIIJJJHJJIJJJJIJHHHHHFFFFFCCC	AS:i:-56	XN:i:0	XM:i:5	XO:i:2	XG:i:7	NM:i:12	MD:Z:32^CCC3G1A0C0C0C2^GCCA58	YT:Z:UU
+37/2	0	chrUn_gl000220	113795	255	31M2I22M1I45M	*	0	0	GCGCAGTCCGCCCGGAGGATTCAACCCGGCGGCGCGCGTCCGGCCGTGCCGGTGGTCCCGGCGGATCTTTCCCGCTCCCCGTTCCTCCCGCCCCCTTCACC	@@CFFFFFHHHHGGJHIJHIJGHHJJJAFEBEE?@BBDDBD6BDD5:B8B<88@@5@C?@B7<B87;CAAAADD<&&&8590<CCCDCBD&&0590<?CCC	AS:i:-40	XN:i:0	XM:i:6	XO:i:2	XG:i:3	NM:i:9	MD:Z:34G11T3C21C14A5C4	YT:Z:UU
+41/1	4	*	0	0	*	*	0	0	GTGTGTTTCCAAAACCTAGGTAATACCAGACATTTGGGGGAAAAGAAAGTGGAGGGTATGGTAAGTCTTTGTTTTGGCCACATCCTATTTACTTAGAAGAT	CCCFFFFFHHHHHJJJJJJJHIIJIJJJJIJJJJJJJJJJHIJJIIJIJFGIGIJJ?EHHFFFFFFFEEEEEEDDDDDDDDDDDDDDDDEEEDDEDDDDC>	YT:Z:UU
+42/2	4	*	0	0	*	*	0	0	GCCGGACGATATTTTATCGGACGACACGATCGTGGATCTTTTGAAGGGCCAGGGCGTTGAGATTGCACGGCGTACGGTTGCCAAGTACCGCGAGGCGATGA	@CCFFFFFGHHHHJJJJIJJJJJIJGGGIIJIDGIEHFHHHHFDEFFDDDDBDDDD@BD?C?CCDCDDDBDDBDDBDCDBBDDDDDDDDDDD>BBDD3@BC	YT:Z:UU
+42/1	4	*	0	0	*	*	0	0	CCCAGATACACGCAGACTCATGATACCTCACCTTCATTGTGACTTTTCATTGCCCGAGCCAATCTACACTGTTGTTGAGAGACGAGCAAAGACATCCCGCA	CCCFFFFFHHHHHJJJJJJIJIJJJJJJJJJJJJJIJJJHIJIIJJJJIJJJIIIJGJJJJJIHHHHFFFFFFFEAEEEDDDDDDDDDDDDDDDDD>ABBB	YT:Z:UU
+44/1	4	*	0	0	*	*	0	0	GGTTTTGCATTTCTCTGTGACCCAATGGGAAGCACAGCTTCTCTCCAGAGTTTGGAGAATTCAGGGTATTCATCTTGCTGCTGCAAAGTTAGTAGTTGGAA	CCBFFFFFHHHHHJJJJJJJJJJJJJJJJIJJJJJJJJJJIJJJJJJJHIIIJJJJJJJJJJIJIJFHHJJJJHHHHHHHFFFFFECECEEDDEEEDDDDC	YT:Z:UU
+41/2	4	*	0	0	*	*	0	0	TGATTGCAAAATAAGAATCCGGTTTAGTCTGATTGTATACAGGTTGCCAAATTCAATGTTCACCTTCTATGATGACAACCGTTTGTCTATACTTTTCATAA	CCCFFFFFHHHHHJJJJJJJJJIJJIJHGIIHJJJHGIJJJJJHHIJJJJJIJJJJIIJJJIJJJJJJJJJJJJJJHHHHFFDEDEEDFEDDEDEDDDEEC	YT:Z:UU
+18/1	4	*	0	0	*	*	0	0	CTTGTCTAACATATAAGGGCCCTGGGTTCGATCCCCAGTACTGAAAGAAAAAAAAAAAAAAAAAAGCAAACAAGAAATACAAAACTGATAATACTAGAAAC	CCCFFFFFHHHHHJJJJJJJJJJJJJHIJJJJJJJJJJHIJJJJJJJJJJJJJJHFBDDDDDDDD>8ACDADD<CAC>C>CDCBBC:A>CC:>>>>:ACD9	YT:Z:UU
+43/1	16	chrUn_gl000220	158975	255	46M3D5M1I5M5D44M	*	0	0	TGGTTCCCTCCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCC	BDBBDDDDDDDDDDDDDDDDDDDEDDDDDBBDDDDDDDDDDFFFFFHHGHHJIHGJJJIJJJJJJJJJJJJJJJJJJJJJJJJJJJJJHHHHHFFFFFCCC	AS:i:-57	XN:i:0	XM:i:3	XO:i:3	XG:i:9	NM:i:12	MD:Z:46^CCC3G2C0C2^CGCCA44	YT:Z:UU
+43/2	4	*	0	0	*	*	0	0	GAGTCTAACGCGTGCGCGAGTCAGGGGCTCGTCCGAAAGCCGCCGTGGCGCAATGAAGGGGAAGGGCCCCGTCCGGGGGCCCGAGGGGGGATCCCGAGGCC	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJIG@DHEEHFDDEDDD@?BDBBBBDDDD:@C&057@DBBDD@>9@9<BDDD@B<<)9>9@BD&0?CA<09>B<	YT:Z:UU
+45/1	4	*	0	0	*	*	0	0	CGGTGCGGCTACCTTGTTACGACTTAGTCCCAATCGCCAATCCCACCTTAGACGGCTCCCTCCCAAAAGGGTTAGGCCACCGGCTTTGGGTGTTACCAACT	CCCDFFFFHHHHHJJIJJJJJJJJJJIJJJJJJIJGIJJJIJJJGJJJJIJJJHHFFFDDEEDDDDDBDDDDDDDCDDDDDDDDD>CDDD8<>BDDDDDBC	YT:Z:UU
+44/2	4	*	0	0	*	*	0	0	TAGTGGTACTAACATCCCTCCTACCATCAATGTCCAAGTGATTTTAATTTATAGTCTCTTAGTTTAACACTGTTCTTAGAAAAAAATTGGCAAACTCAAGC	BCCDFFFFHHHHHJJJJJIJJIJJJJJJJJJJIIIJJJGHIJJJJJJJJJJIJJJJJJJJJJIJJJJIJJJIIJJJJJJIJIGHFDEEEDDDDDBDDD>CA	YT:Z:UU
+46/1	4	*	0	0	*	*	0	0	ACCACCACTGCCCAACAGCAAATTCTTAAAATAAAGTTAGTTATATTCAGTCATTAGCCAGAAACTTAGTATATCTGACTCCTTCCTATTTATGATCAAGC	CCCFFFFFHHHHHJJIJIJJJJJJJJJJJJJJJJJJGHIJJJJJJJJIIJHIJJJIJIJGIGGIJJJJIGIIJJJGHGHGHHHHFFFFFFFFEEDEEDCDC	YT:Z:UU
+46/2	4	*	0	0	*	*	0	0	GCAAAACTGGTTTTAGTTTCTTATAATGAATTGGACAATGTTTTTGAACGTCAGAAGAGGCAAGTGTCAAACTTTTCTCCCTTTAAAAAAAGCCCTAAGTG	CCCFFFFFHHHHHJIJHJJJJJIJJJJJJJJJIJJJJJJJHJJJJJJJJJJIIJJJJJJJJJIJJIJJJJHHHHHHFFFFFFEEEEEEDBBDBDDDDDDDD	YT:Z:UU
+47/2	4	*	0	0	*	*	0	0	CGCTTTGCCATCGATGCCGGTGCCCTGCACATTTTCGACAAGGCGACCGGTAGGCGGCTCTAGCCTTCCAGCAGATAGTCCGGCATGTCCTTCAGCCGGGC	CCCFFFFFHHHHHJJJJJJJGIJJJIJJJJIJJJJIJIJJJJJIJJIJFFDDDDCBDDDDDDDDDDDDDDDDDDDDDCDECDDDDDDDDECCCDDDCB>B5	YT:Z:UU
+48/2	4	*	0	0	*	*	0	0	CACCTGATGGCTTTGCTCATGTCCTTCACTGGAAGAACTTGTGATGTCCCCTTTTGTCCTGATTCCTTTACTGAATCACTCTGTCACTGCCCTTTTCTGCT	CCCFFFFFHHHHGJJJJJJJJJJJJJJJJIJJJJJJJJJJFHIJJJJJJJHIIJIJJJJJJJJJJJJJJJJJJJJJJIJHHHHHHC@DDFFDEEEEEDDDC	YT:Z:UU
+36/2	0	chrUn_gl000220	115012	255	37M3D11M4D53M	*	0	0	CCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCCGAAACGATC	BCCFFFFFHFFHFJJIJIJJIJJJJJIIGHGIJIIHIIDIJJHIIIJIIJBHDBDFFEEEEDDDBBDDDD@BDDDDEDDCCBACACCCDDDDDDBBDDBDD	AS:i:-60	XN:i:0	XM:i:5	XO:i:2	XG:i:7	NM:i:12	MD:Z:37^CCC3G1A0C0C0C2^GCCA53	YT:Z:UU
+45/2	4	*	0	0	*	*	0	0	GTTAAGTCCCGCAACGAGCGCAACCCTCGTTCTATGTTGCCAGCACGTTATGGTGGGGACTCATAGGAGACTGCCGGGGTCAACTCGGAGGAAGGTGGGGA	CCCFFFFFHHHHHJJJJJIJJJJJJJJJJJJJJJJJIJJJJGIJJHHHGFBEFDEDDDBDDDDDDDDDDDDDDDDBDDD9@DDCDDDDDDDDBBD9@BDD?	YT:Z:UU
+47/1	4	*	0	0	*	*	0	0	GCCGAAGGACATTTCCATGCTGGCGCCGCTGACCGAGCGTGAATTGCTCTACCGGCTCCTGACGGGCGAACAGTCCGGCCGCCTGCGCCAGATCGCCCAGT	CCCFFFFFHHHHHJIJJJJJJIJJJJJJJJJJJJIFIJJHHHHHFFFFFFEEEDDDDDDDDADDDDDDBDDDCAC@?B@BBBD9??55<<91>?8?BD084	YT:Z:UU
+50/1	4	*	0	0	*	*	0	0	TTCCTCTCGCCGGGAGGCGCGTCCCGGGGTGTCCTCGTCGCGCGCGTGCGCGCCCGGTGCGCACCAGCGTCGGCGTCCGGAGTCGTGGAAGGAGGGGGTCG	BBBDFFFFHHHHHJJIIJJJJJJJJJJJF<ACDDDDDDDBDDCDDD@BDDDDDDDDD7@BDDDDDDDDDDD>BDDDDDDDDDBCDDDDBCDD<BDDDD55<	YT:Z:UU
+49/1	4	*	0	0	*	*	0	0	CTGGTATGTGGTTGCTGGAACTGTACCCCAGTCCTGCAGAAGAGCAGCCAGTGCTCTTAACCACTGAACCATCTGGCCCTTTAATTCATTTTGAGCTGGTT	CCCFFFFFHHHHHJJJJJJJJJJJGJJJJJJJJJJJJJJJJJJJJJJIJJJIIJJJJJJJJJJJJJIIJJHHHHHFFFFDEEEEEEEDEEEEDDDDDDDDD	YT:Z:UU
+51/2	4	*	0	0	*	*	0	0	AGAGAAAGCAAGAGTCACAGGCTGGGAGAACCCCATCTTGGACCCTCTTTCCTGAGAACATCCACATTCCACCCAGGTTGCCACTCCCGTGGGACCCAGGT	@CCFFFFFHHHHHJIJIJJJIIJJJJIE:C?CGGIIIIIGJIIJJJIIIJJJJJJJJ=8=@CHIGHHHEHFDFFDDD@CCBDDDDDDDDDDDDDDDDBDD>	YT:Z:UU
+52/1	4	*	0	0	*	*	0	0	GCCTTGACATTTGTAACTCTTGTATCAAGTTTGAAGAAGAGGACTTTATAAGTTATTCTTTCTCTGAGATGGATGCTTTTCCAAATTATCACAGCTAATGA	CCCFFFFFHHHHHIJJGIIIJJIIJIJJJIJJJIJJIJJIJJGIJJJIJJIIIJJJJJJJJJJJJJJJJJJJJJJJIJJJJJJJHHHGHHFFEFFECEEDC	YT:Z:UU
+53/1	4	*	0	0	*	*	0	0	GTCGCGCCGTGGGAGGGGCGGCCCGGCCCCCCCCGCCCCGGTTTGACCCGGGGGAAGGGGGGGGGGGCGCGGCGACGGGGATCAGGCTTCCTCGGCCCCGG	CCCFFFFFHFHHHJJJJJJJHFDDDDDDDDDDDD69BDD05)07<>CCD58BDD-5?BBDDDD&5&0)0&5>-5&&09>&&)8<@C38CCC@>@52<<@99	YT:Z:UU
+50/2	4	*	0	0	*	*	0	0	GGGAGCGGCGACGACGGCGGGAGCGGCCTCAAGAGTCTGAGGAGAGCCAGGAGGCTCCCCTACGGGCCCCACCGCCCCCCGACCCCGGGGGCGGAGGGGCG	CCCFFFFFHHHHHJJJIIJDD@8@BBBB3?CC>><3>A>ACDCBA<ADDDAB<@<BBBDDCB?7BB@B@DBBDDDBDDDD99B@BD9<@BD-9@5<@@D##	YT:Z:UU
+48/1	4	*	0	0	*	*	0	0	CCAGCCTGTTGTGATAGGTTTTGAGGGGGTTATGGGTAAGAAGTTAACTGCTTGTTCAATTTTATGTATGACTTACACTGAAGTGAAGGAATGCTGGATAA	CCCFFFFFHHHHHHJJJJDHIJIGHIJJJDHFHIJJBGEGHGGHIJIJJCHIJJGHHHFHFFFFFFFFEEEEEEDCDDDDDDDCDDEDCDDDDDDDDDDDC	YT:Z:UU
+52/2	4	*	0	0	*	*	0	0	CTTAAATGAACATAAATATATAAAGTCTTTTTGTTTGTTTTGTATTTAGTAGAGTTTAAAATCTTGGCTCTTACTGTGGTACAAGCCCAAAATAAGAACAA	@CCDDFFFGHHHHIJIJIJJJJJJJJJIJJJJJIJJJJJJJIGGIJJJJIFGGIHIJJHIJJIJIGIJJJJJJIJJHIIIJHHGHHFFDFECDEEEDDDD;	YT:Z:UU
+51/1	4	*	0	0	*	*	0	0	TGGGCATCGAGGCTTGGCACTCAGATAGGAGGAGGAGGAATGGCTTTCTCCTGTTTTCTCTGGCATCACCCCTGCTGCCAGTCTCCTTTGATCCTGCTGCT	CCCFFFFFHHHHHJJJJJJJJJJJIJIIJGIJGHJDFHGIIIIJJJJJIJJJIIJJJHHHHHHFFFFFEEEDDDDDDDDDDDDDDDDDDDCDDDDDDDDDD	YT:Z:UU
+49/2	4	*	0	0	*	*	0	0	GCATGGTCATGGTATCTCTTCACAGCAATAAAAACCCTACTAATAACTAAGACAATCCAATCTCTTATCTGTCATAAAACCAGCTCAAAATGAATTAAAGG	@C@FFFEEHHHHHFIJJJJJJJJJJJJJJIJJJIJIJJJJJJJJIJJJJJJJJJHJJJJJJJJJJJJJJJIJJJJJIHHHHHFFFFEDEEEEDDEEDCCDD	YT:Z:UU
+53/2	4	*	0	0	*	*	0	0	CCGGAGGATTCAACCCGGCGGCGCGCGTCCGGCCGTGCCGGTGGTCCCGGCGGATCTTTCCCGCTCCCCGTTCCTCCCGACCCCTTCACCCGCGCGTCGTT	@BCFFFFFHHHGGJJJJJJJJJJHFDDDDDDDDDDDDDDDD@BDBDDDDDBBD@DDDDDDDDDDDDDDDD@BDDDDDDDBDDDD?BD:ACDBBBDDDDDD<	YT:Z:UU
+55/1	4	*	0	0	*	*	0	0	GATCCTCTTGCCTCTGGTGCTGCCAGCACCACCACCCAGCAAGAAAACTTTTAAAGACAGTCTCATTAGCTGAGCATGGTCTTGAACTTAACATATAGCTA	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJIJJJJIJJJJJJJJJJJJJJJJJJJJJIIJJJJJEIHHHHHFCEFFCCEEEEDDDDDDDDDDEEDCCE	YT:Z:UU
+56/1	4	*	0	0	*	*	0	0	CTGCGGTTATGAGTACGACCAGGCGTGAGACTAATGCGTTCCGTGGGATTTTCAAGGACTGTCAGGGACGCACCGGACACAACACAAGTGTTGTGCTCTAC	@CCFFFFFHHHHHIHJJJJJIJJJJIJIIHJJJIJJIJIJJIJFHIJJJHHHHHHFFFFFDEEEEEDDDDDBDDBDBBDDCDBDDDD?>ACDBCDDDDDDC	YT:Z:UU
+56/2	4	*	0	0	*	*	0	0	CGCTCAAGCGTGTTACCCATACCTCACCGTTGATGTTAAAGTGATGCATCAACGAGTAGGCAGGCGTGGAGGTTGTGAAGAAGCCTTGGCAGTGATGCCGG	CCCFFFFFHHHHHJJJJJIJIJJJJJJJJJJJJJJIIEIIIGGGIJIIJIHGIIIJGHIJJJJHHHBDDACD@BD@BDEDDDDDDDCCDDDD<CCCDDDDD	YT:Z:UU
+40/1	4	*	0	0	*	*	0	0	CCTCTCTCTCTCTCCCTCTCTCTCTCTCTGCCCTTCTGTGCCTCTACTACCCTCTTAACTCCCCTCTCCGTGCCCTGAATAAACTCTATTCTATACTGTAC	BBBFFFFFHHHHHJJJJJJJJJJJJJJJIIJIGGHIJJGHGIJJIEIJGIGIJIIJHFHIJIIIJJJJJIAHHFFEDDCEECCEEDDDDEEEEEEDECDDD	YT:Z:UU
+57/2	4	*	0	0	*	*	0	0	TGCACCGGCGGCGGGGTGAGGTGTAGGGGGATGGGCGGGACGCTGGGAGGAGCCAATGACCCGTGCGAACGTTCAGGTCGATGGCAGCTGCGACGTGCGTC	#####################################################################################################	YT:Z:UU
+58/1	4	*	0	0	*	*	0	0	AGCAGGAAGTAGCCCAAGAGATACGACGCCCTCATTCCCTTTTCATCATCGGCTTCCCCTTCCCTTTTTCTTTCCTTCTCTTTATAATATAAAAAAGGGAG	??@DADDDDHHHAGHIBB;FACDB@DHIF@??G>DBBEBFG?C=@=F4==F4ADGHCHH;;3@DC>;;>CCCCC(-(5:A>3:;>>;BCDDE;@><<BBBB	YT:Z:UU
+59/1	4	*	0	0	*	*	0	0	GCCCAAGGTACACGCCCCTCAGCTCTGCCCTACATAGAAATGGGGCTTGGACAGGAGGTACTGTCCAGTCTAGGTGCTGGCCAGCGGAAAAAATCCTCAGT	BCCFFFFFFFHHHJJJJJJJJJJJJIJJJJJJJJJIIJJJJJJJJJJJJJIIIJJIIHAEHFFFFFFECCECCECACDDD?ACB?=B9BDCA<8>CCCA>:	YT:Z:UU
+59/2	4	*	0	0	*	*	0	0	GTGTTTATGAATCAAGGCAGGCTAGGTGGAACCTGTTCAAGGCCGCTGGATGCTGAAGGATGGGATGGTCGAGTCACGTGTTCATCGCCTCCCTTTTCTCC	BBBDFFFFHHHHHJHIJJJJJJJJJJCFHIJJJJIGHGGIJJJJJJJIJHHHIJJJJJJHHHHFFFFDDEDDDDDDDDBB?BDDEEDDDDDDDDDDDDDDC	YT:Z:UU
+55/2	4	*	0	0	*	*	0	0	TCTGACCTGTTATTCTCCATATATTTATATTTATATTTCTCCATATATTTATGAATATTATATGAGGTATATAGTGTGGAAGTAAAGGCCTATAATTCCAG	CCCFFFFFHHHHHJJJJJJJJJJJJJIJIJJJIIGIJJJJJIIJJJJJJJJIJIJJJJJJJJJJIJJFHIJJJJIJIIJJJJJJJIIIJJHHHHFGFFFFB	YT:Z:UU
+54/1	16	chrUn_gl000220	158967	255	54M3D11M4D36M	*	0	0	CTAGTAGCTGGTTCCCTCCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTC	DDDDDDDDDDDDBDBDDDDDDDDDDDDDDDDEDDDDDDDDDDDDFFHHJJJJJJJJJIHGIHGJJJJJJJJJJJJJJJJJJJJJJJJJHHHHGFFFFFCCC	AS:i:-58	XN:i:0	XM:i:5	XO:i:2	XG:i:7	NM:i:12	MD:Z:54^CCC3G1A0C0C0C2^GCCA36	YT:Z:UU
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output.fastq	Mon Mar 21 17:33:26 2016 -0400
@@ -0,0 +1,20 @@
+@34/1
+GTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCCGAAACGATCTCAAC
++
+DBDBDDDDDEDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD?FFHHHJJJJJJIHJJJJIJJJJIJJJJJJIIJIIJJJHJJIJJJJIJHHHHHFFFFFCCC
+@37/2
+GCGCAGTCCGCCCGGAGGATTCAACCCGGCGGCGCGCGTCCGGCCGTGCCGGTGGTCCCGGCGGATCTTTCCCGCTCCCCGTTCCTCCCGCCCCCTTCACC
++
+@@CFFFFFHHHHGGJHIJHIJGHHJJJAFEBEE?@BBDDBD6BDD5:B8B<88@@5@C?@B7<B87;CAAAADD<&&&8590<CCCDCBD&&0590<?CCC
+@43/1
+TGGTTCCCTCCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCC
++
+BDBBDDDDDDDDDDDDDDDDDDDEDDDDDBBDDDDDDDDDDFFFFFHHGHHJIHGJJJIJJJJJJJJJJJJJJJJJJJJJJJJJJJJJHHHHHFFFFFCCC
+@36/2
+CCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTCTTGGGGCCGAAACGATC
++
+BCCFFFFFHFFHFJJIJIJJIJJJJJIIGHGIJIIHIIDIJJHIIIJIIJBHDBDFFEEEEDDDBBDDDD@BDDDDEDDCCBACACCCDDDDDDBBDDBDD
+@54/1
+CTAGTAGCTGGTTCCCTCCGAAGTTTCCCTCAGGATAGCTGGCGCTCTCGCAGAGACTCGAAACCCGCAGTTTTATCCGGTAAAGCGAATGATTAGAGGTC
++
+DDDDDDDDDDDDBDBDDDDDDDDDDDDDDDDEDDDDDDDDDDDDFFHHJJJJJJJJJIHGIHGJJJJJJJJJJJJJJJJJJJJJJJJJHHHHGFFFFFCCC