| Next changeset 1:9cb93e3df272 (2020-09-18) |
|
Commit message:
"planemo upload for repository https://github.com/QFAB-Bioinformatics/metaDEGalaxy/tree/master/reheader commit 6783cd68521863b34f8e77cbb7ba404700c72313-dirty" |
|
added:
reheader.py reheader.xml test-data/F3D0_R1.fastq test-data/F3D0_R2.fastq test-data/test_output.txt |
| b |
| diff -r 000000000000 -r 331fd79a9341 reheader.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/reheader.py Thu Sep 03 01:32:31 2020 +0000 |
| [ |
| @@ -0,0 +1,85 @@ +#!/usr/bin/env python +import sys +from Bio.Seq import Seq +from Bio import SeqIO +from Bio.SeqRecord import SeqRecord +from os.path import basename +import os +import re +import argparse + +parser = argparse.ArgumentParser( + description="reformat the read name (header) by appending the sample name to the read name\n" + + "Example:\n python reheader.py -n F3D0_R1.fastq -i test-data/F3D0_R1.fastq -o test-data/test -l mylog -d test-data/") +parser.add_argument("-v","--version",action="version",version="%(prog)s 1.0") +parser.add_argument("-n","--samplename",dest="samplename",default=False,help="input sample name") +parser.add_argument("-i","--input",dest="inputfile",default=False,help="input filename in FASTQ format") +parser.add_argument("-l","--log", dest="logfile",default=False,help="output log file") +parser.add_argument("-o","--outfile",dest="outputfile",default=False,help="output filename") +parser.add_argument("-d","--outdir",dest="outputdir",default=False,help="output directory") + + +if(len(sys.argv) == 1): + parser.print_help(sys.stderr) + sys.exit() + +args = parser.parse_args() + + +filename = args.samplename +infile = args.inputfile +str_to_add = os.path.splitext(basename(filename))[0] +outfile = args.outputfile +outdir = args.outputdir +logfile = args.logfile + + +rdict = { + '_R1': '/1', + '_R2': '/2', + '_1': '/1', + '_2': '/2', +} + +rdict_remove = { + '_R1': '', + '_R2': '', + '_1': '', + '_2': '', +} + +def makesubs(s): + for pattern, repl in rdict.items(): + pat1 = pattern +'_?[A-Za-z0-9]+$' + pat2 = pattern + combined_pat = r'|'.join((pat1, pat2)) + s = re.sub(combined_pat, repl,s) + return s + +def makesubs_remove(s): + for pattern, repl in rdict_remove.items(): + pat1 = pattern +'_?[A-Za-z0-9]+$' + pat2 = pattern + combined_pat = r'|'.join((pat1, pat2)) + s = re.sub(combined_pat, repl,s) + return s + +def appendStringToSequenceHeader(inputfile,header_to_add): + records=[] + for seq_record in SeqIO.parse(inputfile, "fastq"): + header =seq_record.id + header = "{0}".format(header) + "_" +header_to_add + record = SeqRecord(seq_record.seq,id=header,description="") + record.letter_annotations["phred_quality"]=seq_record.letter_annotations["phred_quality"] + records.append(record) + return records + +str_to_search = makesubs_remove(str_to_add) +str_to_add = makesubs(str_to_add) +final_records=[] +outlogfile=open(os.path.join(outdir,logfile),"w") + +final_records=appendStringToSequenceHeader(infile,str_to_add) +outlogfile.write(str_to_search) +SeqIO.write(final_records, outfile , "fastq") +outlogfile.close() |
| b |
| diff -r 000000000000 -r 331fd79a9341 reheader.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/reheader.xml Thu Sep 03 01:32:31 2020 +0000 |
| b |
| @@ -0,0 +1,69 @@ +<tool id="rename" name="reheader" version="1.0.0" force_history_refresh="True"> + <description>Rename sequence header in FASTQ file</description> + <requirements> + <requirement type="package" version="1.74">biopython</requirement> + </requirements> + <version_command> + python ${__tool_directory__}/reheader.py --version + </version_command> + <command detect_errors="aggressive"> + python ${__tool_directory__}/reheader.py -n $inputfile.display_name -i $inputfile -o $output -d $log.files_path -l $log + </command> + + <inputs> + <param format="fastq,fastq.gz,fastq.bz2" name="inputfile" type="data" label="Short read data from your current history"/> + </inputs> + <outputs> + <data name="output" format="fastq" label="${tool.name}.${inputfile.display_name}"/> + <data name="log" format="txt" label="${tool.name}.${inputfile.display_name}.log"/> + </outputs> + <tests> + <test> + <param name="inputfile" value="F3D0_R1.fastq"/> + <output name="outout" file="test_output.txt"/> + </test> + </tests> + + <help> +**what it does** + +Renames sequence header in FASTQ file by appending SAMPLEFILENAME/1 to the end of the sequence header for the forward read or SAMPLEFILENAME_2 to the end of the sequence header for the reverse read using python_ and Biopython_ + +.. _python : https://www.python.org/downloads/release/python-370/ +.. _Biopython : https://biopython.org/ + +----- + +Forward Read +------------ +**Input** + +- **sequence header** : @M03033:102:000000000-ATN0P:1:1101:9793:1042 + +**Output** + +- **sequence header** : @M03033:102:000000000-ATN0P:1:1101:9793:1042_Sample1/1 + +Reverse Read +------------ + +**Input** + +- **sequence header** : @M03033:102:000000000-ATN0P:1:1101:9793:1042 + +**Output** + +- **sequence header** : @M03033:102:000000000-ATN0P:1:1101:9793:1042_Sample1/2 + +----- + +========= +Resources +========= + +**Wrapper Authors** + +QFAB Bioinformatics (support@qfab.org) + </help> + +</tool> |
| b |
| diff -r 000000000000 -r 331fd79a9341 test-data/F3D0_R1.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/F3D0_R1.fastq Thu Sep 03 01:32:31 2020 +0000 |
| b |
| b'@@ -0,0 +1,4000 @@\n+@M00967:43:000000000-A3JHG:1:1101:18327:1699 1:N:0:188\n+NACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+#>>AABABBFFFGGGGGGGGGGGGGGGGHHHHHHHGGGHHHHHGHGGGGGGGHGGGGGGHHHHHHHHHHGGGGGHHHHGHGGGGGGHHBGHGDGGGGGHHHGGGGHHHHHHHHGGGGGHG@DHHGHEGGGGGGBFGGEGGGGGGGG.DFEFFFFFFFDCFFFFFFFFFFFFFFFFFFFFFFFFFFDFDFFFEFFCFF?FDFFFFFFFFAFFFFFFFFFFFBDDFFFFFEFADFFFFFBAFFFA?EFFFBFF\n+@M00967:43:000000000-A3JHG:1:1101:14069:1827 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+3AA?ABBDBFFBEGGEGGGGAFFGGGGGHHHCGGGGGGHFGHGGCFDEFGGGHGGGEGF1GGFGHHHHHGGEGGHHHHHFGGGGGGHHHHHGGGGCDDGHHGGGFHHHHHHHHCD@CCHGGGGHEHGGG@GFGGGGGGG@BGGGEGCEBFFFBFFB;9@EFFFEFFFFFFFFFFFFAFBBBFFFFFBBBFFFFBBBFFFFFFFFFFFBBBBBBBFFFFFFFFFDDFAFFFFF.AF9/FBBBBB.EAFFE?F\n+@M00967:43:000000000-A3JHG:1:1101:18044:1900 1:N:0:188\n+TACGGAGGATGCGAGCGTTGTCCGGAATCACTGGGCGTAAAGGGCGCGTAGGCGGTTTAATAAGTCAGTGGTGAAAACTGAGGGCTCAACCCTCAGCCTGCCACTGATACTGTTAGACTTGAGTATGGAAGAGGAGAATGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGATTCTCTGGGCCAAGACTGACGCTGAGGCGCGAAAGCGTGGGGAGCAAACA\n++\n+BA@BBBABBFFFGGGGGGGGGGHGGGGGHHHHGHHHGGGGHHHGGGGGGGGGHGGGGGGFHFHHHFHHHHHHGFHHHHHHFHGGFGHHGGHGGHHFHHHHHHHHHHHHHHGHHHGHDHHHHHHHHHGHDGFFHHHDHHHHHHHHHHHHFHHHDGHHGGGGGHHHHHGGGGEGGGGGBF0BGGGGGGGGFFFGGGDGFFFFBCFBFFBFFBEEFFEFFFFFFFFBDA9EEFFF>DFFFFDDFFFADDFFFF.\n+@M00967:43:000000000-A3JHG:1:1101:13234:1983 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGAAGATCAAGTCAGCGGTAAAATTGAGAGGCTCAACCTCTTCGAGCCGTTGAAACTGGTTTTCTTGAGTGAGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCTCAACTGACGCTCATGCACGAAAGTGTGGGTATCGAACAG\n++\n+BCBCCCCCCFFFGGGGGGGGGGHGGGGGHHHHHHHGGGHHHHHGHGGGGGGGGGGGGGHHHHHHHHHHHGGGGGHHHHHHHHHGGHHHHHHHHGHHHGHGGGGGGHHHHHHHGHHGHHHGHHHHHHHGGGGGHHHHHHGGGGGHHGGGGHGHHHHGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFAFFFFFFFFAFFFFFFFFFFFFFFBFFFFFFFFFFEFFFF\n+@M00967:43:000000000-A3JHG:1:1101:16780:2259 1:N:0:188\n+TACGTAGGTGGCAAGCGTTATCCGGATTTACTGGGTGTAAAGGGCGTGTAGGCGGGGACGCAAGTCAGATGTGAAAACCACGGGCTCAACCTGTGGCCTGCATTTGAAACTGTGTTTCTTGAGTACTGGAGAGGCAGACGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGTCTGCTGGACAGCAACTGACGCTGAGGCGCGAAAGCGTGGGGAGCAAACG\n++\n+BBCCBFFCDFCFGGGGGGGGGGHGGGGGHHHHHHHHGGHFHHHGGGEGGHHHGGGGGGGEFGGGHHHHFHHHHGHHHHG/FGGGFGHFHHHGA?CGEHGGGHGHHFBGHHHGGHHHHDGHHFHFHHHHHEC0DDDEDAFFDDEGHHHHHEHHHHHEGGAFFGGGGBDGGGGGGFFFGGFFG.FFFFFFFFFFFFFFFDEFF@=FFFFFF.AFFEFFFFFFFDFFFBFF.:;F-@FFBDAFAAC.;9?FFE.\n+@M00967:43:000000000-A3JHG:1:1101:19378:2540 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGAAGATCAAGTCAGCGGTAAAATTGAGAGGCTCAACCTCTTCGAGCCGTTGAAACTGGTTTTCTTGAGTGAGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCTCAACTGACGCTCATGCACGAAAGTGTGGGTATCGAACAG\n++\n+CCCCCCCCCFFFGGGGGGGGGGHGGGGGHHHHHHGGGGHHHHHGFFGGGGGGGGGGGGHHHGHHHHHHHGGGGGHHHHHHHHHGGHHGGHHHHHHHHHHGGGGGGGHHHHHHHHHGHHHHHHHHHHHGGGGGGHHHHHGGGGGHHGGGGHGHHHHGCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFDFFFFFFADFFFFFFFFAFFFFFFFFF\n+@M00967:43:000000000-A3JHG:1:1101:17674:2779 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGACTCTCAAGTCAGCGGTCAAATCGCGGGGCTCAACCCCGTTCCGCCGTTGAAACTGGGAGCCTTGAGTGCGCGAGAAGTAGGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCCTACCGGCGCGCAACTGACGCTCATGCACGAAAGCGTGGGTATCGAACAG\n++\n+ABB@BABDAFFFGGGGGGGGGGHGGGGGHHGHHGHGGGHFDHHGHEGCGGGGGGGFEECGHHHHHHHHHGGGGGHHHGGGGGGGGGHGHGHGGGGHHHGGGGDGGEHHHGFHHGGGFHHGHFHFHCGGGGGGGFGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF?FFFFFFFFFFFFFAFEFFFFFFFFFFBFFDFFFFFFFFFDBF/AEFFFF\n+@M00967:43:000000000-A3J'..b'FFF/;B?BBFFFFFFBDFFFBD9ABBBEDEDFFFBF\n+@M00967:43:000000000-A3JHG:1:1104:14390:26638 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+BABBBBBDBFFFGGGGGGGGGGHFGGGGHHHHHHHGGGHHHHHGFGGGGGGGGGGGGGGHHHHHHHHHHGGEGGHHHHHHGGGGGGHHHGHGGGGGGGHHHGGGGHHHHHHHHGGGGGHGGGHHGGGGGACEFGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFF9.:BFFFFFDFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF.:AFBFFADAFA?EFFFFFF\n+@M00967:43:000000000-A3JHG:1:1104:16844:26739 1:N:0:188\n+TACGTAGGGAGCGAGCGTTATCCGGATTTATTGGGTGTAAAGGGTGCGTAGGCGGTAATACAGGTCTTTGGTATAAGCCCGAAGCTTAACTTCGGTAAGCCAGAGAAACCGTATAACTAGAGTATGGTAGAGGCAAGTGGAATTTCTAGTGTAGCGGTAAAATGCGTAGATATTAGAAGGAACACCAGTGGCGAAGGCGACTTGCTGGGCCAAAACTGACGCTGAGGCACGAAAGCGTGGGGAGCAAATAG\n++\n+BBBBCFFCCFCCGGGGGGGGGGHGGGGGHHHHHHHGGGHHHHHGHGGGGGGGHGGGGGHHHHHHHHHHHHHHHHHHHHHHGGGGGHHHHHHHHGGGHGHHHHHHHHHHGHGGHHHHGHHHHHHHHHHHHHHHHGGHGHHHHHHHHHHHHFHHHHHGGGGGHHHHHGCADGHHHHHHHHHGHHHGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFFBDADFFF?DFFFFFFF\n+@M00967:43:000000000-A3JHG:1:1104:17020:27386 1:N:0:188\n+TACGTAGGTGGCAAGCGTTATCCGGATTTATTGGGCGTAAAGGGAACGCAGGCGGTCTTTTAAGTCTGATGTGAAAGCCTTCGGCTTAACCGGAGTAGTGCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAACTCCATGTGTAGCGGTGAAATGCGTAGATATATGGAAGAACACCAGTGGCGAAAGCGGCTCTCTGGTCTGTAACTGACGCTGAGGTTCGAAAGCGTGGGTATCGAACA\n++\n+ABBBAFFBFFBFGGGEEGGGGGHGGGGGHFFHHHHCEEFGGHHGHEGGGGGGCGCEEGHHHFGHHHHHHHGGFEHHHHHHHHHGGGGHHHGEGECFFHEGFHGHHH?GHHHHHHHGAFHHHHHGHHHHHHHHHHFGHGFFFHHHHHHHHHHHHHHGGGGGGHHGHGGGGFEHGHFCGHFFHHHHHHGGFFGGGGGGFGBEGGFFFFFFF/BFFBFFFFFBFFDFDAFFFFBCEFFFFDDDCFAFFBDFFF9\n+@M00967:43:000000000-A3JHG:1:1104:14705:27440 1:N:0:188\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGGATGCCAAGTCAGCGGTAAAAAAGCGGTGCTCAACGCCGGCGAGCCGTTGAAACTGGCGTTCTTGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGGGGCACGAAAGCGTGGGTATCGAACAG\n++\n+AAAAAAAD11BFCGGGGGG?GGFGCGGGHHHHHGH//BGHHHBEBC?EAEEGHGG//?>CCGGFHHFHEE@E?FHHBE/C/<A@CHHHHHGCCC/<<--<ACGGGGHHFCBF.-:;CEG0BB0BF?EGGGGGG0BBBFB@@@AFBBBF-:A-9FF=@A@?/9BBBB/B/B//F/99A-99=F-;BFA@-A?/;9@@BEFBFFBFBF-99-;9-9B-//99-9B---9999-;:B-/;AB@@-/;/BFFE@/\n+@M00967:43:000000000-A3JHG:1:1104:19959:27794 1:N:0:188\n+TACGTAGGGAGCAAGCGTTGTCCGGAATTACTGGGTGTAAAGGGAGCGTAGGCGGGATAGTAAGTTGATTGTTAAAACTATCGGCTCAACTGGGAGATGCCATCAAAACTGCAATTCTTGAGTGAAGTAGAGGTAGGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAATGGCGAAGGCGGCCTGCTGGGGCTCTACTGACGCTGAGGCACGAAAGTGTGGGGAGCAAACAG\n++\n+AAAAAFF?11AAAGGGGGGCGGHGGG?FHH1FHHHHAEHHHHHGEGGEEGCGGGGG/E01BGHHHB122@BG2BGH11F12>EEGGGHHHB0</0//BFG11B11FFHGGF11121@G2>GHHF11<?<GF1G1?1FGCGGGF0D0D0=GFGHGHGG.:CGCGFB.ECGEGGGBF0F0/A?BF9FF..9/9B@-A@B?---9A:F//B---;-///9//9BAA--AE9-A-9ABA/:/:9B-----;FFFF\n+@M00967:43:000000000-A3JHG:1:1104:14252:27883 1:N:0:188\n+TACGTAGGTGGCAAGCGTTGTCCGGATTTACTGGGTGTAAAGGGCGTGTAGCCGGGAAGGCAAGTCAGATGTGAAATCTGGAGGCTCAACCTCCAAATTGCATTTGAAACTGTTTTTCTTGAGTACTGGAGAGGCAGACGGAATTCCTTGTGTAGCGGTGAAATGCGTAGATATAAGGAGGAACACCAGTGGCGAAGGCGGTCTGCTGGACAGCAACTGACGGTGAGGCGCGAAAGCGTGGGGAGCAACCA\n++\n+BBBABFF>FFB@FFFEGGGGGGHGGGGGHHHHHHGHG2FHGHHDGGGGGHHHHGGGGGGGHHHGFHHHHHHHHGHHHHHHHHHGGHHHHHHHHHHHHHHHHGHHHHHHHHHHHHHHGHHHHFHHHHHHHHGHGGGEFGHGGGGGHHHHHEHHHHHGGGGGEDHHHFGGGGGHHHDFHGFEGGGHHHGHHHGHGGGGFGE?-ACGFFFFFFFFFFFFFFFFFFF.AFFFFFCCCFA?FAADFFF?DFFB...\n+@M00967:43:000000000-A3JHG:1:1104:16540:28173 1:N:0:188\n+TACGTAGGGGGCGAGCGTTATCCGGATTCACTGGGCGTAAAGGGAGCGCAGGCGGCACGGCAAGTCCGGTGTGAAAGCCCGGGGCCCAACCCCGGGACTGCACTGGAAACTGCCGGGCTGGAGTGCCGGAGGGGTAAGCGGAATTCCTGGTGTAGCGGTGAAATGCGCAGATATCAGGAGGAACACCAGTGGCGGAGGCGGCTTACTGGACGGTCACTGACGCTGAGGCTCGAAAGCGTGGGGAGCAAACA\n++\n+AAAABFFBBBBBGGGGGGGGGGHGGGGGHHHHHHHGFEEGGHHGGGGGGGGGGGGGGGGGGGGGHHHGGCGGGHHHHHHHGGGGGGGGHGGGGGGGGGHHHHGHHGFGGGGGGGGGGFGGGGGGGGGGFCBFFFFFFFFFFFFFFFFFFFFFFFFFFF9DAFFFFFFFFFFFFFEFFFFFFFFFFFFFFBFFFFF>BFDFD;@@BFFFF9.=@.;A9/FEFFB?D.FE?EFFBFDD/..9@DF.;.;BB9.\n' |
| b |
| diff -r 000000000000 -r 331fd79a9341 test-data/F3D0_R2.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/F3D0_R2.fastq Thu Sep 03 01:32:31 2020 +0000 |
| b |
| b'@@ -0,0 +1,4000 @@\n+@M00967:43:000000000-A3JHG:1:1101:18327:1699 2:N:0:188\n+CCTGTTTGATCCCCGCACTTTCGTGCCTCAGCGTCAGTAGGGCGCCGGTATGCTGCCTTCGCAATCGGGGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCGCCCACTCGAGCCCGGCAGTTTCAACGGCTGTACGGGGTTGAGCCCCGCAATTTTACCGCTGACTTGGCAGGCCGCCTACGCACCCTTTAAACCCAATAAATCCGGATAACGCCCGCCTCCTCCGT\n++\n+BABBBFFFFFFFFGEGGGGGGGGHGGHHHHHGGGGGGHHHGGHGGGGGGEGHHGHHHHHHHG11EHEGCGGGHHHHGFGEHHHHGHHHHHHHHHFHHGGGCGHHHHHGGGGGGFHG@EGHHHHHHGFGGGGGGAGDGGGGGGFGHHGHHECDGGGGFGGGDFG:D-FFGGGGFCFBDFBFFFF.;--.BBFFFF/9B.;>-@DDF.-9B@DFFF/FFF/AFAB/B/B/FFF.AAFBA--@=F-.AA.;B..\n+@M00967:43:000000000-A3JHG:1:1101:14069:1827 2:N:0:188\n+CCTGTTTGATCCCCGCACTTTCGTGCCTCAGCGTCAGTAGGGCGCCGGTATGCTGCCTTCGCAATCGGGGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCGCCCACTCGAGCCCGGCAGTTTCAACGGCTGTACGGGGTTGAGCCCCGCAATTTTACCGCTGACTTGGCAGGCCGCCTACGCACCCTTTAAACCCAATAAATCCGGGTAAAGCTCGCATCCTCCGG\n++\n+BABBBFFFFFFFGGGGGGGGGGGHGGHHHHHHGGGGGHHHHFHGGGGGGFGHHGHHHHGHGGGGGHGGCGGGHHHHGGGDBGHHHFHHHHHHGHHHHGGGGGHHHEHGGGGGHHHGGGGHGHHHHGGGGGGHHEHGGGGG@FCCGHHHHGHGGGAGFBCEFGGGADGFGFGGFF;ADFFFFFBDFA>.BFFEFFFFFAFF;D..FFAA;BFEFFBBFFFF?EFFFFB/FDC-99B9FB/;AF-ABFFFF.-\n+@M00967:43:000000000-A3JHG:1:1101:18044:1900 2:N:0:188\n+CCTGTTTGCTCCCCACGCTTTCGCGCCTCAGCGTCAGTCTTGGCCCAGAGAATCGCCTTCGCCACTGGTGTTCCTCCTAATATCTACGCATTTCACCGCTACACTAGGAATTCCATTCTCCTCTTCCATACTCAAGTCTAACAGTATCAGTGGCAGGCTGAGGGTTGAGCCCTCAGTTTTCACCACTGACTTATTAAACCGCCTACGCGCCCTTTAAGCCCAGTGATTCCGGACAAAGCCCGCCTCCCCCG\n++\n+>AAABFFFFFFFGGGGGGGGGGGGGGGGGGHGGGGGGHHGHGGCEGEGGGHHFHFGGGGHHGGGGBFHHGHGGHGHGHEHH4FHHHGGGGGGHHHFHGGGGGHH3FHGFGGHHHFFHDFHFHHHHHHHFHHHHDFG2FFFHHHHFHHHHHGHHGHGG/000<CDDHC01<CE/<.C=DDGDDGHBC/CCGHBCGF00CGB-;:E.CFGAAFG?CF009:9E..;//BB/:/9.9-.::B.9---9:..9D-\n+@M00967:43:000000000-A3JHG:1:1101:13234:1983 2:N:0:188\n+CCTGTTCGATACCCACACTTTCGTGCATGAGCGTCAGTTGAGCGCCGGTATGCTGCCTTCGCAATCGGAGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCGCTCACTCAAGAAAACCAGTTTCAACGGCTCGAAGAGGTTGAGCCTCTCAATTTTACCGCTGCCTTGATCTTCCGCCTGCGCCCCCTTTAAACCCAATAAATCCGGATAACGCTCGCATCCTCCGT\n++\n+ABBBBFFBBFBBGGGGGGGGGGGHGGHGGHHHGGGGGHHHHHHGGGGGGGGHHHHHHHHHHGGGGHGGGGGHHGHHGGGGHHHHHGHHHHHHHHHHHGGGGGHHHHGGGGGGHHHGGGGHHHHFHGGGHGHHGHHHHHHHHDGHGHHHFHHGGGFDC<C.AGHC/C/C0GG//:GBB0CF9009??-..9;/0;;B90;.;D.C.B9.:-@-:F/BF//;EE./;::/9.@.A.B:.9.9;-999B99/.;\n+@M00967:43:000000000-A3JHG:1:1101:16780:2259 2:N:0:188\n+CCTGTTTGCTCCCCACGCTTTCGCGCCTCAGCGTCAGTTGCTGTCCAGCAGACCGCCTTCGCCACTGGTGTTCCTCCTAATATCTACGCATTTCACCGCTACACTAGGAATTCCGTCTGCCTCTCCAGTACTCAAGAAACACAGTTTCAAATGCAGGCCACAGGTTGAGCCCGTGGTTTTCACATCTGACTTGCGTCCCCGCCTACACGCCCTTTACACCCAGTAAATCCGGATAACGCTTGCCACCTACG\n++\n+>AABAFFFFFFFGGGGGGGGGGGGGGGGGHHHGGGGGGHHHGFHHHHHHHHBGHGGGGGHHGGGGFHHFHHEFHHHHHEHHFHHHFFGGGGGHHFHHGGGGGFHHHHHHHHHHHHHHHGHHHHHHHFHHHHHHHHHEHHHHHHHHHHHGGGHHHHGGGGHGFHHHFG<CFHGGCDC.EHGHHGFHHGG0CGHHHCD@BFGF?B?DGGGDCEFFFFFFFFB?EFFFF/BFFD;?FFFBD=AEEFB/FFFFF?\n+@M00967:43:000000000-A3JHG:1:1101:19378:2540 2:N:0:188\n+CCTGTTCGATACCCACACTTTCGTGCATGAGCGTCAGTTGAGCGCCGGTATGCTGCCTTCGCAATCGGAGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCGCTCACTCAAGAAAACCAGTTTCAACGGCTCGAAGAGGTTGAGCCTCTCAATTTTACCGCTGACTTGATCTTTCGCCTGCGCACCCTTTAAACCCAATAAATCCGGATAACGCTCGCATCCTCCG\n++\n+BCCCCFFCCFCCGGGGGGGGGGGHGGHGHHHHGGGGGHHHHHHGGGGGCEFHHHHHHHHHHGGGGHGGGGGHHHHHGGGGHEHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHGGGGHHHHHHGGGHGHHHHHHBCGHHGHHFGHHHHHGGGGGC-C...;;CC/;0C/:0:0<0;BFB0:A??-9/0;90;9CF0;9?D...--@-@9.BB0BBB9AEEF0FF009.=...0B.-.9--;.@B0:B;\n+@M00967:43:000000000-A3JHG:1:1101:17674:2779 2:N:0:188\n+CCTGTTCGATACCCACGCTTCCGTGCATGAGCGTCAGTTGCGCGCCGGTAGGCTGCATTCGCAATCGGAGTTCTGCGCGCTATCTATTCATCTTACCGCTACACCACGCCTTCCGCCTACTTCTCGCGCACTCAAGGCTCCCAGTTCCAACCACGGGACGGGGCTGAGCCCCGCGAATTTCCCCCTGACTTGACAGTCCGCCCGCGACCCCTTTAAAACTAAACACTCCGCATTACGCTAGCACCCTCCT\n++\n+1>AA>1BA1D>A1BGEAAE00B0BCEC221BFEAEE?D211AAEAE?/>///B@0>0@FG2/EFC1A//>/BB21B//>/>//0B2B2BGF222B>CBBC/><FGG/<C/?EHCEGC/.><<111F.C@-<</<C0/./;C:/0/:0<CGG....--.99A-A-;--;B/AB--;--9;///9--;-;B///9//99//9---9----;---9///9/B/////:9//------;9;----9//99-A9;\n+@M00967:43:000000000-A3JHG:1'..b'9-99:F/B//BA....;9/B/.9---;--:/9.--9-..;..\n+@M00967:43:000000000-A3JHG:1:1104:14390:26638 2:N:0:188\n+CCTGTTTGATCCCCGCACTTTCGTGCCTCAGCGTCAGTAGGGCGCCGGTATGCTGCCTTCGCAATCGGGGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCGCCCACTCGAGCCCGGCAGTTTCAACGGCTGTACGGGGTTGAGCCCCGCAATTTTCCCGCTGACTTGGCAGGCCGCCTACGCCCCCTTTAAACCCAATAAATCCGGATAACGCCCGCATCCCCCG\n++\n+111AAFFFFFFFGGGGGEGGGGGHGGHHHHHCEFGGFHHHCFHGGGGGGGFHHHHHGHGHHGCGGHGGGGGGHHHHGGGGGBGHHHHHGFHGHGHHHGGGGGHHHHGFCGGGHHHGCCFHHHHHHGGGGGGHGEF<EGGGGGGGGHHHGHFCCG@@9CFGB@GGG@FFEGGG---@AB9///FA?@=-B/FFFFBB@?E-99999A--9@@-AB/;BFBBFFBF9FBBBB@-:-B/9--99@-;@F--;9\n+@M00967:43:000000000-A3JHG:1:1104:16844:26739 2:N:0:188\n+CCTATTTGCTCCCCACGCTTTCGTGCCTCAGCGTCAGTTTTGGCCCAGCAAGTCGCCTTCGCCACTGGTGTTCCTTCTAATATCTACGCATTTTACCGCTACACTAGAAATTCCACTTGCCTCTACCATACTCTAGTTATACGGTTTCTCTGGCTTACCGAAGTTAAGCTTCGGGCTTATACCAAAGACCTGTATTACCGCCTACGCACCCTTTACACCCAATAAATCCGGATAACGCTCGCTCCCTACGT\n++\n+AABBBFFFFFFFGGGGGGGGGGGHGGHHGHGGGGGGGHHHHHGHGHGHHHHHHHGGGGGHHGGGGHHHBFHGHHHHHHFHHHHHHHHGGGGGHHHHGGGCGGGHFFHFHHHHHHGGHHGHHHHHGHHHHHHHHHGHHHHHHHHGHGGFFEHHHHHHGHHFEGGGHHHHHHFFCGGCGHHHHHGCGHHCBGFHB00CFGHGCFDGBA?DGGF?FGGGGGGGGGGGGFFFBABDFDFFDF.;-C.ADE/FFF9\n+@M00967:43:000000000-A3JHG:1:1104:17020:27386 2:N:0:188\n+CCTGTTCGATACCCACGCTTTCGAACCTCAGCGTCAGTTACAGACCAGAGAGCCGCTTTCGCCACTGGTGTTCTTCCATATATCTACGCATTTCACCGCTACACATGGAGTTCCACTCTCCTCTTCTGCACTCAAGTCTCCCAGTTTCCAATGCACTACTCCGGTTAAGCCGAAGGCTTTCACATCAGACTTAAAAGACCGCCTGCGTTCCCTTTACGCCCAATAAATCCGGATAACGCTTGCCACCTACG\n++\n+3A?AAFFBBFBBGGFGGGGGGGGGGHGHHGHFGGGGGHHHHHHGHHHHHHHGHHGGGGGHGGGGGHHHHHHGHHHHHEHHHHHHHHGGGGGGHHHFHGGGGGFHHGHHFHHHHHHGHFHHHHGHHHHGGHHHHGHGHFHH<FGHHHGHGHGHHHHHHHHHHGGGCFG<FHHGCCDGHHHHHHHHHHHHFHHHHHHFFGHGGDCBGFDFDEGGGGGGFGGGFF..FFFFFFFADFFFBACFFEFFFFFFEF9\n+@M00967:43:000000000-A3JHG:1:1104:14705:27440 2:N:0:188\n+CCTGTTCGATACCCACGCTTTCGTGCCTCAGCGTCAGTAGGACGCCGGTATGCCGCCTTCGCAATCGGAGTTCTGCGTGATATCTATGCATTTCACCGCTACACCACGCATTCCGCATACTTCTCTCCCACTCAAGAACGCCAGTTTCTACGGCTCGACGGCGTTGAGCACCGCTGTTTTACCGCTTGCCTGGCTTCCCCCCTACCCCCCCCTTAAACCCAATAAACTCGGGCAACCGCCTCATCCTCCC\n++\n+AAAAAF3AAFAAGGGAGGGGGG?0FFHHHHHHGGGGGHCG111/AFAA///FG1EFEGGHGGCF/@FGEEGHHG1B/>/?GHHHGH2FGHHH2B>GHEGGGGEHF//?A/ECHH0C///?1@GF111<FHGHHGHH10F.>CCGHHHH000<<C-A--;--9:9A-.;00;.-;--9;FFBFF------/;/9--:-BB--;--9-9--;9-@F/BBBBB?E-/9B//;-9;--;A------//;B//9-\n+@M00967:43:000000000-A3JHG:1:1104:19959:27794 2:N:0:188\n+CCTGTTCGATACCCACACTTTCGAGCCTCAGCGTCAGTTAAAGCCCAGCAGGCCGCCTTCGCCACTGGTGTTCCTCCTAATATCTACGCATTTCACCGCTACACCAGGAATTCCGCCTGCCTCTACTTCACTCAAGAACTGCAGTTTCGATCGCATCTACCAGTTGAGCCCCTAGATTTAACCACCAACTTACTATCCCGCCTACGCCCCCCTTACACCCAGTAATTCCGGACAACGCTTGCTCACTCCG\n++\n+1AA@AF3>1BADBCGC1EBGGFG00EEFGHHHGGFGCHFAD21ABF/FGF0FCAEEGFGHHG/E/1BG1BG?G22B0B1FHFHHFB1FEFGEHHGHHGGGGGHH0/0?/BCGHHHGCC/F0FCFB12@@GGHH0GF1<011?GG1<G1/?10>AF.-1>=11<DGHDCCG...</0<;C0;C00..::.CFF00C0BFGBGG?ABC@--AA-9F;/9/;AE-;BF/B/BB@---;-;9-99A9//://:-\n+@M00967:43:000000000-A3JHG:1:1104:14252:27883 2:N:0:188\n+CCTGTTTGCTCCCCACGCTTTCGCGCCTCACCGTCAGTTGCTGTCCAGCAGACCGCCTTCGCCACTGGTGTTCCTCCTTATATCTACGCATTTCACCGCTACACAAGGAATTCCGTCTGCCTCTCCAGTACTCAAGAAAAACAGTTTCAAATGCAATTTGGAGGTTGAGCCTCCAGATTTCACATCTGACTTGCCTTCCCGGCTACACGCCCTTTACACCCAGTAAATCCGGACAACGCTTGCCACCTACG\n++\n+BBABBFFFFFFFG?E?EGGGGGGGGGGGGFHHGHGGGHGHHFGGFHHFHHHHHHGGGGGHHGGGGHGHGGHGHHHHFHHFHHHHHFHGGFFGHHHHHGGGEEGHHFHHFHHHHFHHHFGHHHHFHHEHHHHGHHGGHHHHHGHHHHHDHHHHFHHFFHHHHHHGDHHHHHH/<GFHFGDGCB;GHBFHGGHHHHHHH0CGDF?DAFGGGGGG?FGGGGBGFGGFFFBFF/ADABC..@CDEEEFBFF.9B.\n+@M00967:43:000000000-A3JHG:1:1104:16540:28173 2:N:0:188\n+CCTGTTTGCTCCCCACGCTTTCGAGCCTCAGCGTCAGTGACCGTCCAGTAAGCCGCCTCCGCCACTGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCGCTTACCCCTCCGGCACTCCAGCCCGGCAGTTTCCAGTGCAGTCCCGGGGTTGGGCCCCGGGCTTTCACACCGGACTTGCCGTGCCGCCTGCGCTCCCTTTACGCCCAGTGAATCCGGATAACGCTCGCCCCCCACG\n++\n+3AABBF5CFFFFGGCFGFGGGGGGHGGHHHHCEFGGGHBHHHGGGGEHHHHFHHGGGGGHGGGGGHHH@GHGHHHEHHHHHHHHHEFEFGGGHHHHHGGGGGHHHFHGHHHHHHEDDFDHGFHCDD/<<?@FHFHBGGDD?CCDHFHFHHGFHHFFHHFHF@@CBFGG?FEBC?CCCFEFF/FBFF-@>B;BFFFCD-.;@AC..9---;AFFBBFFEDFFDF/;///BBAFA..BFD=.;--@D=B--;-\n' |
| b |
| diff -r 000000000000 -r 331fd79a9341 test-data/test_output.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_output.txt Thu Sep 03 01:32:31 2020 +0000 |
| b |
| b'@@ -0,0 +1,4000 @@\n+@M00967:43:000000000-A3JHG:1:1101:18327:1699_F3D0/1\n+NACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+#>>AABABBFFFGGGGGGGGGGGGGGGGHHHHHHHGGGHHHHHGHGGGGGGGHGGGGGGHHHHHHHHHHGGGGGHHHHGHGGGGGGHHBGHGDGGGGGHHHGGGGHHHHHHHHGGGGGHG@DHHGHEGGGGGGBFGGEGGGGGGGG.DFEFFFFFFFDCFFFFFFFFFFFFFFFFFFFFFFFFFFDFDFFFEFFCFF?FDFFFFFFFFAFFFFFFFFFFFBDDFFFFFEFADFFFFFBAFFFA?EFFFBFF\n+@M00967:43:000000000-A3JHG:1:1101:14069:1827_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+3AA?ABBDBFFBEGGEGGGGAFFGGGGGHHHCGGGGGGHFGHGGCFDEFGGGHGGGEGF1GGFGHHHHHGGEGGHHHHHFGGGGGGHHHHHGGGGCDDGHHGGGFHHHHHHHHCD@CCHGGGGHEHGGG@GFGGGGGGG@BGGGEGCEBFFFBFFB;9@EFFFEFFFFFFFFFFFFAFBBBFFFFFBBBFFFFBBBFFFFFFFFFFFBBBBBBBFFFFFFFFFDDFAFFFFF.AF9/FBBBBB.EAFFE?F\n+@M00967:43:000000000-A3JHG:1:1101:18044:1900_F3D0/1\n+TACGGAGGATGCGAGCGTTGTCCGGAATCACTGGGCGTAAAGGGCGCGTAGGCGGTTTAATAAGTCAGTGGTGAAAACTGAGGGCTCAACCCTCAGCCTGCCACTGATACTGTTAGACTTGAGTATGGAAGAGGAGAATGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGATTCTCTGGGCCAAGACTGACGCTGAGGCGCGAAAGCGTGGGGAGCAAACA\n++\n+BA@BBBABBFFFGGGGGGGGGGHGGGGGHHHHGHHHGGGGHHHGGGGGGGGGHGGGGGGFHFHHHFHHHHHHGFHHHHHHFHGGFGHHGGHGGHHFHHHHHHHHHHHHHHGHHHGHDHHHHHHHHHGHDGFFHHHDHHHHHHHHHHHHFHHHDGHHGGGGGHHHHHGGGGEGGGGGBF0BGGGGGGGGFFFGGGDGFFFFBCFBFFBFFBEEFFEFFFFFFFFBDA9EEFFF>DFFFFDDFFFADDFFFF.\n+@M00967:43:000000000-A3JHG:1:1101:13234:1983_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGAAGATCAAGTCAGCGGTAAAATTGAGAGGCTCAACCTCTTCGAGCCGTTGAAACTGGTTTTCTTGAGTGAGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCTCAACTGACGCTCATGCACGAAAGTGTGGGTATCGAACAG\n++\n+BCBCCCCCCFFFGGGGGGGGGGHGGGGGHHHHHHHGGGHHHHHGHGGGGGGGGGGGGGHHHHHHHHHHHGGGGGHHHHHHHHHGGHHHHHHHHGHHHGHGGGGGGHHHHHHHGHHGHHHGHHHHHHHGGGGGHHHHHHGGGGGHHGGGGHGHHHHGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFAFFFFFFFFAFFFFFFFFFFFFFFBFFFFFFFFFFEFFFF\n+@M00967:43:000000000-A3JHG:1:1101:16780:2259_F3D0/1\n+TACGTAGGTGGCAAGCGTTATCCGGATTTACTGGGTGTAAAGGGCGTGTAGGCGGGGACGCAAGTCAGATGTGAAAACCACGGGCTCAACCTGTGGCCTGCATTTGAAACTGTGTTTCTTGAGTACTGGAGAGGCAGACGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGTCTGCTGGACAGCAACTGACGCTGAGGCGCGAAAGCGTGGGGAGCAAACG\n++\n+BBCCBFFCDFCFGGGGGGGGGGHGGGGGHHHHHHHHGGHFHHHGGGEGGHHHGGGGGGGEFGGGHHHHFHHHHGHHHHG/FGGGFGHFHHHGA?CGEHGGGHGHHFBGHHHGGHHHHDGHHFHFHHHHHEC0DDDEDAFFDDEGHHHHHEHHHHHEGGAFFGGGGBDGGGGGGFFFGGFFG.FFFFFFFFFFFFFFFDEFF@=FFFFFF.AFFEFFFFFFFDFFFBFF.:;F-@FFBDAFAAC.;9?FFE.\n+@M00967:43:000000000-A3JHG:1:1101:19378:2540_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGAAGATCAAGTCAGCGGTAAAATTGAGAGGCTCAACCTCTTCGAGCCGTTGAAACTGGTTTTCTTGAGTGAGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCTCAACTGACGCTCATGCACGAAAGTGTGGGTATCGAACAG\n++\n+CCCCCCCCCFFFGGGGGGGGGGHGGGGGHHHHHHGGGGHHHHHGFFGGGGGGGGGGGGHHHGHHHHHHHGGGGGHHHHHHHHHGGHHGGHHHHHHHHHHGGGGGGGHHHHHHHHHGHHHHHHHHHHHGGGGGGHHHHHGGGGGHHGGGGHGHHHHGCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFDFFFFFFADFFFFFFFFAFFFFFFFFF\n+@M00967:43:000000000-A3JHG:1:1101:17674:2779_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGCAGGCGGACTCTCAAGTCAGCGGTCAAATCGCGGGGCTCAACCCCGTTCCGCCGTTGAAACTGGGAGCCTTGAGTGCGCGAGAAGTAGGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCCTACCGGCGCGCAACTGACGCTCATGCACGAAAGCGTGGGTATCGAACAG\n++\n+ABB@BABDAFFFGGGGGGGGGGHGGGGGHHGHHGHGGGHFDHHGHEGCGGGGGGGFEECGHHHHHHHHHGGGGGHHHGGGGGGGGGHGHGHGGGGHHHGGGGDGGEHHHGFHHGGGFHHGHFHFHCGGGGGGGFGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF?FFFFFFFFFFFFFAFEFFFFFFFFFFBFFDFFFFFFFFFDBF/AEFFFF\n+@M00967:43:000000000-A3JHG:1:1101:18089:2781_'..b'BBFFADBFFFFF?BBBBFFFFFFF/;B?BBFFFFFFBDFFFBD9ABBBEDEDFFFBF\n+@M00967:43:000000000-A3JHG:1:1104:14390:26638_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCCTGCCAAGTCAGCGGTAAAATTGCGGGGCTCAACCCCGTACAGCCGTTGAAACTGCCGGGCTCGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACCCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGAGGCACGAAAGTGCGGGGATCAAACAG\n++\n+BABBBBBDBFFFGGGGGGGGGGHFGGGGHHHHHHHGGGHHHHHGFGGGGGGGGGGGGGGHHHHHHHHHHGGEGGHHHHHHGGGGGGHHHGHGGGGGGGHHHGGGGHHHHHHHHGGGGGHGGGHHGGGGGACEFGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFF9.:BFFFFFDFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF.:AFBFFADAFA?EFFFFFF\n+@M00967:43:000000000-A3JHG:1:1104:16844:26739_F3D0/1\n+TACGTAGGGAGCGAGCGTTATCCGGATTTATTGGGTGTAAAGGGTGCGTAGGCGGTAATACAGGTCTTTGGTATAAGCCCGAAGCTTAACTTCGGTAAGCCAGAGAAACCGTATAACTAGAGTATGGTAGAGGCAAGTGGAATTTCTAGTGTAGCGGTAAAATGCGTAGATATTAGAAGGAACACCAGTGGCGAAGGCGACTTGCTGGGCCAAAACTGACGCTGAGGCACGAAAGCGTGGGGAGCAAATAG\n++\n+BBBBCFFCCFCCGGGGGGGGGGHGGGGGHHHHHHHGGGHHHHHGHGGGGGGGHGGGGGHHHHHHHHHHHHHHHHHHHHHHGGGGGHHHHHHHHGGGHGHHHHHHHHHHGHGGHHHHGHHHHHHHHHHHHHHHHGGHGHHHHHHHHHHHHFHHHHHGGGGGHHHHHGCADGHHHHHHHHHGHHHGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFFBDADFFF?DFFFFFFF\n+@M00967:43:000000000-A3JHG:1:1104:17020:27386_F3D0/1\n+TACGTAGGTGGCAAGCGTTATCCGGATTTATTGGGCGTAAAGGGAACGCAGGCGGTCTTTTAAGTCTGATGTGAAAGCCTTCGGCTTAACCGGAGTAGTGCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAACTCCATGTGTAGCGGTGAAATGCGTAGATATATGGAAGAACACCAGTGGCGAAAGCGGCTCTCTGGTCTGTAACTGACGCTGAGGTTCGAAAGCGTGGGTATCGAACA\n++\n+ABBBAFFBFFBFGGGEEGGGGGHGGGGGHFFHHHHCEEFGGHHGHEGGGGGGCGCEEGHHHFGHHHHHHHGGFEHHHHHHHHHGGGGHHHGEGECFFHEGFHGHHH?GHHHHHHHGAFHHHHHGHHHHHHHHHHFGHGFFFHHHHHHHHHHHHHHGGGGGGHHGHGGGGFEHGHFCGHFFHHHHHHGGFFGGGGGGFGBEGGFFFFFFF/BFFBFFFFFBFFDFDAFFFFBCEFFFFDDDCFAFFBDFFF9\n+@M00967:43:000000000-A3JHG:1:1104:14705:27440_F3D0/1\n+TACGGAGGATGCGAGCGTTATCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGGATGCCAAGTCAGCGGTAAAAAAGCGGTGCTCAACGCCGGCGAGCCGTTGAAACTGGCGTTCTTGAGTGGGCGAGAAGTATGCGGAATGCGTGGTGTAGCGGTGAAATGCATAGATATCACGCAGAACTCCGATTGCGAAGGCAGCATACCGGCGCCCTACTGACGCTGGGGCACGAAAGCGTGGGTATCGAACAG\n++\n+AAAAAAAD11BFCGGGGGG?GGFGCGGGHHHHHGH//BGHHHBEBC?EAEEGHGG//?>CCGGFHHFHEE@E?FHHBE/C/<A@CHHHHHGCCC/<<--<ACGGGGHHFCBF.-:;CEG0BB0BF?EGGGGGG0BBBFB@@@AFBBBF-:A-9FF=@A@?/9BBBB/B/B//F/99A-99=F-;BFA@-A?/;9@@BEFBFFBFBF-99-;9-9B-//99-9B---9999-;:B-/;AB@@-/;/BFFE@/\n+@M00967:43:000000000-A3JHG:1:1104:19959:27794_F3D0/1\n+TACGTAGGGAGCAAGCGTTGTCCGGAATTACTGGGTGTAAAGGGAGCGTAGGCGGGATAGTAAGTTGATTGTTAAAACTATCGGCTCAACTGGGAGATGCCATCAAAACTGCAATTCTTGAGTGAAGTAGAGGTAGGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAATGGCGAAGGCGGCCTGCTGGGGCTCTACTGACGCTGAGGCACGAAAGTGTGGGGAGCAAACAG\n++\n+AAAAAFF?11AAAGGGGGGCGGHGGG?FHH1FHHHHAEHHHHHGEGGEEGCGGGGG/E01BGHHHB122@BG2BGH11F12>EEGGGHHHB0</0//BFG11B11FFHGGF11121@G2>GHHF11<?<GF1G1?1FGCGGGF0D0D0=GFGHGHGG.:CGCGFB.ECGEGGGBF0F0/A?BF9FF..9/9B@-A@B?---9A:F//B---;-///9//9BAA--AE9-A-9ABA/:/:9B-----;FFFF\n+@M00967:43:000000000-A3JHG:1:1104:14252:27883_F3D0/1\n+TACGTAGGTGGCAAGCGTTGTCCGGATTTACTGGGTGTAAAGGGCGTGTAGCCGGGAAGGCAAGTCAGATGTGAAATCTGGAGGCTCAACCTCCAAATTGCATTTGAAACTGTTTTTCTTGAGTACTGGAGAGGCAGACGGAATTCCTTGTGTAGCGGTGAAATGCGTAGATATAAGGAGGAACACCAGTGGCGAAGGCGGTCTGCTGGACAGCAACTGACGGTGAGGCGCGAAAGCGTGGGGAGCAACCA\n++\n+BBBABFF>FFB@FFFEGGGGGGHGGGGGHHHHHHGHG2FHGHHDGGGGGHHHHGGGGGGGHHHGFHHHHHHHHGHHHHHHHHHGGHHHHHHHHHHHHHHHHGHHHHHHHHHHHHHHGHHHHFHHHHHHHHGHGGGEFGHGGGGGHHHHHEHHHHHGGGGGEDHHHFGGGGGHHHDFHGFEGGGHHHGHHHGHGGGGFGE?-ACGFFFFFFFFFFFFFFFFFFF.AFFFFFCCCFA?FAADFFF?DFFB...\n+@M00967:43:000000000-A3JHG:1:1104:16540:28173_F3D0/1\n+TACGTAGGGGGCGAGCGTTATCCGGATTCACTGGGCGTAAAGGGAGCGCAGGCGGCACGGCAAGTCCGGTGTGAAAGCCCGGGGCCCAACCCCGGGACTGCACTGGAAACTGCCGGGCTGGAGTGCCGGAGGGGTAAGCGGAATTCCTGGTGTAGCGGTGAAATGCGCAGATATCAGGAGGAACACCAGTGGCGGAGGCGGCTTACTGGACGGTCACTGACGCTGAGGCTCGAAAGCGTGGGGAGCAAACA\n++\n+AAAABFFBBBBBGGGGGGGGGGHGGGGGHHHHHHHGFEEGGHHGGGGGGGGGGGGGGGGGGGGGHHHGGCGGGHHHHHHHGGGGGGGGHGGGGGGGGGHHHHGHHGFGGGGGGGGGGFGGGGGGGGGGFCBFFFFFFFFFFFFFFFFFFFFFFFFFFF9DAFFFFFFFFFFFFFEFFFFFFFFFFFFFFBFFFFF>BFDFD;@@BFFFF9.=@.;A9/FEFFB?D.FE?EFFBFDD/..9@DF.;.;BB9.\n' |