Merge pull request #238 from KevinMenden/splitcigar

Splitcigar

.github/filters.yml

@@ -140,6 +140,10 @@ gatk4_samtofastq:
   - software/gatk4/samtofastq/**
   - tests/software/gatk4/samtofastq/**
 
+gatk4_splitncigarreads:
+  - software/gatk4/splitncigarreads/**
+  - tests/software/gatk4/splitncigarreads/**
+
 gffread:
   - software/gffread/**
   - tests/software/gffread/**

software/gatk4/splitncigarreads/functions.nf

@@ -0,0 +1,59 @@
+/*
+ * -----------------------------------------------------
+ *  Utility functions used in nf-core DSL2 module files
+ * -----------------------------------------------------
+ */
+
+/*
+ * Extract name of software tool from process name using $task.process
+ */
+def getSoftwareName(task_process) {
+    return task_process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()
+}
+
+/*
+ * Function to initialise default values and to generate a Groovy Map of available options for nf-core modules
+ */
+def initOptions(Map args) {
+    def Map options = [:]
+    options.args          = args.args ?: ''
+    options.args2         = args.args2 ?: ''
+    options.publish_by_id = args.publish_by_id ?: false
+    options.publish_dir   = args.publish_dir ?: ''
+    options.publish_files = args.publish_files
+    options.suffix        = args.suffix ?: ''
+    return options
+}
+
+/*
+ * Tidy up and join elements of a list to return a path string
+ */
+def getPathFromList(path_list) {
+    def paths = path_list.findAll { item -> !item?.trim().isEmpty() }  // Remove empty entries
+    paths = paths.collect { it.trim().replaceAll("^[/]+|[/]+\$", "") } // Trim whitespace and trailing slashes
+    return paths.join('/')
+}
+
+/*
+ * Function to save/publish module results
+ */
+def saveFiles(Map args) {
+    if (!args.filename.endsWith('.version.txt')) {
+        def ioptions = initOptions(args.options)
+        def path_list = [ ioptions.publish_dir ?: args.publish_dir ]
+        if (ioptions.publish_by_id) {
+            path_list.add(args.publish_id)
+        }
+        if (ioptions.publish_files instanceof Map) {
+            for (ext in ioptions.publish_files) {
+                if (args.filename.endsWith(ext.key)) {
+                    def ext_list = path_list.collect()
+                    ext_list.add(ext.value)
+                    return "${getPathFromList(ext_list)}/$args.filename"
+                }
+            }
+        } else if (ioptions.publish_files == null) {
+            return "${getPathFromList(path_list)}/$args.filename"
+        }
+    }
+}

software/gatk4/splitncigarreads/main.nf

@@ -0,0 +1,41 @@
+// Import generic module functions
+include { initOptions; saveFiles; getSoftwareName } from './functions'
+
+params.options = [:]
+def options    = initOptions(params.options)
+
+process GATK4_SPLITNCIGARREADS {
+    tag "$meta.id"
+    label 'process_medium'
+    publishDir "${params.outdir}",
+        mode: params.publish_dir_mode,
+        saveAs: { filename -> saveFiles(filename:filename, options:params.options, publish_dir:getSoftwareName(task.process), publish_id:meta.id) }
+
+    conda (params.enable_conda ? 'bioconda::gatk4:4.1.9.0' : null)
+    if (workflow.containerEngine == 'singularity' && !params.singularity_pull_docker_container) {
+        container 'https://depot.galaxyproject.org/singularity/gatk4:4.1.9.0--py39_0'
+    } else {
+        container 'quay.io/biocontainers/gatk4:4.1.9.0--py39_0'
+    }
+
+    input:
+    tuple val(meta), path(bam)
+    tuple path(fasta), path(fai), path(dict)
+
+    output:
+    tuple val(meta), path('*.split_cigar.bam'), emit: bam
+    path  '*.version.txt'                     , emit: version
+
+    script:
+    def software = getSoftwareName(task.process)
+    def prefix   = options.suffix ? "${meta.id}.${options.suffix}" : "${meta.id}"
+    """
+    gatk SplitNCigarReads \\
+        -R $fasta \\
+        -I $bam \\
+        -O ${prefix}.split_cigar.bam \\
+        $options.args
+
+    gatk --version | grep Picard | sed "s/Picard Version: //g" > ${software}.version.txt
+    """
+}

software/gatk4/splitncigarreads/meta.yml

@@ -0,0 +1,61 @@
+name: gatk4_splitncigarreads
+description: Splits reads that contain Ns in their cigar string
+keywords:
+  - vcf
+  - merge
+tools:
+  - gatk4:
+      description: |
+        Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools
+        with a primary focus on variant discovery and genotyping. Its powerful processing engine
+        and high-performance computing features make it capable of taking on projects of any size.
+      homepage: https://gatk.broadinstitute.org/hc/en-us
+      documentation: https://gatk.broadinstitute.org/hc/en-us/categories/360002369672s
+      doi: 10.1158/1538-7445.AM2017-3590
+params:
+  - outdir:
+      type: string
+      description: |
+        The pipeline's output directory. By default, the module will
+        output files into `$params.outdir/<SOFTWARE>`
+  - publish_dir_mode:
+      type: string
+      description: |
+        Value for the Nextflow `publishDir` mode parameter.
+        Available: symlink, rellink, link, copy, copyNoFollow, move.
+  - enable_conda:
+      type: boolean
+      description: |
+        Run the module with Conda using the software specified
+        via the `conda` directive
+  - singularity_pull_docker_container:
+      type: boolean
+      description: |
+        Instead of directly downloading Singularity images for use with Singularity,
+        force the workflow to pull and convert Docker containers instead.
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test']
+  - bam:
+      type: list
+      description: BAM/SAM/CRAM file containing reads
+      pattern: "*.{bam,sam,cram}"
+  - fasta:
+      type: tuple of files
+      description: |
+        Tuple of fasta file (first), sequence dict (second) and fasta index (third)
+      pattern: ["*.fasta", "*.dict", "*.fai"]
+output:
+  - bam:
+      type: file
+      description: Output file with split reads (BAM/SAM/CRAM)
+      pattern: "*.{bam,sam,cram}"
+  - version:
+      type: file
+      description: File containing software version
+      pattern: "*.version.txt"
+authors:
+  - "@kevinmenden"

tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.fna.fai

@@ -0,0 +1 @@
+MT192765.1	29829	120	80	81

tests/data/fastq/dna/sarscov2_1.fastq.gz

@@ -0,0 +1,400 @@
+@ERR5069949.2151832 NS500628:121:HK3MMAFX2:2:21208:10793:15304/1
+TCATAAACCAAAGCACTCACAGTGTCAACAATTTCAGCAGGACAACGCCGACAAGTTCCGAGGAACATGTCTGGACCTATAGTTTTCATAAGTCTACACACTGAATTGAAATATTCTGGTTCTAGTGTGCCCTTAGTTAGCAATGTGCGT
++
+AAAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEE<EEAAAEEEEEEEEEAAAAEAEEEAEEEEEE<AAAA
+@ERR5069949.576388 NS500628:121:HK3MMAFX2:4:11501:11167:14939/1
+ACTGTTTTCTTTGTAGAAAACATCCGTAATAGGACCTTTGTATTCTGAGGACTTTGTAAGTAAAGCACCGTCTATGC
++
+AAA6AEEEEEEEEEAEEE/6EEAEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEE<AAEEEEEEEEEEE</EEEA/AE
+@ERR5069949.501486 NS500628:121:HK3MMAFX2:1:21305:3080:7962/1
+GTACATAACAGCTTCTACACCCTTAAGGGTTGTCTGCTGTTGTCCACAAGTTTTACACACCACGTTCAAGACTCTTTTGCAAGAATCTAAATTGGCATGTTGAAACAAGTAACTCATTGTTTCTCTAACATCACCTAACTCACCTA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEE</EEEEEE/EEAAEEEEEE<EEEEEEEEA<E6EEEEAEAAAAE
+@ERR5069949.1331889 NS500628:121:HK3MMAFX2:3:11605:17658:18222/1
+GTCTACAAGCTGGTAATGCAACAGAAGTGCCTGCCAATTCAACTGTATTATCTTTCTGTGCTTTTGCTGTAGATGCTGCTAAAGCTTACAAAGATTATCTAGCTAGTGGGGGACAACCAATCACTAATTGTG
++
+A/AAAEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEEEEEEEAEEEEEEEEEEEEEEE/EEEEE<AEAEEEEE/EAEAEEE/AEEEEEEEEEEEEEEEEEEEEAE/EEEEEEEEEEEEEEEEEEEEEEEA<EE
+@ERR5069949.2161340 NS500628:121:HK3MMAFX2:2:21107:15810:18317/1
+AACCAGAATATTTCAATTCAGTGTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGTCGGCGTTGT
++
+A/AA//EEAEA/E/AEEEE6EE/EEEA/6AEEEEEEEEE6EEEAEAEE//A/EEEEEE//E/E/A//E/E/<<EE</E/E
+@ERR5069949.973930 NS500628:121:HK3MMAFX2:2:21105:5548:5659/1
+CTAAAAGTGAAGTCAAAATTGTGAGTAACAACCAGTGGTGTGTACCCGTGATTGTTCTTTTGACTGCACTTTGGAAAGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEAE/E///EE/EE/AE6/A//E/E//</E//E/EE/E////<
+@ERR5069949.2417063 NS500628:121:HK3MMAFX2:1:11310:4428:10115/1
+CCTGTTGTCTATGTGATAGACGTGCCACATGCTTTTCCACTGCTTCAGACACTTATGCCTGTTGGCATCATTCTATTGGATTTGATTACGTCTATAATCCGTTTATGATTGATGTTCAACAATGGGGTTTTACAGGTAACCTACAAAGCA
++
+AAAAAEEEAEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAEEEE/EEEEAEEEEE<EEEEAEE<AAEAEEAAEEEEAEEEEEEEEEEEEEEEEEEEEEAEAEEE<AEE<EAE/EEA<EEEEAEEEAEEEEEEEA/EEEAE/AEEEA
+@ERR5069949.376959 NS500628:121:HK3MMAFX2:1:21205:12113:10369/1
+GCTCCATATATAGTGGGTGATGTTGTTCAAGAGGGTGTTTTAACTGCTGTGGTTATACCTACTAAAAAGGCTGGTGGCACTACTGAAATGCTAGCGAAAGCTTTGAGAAAAGTGCCAACAGACAATTATATAACCACTTACCCGGGTCAGG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAAAAAEEEEEEEEEAEEAEEEEEEEEEEEEAAAEEAEA
+@ERR5069949.1088785 NS500628:121:HK3MMAFX2:1:21109:26270:17779/1
+CAGTAGTCTTACTCTCAGTTTTGCAACAACTCAGAGTAGAATCATCATCTAAATTGTGGGCTCAATGTGTCCAGTTACACAATGACATTCTCTTAGCTAAAGATACTACTGAATCCTTTGAAAAAAAGGTTTCACTACTTTCGGTTTTG
++
+AAAAAE/EAEE<EEA///<AEEE/EE<AEEE<EA/EEEEEEE/EAAAEEEEEE<E/E6AE<<E/EEA//</E/EEE/EEE/EE/E/<<EEAAAE<EEEEEE/EAEA//<//AA/E</A<<E/EEEE/AEE<E/<EAE</A6///AEEAA
+@ERR5069949.1066259 NS500628:121:HK3MMAFX2:1:11312:18369:8333/1
+CCTTATGACAGCAAGAACTGTGTATGATGATGGTGCTAGGAGAGTGTGGACACTTATGAATGTCTTGACACTCGTTTATAAAGTTTATTATGGTAATGCTTTAGATCAAGCCATTTCCATGTGGGCTCTTATAATCTCTGTTACTTC
++
+AAAAAEAEEAEEEEEEEEEEEEEEEEAEEEEAEEEEEEEEAEEEEEEEEEEEEEEEEE/EAEEEEEE/6EEEEEEEEEEAEEAEEE/EE/AEEAEEEEEAEEEA/EEAAEAE<AEEAEEEAEAEEEAEAEEAE/AEEEEAEEEEAEA
+@ERR5069949.2832676 NS500628:121:HK3MMAFX2:4:11602:20035:16369/1
+GTCAACCAAAATGCACAAGCTTTAAACACGCTTGTTAAACAACTTAGCTCCAATTTTGGTGCAATTTCAAGTGTTTTAAATGATATCCTTTCACGTCTTGACAAAGTTGAGGCTGAAGTGCAAATTGATAGGTTGATCA
++
+AAAA6EEEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<E/EAEEAEEEAEEAEEEEEAEEEEEEEEEEEEEEAEEAEEEEEAAEEEEEEA<AEEEAAAAEEEEE<EEAAAEEAEEAAEEEEA
+@ERR5069949.2953930 NS500628:121:HK3MMAFX2:3:21504:10316:13537/1
+GTGCTCAAAGGAGTCAAATTACATTACACATAAACGAACTTATGGATTTGTTTATGAGAATCTTCACAATTGGAACTGTAACTTTGAAGCAAGGTGAAATCAAGGATGCTACTCCTTCAGATTTTGTTCGCGCTACTGCAACGATACCGAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEE/EEEEEE<AA/EA<A<AAAA<AAEAAEAEE/E<<E<EAAEA
+@ERR5069949.324865 NS500628:121:HK3MMAFX2:1:11102:17526:14721/1
+CACAAACTCTAAAAGAATGTATAGGGTCAGCACCAAAAATACCAGCAGATAATAATGTTGCAAGTAGAACTTCGTGCAGATTAAAATTTTCATAAGCACTCTAAAGAAGTTGAATGTCTTCAAATTTCTTAACATTAGGGCCCACAACAAG
++
+AAAAAEA/A<EAA/AE/EE/EE//////EA/EEE/E/EEEE//E/6//EA//<AA/A/EEEAA/EEEE/EEEA/E/</AEE////AEEEE//<E//EAE/A///<EEE//<E<<EEE<///A//E/E/EEEAA/<A////<A/AEAAA//E
+@ERR5069949.2185111 NS500628:121:HK3MMAFX2:2:11210:1473:2754/1
+GATGAAACATCATGCGTGATAACACCCTTATAAAACATTTTAAAGCATTGAGCTGATTTGTCTTTATGTGCTTAAAGCTTATTATCATAAACCAAAGCACTCACAGTGTCAACAATTTCAGCAGGACAACGCCGACAAGTTCCGAGGAAC
++
+AAAA6EEAEAAAAAEEEAEEEEEEA<<EEE/AEEE</EEEE/EEE<AA//EAA</A/E<A</EEAE/E/<EEE/AEAA/AE//EEEEEEAEAEE/EEAEA/AEEA/A//A<EEAAAEEEEEEEEE/A6<6//AE<<A<E/AEE<EE/AEA
+@ERR5069949.937422 NS500628:121:HK3MMAFX2:2:11104:8644:10925/1
+TTACCAATGTGCTATGAGGCCCAATTTCACTATTAAGGGTTCATTCCTTAATGGTTCATGTGGTAGTGTTGGTTTTAACATAGATTATGGAATTACCAACTGGAGTTCATGCTGGCACAGACTTAGAAGGTAACTTTTATGGACCTTTTGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEA<E<EE/EAE/E//EE/EEEEEEEAEEE<E/EEE/EEA<EEE</66EE<A/AA<EEEA<EAE/A//AEEAEE///<A<EEEEEEA
+@ERR5069949.2431709 NS500628:121:HK3MMAFX2:3:21609:20858:5879/1
+GTTTATGATTGATGTTCAACAATGGGGTTTTACAGGTAACCTACAAAGCAACCATGATCTGTATTGTCAAGTCATGGTAATGCACATGTAGCTAGTTGTGATGCAATCATGACTAGGTGTCTAGCTGTCCACGAGTGCTTTGTTAAGCGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAAAEAEAEE<AEEEEAEEAA<AEA
+@ERR5069949.1246538 NS500628:121:HK3MMAFX2:4:11505:20131:7933/1
+AGTATGGACAATTCACCTAATTTAGCATGGCCTCTTATTGTAACAGCTTTAAGGGCCAATTCTGCTGTCAAATTACAGAATAATGAGCTTAGTCCTGTTGCACTACGACAGATGTCTTGTGCTGCCGGTACTACACAAACTGCTTGCA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEEEA/EEEEEAEEEEEE/EEEEEEEEEEEEAAAAEEAEEEEEEEEEEEEEEEEE
+@ERR5069949.1189252 NS500628:121:HK3MMAFX2:2:21308:10087:1427/1
+CTATAACACATATAAAAATACGTGTGATGGTACAACATTTACTTATGCATCAGCATTGTGGGAAATCCAACAGGTTGTAGATGCAGATAGTAAAATTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEAEEE
+@ERR5069949.2216307 NS500628:121:HK3MMAFX2:1:21302:10479:19538/1
+CCAAGCAGGGTTACGTGTAAGGAATTCTCTTACCACGCCTATTTGTGGCCTGATAATTGCAGATGAAACATCATGCGTGATAACACCCTTATAAAACATTTTAAAGCATTGAGCTGATTTGTCTTTATGTGCTTTAAGCTTATTATC
++
+AAAE/AEAA/EEE/AE6EEEA/EEEEEEEEEE/<EAE<<6E/EE6E<A/AEE/EAE<</<EAEE/EEEEEEAEA/EAEEEEEEAEA//EAE<<EEA<EEE</<AA<</A<<A<AAAEEE/A6A6A//<AAAAEA6<A<<EAAA/AEA
+@ERR5069949.3273002 NS500628:121:HK3MMAFX2:4:11506:23944:16676/1
+GGTGTATTCAAGGCTCCCTCAGTTGCAACCCATATGATGCCGTCTTTGTTAGCACCATAGGGAAGTCCAGCTTCTGGCCCAGTTCCGAGGTAGTAGAAATACCATCTTGGACTGAGATCTTTCATTTTACCGTCACCACCACCCATTC
++
+AAAAAAEEEEEEEEEEEEEEEEEEEE/EEEE/EAEEEEEEEEEAEEEAEAEEEEEEEEEEEE<EEEAE/AEAEE/<EEEE/A/E/E/<AE<EEEAE//<EEAAAEA/E<<<EE<AA/E/6E6A//AEEE<EAEA/EE/<E/A//A/EA
+@ERR5069949.3277445 NS500628:121:HK3MMAFX2:4:11601:26589:14031/1
+TGACCAAATTGGCTACTACCGAAGAGCTACCAGACGAATTCGTGGTGGTGACGGTAAAATGAAAGATCTCAGTCCAAGATGGTATTTCTACTACCTAGGAACTGGGCCAGAAGCTGGACTTCCCTATGGTGCTAACAAAGACGGCATCATA
++
+AA/AAEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEAEEEEEEEEEEEAEEEEEEEEEEEEEE/AEEEAEEEAAEEAEEE<E/EEEEEEEEAEEEEA/EEEEEEEAAA//EEEEEAEA<AEEAAEAEAEAEAEEEEAEEEEEEA/EEAA
+@ERR5069949.3022231 NS500628:121:HK3MMAFX2:1:21102:10897:17190/1
+GTACGAACTTATGTACTCATTCGTTTCGGAAGAGACAGGTACGTTAATAGTTAATAGCGTACTTCTTTTTCTTGCTTTCGTGGTATTCTTGCTAGTTACACTAGCCATCCTTACTGCGCTTCGATTGTGTGCGTACTGCTGCAATAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEAAAAEEEEEEEEAEEE
+@ERR5069949.184542 NS500628:121:HK3MMAFX2:2:21205:8408:1786/1
+TCAATGGGGAATGTCCAAATTTTGTATTTCCCTTAAATTCCATAATCAAGACTATTCAACAAAGGGTTGAAAAGAAAAAGCTTGATGGCTTTATGGGTAGAATTCGATCTGTCTATCCAGTTGCGTCTCCAAATGAATGCAACCAAATGTG
++
+AAAAAEEEEE/EA/E/EEE/EEAEE/E/EE/EEEE//AEE/E/EEE//EEE/</E/EE<<//EE/EE<EEEEEAE/E/EAAEEEAEAEE<E</EEEE/E//E<<<///E<//A<AE</<AEEEAAE///EE</EE//AA///<E</A<AEA
+@ERR5069949.540529 NS500628:121:HK3MMAFX2:3:21609:5345:11242/1
+CCGTCTATGCAATACAAAGTTTCTTTAGAAGTTATATGTTTATAGTGACCACACTGGTAATTACCAGTGTACTCACTAGCACAAGTAAATGTACCATGCTTAAGTTCATACTGAGCAGGTGGTGCTGACATCATAACAAAAGGTGACTC
++
+AAAAAE/A/EEEEAA/EEAEEEAEEEEEEE6EEAEEEEEEAE/E<AE/EEEEEEEEEEAEEEEEEEEEEEEEEAAEAEAEEEAEEE<EEEEEEEAE</E<EEEA<E/EEE<EEEE<EE<EEEA<A6EE</EEE<AEAA<EE<EEAEAEA
+@ERR5069949.686090 NS500628:121:HK3MMAFX2:1:11212:8601:4499/1
+ATGTCAGATTTATGTGACAATTTAAGACATTCAACAACATCTTTAGTTTCTACATCTGAATCAACAAACCCTTGCCGAGCTGCTGAAATAAAAGTAGATAAGACATTGTCTAAGGACACATTCTTTGCAAGTTCAGCTTCTGCAGTTGCA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEE/EEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEE6EEEAEEEEEEEE/EEAE<EEEAEEEEEAEEEEEEEEEEEEEAEAEEE
+@ERR5069949.2787556 NS500628:121:HK3MMAFX2:4:11609:10602:17177/1
+GCTGCTAGAGACCTCGTTTGTGCACAAAAGTTTAACGGCCTTACTGTTTTGCCACCTTTGCTCACAGATGAAATGATTGCTCAATACACTTCTGCACTGTTAGCGG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAE
+@ERR5069949.2650879 NS500628:121:HK3MMAFX2:1:21111:16524:5801/1
+GCTATAACGCAGCCTGTAAAATCATCTGGTAATTTATAATTATAATCAGCAATCTTTCCAGTTTGCCCTGGAGCGATTTGTCTGACTTCATCACCTCTAATTACAAATGAATCTGCATAGACATTAGTAAAGCAGAGATCATTTAATTTA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEAA<AEEEEEEEEEEEAEEEEEE/EEE<EEAEEEAE
+@ERR5069949.2064910 NS500628:121:HK3MMAFX2:3:21608:24504:13289/1
+CTACCCTTCTGCTCGCATAGTGTATACAGCTTGCTCTCATGCCGCTGTTGATGCACTATGTGAGAAGGCATTAAAATATTTGCATATAGATAAATGTAGTAGAATTATACCTGCACGTGCTCGTGTAGAGTGTTTTGATAAATTCAAAG
++
+AAAAAEA<AEEE/EE/EE/EEE/EEE//E/EA/EEEEEEEEEEE<E/A</</A/AEAEEEE/EAE/AEEE<AAE/E//EA</<//E<A/AE<EAEEE///<EA/E<AEAAAA/A/</EEEEEAE/A/<A/AA/EAE<AE/6/<<A</<A
+@ERR5069949.2328704 NS500628:121:HK3MMAFX2:4:21612:9379:15958/1
+GATTTTTAAGTGTGTCACTTAACATTTGTACAATCTTTATACGCACTACATTCCAAGGAAGTCCTTTGTACATAAGTGGTATGAGGTGTTTAAATTGATCTCCAGGCGGTGGTTTAGCACTAACTCTGGAAAAATCTGTATTATTAGGTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/E/EEEEEEEAAEEEEEAAEAEEEEE<AEA/EEEEEEAEAEEAEAAAE/<EEAAE/EEAEEEEEEEAAEEEEEEAEAAEA
+@ERR5069949.1067032 NS500628:121:HK3MMAFX2:4:21403:19263:10183/1
+GTACAAAAATAGCCTAAGAAACAATAAACTAGCATTATACACTGAAGTGTATTTCCAGTTATGAAGAAAATAGGGCAATACTCAACACACATAAAAACAATACCTCTGGCCAAAAACATGACAGTTGTAACTACACCTGAGTAGTTAGAA
++
+AAAAA/EEEAE/EAEEEEEEEEEEEE/E/EEEEEEEE/EEEE<EE//EAEEEE/EAAAEEEEAEEEEEE6AEEEA6EEE//EAE6A/EE<EE/EEEEAA///A<EEEEEEAAEEE/EE/AAA/E/<6<AAEEEA<</EEAEE/EAAE<EE
+@ERR5069949.3338256 NS500628:121:HK3MMAFX2:1:11108:15722:2450/1
+AGAGTAGACTATATATCGTAAACGGAAAAGCGAAAACGTTTATATAGCCCATCTGCCTTGTGTGGTCTGCATGAGTTTAGGCCTGAGTTGAGTCAGCACTGCTCATGGATTGTTGCAATTGTTTGGAGAAATCATCCAAATCTGCAGCAGG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/AEEAEEEEEEEEEEAEEEEE<EAAEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEE<AAEEEEEEAEEA<AEEEEEEEEEEEA
+@ERR5069949.1412839 NS500628:121:HK3MMAFX2:3:11609:15519:12102/1
+CATTAGCACAAGTTGTAGGTATTTGTACATACTTACCTTTTAAGTCACAAAATCCTTTAGGATTTGGATGATCTATGTGGCAACGGCAGTACAGACAACACGATGCACCACCAAAGGATTCTTGATCCATATTGGCTTCCGGTGTAA
++
+AAAAAEEAEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEE<AEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEE/EEEE/EEEAAAEEEE/EEEEE<E/<EAE<AE<<A<EAAAAAAEAA<AEE
+@ERR5069949.1538968 NS500628:121:HK3MMAFX2:2:11304:14669:15804/1
+AAGTTGGCGTATACGCGTAATATATCTGGGTTTTCTACAAAATCATACCAGTCCTTTTTATTGAAATAATCATCATCACAACAATTGTATGTGACAAGTATTTCTTTTAATGTGTCACAATTACCTTCATCAAAATGCCTTAAAGCATAG
++
+AAAAAEEEEEEEAEEAEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEAEAEEEEAEEEEA<EEEEEEAE/EEEE/EEAEEE/AEEE/EEEAEEE//EEEE/EEEEEEE<EEEEA/<6/6AAAEEEAA6AA6<AEA<EA</AE/E<AA6EEA
+@ERR5069949.147998 NS500628:121:HK3MMAFX2:3:21611:12530:12698/1
+GCAATTCATAGCTCTTTTCAGAACGTTCCGTGTACCAAGCAATTTCATGCTCATGTTCACGGCAGCAGTATACACCCCTCTTAGTGTCAATAAA
++
+AAAAAEAEEEEEEEEEEEEEEAAEE/EEE/EE/E/EEE6EEE/EEEE//<EAEAE/A/EEE/EE6<AEEE<EAEA<E/<A/EEE/A//EEE<EE
+@ERR5069949.366975 NS500628:121:HK3MMAFX2:1:11105:12446:15727/1
+GTGTAACCATTTAAACCCTGACCCGGGTAAGTGGTTATATAATTGTCTGTTGGCACTTTTCTCAAAGCTTTCGCTAGCATTTCAGTAGTGCCACCAGCCTTTTTAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEAEA<EEEEEEEEEEAEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEE
+@ERR5069949.1372331 NS500628:121:HK3MMAFX2:3:11601:26764:8562/1
+ATCCTTTAGGATTTGGATGATCTATGTGGCAACGGCAGTACAGACAACACGATGCACCACCAAAGGATTCTTGATCCATATTGGCTTCCGGTGTAACTGTTATTGCCTGACCAGTACCAGTGTGTGTACACAACATCTTAACACAATTAGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEAEEEEEEEEAEEEEEEEEEAEEEEEAAAEEEEEEE<<EE
+@ERR5069949.1709367 NS500628:121:HK3MMAFX2:1:11301:9990:2022/1
+AAGTGTATCTTGATCCTCATAACTCATTGAATCATAATAAAGTCTAGCCTTACCCCATTTATTAAATGGAAAACCAGCTGATTTGTCTAGGTTGTTGACGATGACTTGGTTAGCATTAATACAGCCACC
++
+AAAAAEEEAEEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAAEAAE<AEA<EEE/EEA<AAEE/AEEEEEEEE/EAEAEEA<EAEA/EEEAEEAAAAAAAAEAAEAAE/A6EEAEEEEAAAEE/AA
+@ERR5069949.2388984 NS500628:121:HK3MMAFX2:2:21306:17580:3896/1
+TGTCTATGTGATAGACGTGCCACATGCTTTTCCACTGCTTCAGACACTTATGCCTGTTGGCATCATTCTATTGGATTTGATTACGTCTATAATCCGTTTATGATTGATGTTCAACAATGGGGTTTTACAGGTAACCTACAAAGCAACCAT
++
+AAAAAEEEEEEEEE/EEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEAEEEEAE/EAEEEEEEEEEEAEEEE/EEEEEEEEAEEAEEAEEEEEAEEEEAEEEE<A/EA<E
+@ERR5069949.1132353 NS500628:121:HK3MMAFX2:2:21111:21744:6806/1
+GATTTAGCCACATTCAAAGACTTCTTCAACTTTTTAAGAACAACTTCAGAATCACCATTAGCAACAGCCTGCTCATAAGCTTCTTGAGCAGTAGCAAAAGCTGCATATGATGGAAGGGAACTAAACTCTGAGGCTATAGCTTGTAAGGTT
++
+AAAAAEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEEEEEAEAAEEAE<EEAAEEEEEEEEEEE<A<EEEEEEEEEEEEEEAEE
+@ERR5069949.1151736 NS500628:121:HK3MMAFX2:4:21406:11928:13664/1
+TCTAAGCATAGTGAAAAGCATTGTCTGCATAGCACTAGTAACTTTTGCCCTCTTGTCCTCAGATCTAGCCTGTTTATACATTTGGGTCATAGCTTGATCAGCCATCTTTTCCAACTTACGTTGCATGGCTGCATCACGGTCAAATTCAGAT
++
+AAAAAEEAEEEEEAEE/EEEEAEEEEEEEEEEE/EEEEEEEEEE/EEAEEEEEE/EE<EEEEEEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEEEE<AEEEEEAAEAAEEEEEEA/<AAEEAAEEEEAEAAAE/EE//AEAAAAAA
+@ERR5069949.479807 NS500628:121:HK3MMAFX2:3:11402:16906:16736/1
+TATCTGCCCATTTAATAGAAGTTAAACCATTAACTTGTGGGTATTTCCACTTTTTAGTGTGATTTAATGCTGACATGTACCTACCCAGAAAACTAGGATCAGTTGTGTGGTAGTACTCAAAAGCCTCAACACGTAGAGTGTCATCATTAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEAEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEEEEEEEAEEEEEEEEEEEAAAAEAAA/EEAAEAA/AEEAAAEEEAEAAAEEEE/AA
+@ERR5069949.2176303 NS500628:121:HK3MMAFX2:4:11608:21150:1680/1
+TTACCTGCACCACGCACATTGCTAACTAAGGGCACACTAGAACCAGAATATTTCAATTCAGTGTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGTCGGCGTTGTCCTGCTGAAATTGTTGACACTGTGAGTGCTT
++
+AAAAAEEEEEEAEEEEEEAEEEEEEEEEEEEEEAEAEEEEEAEEEEEEEEEAEEEEEEEEEEE<EEEEAEEEEEEEEEEEEEEEEEEEEEEAEEEE<A<AEEEAE/EE<EEEEAAAAAAEEAA<AAAEEEEEE<EEEEEEEEAEE<EEEEA
+@ERR5069949.2772897 NS500628:121:HK3MMAFX2:1:11106:12629:15009/1
+GTTGAAAGTAGATCTTCAATAAATGACCTCTTGCTTGGTTTTGATGGATCTGGTAATATTTGTGAAAAATTAAAACCACCAAAATCTTTAATTGGTGGTGTTTTGTAAATTTGTTTGACTTGTGCAAAAACTTCTTGGGTGTTTTTGTCTT
++
+AAAAAEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAAEEEEEEEEAEEEEEEEEEEEEEEEEEEEEAAEEAA<EEEEEEA
+@ERR5069949.1020777 NS500628:121:HK3MMAFX2:3:11408:13234:17983/1
+GGATTAGTAACACTACAGCTGATGCATACATAACACAGTCTTTTAGCTTAAAACCAGACAAACTAGTATCAACCATATCCAACCATGTCATAATACGCATCACCCAACTAGCAGGCATATAG
++
+AAAAAEEEE6E6EEAEEE<EEEEEE/EEEEEE6EEEEEEEEEEEEE/EEEEE/EEEEEEEEAEEEEAAEEA/EEEEE<EAEAEAEAE<EAEEEA<EE<AEE/EAAA6AA/E6AAAA6AEEEE
+@ERR5069949.465452 NS500628:121:HK3MMAFX2:3:11503:20241:8941/1
+ACCTGATGCTGTTACAGCGTATAATGGTTATCTTACTTCTTCTTCTAAAACACCTGAAGAACATTTTATTGAAACCATCTCACTTGCTGGTTCTTATAAAGATTGGTCCTATTCTGGACAATCTACACAACTAGGTATAGAATTTCTTAAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEE/EEAEEEE/EEEEEEEEEEAEEEEEEEEEEEEEEE<EEEEAAAEAEEEEEEAA6AAEEEEEA<EEEEE</EEAEE/EE
+@ERR5069949.1704586 NS500628:121:HK3MMAFX2:3:11408:18502:19140/1
+GATAAACGCACTACGTGCTTTTCAGTAGCTGCACTTACTAACAATGTTGCTTTTCAAACTGTCAAACCCGGTAATTTTAACAAAGACTTCTATGACTTTGCTGTGTCTAAGGGTTTCTTTAAGGAAGGAAGTTCTGTTGAATTAAAACA
++
+AAAA6EEEE/EE6EEEEEEEEEEEEEEEEEE<EEEEEEEE6EEAEEEEEA<EEEEE66EEEEE///EEEAEEEE<EEEEEEA/EE/EEEEEEEAE<E<AA<AAAEEAE/AEE<E<AA<EAAEEAE/AEE/E/EAEAAAEE/EA/A//EE
+@ERR5069949.1258508 NS500628:121:HK3MMAFX2:2:21106:23201:10887/1
+CCCTTTGTTGTGTTGTAGTAAGCTAACGCATTGTCATCAGTGCAAGCAGTTTGTGTAGTACCGGCAGCACAAGACATCTGTCGTAGTGCAACAGGACTAAGCTCATTATTCTGTAATTTGACAGCAGAATTGGCCCTTAAAGCTGTTACAA
++
+AAAAAEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAAAAAAEEEAAE/AE</EEEEAA/AAEEEEEEAEEEAAAEEE/AEEAAE
+@ERR5069949.986441 NS500628:121:HK3MMAFX2:2:11310:5407:13956/1
+GTTACTCACAATTTTGACTTCACTTTTAGTCTTAGTCCAGAGTACTCAATGGTCTTTGTTCTTTTTTTTGTATGAAAATGCCTTTTTACCTTTTGCTATGGGTATTATTGCTATGTCTG
++
+AAAAAEAEEEEEEE/EEE/EEEEEAEEEEEEEEEEEEEEEEEEEEE</EAAEA/EEEEEEEEAEAAEEEEEEEEEEEEE/E//<EAE/6///EE//E/EEE///E<EEEEA</A<<//<
+@ERR5069949.2674295 NS500628:121:HK3MMAFX2:1:21303:10381:8722/1
+GTTATAACACTGACACCACCAAAAGAACATGGTGTAATGTCAAGAATCTCAAGTGTCTGTGGATCACGGACAGCATCAGTAGTGTCAGCAATGTCTCTGCCAAATTGTTGGAAAGGCAGAAACTTTTTGTTAGACTCAGTAAGAACAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEE<EEEEEE<EEEEEEE<AEAAA<EAAE/EEEAA<AAEEE<AAEEEE<EEEEEEEEEEEEAEEE
+@ERR5069949.885966 NS500628:121:HK3MMAFX2:4:11610:19682:20132/1
+GTCTAATCATAATTTCTTGGTACAGGCTGGTATTGTTCATCTCAGGGTTATTGGACATTCTATGCAAAATTGTGTACTT
++
+AAA//E/EAA/E//E//E//E/E//AE/A/E//EAEA///AE//E///E/EEE6EEEAEEA///E/AEE/EAEE/E//E
+@ERR5069949.2342766 NS500628:121:HK3MMAFX2:4:21501:7762:9816/1
+CAGGTTATGTTGATACACCTAATAATACAGATTTTTCCAGAGTTAGTGCTAAACCACCGCCTGGAGATCAATTTAAACACCTCATACCACTTATGTACAAAGGACTTCCTTGGAATGTAGTGCGTATAAAGATTGTACAAATGTTAAGTGA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEAEAEAAAAEAEEEAEEAEEE/EEEE<AAEEAEA
+@ERR5069949.3122970 NS500628:121:HK3MMAFX2:2:21304:20157:2460/1
+TGCCAATCCTGTAGCGACTGTATGCAGCAAAACCTGAGTCACCTGCTACACGCTGCGAAGCTCCCAATTTGTAATAAGAAAGCGTTCGTGATGTAGCAACAGTGATTTCTTTAGGCAGGTCCTTGAT
++
+AAA/AEEEEEEEAEEEAEEEEEEEEEEEEE6AEEAEAEAEEAE<AEEEEEEEEAEEAEEAEEEEEAEE</EEEEAEEAAEE/EEEAEEEEEEEEE<A//EEE/EAAEEEEEEAA/AEEAEAAA6//A
+@ERR5069949.3279513 NS500628:121:HK3MMAFX2:1:21110:4328:5711/1
+TGCGTCTTGGGTCCCCGCTCTCACTCAACATGGCAAGGAAGACCTTAAATTCCCTCGAGGACAAGGCGGTCC
++
+AAAAAEEEEA/E//AE/E/EEEE6EEEE/E6/6A/E6E/6AEA/E<EEEE//E/EEEA//AEE///A//E/A
+@ERR5069949.309410 NS500628:121:HK3MMAFX2:1:11107:4644:2079/1
+GAAGAAGATTGGTTAGATGATGATAGTCAACAAACTGTTGGTCAACAAGACGGCAGTGAGGACAATCAGACAACTACTATTCAAACAATTGTTGAGGTTCAACCTCAATTAGAGATGGAACTTACACCAGTTGTTCAGACTATTGAAGTGA
++
+AAAAA//EEEEA6EEEAE</EEE/EEEEE/EE6EEEEEEEEEEEEEEEAEEAAEEEEEEEEEAEEEEEE/EEAEEEEEAEAEEE/EEAEEE<AEEEAA////EEEEEEEEA//A/EE/EAAEA/AE<EE/E//E/</AEAEAE/AEA/AEA
+@ERR5069949.532979 NS500628:121:HK3MMAFX2:2:21104:15575:19089/1
+CATGGGCACACTTTCTTATGAACAATTTAAGAAAGGTGTTCAGATACCTTGTACGTGTGGTAAACAAGCTACAAAATATCTAGTACAACAGGAGTCACCTTTTGTTATGATGTCAGCACCACCTGCTCAGTATGAACTTAAGCATGGTA
++
+AAAAAEEEEEEEEEEEEEEEEEEE/EEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/A/</EEE<EEE<<AEEEEEA/AA/AE/EE/EEEEAEA</EAEEE<<AEEEEE
+@ERR5069949.2888794 NS500628:121:HK3MMAFX2:3:11604:17798:16811/1
+CCTTGAATGAGTCTAATTCAGGTTGCAAAGGATCATAAACTGTGTTGTTGACAATTCCTATTACAACATCACAGTTACCAGACACAAATGTGTTGTCTGTAGTAATGATTTGTGGTTCATAAAAATTCCTTTGTGTTACAAACCAGTGTGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEE/EEEEEEEEEAEEEEAAAAEEAAEEEEAAAEEEEEEEEEEEEEEEAAAEAA
+@ERR5069949.2205229 NS500628:121:HK3MMAFX2:1:11212:26331:4199/1
+CACTAGAACCAGAATATTTCAATTCAGTGTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGTCGGCGTTGTCCTGCTGAAATTGTTGACACTGTGAGTGCTTTGGTTTATGATAATAAGCTTAAGCACATAAAGA
++
+AAAAAEEEEEEA/EEEEEEEAEEEEEEEEEAEAEAEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEAEEEEEEAEEEEEE6EE<E<<EEEAEEEEEEEEEEEEEEEEA<AEEEEEEAAAEEEEEEEEEEEEEEEE<EAAEEEEE
+@ERR5069949.786562 NS500628:121:HK3MMAFX2:2:11109:4461:9003/1
+AACACCAGCTTCTGATCTTTCACAAGTGCCGTGCCTACAGTACTCAGAATCAAAAGTTGTTACCACTCTAACAGAACCTTCAAGGTAGGTGTTAGGAAATTGAATAATAGAGCCATCCATGAGCACATAACGTGTGTCAGGGCGTAAACTT
++
+A6AAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEAAEEEEEEEEAEAEAEEEA/AAAEEEEAEEEEAEE<EEAAEEAAA/EA<EAEA
+@ERR5069949.919671 NS500628:121:HK3MMAFX2:4:21608:23884:15511/1
+CTGTAATAGTTGTGTCCGTACCAGCTGCTTGTGCTGTTTGCCTGTCAACAAAAGGTCCATAAAAGTTACCTTCTAAGTCTGTGCCAGCATGAACTCCAGTTGGTAATTCCATATGGTGCATGTAACAAAAAGAGACACAGTCATAATCTAT
++
+AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEAAEEEEEEEEEEEEEEEEEEEEAEEEAEEEAEEEA/EAAAAEEAEEEEEEEEA6AAEEAAAEAAEEEEEEEE
+@ERR5069949.1328186 NS500628:121:HK3MMAFX2:3:21510:25432:14327/1
+TCTTTGTAAGCTTTAGCAGCATCTACAGCAAAAGCACAGAAAGATAATACAGTTGAATTGGCAGGCACTTCTGTTGCATTACCAGCTTGTAGACGTACTGTGGCAGCTAAACTACCAAGTACCATACCTCTATTTAGGTTGTTTAATCCTT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAE/EEEEEEEEAEEE/EEAEEEEEEAEEEEEEEEAAA<EEEEEAEE/<6A<AEEEEE//EEAA<EEAEE<E<EAA<E/6E</EE
+@ERR5069949.870926 NS500628:121:HK3MMAFX2:4:11510:4100:9432/1
+CAACTACACTTAACGGTCTTTGGCTTGATGACGTAGTTTACTGTCCAAGACATGTGATCTGCACCTCTGAAGACATGCTTAACCCTAATTATGAAGATTTACTCATTCGTAAGTCTAATCATAATTTCTTGGTACAGGCTGGTAATGTT
++
+AAAAAEEAEEEEEEEEEAEEEEEEAEEEEEEEEE/EEEEEEEAEEEEEEEAEEA<AE/AEEEEEEEEE<EEEEEEE/AEAEEEEAAEEEEAAEEEEEAEEEEEEEEAEEAAEEEEEAEAA/EEEAAEEEAAAA/A/EA/E/AEEEE/EE
+@ERR5069949.2257580 NS500628:121:HK3MMAFX2:3:11512:13762:8305/1
+AGTTGCAATTTACAAGTCTTGAAATTCCACGTAGGAATGTGGCAACTTTACAAGCTGAAAATGTAACAGGACTCTTTAAAGATTGTAGTAAGGTAATCACTGGGTTACATCCTACACAGGCACCTACACACCTCAGTGTTGACACTAAATT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEEEEEEEAEEEEEEAEEEEEEEEEEEEAEEEEEEEEEAEEEEEEAEE/EEAEA<EAEEEAEEEEEEEEE<EEAAAEAEEE<EA
+@ERR5069949.3249622 NS500628:121:HK3MMAFX2:2:21101:6979:17698/1
+TTGCCATGTTGAGTGAGAGCGGTGAACCAAGACGCAGTATTAATGGGTAAACCTTGGGGCCGACCTTGTTTTTATCG
++
+AAAAAE/EA/EEAEEEAEEAE/EE/EEE/EEEEEEEE/E//A//EEEEE/EEEEEE/AEE/EA////<<E<///E/E
+@ERR5069949.611123 NS500628:121:HK3MMAFX2:4:21402:10868:14846/1
+ACTCTAGATAATTCATTAGGTTTCTTAATAGTAAGACTAGAATTGTCTACATAAGCAGCCATTAGATCTGTGTGGCCAACCTCTTCTGTAATTTTTAAACTATTATTTGCTGGTTTAAGTATAAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEAEEEEEEEEE/AEEEEAEEE<AEEEEEEEEEEEAEEE
+@ERR5069949.651338 NS500628:121:HK3MMAFX2:2:11307:9572:4021/1
+CTATAACATTAATAGGCAATGAACCTTTAGTGTTATTAGCTCTCAGGTTGTCTAAGTTAACAAAATGAGAGAGAGAATGTCTTTCATAAGTCTTTTGACCAGCTTTAACAAAGTAAAGATGGATGGAACCATTCTTCAGAGT
++
+AAAAAEEEAEAEEEEEEEEEAEEEEA/EEEEA/EE/<AEEAEEE/EA<//EE/A/E/EEEEEEEEEEEEEE/EAEEEEAAEE/E<AEEE/<<EEEEE<E/EEE<A/A/<<EAE/A</EEEA<EE/EE//AEEE</6/A///A
+@ERR5069949.169513 NS500628:121:HK3MMAFX2:2:21303:6012:2171/1
+AATCAAGACTATTCAACCAAGGGTTGAAAAGAAAAAGCTTGATGGCTTTATGGGTAGAATTCGATCTGTCTATCCAGTTGCGTCACCAAATG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EE6EEEE
+@ERR5069949.155944 NS500628:121:HK3MMAFX2:2:21111:24075:19469/1
+GGATAGACAGATCGAATTCTACCCATAAAGCCATCAAGCTTTTTCTTTTCAACCCTTGGTTGAATAGTCTTGATTATGGAATTTAAGTGAAATACAAAATTTGGACATTCCCCATTGAAGGTGTCAAATTTCTTTGCCAATTTAATTTCA
++
+AAAAAEEAEEEEEEE6EEEEEEEEEAEEEA<EEEEEEEEEA66AEA/EAEEEAEE/EEEEEAE/A/EEEEA6EEE/EA/EEEEE/A//EEE/AAEEEAAEEEEAA/<EEAAAAAA6/AEE/EEEAEAAA/66/6/AEA/A6/A</AE<AE
+@ERR5069949.2033605 NS500628:121:HK3MMAFX2:4:21508:23814:9228/1
+CTACACGAGCACGTGCAGGTATAATTCTACTACATTTATCTATAGGCAAATATTTTAATGCCTTCTCACATAGTGCATCAACAGCGGCATGAGAGCAAGCTGTATACACTATGCGAGCAGAAGGGTAGTAGAGAGCTAGGCCAATAGCAA
++
+AAAAAEE/EEEEEEEEEEEEEEEEEAEEEEEEEEEAEEAEEEEEAEEEEEEEEAAEEEEE<EEE<EEEEEAEEAEEEAAE/EEEA/EEA<EEEEEEEEEEEE<EE/AEEAEEAA<AEA/EE/EEAE6E/AE/<E/AAAAEEE<AAE<AAA
+@ERR5069949.2730382 NS500628:121:HK3MMAFX2:4:21605:26351:3903/1
+ACTGAATTTTCTGCACCAAGTGACATAGTGTAGGCAATGATGGATTGACTAGCTACACTACGTGCCCGCCGAGGAGAATTAGTCTGAGTCTGATAACTAGCGCATATACCTGCACCAATGGGTATGTCACACTCATATGAGT
++
+AAAAAE/EEEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEE<EEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEE<EEAAAEEE<AEA<<AAE<AEAEAEA<A/AAAEEE<A<AA<A
+@ERR5069949.2125592 NS500628:121:HK3MMAFX2:1:11307:14726:9430/1
+ATGTGAGAAGGCATTAAAATATTTGCCTATAGATAAATGTAGTAGAATTATACCTGCACGTGCTCGTGTAGAGTGTTTTGATAAATTCAAAGTGAATTCAACATTAGAACAGTATGTCTTTTGTACTGTAAATGCATTGCCTGAGACGAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEAEEEEEEEAEEEEA<EEEEEEEAEEEEEEAEEEEEEAAEEEEEE<EE/A/EEEA
+@ERR5069949.1062611 NS500628:121:HK3MMAFX2:1:11304:13814:5808/1
+AACTAGCATTATACACTGAAGTGTATTACCAGTTATGAAGAAAATAGGGCAATACTCAACACACATAAAAACAATACCTCTGGCCAAAAACATGACAGTTGTAACTACACCTGAGTAGTTAGAAGTAACAGAGATTATAAGAGCCCACATG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEE/AEAEEEEEEEE6AEEAEEEEEEEEEEAEEEEAEEEEEEEEEEEEEEEEEE<AEEEEEEEEEEEEEAEEEEEEAEEEAEAEEEEAAEEEAAEEA
+@ERR5069949.1778133 NS500628:121:HK3MMAFX2:1:11203:21125:7007/1
+CATCTTATTTCTATAGAGACACTCATAAAGTCTGTGTTGTAAATTGCGGACATACTTATCGGCAATTTTGTTACCATCAGTAGATAAAAGTGCATTAACATTGGCCGTGACAGCTTGACAAATGTTAAAAACACTATTAGCATAAGCAGTT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEAAAEEEEA<EEEEEEEEEEEEAEAEAEEEEAEEEEEEEEEEEEEEAEAEEEEEEEAEEEEEEEEEEEEEEEEEEA<AAAAAAEEEEEAEEA
+@ERR5069949.3057020 NS500628:121:HK3MMAFX2:2:21204:6342:14557/1
+CAATTTGCCTATGCCAACAGGAATAGGTTTTTGTATATAATTAAGTTAATTTTCCTCTGGCTGTTATGGCCAGTAACTTTAGCTTGGTTGTACGC
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE
+@ERR5069949.2972968 NS500628:121:HK3MMAFX2:4:21405:4877:14897/1
+AGTAAACTGTTACAAACAACAACAGCAAGTTGCAAACAAAGTGAACACCCTTGGAGAGTGCTAGTTGCCATCTCTTTTTGAGGGTTATGATTTTGGAATCGCTCTGAAAAACAGCAAGAAGTGCAACGCCAACAATAAGAC
++
+AAAAAEEEE//EEEE6AE/AE/E6EEEEEEEE//E/EEEEE//EEA/A<EAEEEEEEEAEEEE/<E</EE/E<E<EE/E/E/E<///A/AE//EEEE//EA/EEAEE//EAA66EEEEEE/E/6//A/EE/EE/A//A//6
+@ERR5069949.2734474 NS500628:121:HK3MMAFX2:4:11507:10751:13427/1
+TGGCAATAGAGTTATTAGAGTAAGCAACTGAATTTTCTGCACCAAGTGACATAGTGTAGGCAATGATGGATTGACTAGCTACACTACGTGCCCGCCGAGGAGAATTAGTCTGAGTCTGATAACTAGCGCATATACCTGCACCAATGGGT
++
+AA/AAAEEEEEEEAEEEE//EEEEEEEEEEEEEEEEEAEEEEE/EEEE/EEAEEEEEEEE/EEEEEEEEE/E/EEE/AEA/EEEAAEE/EEAAEAAEAEEEEAEEEEEEEEA<EEAEE/A<EEEEEEEAAEEEEAA<AAAAE/AEE/AA
+@ERR5069949.856527 NS500628:121:HK3MMAFX2:2:11301:22893:2784/1
+CAACTACACTTAACGGTCTTTGGCTTGATGACGTAGTTTACTGTCCAAGACATGTGATCTGCACCTCTGAAGACATGCTTAACCCTAATTATGAAGATTACTCATTCGTAAGTCTAATCATAATTTCTTGGTACAGGCTGGTAATGTTCAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAAEEEEEEEAEEEE<AAAEEAEEEAEEEE/EEEEEAEEEAEAEE//<//AE//E<E//<//EA/A<EA6EE//E/AA/6A//6/AEE<EE6AE/AE
+@ERR5069949.2098070 NS500628:121:HK3MMAFX2:3:11406:20552:1667/1
+GCTCTCTACTACCCTTCTGCTCGCATAGTGTATACAGCTTGCTCTCATGCCGCTGTTGATGCACTATGTGAGAAGGCATTAAAATATTTGCCTATAGATAAATGTAGTAGAATTATACCTGCACGTGCTCGTGTAGAGTGTTTTGATAAAT
++
+AAAAAEEEEEEE/EE/EA/EEA/EEEEA/<6EEEEEEA<AE/EAEEEE/AAA<A<E<EEA<EEA/EEAEEEEAEEAEE/EAEEEEE//6/6AEEA</EEAEEEEA<EE/AE66AEAEAAAA/<AA<<<<EAAEAEAEE/EA6EEEEEAEEE
+@ERR5069949.1552198 NS500628:121:HK3MMAFX2:1:11108:11797:17621/1
+ATATTACGCGTATACGCCAACTTAGGTGAACGTGTACGCCAAGCTTTGTTAAAAACAGTACAATTCTGTGATGCCATGCGAAATGCTGGTATTGTTGGTGTACTGACATTAGATAATCAAGATCTCAATGGTAACTGGTATGATTTCGGT
++
+AAAAAEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEAEAEEEEEEEEEEEAEEEEEEEEEAEEEEEEEEEEEEEEA
+@ERR5069949.2385514 NS500628:121:HK3MMAFX2:3:21611:13560:16508/1
+GAGTTGACATCTATGAAGTATTTTGTGAAAATAGGACCTGAGCGCACCTGTTGTCTATGTGATAGACGTGCCACATGCTTTTCCACTGCTTCAGACACTTATGCCTGTTGGCATCATTCTATTGGATTTGATTACGTCTATAATCCGTTT
++
+A/AAAE6EEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EE<EEEEE/A/AE6EEA/E/EEEEEA/AEEEEEEE<EEEAAEEEEEEE/EAEEEEEE/EEAAAEEE/<<AEEEEEEAEAEA<E/AAAAAAAA/<EE/EEAEEAE
+@ERR5069949.2270078 NS500628:121:HK3MMAFX2:4:21510:16801:19554/1
+ATAAACATGTTAGGGTAACCATTAACTTGATAATTCATTTTAAAACCCATCATAGAGATGAGTCTTCTATAGGTCATGTCCTTAGGTATGCCAGGTATGTCAACACATAAACCTTCAGTTTTGAATTTAGTGTCAACACTGAGGTGTGTAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEE/EEEAEEEEEE<EEEEEAEEEEE/<<EEA<EEEA<A<A<6AEE<AEEEEEEEEEAEEEE//AEEAEEEEEEA
+@ERR5069949.114870 NS500628:121:HK3MMAFX2:4:11410:23953:14505/1
+AAAGGAGCTGGTGGCCATAGTTACGGCGCCGATCTAAAGTCATTTGACTTAGGCGACGAGCTTGGCACTGATCCTTATGAAGATTTTCAAGAAAACTGGAACACTAAACATAGCAGTGGTGTTACCCGTGAACTCATGCGTGAGCTTAAC
++
+AAAAA/EAEEEEEEAEEEEEEEEE/EEEAEEEEEEEEEEEEEEEEAEAEAEEAEEEEEEEEEEEEEEEEAEAA<EE</AAEEEEEEEAEEEAEEEEEEEEEEEEEEEEEEEEEAEEA<AEEEEEEA<<AAE<EEEEEAE<AAEAAAEAEE
+@ERR5069949.2668880 NS500628:121:HK3MMAFX2:1:21104:5191:15623/1
+GTTTGTGGACCTAAAAAGTCTACTAATTTGGTTAAAAACAAATGTGTCAATTTCAACTTCAATGGTTTAACAGGCACAGGTGTTCTTACTGAGTCTAACAAAAAGTTTCTGCCTTTCCAACAATTTGGCAGAGACATTGCTGACACT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAE<EEAEEEEA<EAA<AAAEEEEEE<AEA<EAAEE
+@ERR5069949.257821 NS500628:121:HK3MMAFX2:2:21307:14790:19018/1
+TATGTAGCCATACTCCACTCATCTAAATCAATGCCCAGTGGTGTAAGTAATTCAGATACTGGTTGCAAAGTTTTTATGACAGCATCTGCCACAACACAGGCGAACTCATTTACTTCTGTACCGAGTTCAACTGTATAGG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEE6EEEAE/EE/EEEE<EEEE/EEEEEEEEEEEA/EEEEEEEAEEE/E<EE/EAAEEE/<E/<E<<E/A/EEAEEA<AAAEAE/<AE<EE<EA/A
+@ERR5069949.2243023 NS500628:121:HK3MMAFX2:2:21301:7521:3094/1
+TTTCAAGACTTGTAAATTGCAACTTGTCATAAAGGTCTCTATCAGACATTATGCAAAGTATGCCTACTTTTGCTCTGGTAATAGCAACATTAAATCTGTTTACATTACAAGAGTGAGCTGTTTCAGTGGTTTGAGTGAATATGACATAGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEEEAEEEEEEEEEEEEEEEEEAEEEE<EEEAEAAEEEEEEEEEA<AAEEEE<EEEEAEAEE/E<AE<EAEEAE<EE
+@ERR5069949.2605155 NS500628:121:HK3MMAFX2:1:21311:18812:15592/1
+GTTCTTACCTTTCTTTTCCAATGTTACTTGGTTCCATGCTATACATGTCTCTGGGACCAATGGTACTAAGAGGTTTGATAACCCTGTCCTACCATTTAATGATGGTGTTTATTTTGCTTCCACTGAGAAGTCTAACATAATAAGAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEE/EEEEEEEEEEEEEE<EEAEEEAEAEAEEEEEEEEAAEEEEE<EEAEAEEEAA<E<EAAE</E/AA
+@ERR5069949.1340552 NS500628:121:HK3MMAFX2:2:11309:7320:9079/1
+TGCCTGACCAGTACCAGTGTGTGTACACAACATCTTAACACAATTAGTGATTGGTTGTCCCCCACTAGCTAGATAATCTTTGTAAGCTTTAGCAGCATCTACAGCAAAAGCACAGAAAGATAATACAGTTGAATTGGCAGGCACTTCTGTT
++
+AAAAAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEAEEE<EEEEEEEEEAAEEEEEEEEAAEEEEEEEEEE/EEEAAEE<A/EEAAAAEA
+@ERR5069949.1561137 NS500628:121:HK3MMAFX2:2:11303:23159:8428/1
+ATGTAAGGCTTTGTTAAGTCAGTGTCAACATGTGACTCTGCAGTTAAAGCCCTGGTCAAGGTTAATATAGGCATTAACAATGAATAATAAGAATCTACAACAGGAACTCCACTACCTGGCGTGGTTTGTATGAAATCACCGAAATCATAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEE/EEEEEEEEEEEEEEEEEEAEEAEEAEEEEEEEEEEEEEEEE/EEEEEEEAEEEEEAAAEE/EE<EEAA/A<A<EEAEE/AAAE<EA
+@ERR5069949.2361683 NS500628:121:HK3MMAFX2:2:21103:8659:7427/1
+GATTTTTCCAGAGTTAGTGCTAAACCACCGCCTGGAGATCAATTTAAACACCTCATACCACTTATGTACAAAGGACTTCCTTGGAATGTAGTGCGTATAAAGATTGTACAAATGTTAAGTGACACACTTAAAAATCTCTCTGACAGAGT
++
+AAAAA/AA/EAEEAAEEEEAEE/EAEAAA<EEAAEEEEEE/EEEEEEEEEAAEEE/AEEE/EEEA<EEEAEEEEE<AAA<EEEEEEEEEAAEEA<AEEEEEEEEE<EEE/<AEEEEEAEEA<//A<EEAAEE6EAEEAAAEEAA6AEEA
+@ERR5069949.2521353 NS500628:121:HK3MMAFX2:3:11508:26680:5334/1
+CTTTTACAAGACTTCAGAGTTTAGAAATTGTGGCTTATAATGTTGTAATTAAGGGACACTTTGATGGACAACAGGGTGAAGTACCAGTTTCTATCATTAATAACTCTGTTTACACAAAAGTTGATGGTGTTGATGTAGAATTGTTTGAAA
++
+AAA/AE/6E6EEEEAEE/EE/EEE/EE/EA/EAEA//EEEEE6EAEAE/EEEEEE/EAE////EEA/EEEEEEEEEEEEEE///A/EEAEEEEEEEE<AEAEEE/AE/E<E/EEEEEA/E///AE/66AEEAEEE<E//E/EA/A<6AEE
+@ERR5069949.1261808 NS500628:121:HK3MMAFX2:3:11601:6846:14627/1
+GTGAAATTAGTATGGACAAATCACCTAATTTAGCATGGCCTCTTATTGTAACAGCTTTAAGGGCCAATTCTGCTGTCAAATAACAGAATAATGAGCTTAGTCCTGTTGCACTACGACAGATGTCTTGTGCTGCCGGAACTACACAAACT
++
+AAAA/AAA/EEEE/AE<EE//EE/EEEEEA<EAEEE66//EE/E/A6EEEA/E<<E/E//A/EE/AEEE6/E6EA<EEE/E//EAAEEEA/E//EEE/EA/A/A</<<E/AA<EEE<E<EEEE<A6</<A//A/AA/<<<AA<<A<AEA
+@ERR5069949.2734873 NS500628:121:HK3MMAFX2:4:21611:11682:9300/1
+ACATAGTGTAGGCAATGATGGATTGACTAGCTACACTACGTGCCCGCCGAGGAGAATTAGTCTGAGTCTGATAACTAGCGCATATACCTGCACCAATG
++
+AAAAAEEEEEEEEEEEEEEEEEEEE/EEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@ERR5069949.3017828 NS500628:121:HK3MMAFX2:1:21111:18284:17510/1
+ATGATGAACCGACGACGACTACTAGCGTGCCTTTGTAAGCACAAGCTGATGAGTACGAACTTATGTACTCATTCGTTTCGGAAGAGACAGGTACGTTAATAGTTAAT
++
+AAAAAE6EEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEE
+@ERR5069949.573706 NS500628:121:HK3MMAFX2:3:21407:23793:9689/1
+GTACGAACTTAAGCATGGTACATTTACTTGTGCTAGTGAGTACACTGGTAATTACCAGTGTGGTCACTATAAACATATATCTTCTAAAGAAACTTTGTATTGCATAGACGGTGCTTTACTTACAAAGTCCTCAGAATACAAAGGTCCTAT
++
+AAAAA6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEAAEEEAEEEEEEEEEEE
+@ERR5069949.1980512 NS500628:121:HK3MMAFX2:1:21303:8321:13259/1
+ATATTGAGTGTTGGGTATAAGCCAGTAATTCTAACATAGTGCTCTTGTGGCACTAGTGTAGGTGCACTTAATGGCATTACTGTATGTGATGTCAGCACAAAATAATCACCAACATTTAATTTGTAAGTTGTTGTACCTCGGTAAACAACAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEAEEEEEEEAEAEEEEAAEAEEEEEAEEEAEEEAEAEEEAEAAAEAEAAEEEEEAAEEEEEEAEEEEEEEEEE
+@ERR5069949.1014693 NS500628:121:HK3MMAFX2:1:11109:13494:20222/1
+GTCTATATGCCTGCTAGTTGGGTGATGCGTATTATGACATGGTTGGATATGGTTGATACTAGTTTGTCTGGTTTTAAGCTAAAAGACTGTGTTATGTATGCATCAGCTGTAGTGTTACTAATCCTTATGACAGCAAGAACTGTGTATGAT
++
+AAAAAAEEAEEE6EAE//E/EEE6AEAA/EAEAEE6/E//EAE/EEEEAEE/EEE/EAEEEEEAE/EEEEEAEEEEEAAEEAEEE/AE/EAEAEEEEEEEEEEEEEE/AE/E/E/<<<AA<E<AEE</EEEEA6<AEEAAAA//A//EEE
+@ERR5069949.3184655 NS500628:121:HK3MMAFX2:3:21606:7576:13093/1
+TTCCAGAAGAGCAAGGTTCTTTTAAAAGTACTGTTGTACCTCTAACACACTCTTGGTAGTGATAAAGCTCACAAGTAGCGAGTGTTATCAGTGCCAAGAAAAGAATAATTTTCATGTTCGTTTAATCAATCTCCATTGGTTGCTCTTCAT
++
+AAAAAEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEA/EEEEEEEEEEEEEAEEEEEEE/EAEAEEA<A/EE/E<E/EEEEEEEEEE<E<EEEEA/<<<A/E///AEE<EEE/<A6E/A/6A<6AE/<E6EAAA
+@ERR5069949.29668 NS500628:121:HK3MMAFX2:2:21202:17958:8710/1
+CTCCAAAGCCACGTACGAGCACGTCGCGAACCTGTAAAACAGGCAAACTGAGTTGGACGTGTGGTTTCTAGTTGTAACCAGGGACAAGG
++
+AA/AAEE6////EEEAEEE<EAEA/AEEEEEAEEEAEEEAE/A6EA/E/EA6A<A<AAA</////EEE6///<</EE/EE/E/6////E
+@ERR5069949.3258358 NS500628:121:HK3MMAFX2:3:21508:15641:1056/1
+GTTCCTAGGTAGTAGAAATACCATCTAGGACTGAGATCTTTCATTTTACAGTCACCACCACGAATTCGTCTTGTAGCTCTTCGGTAGTAGCCAATTTGGTCATCTGGACTGCTATTGGTGTTAATTGGAACGCCTTCTCCTCGAGGGAATT
++
+AAA/AAEAEEAEEEEEEA/EEEEEEE//EE/EAEE/EE/A6EEEAEAEE/EE//EEEE/AE/AAAE/AAEE/A/E//AEE/AA/EE/EE6EE<///E//AEE<EEA<EE//EE<<</<EAEEEEEE//</A<EE<</<EA/A/<A<A/EEE
+@ERR5069949.1476386 NS500628:121:HK3MMAFX2:2:21303:26270:9400/1
+TAATGACCCTGTGGGTTTTACACTTAAAAACACAGTCTGTACCGTCTGCGGTATGTGGAAAGGTTATGGCTGTAGTTGTGATCAACTCCGCGAACCCATGCTTCAGTCAGCTGATGCACAATCGTTTTTAAACGGGTTTGCGGTGTAAGTG
++
+AAAAA/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEAEEEEEEE/AEEE/EEEEEE/AEE/EEAE/EEE<EA/<EEA/EEEEE/EEEEAAEEEAAAAEEAEEE
+@ERR5069949.2415814 NS500628:121:HK3MMAFX2:1:11107:11754:7444/1
+GCTTTGAGTTGACATCTATGAAGTATTTTGTGAAAATAGGACCTGAGCGCACCTGTTGTCTATGTGATAGACGTGCCACATGCTTTTCCACTGCTTCAGACACTTATGCCTGTTGGCATCATTCTATTGGATTTGATTACGTCTATAATC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEAEEEEEEAEEEEEEEAEEEAEEEEEAEEEEEAEAAAE<AEA

tests/data/fastq/dna/sarscov2_2.fastq.gz

@@ -0,0 +1,400 @@
+@ERR5069949.2151832 NS500628:121:HK3MMAFX2:2:21208:10793:15304/2
+ATGTGTACATTGGCGACCCTGCTCAATTACCTGCACCACGCACATTGCTAACTAAGGGCACACTAGAACCAGAATATTTCAATTCAGTGTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGTCGGCGTTGTCCTG
++
+AAAAAEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEAEEEEEAAEEEEEEEEEAAEAAA<<EAAEEEEEEEAAA<<<AE
+@ERR5069949.576388 NS500628:121:HK3MMAFX2:4:11501:11167:14939/2
+GCATAGACGGTGCTTTACTTACAAAGTCCTCAGAATACAAAGGTCCTATTACGGATGTTTTCTACAAAGAAAACAGT
++
+AAAAA6EEAEEEEEAEEAEEAEEEEEEA6EEEEAEEAEEEEE6EEEEEEAEEEEA///A<<EEEEEEEEEAEEEEEE
+@ERR5069949.501486 NS500628:121:HK3MMAFX2:1:21305:3080:7962/2
+TTACAGAGCAAGGGCTGGTGAAGCTGCTAACTTTTGTGCACTTATCTTAGCCTACTGTAATAAGACAGTAGGTGAGTTAGGTGATGTTAGAGAAACAATGAGTTACTTGTTTCAACATGCCAATTTAGATTCTTGCAAAAGAGTCTTGAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEAA/EEE/<AEEAAEAEA</EEEAEAAAAAEE
+@ERR5069949.1331889 NS500628:121:HK3MMAFX2:3:11605:17658:18222/2
+CACAATTAGTGATTGGTTGTCCCCCACTAGCTAGATAATCTTTGTAAGCTTTAGCAGCATCTACAGCAAAAGCACAGAAAGATAATACAGTTGAATTGGCAGGCACTTCTGTTGCATTACCAGCTTGTAGAC
++
+AAAAAEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEE/EEEEEEAAEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<AEEEEEEEEEAEEEEEEEEAEEEEE/A
+@ERR5069949.2161340 NS500628:121:HK3MMAFX2:2:21107:15810:18317/2
+GGACAACGCCGACAAGTTCCGAGGAACATGTCTGGACCTATAGTTTTCATAAGTCTACACACTGAATTGAAATATTCTGGTT
++
+AA6AAAE6/EA/<E/E/<//EE///EAEAEEEE////6EEE//EEE/EAEEE/AEA/EA<EA</AEEE/EE/AEAE</E//A
+@ERR5069949.973930 NS500628:121:HK3MMAFX2:2:21105:5548:5659/2
+CCTTTTGATGTTGTTAGACAATGCTCAGGTGTTACTTTCCAAAGTGCAGTGAAAAGAACAATCAAGGGTACACACCACTGGTTGTTACTCACAATTTTGACTTCACTTTTAG
++
+AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEAE<EAEEEEEEAEEAEE
+@ERR5069949.2417063 NS500628:121:HK3MMAFX2:1:11310:4428:10115/2
+CAATAGTCCAGTCAACACGCTTAACAAAGCACTCGTGGACAGCTAGACACCTAGTCATGATTGCATCACAACTAGCTACATGTGCATTACCATGGACTTGACAATACGGATCATGGTTGCTTTGTAGGTTACCTGTAAAACCCCATTGTT
++
+AAAAAEEEEEEEEE/E6EEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEEEAEEEEEEEAEEEEEAEEEEEAAEAEE/EEAE/EEEAAAA6A<E/EAEAE/EEAEAE<//EEEEE<AEEA
+@ERR5069949.376959 NS500628:121:HK3MMAFX2:1:21205:12113:10369/2
+AGATAATAGATGGTAGAATGTAAAAGGCACTTTTACACTTTTTAAGCACTGTCTTTGCCTCCTCTACAGTGTAACCATTTAAACCCTGACCCGGGTAAGTGGTTATATAATTGTCTGTTGGCACTTTTCTCAAAGCTTTCGCTAGCATTT
++
+AAAAAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEE<EEEEEA<A<AAAAAAA<<<EAEEEE<AEEEAAAAA<AAAEAAAEE
+@ERR5069949.1088785 NS500628:121:HK3MMAFX2:1:21109:26270:17779/2
+GCATTTCTTCACAAAGCTTGTTTATGTCTACAGCACCCTGCATGGAAAGCAAAACAGAAAGTAGTGAAACCATTTTTTCAAAGGCTTCAGTAGTATCTTTAGCTAAGAGAATGTCATTGTGTAACTGGACACATTGAGCCCACAATTTAG
++
+AAAAAEEEEAEEEEEEEEEE/EEEEEEEEEEEEEEEEAEE/EEAEEE/AAEEEEEEEEEEEEEEEEEEEE/EEEEEEEAEEEEEEEEEAEEEE</EEE/EA/<EAEAEEEEEEAEEA<EEEEEAEEEA<AAAA/EEAE<E//E<EEEEEA
+@ERR5069949.1066259 NS500628:121:HK3MMAFX2:1:11312:18369:8333/2
+GTACAAAAATAGCCTAAGAAACAATAAACTAGCATTATACACTGAAGTGTATTACCAGTTATGAAGAAAATAGGGCAATACTCAACACACATAAAAACAATACCTCTGGCCAAAAACATGACAGTTGTAACTACACCTGAGTAGTTAGAAG
++
+AAAAAEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEA/EAA</EEEEEEEEEE/AE//A/E<AE<AA<AEEE/AAEAAAEEAEEA<A6AEEA<EEAEEEEEEEAA//EE
+@ERR5069949.2832676 NS500628:121:HK3MMAFX2:4:11602:20035:16369/2
+TGATCAACCTATCAATTTGCACTTCAGCCTCAACTTTGTCAAGACGTGAAAGGATATCATTTAAAACACTTGAAATTGCACCAAAATTGGAGCTAAGTTGTTTAACAAGCGTGTTTAAAGCTTGTGCATTTTGGTTGAC
++
+AA6AAEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEAEEE<AEEEEEEE<AAEEAEEAEEEEEEEE/EEEAAEEEEAAA<A<EEAEE/EAAEEEAEEEEA/<AEE6AEAE/<EEEE<A
+@ERR5069949.2953930 NS500628:121:HK3MMAFX2:3:21504:10316:13537/2
+CCCTTGGAGAGTGCTAGTTGCCATCTCTTTTTGAGGGTTATGATTTTGGAAGCGCTCTGAAAAACAGCAAGAAGTGCAACGCCAACAATAAGCCATCCGAAAGGGAGTGAGGCTTGTATCGGTATCGTTGCAGTAGCGCGAACAAAAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEE/EEEEAEEEEEEAAEEEEEEEEEEAEE/EEEEAEEEEEEEEEAEEEEEEEAEEEEAAEEEEEE/EEEEEEAAAAAEEEEAAEEEEEEEEEE/EEEEEE
+@ERR5069949.324865 NS500628:121:HK3MMAFX2:1:11102:17526:14721/2
+ATGAATCTGATGAATACATAGCTACTAATGGACCTCTTAAAGTGCGTGGTAGTTGTGATTAAAGCGGACACATACTTGCTAAACACTCTCTTCATGATGTC
++
+A/AAAEEEEA6AA6EE//EEA/EEEAE/EA/A////E</EEAA//EEA////EAE<///E/AEA</AAE/EA//E<EAAAE/AA//AEE//A/AE//</EE
+@ERR5069949.2185111 NS500628:121:HK3MMAFX2:2:11210:1473:2754/2
+GCTAAGCACTATGTGTACATTGGCGACCCTGCTCAATTACCTGCACCACGCACATTGCTAACTAAGGGCACACTAGAACAAGAATATATCAATTCAGTGTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGT
++
+AAAAA6/E//AEEAEEEEE/EEE/E/EEEE/E/E6E///E<E/<EE</A/<EEE//E6/<EEEAE<//<E/E/6EEE/EAAA<//EE//<EEE/AEEAA/A<EAE/EEA<//AEAEEE/<A/E</AA<EEA<<EAAEA<<</E/EEE
+@ERR5069949.937422 NS500628:121:HK3MMAFX2:2:11104:8644:10925/2
+GGTTAAAGTCATTAAGAGTTGTGGTAAATCGATTGAGAAACCACCTGTCTCCATTTATAACAGCAGCGTACAACCAAGCTAAAACATTAACTGTAATAGTTGTGTCCGTACCAGCTGCTTGTGCTGTTTGCCTGTCAACAAAAGGTCCATA
++
+AAAAAEEEAEEEEAEEEEEEEEEEEEEEEEEEEEEEEE6/EE/EEEEEEAAAAAEEEEEEAEEE<EEEAEE//EEEEEEEEEEAEEEEEEEAEAAEEAEA<EA<AEEEA/EEEAEAAAA<AAAAAAAAA6EA/AAE<EAAEEE/EE<A/AA
+@ERR5069949.2431709 NS500628:121:HK3MMAFX2:3:21609:20858:5879/2
+ATTATAGGATATTCAATAGTCCAGTCAACACGCTTAACAAAGCACTCGTGGACAGCTAGACACCTAGTCATGATTGCATCACAACTAGCTACATGTGCATTACCATGACTTGACAATACAGATCATGGTTGCTTTGTAGGTTACCTGTAAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAEEEAEEEEEEEEEEEEEEEEEEAEEAAAEAEEEEAAAA<EAAAAAAA
+@ERR5069949.1246538 NS500628:121:HK3MMAFX2:4:11505:20131:7933/2
+GTTGTAGTAAGCTAACGCATTGTCATCAGTGCAAGCAGTTTGTGTAGTACCGGCAGCACAAGACATCTGTCGTAGTGCAACAGGACTAAGCTCATTATTCTGTAATTTGACAGCAGAATTGGCCCTTAAAGCTGTTACAATAAGAGGCCAT
++
+AAAAAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAAEEEAEEAAA<AEEAAAEAAEEEEAAEEEAAAAAA
+@ERR5069949.1189252 NS500628:121:HK3MMAFX2:2:21308:10087:1427/2
+CAATTTTACTATCTGCATCTACAACCTGTTGGATTTCCCACAATGCTGATGCATAAGTAAATGTTGTACCATCACACGTATTTTTATATGTGTTATAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEE
+@ERR5069949.2216307 NS500628:121:HK3MMAFX2:1:21302:10479:19538/2
+GTGTAGACTTATGAAAACTATAGGTCCAGACATGTTCCTCGGAACTTGTCGGCGTTGTCCTGCTGAAATTGATGACACTGTGAGTGCTTTGGTTTATGATAATAAGCTTAAAGCACATAAAGACAAATCAGCTCAATGCTTTAAA
++
+A//AA6EEE/EEEEAEEEEE/EE/A/6E/EAE</EE/AEE/AEEAEAE/E//AA<6AA6<E/EE/EEE/AE/6/EEE<//E/E6A/<E/////EAE<E<</AAAEE//AAEE/6AAA6/AEA/AAEEAAAA////AA/6E<6AAE
+@ERR5069949.3273002 NS500628:121:HK3MMAFX2:4:11506:23944:16676/2
+TGGCAAGGAAGACCTTAAATTCCCTCGAGGACAAGGCGTTCCAATTAACACCAATAGCAGTCCAGATGACCAAATTGGCTACTACCGAAGAGCTACCAGACGAATTCGTGGTGGTGACGGTAAAATGAAAGATCTCAGTCCAAGATGGTA
++
+AAAAAEEEEEEEEEEEEEEEEEEAEEAEEEEEEEAEEE/EEA/EEEEEE6EEEAEEEEEEEEEAEEEEAEEEEE<EEEE<EEE/EE//AE/EEAEAEAE/EAAEEA6AEEE/<E<</EEE/E</AEA//A/EA6<AAEEEA/AEE6AE/E
+@ERR5069949.3277445 NS500628:121:HK3MMAFX2:4:11601:26589:14031/2
+CTCAGTTGCAACCCATATGATGCCGTCTTTGTTAGCACCATAGGGAAGTCCCGCTTCTGGCCCCGTTCCTAGGTAGTAGAAATACCATCTTGGACTGAGAT
++
+AAAAAEEEEE6EEEEEEEEEEEEEEE//6EEAE//A/AEE<EAAE/A/E/E//E/E///E/EE///</A/E/<//E<EA/AE/EA</<//<AE///A/</<
+@ERR5069949.3022231 NS500628:121:HK3MMAFX2:1:21102:10897:17190/2
+ATATTGCAGCAGTACGCACACAATCGAAGCGCAGTAAGGATGGCTAGTGTAACTAGCAAGAATACCACGAAAGCAAGAAAAAGAAGTACGCTATTAACTATTAACGTACCTGTCTCTTCCGAAACGAATGAGTACATAAGTTCGTAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEAEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEEEEEEEEEEEEEAAEAAEE<EAEEEE<EEEEEAEEEEEEAAAE
+@ERR5069949.184542 NS500628:121:HK3MMAFX2:2:21205:8408:1786/2
+GGCACCTTCTTTAGTCAAATTCTCAGTGCCACAAAATTCGCAAGTGGCTTTAACAAAATCGCACGT
++
+AAA/AAEAEE/AAAE6/A/EE6AAAAEEE/E//AA6EEE/6/A/A///EEA<EEEAE6A/E////E
+@ERR5069949.540529 NS500628:121:HK3MMAFX2:3:21609:5345:11242/2
+TGGGCACACTTTCTTATGAACAATTTAAGAAAGGTGTTCAGATACCTTGTACGTGTGGTAAACAAGCTACAAAATATCTAGTACAACAGGAGTCACCTTTTGTTATGATGTCAGCACCACCTGCTCAGTATGAACTTAAGCATGGTACAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEAEAEEEEEEEEEEAAEEEEEEEEEEAEEEEEEE<EA<AEEEAEEEAEEEEAEEAEEEAEEEEEEEAAAEEE
+@ERR5069949.686090 NS500628:121:HK3MMAFX2:1:11212:8601:4499/2
+GATCAGGCATTAGTGTCTGATGTTGGTGATAGTGCGGAAGTTGCAGTTAAAATGTTTGATGCTTACGTTAATACGTTTTCATCAACTTTTAACGTACCAATGGAAAAACTCAAAACACTAGTTGCAACTGCAGAAGCTGAACTTGCAAAGA
++
+AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEAEEEEEEEE/EEEEEEEEEEEEEEAEEEEEAEEEE<AEE/EEEEEEEAAAAEEEEEEEEEEEAEAEEEEEAEEEEA
+@ERR5069949.2787556 NS500628:121:HK3MMAFX2:4:11609:10602:17177/2
+CCGCTAACAGTGCAGAAGTGTATTGAGCAATCATTTCATCTGTGAGCAAAGGTGGCAAAACAGTAAGGCCGTTAAACTTTTGTGCACAAACGAGGTCTCTAGCAGC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAEEEEEEAEEEEE<AAA<EEEE<AEEEEEEEEAAEE
+@ERR5069949.2650879 NS500628:121:HK3MMAFX2:1:21111:16524:5801/2
+ATATAATTCCGCATCATTTTCCACTTTTAAGTGTTATGGAGTGTCTCCTACTAAATTAAATGATCTCTGCTTTACTAATGTCTATGCAGATTCATTTGTAATTAGAGGTGATGAAGTCAGACAAATCGCTCCAGGGCAAACTGGAAAGATT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEAEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEEEEEEEEEAEEE<A<EEEEAAAAEEEEEEEEEEEE
+@ERR5069949.2064910 NS500628:121:HK3MMAFX2:3:21608:24504:13289/2
+ACATACTGTTCTAATGTTGAAGTCACTTTGAATTTATCAAAACACTCTACACGAGCACGTGCAGGTATAATTCTACTACATTTATCTATAGGCAAATATTTTAATGCCTACTCACAT
++
+A//AAE/EEEA/EEE/A/EEE6EE/EEEEAE6AA<EAE<EEEEEE/E</<EA/6EEA<EAAEEEEEAEAE/AEAEE/EA<E//EEEEEEEE/EAAE<AA///AE<EAAA/A////<<
+@ERR5069949.2328704 NS500628:121:HK3MMAFX2:4:21612:9379:15958/2
+CATGGATTGGCTTCGATGTCGAGGGGTGTCATGCTACTAGAGAAGCTGTTGGTACCAATTTACCTTTACAGCTAGGTTTTTCTACAGGTGTTAACCTAGTTGCTGTACCTACAGGTTATGTTGATACACCTAATAATACAGATTTTTCCA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEAEEEEEEEEEEEEAEEEEEEAEEEEEEEEEEEAEEEAEAEEEA<AEEEEEEEE/EEEEEEEAEEEEEEE/EEEEAEE
+@ERR5069949.1067032 NS500628:121:HK3MMAFX2:4:21403:19263:10183/2
+GCTTTAGATCAAGCCATTTCCATGTGGGCTCTTATAATCTCTGTTACTTCTAACTACTCAGGTGTAGTTACAACTGTCATGTTTTTGGCCAGAGGTATTGTTTTTATGTGTGTTGAGTATTGCCCTATTTTCTTCATAACTGGTAATACA
++
+AAAAAAEEEE666<EEEE/E/EEEAEAEEEE/EEEEEEEE/AEEEAEE/<AAEE<EAEEEA/AEEE/EAEEEE<E<AEEAEEEE<<//EEE/EEE<EEE<A/A/A<EEE///A/6<A<AE<//<EEEE6</<AAAAAAE<A//<</<A/E
+@ERR5069949.3338256 NS500628:121:HK3MMAFX2:1:11108:15722:2450/2
+CAGCAAACTGTGACTCTTCTTCCTGCTGCAGATTTGGATGATTTCTCCAAACAATTGCAACAATCCATGAGCAGTGCTGACTCAACTCAGGCCTAAACTCATGCAGACCACACAAGGCAGATGGGCTATATAAACGTTTTCGCTTTTCCG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEAEEEEEEEEAAAAAEA<AAAEA<AA
+@ERR5069949.1412839 NS500628:121:HK3MMAFX2:3:11609:15519:12102/2
+GTACACACACTGGTACTGGTCAGGCAATAACAGTTACACCGGAAGCCAATATGGATCAAGAATCCTTTGGTGGTGCATCGTGTTGTCTGTACTGCCGTTGCCACATAGATCATCCAAATCCTAAAGGATTTTGTGACTTAAAAGGTAAGT
++
+AAAA6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEAEEEAAAAEEEEEEEEEEAAEEAEAE<EEEAEAEEE/<AAAEAEAA/EAEEEEAEEAAE/AEA/EEEAEEAEAA
+@ERR5069949.1538968 NS500628:121:HK3MMAFX2:2:11304:14669:15804/2
+CACGATGGCAGACCTCGTCTATGCTTTAAGGCATTTTGATGAAGGTAATTGTGACACATTAAAAGAAATACTTGTCACATACAATTGTTGTGATGATGATTATTTCAATAAAAAGGACTGGTATGATTTTGTAGAAAACCCAGATATATTA
++
+AAAAAEEEAEEAEEEAEEEAEEEAE<EEE6EAEA<EAAAEEEEEEEEEEEEEA/</EEEEEEEEEEEEEEEEEEEEEEAEEEEE/AEEEEEEEEAEEEEEEEEEEEEEEEEAEEEAA<AEAEE<AAE<A<AEEEEE/EA6AAA/EE/EEEA
+@ERR5069949.147998 NS500628:121:HK3MMAFX2:3:21611:12530:12698/2
+ATCTAAAGTCATTTGACTTAGGCCACGAGCTTGGCACTGATCCTTATGAAGATTTTCAAGAAAACTGGAACACTAAACATAGCAGTGGTGTTACCCGTGAACTCATGCGTGAGCTTAACGGAGGGGCATACACTCGCTATGTCGATAACAA
++
+AAAAAAE6EEE/EEE/E/EEA6E/EEE/AE/E//EEEEAA/E/E/EAEEEE/EEEEEEE<EE/E/A/A</<E</<AE<///A<AA<//E/AE/E/EEE/EEEEA//E/A</<AE/////E<AAE<EE//EA/<6/A</A//<AAAA<EE/A
+@ERR5069949.366975 NS500628:121:HK3MMAFX2:1:11105:12446:15727/2
+CTAAAAAGGCTGGTGGCACTACTGAAATGCTAGCGAAAGCTTTGAGAAAAGTGCCAACAGACAATTATATAACCACTTACCCGGGTCAGGGTTTAAATGGTTACAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@ERR5069949.1372331 NS500628:121:HK3MMAFX2:3:11601:26764:8562/2
+TCTACAAGCTGGTAATGCAACAGAAGTGCCTGCCAATTCAACTGTATTATCTTTCTGTGCTTTTGCTGTAGATGCTGCTAAAGCTTACAAAGATTATCTAGCTAGTGGGGGACAACCAATCACTAATTGTGTTAAGATGTTGTGTACACA
++
+AAAAAEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEE/EEEEEEE<EEAEEEEEEA/EEA<EAEEEEEEAEAAAEEAAAEEEEEAAA<AAAEEEA
+@ERR5069949.1709367 NS500628:121:HK3MMAFX2:1:11301:9990:2022/2
+GGTGGCTGTATTAATGCTAACCAAGTCATCGTCAACAACCTAGACAAATCAGCTGGTTTTCCATTTAATAAATGGGGTAAGGCTAGACTTTATTATGATTCAATGAGTTATGAGGATCAAGATACACTT
++
+AAAAAEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEEEEEEEEEEEAEEEEEEEEEEEEEEE
+@ERR5069949.2388984 NS500628:121:HK3MMAFX2:2:21306:17580:3896/2
+AGCTACATGTGCATTACCATGGACTTGACAATACAGATCATGGTTGCTTTGTAGGTTACCTGTAAAACCCCATTGTTGAACATCAATCATAAACGGATTATAGACGTAATCAAATCCAATAGAATGATGCCAACAGGCATAAGTGTCTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEEEEEEAEEEEE/<EE<EE<AEEEAEEEEAAAEEEEEE<<AAEAEEEEEEEEAEEEEEAAEEEAEAAAA<EA/EEAAAAEEEAEE<EEE6AAAAE
+@ERR5069949.1132353 NS500628:121:HK3MMAFX2:2:21111:21744:6806/2
+CAACAGGGCAACCTTACAAGCTATAGCCTCAGAGTTTAGTTCCCTTCCATCATATGCAGCTTTTGCTACTGCTCAAGAAGCTTATGAGCAGGCTGTTGCTAATGGTGATTCTGAAGTTGTTCTTAAAAAGTTGAAGAAGTCTTTGAATGTG
++
+AAAAAEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEEEE<EEEEEEEEAAEEAAEEEEEEEEEEEEEAE<AAAAAAE/AEEAEEEEEEEEEEEEEAAAEAAEEA
+@ERR5069949.1151736 NS500628:121:HK3MMAFX2:4:21406:11928:13664/2
+TTTTGCTACTGCTCAAGAAGCTTATGAGCAGGCTGTTGCTAATGGTGATTCTGAAGTTGTTCTTAAAAAGTTGAAGAAGTCTTTGAATGTGGCTAAATCTGAATTTGACCGTGATGCAGCCATGCAACGTAAGTTGGAAAAGATGGCTGAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEEEE<EEEEAEEEEEEEEEEEEEEEEAEEEEEAEEAAEEE<AEAEEE<A/AAEEEEEEEAAAAA<AAAE<EEEEAEEEAEEEEEEAEEAEA/A
+@ERR5069949.479807 NS500628:121:HK3MMAFX2:3:11402:16906:16736/2
+GTTTACAACAGTAGACAACATTAACCTCCACACGCAAGTTGTGGACATGTCAATGACATATGGACAACAGTTTGGTCCAACTTATTTGGATGGAGCTGATGTTACTAAAATAAAACCTCATAATTCACATGAAGGTAAAACATTTTATGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<AEEEEEEEEEE<EEE<AEE/EEEEEEEEEAEEE<AA/EAA<AEEEEEEEAEEAAA
+@ERR5069949.2176303 NS500628:121:HK3MMAFX2:4:11608:21150:1680/2
+ATAACACCCTTATAAAACATTTTAAAGCATTGAGCTGATTTGTCTTTATGTGCTTTAAGCTTATTATCATAAACCAAAGCACTCACAGTGTCAACAATTTCAGCAGGACAACGCCGACAAGTTCCGAGGAACATGTCTGGACCTATAGTT
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAAEEEEEAEAEEAEAEAEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAEEEEEEAEEEEEEEAAEEAAE/<A/EEAEE/A<AAAA<AAAEEEEEEAAAEAEAEEEAEAE
+@ERR5069949.2772897 NS500628:121:HK3MMAFX2:1:11106:12629:15009/2
+CTTTCGTTGCAATATGGCAGTTTTTGTACACAATTAAACCGTGCTTTAACTGGAATAGCTGTTGAACAAGACAAAAACACCCAAGAAGTTTTTGCACAAGTCAAACAAATTTACAAAACACCACCAATTAAAGATTTTGGTGGTTTTAAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEA<EAEEEE<EEEEAEEAAEEEEEEEEEEE
+@ERR5069949.1020777 NS500628:121:HK3MMAFX2:3:11408:13234:17983/2
+CTATATGCCTGCTAGTTGGGTGATGCGTATTATGACATGGTTGGATATGGTTGATACTAGTTTGTCTGGTTTTAAGCTAAAAGACTGTGTTATGTATGCATCAGCTGTAGTGTTACTAATCC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEEEEEEAEEEAAEEEEEEEEEAEEEEA
+@ERR5069949.465452 NS500628:121:HK3MMAFX2:3:11503:20241:8941/2
+GTTGTAAACACCTTAATAGTCCTCACTTCTCTCAAAGAAAGAAGTGTCTTAAGATTGTCAAAGGTGATAACTTCACCATCTAGGTGGAATGTGGTAGGATTACTAGTGTAATATACACTTTTATCACCTCTCTTAAGAAATTCTATACCT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEE/AEE/EEEEEEEAEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAE/AA<AE</AAEEEEEEEEEEEEAAA
+@ERR5069949.1704586 NS500628:121:HK3MMAFX2:3:11408:18502:19140/2
+CTTGGTTAGCATTAATACAGCCACCATCGTAACAATCAAAGTACTTATCAACAACTTCAACTACAAATAGTAGTTGTCTGATATCACACATTGTTGGTAGATTATAACGATAGTAGTCATAATCGCTGATAGCAGCATTACCATCCTGAG
++
+AAA/AEEE<EEEEE/EEEEEE/EEEEEEEEEEEAEEEEEE/EEEEEEEEEEEEEE<EAEEEEEE/EEEAEEAEAEEEAEEAEEEE/EA<AEA<AAEAEEEE<EEEEEA//E<EA<<A<A//EE6AAA</<EEE6AAEA<AAA<EAEE//A
+@ERR5069949.1258508 NS500628:121:HK3MMAFX2:2:21106:23201:10887/2
+GTGTTCCCTTGAACATAATACCTCTTACAACAGCAGCCAAACTAATGGTTGTCATACCAGACTATAACACATATAAAAATACGTGTGATGGTACAACATTTACTTATGCATCAGCATTGTGGGAAATCCAACAGGTTGTAGATGCAGATAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEE/EEEEEEEEEAEEEEEEEEA/AEAEAAAEAEEEEEAEEEEE<<EEEEEEE<AEAEAAA<EEE
+@ERR5069949.986441 NS500628:121:HK3MMAFX2:2:11310:5407:13956/2
+TGTCATAATACGCATCACCCAACTAGCAGGCATATAGACCATATTAAAATAAGCTACAGTGGCAAGAGAAGGTAACAAAAACAAACAGAGAAATGCATGCTTATGTTTGACAAACATCATTGCAAAAGCAGACATAGCAATAATACCCAT
++
+AAAAAEEEEEEEEAEEEEEEEEAEEEEEAEEEEEEAE/EEEEEAEAEEEEEEEAEAEEEEEEAE/AAAEA/EAEAE/EEEAEEEEEEE/EEEAEAAAE</AAEAEEEE/E/E/E<AEEAE/<AEAEEAA<AEAAEAE<E<EEEEEA<A6A
+@ERR5069949.2674295 NS500628:121:HK3MMAFX2:1:21303:10381:8722/2
+CCTAAAAAGTCTACTAATTTGGTTAAAAACAAATGTGTCAATTTCAACTTCAATGGTTTAACAGGCACAGGTGTTCTTACTGAGTCTAACAAAAAGTTTCTGCCTTTCCAACAATTTGGCAGAGACATTGCTGACACTACTGATGCTGTC
++
+AAAAAEEEEEEEEEEAEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEE<EEEEEEAEEE/EEEEEEE<EEAEEAEEEEEAEEAEE<EEAEE<E/AAAAAAE<AEAAAEEEEAEEAEAEE/EEAAAEE
+@ERR5069949.885966 NS500628:121:HK3MMAFX2:4:11610:19682:20132/2
+CTTAGGTCTTAGGATTGGCTGTATCAACCTTAAGCTTAAGTACACAATTTTGCATAGAATGTCCAATAA
++
+A//AA6EEAEEEEE6EEE/EEA/EA///AAE/EAEEEAE6AE/E/E/EEAAE/EAA/E/E/<EA//E/6
+@ERR5069949.2342766 NS500628:121:HK3MMAFX2:4:21501:7762:9816/2
+TCATAGATGTCAACTCAAAGCCATGTGCCCATAAGACAAATACGACTCTGTCAGAGAGATTTTTAAGTGTGTCACTTAACATTTGTACAATCTTTATACGCACTACATTCCAAGGAAGTCCTTTGTACATAAGTGGTATGAGGTGTTTAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAAAEEEEEEEEEEEAEEEEEAEEEAAAA<<AEEEEAEAEAE
+@ERR5069949.3122970 NS500628:121:HK3MMAFX2:2:21304:20157:2460/2
+ATCAAGGACCTGCCTAAAGAAATCACTGTTGCTACATCACGAACGCTTTCTTATTACAAATTGGGAGCTTCGCAGCGTGTAGCAGGTGACTCAGGTTTTGCTGCATACAGTCGCTACAGGATTGGC
++
+AAAAAEE6EEEEEEEAEEEEEEEEEEEAEEEEEEEEEEEEEE/EEEEEAE<EEAEAEEEEEEEEEAAEEEEAEAEEE/AEEE<A<A/AAAAAE/E<A66AEEEEEEEEEEEAE<</6AA<A/6/EA
+@ERR5069949.3279513 NS500628:121:HK3MMAFX2:1:21110:4328:5711/2
+TCTTTGTTAGCACCATAGGGAAGTCCAGATTCTGGCCCAGGTCCTAGGTAGTATAAATACCATCTTTG
++
+AAAA/AEEEEEE/EEAEEEEE/EEEEAE/EEEEAA/<EAE6EA//AEEA/A/</EE6<E//EE//E/E
+@ERR5069949.309410 NS500628:121:HK3MMAFX2:1:11107:4644:2079/2
+TTAGTAGCCGTATTTAAGGCTCCTGCAACACCTCCTCCATGTTTAAGGTAAACATTGGCTGCATTAACAACCACTGTTGGTTTTACCTTTTTAGCTTCTTCCACAATGTCTGCATTTTTAATGTATACATTGTCAGTAAGTTTTAAATAA
++
+AAA/AEE/E/EE66E//EEEEEEEE/AE///EEEEEAEEEE/AAE/EEAAEE/EE/</<E//EE</EAEE6EEE/E</EEEEEEAEE/EEAEEEEEEAEEE<EAAEE/EA//<A</E/6</A/<AEAAAAEAAAA/EAAE<<AEAEE//E
+@ERR5069949.532979 NS500628:121:HK3MMAFX2:2:21104:15575:19089/2
+TGTTTATAGTGACCACACTGGTAATTACCAGTGTACTCACTAGCACAAGTAAATGTACCATGCTTAAGTTCATACTGAGCAGGTGGTGCTGACATCATAACAAAAGGTGACTCCTGTTGTACTAGATATTTTGTAGCTTGTTTACCACACG
++
+AAAAAEEEEEEEEEEEEEEEEAEEEEAEEEEEEE/EEE6EAEEE/EEEEEEEEEEEAEE/EEEEEEEEEEE/EAEEEEEE/EEEEEEEE<EA<AEA/EE/EEEAEEA<EEEE</<AAAEA/EEAEA<AEAEEEE<AE<AE/AA/E</E/EA
+@ERR5069949.2888794 NS500628:121:HK3MMAFX2:3:11604:17798:16811/2
+TCCCTGCACAAGAAAAGAACTTCACAACTGCTCCTGCCATTTGTCATGATGGAAAAGCACACTTTCCTCGTGAAGGTGTCTTTGTTTCAAATGGCACACACTGGTTTGTAACACAAAGGAATTTTTATGAACCACAAATCATTACTACAG
++
+AAAAAEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEAEE<<<6AEE</AAAEEEEEEEAA<EEAAEA
+@ERR5069949.2205229 NS500628:121:HK3MMAFX2:1:11212:26331:4199/2
+CGTGTAAGGAATTCTCTTACCACGCCTATTTGTGGCCTGTTAATTGCAGATGAAACATCATGCGTGATAACACCCTTATAAAACATTTTAAAGCATTGAGCTGATTTGTCTTTATGTGCTTAAGCTTATTATCATAAACCAAAGCACTC
++
+AAAAAEEEEAEEEEEEEAEEA/EEEEEEEEEE<EEEEEEA/EEEEEEE/EEEAEEEAEEEEEEEEAEEE/E/E/EAEEEAEEEEEEEEEEE/EEE/E<EEEAEE<<E<//<AAEEE<AEEAEEEAAAA/AEAEEAEE/EEE/EE<<A<A
+@ERR5069949.786562 NS500628:121:HK3MMAFX2:2:11109:4461:9003/2
+GCATTTCTTACCTAGAGTTTTTAGTGCAGTTGGTAACATCTGTTACACACCATCAAAACTTATAGAGTACACTGACTTTGCAACATCAGCTTGTGTTTTGGCTGCTGAATGTACAATTTTTAAAGATGCTTCTGGTAAGCCAGTACCATA
++
+AAAAAEEEEEEEEAEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEAEEEEEEEEEEEEEEEAEEEEEEE<EEAEEEEE<EEEEEEEEEAEEEAEEEAEEAA6A<EEEEAAEEEEAA/AEEEEEE/EEEEEEE
+@ERR5069949.919671 NS500628:121:HK3MMAFX2:4:21608:23884:15511/2
+CCTTAATGGTTCATGTGGTAGTGTTGGTTTTAACATAGATTATGACTGTGTCTCTTTTTGTTACATGCACCATATGGAATTACCAACTGGAGTTCATGCTGGCACAGACTTAGAAGGTAACTTTTATGGACCTTTTGTTGACAGGCAAAC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAE/EEEEEEAEEEEEEEEAEEEEEEEEEEEAEEEEEEEAEEEEAEEEEAEEEEE6EEEEEEEAAEAEEEEEEE<EEEEEEE6AAEEAEEEAA6AEEAAAAAEEAAEEEAEAEEE
+@ERR5069949.1328186 NS500628:121:HK3MMAFX2:3:21510:25432:14327/2
+GTACTATCTATACAGAACTGGAACCACCTTGTAGGTTTGTTACAGACACACCTAAAGGTCCTAAAGTGAAGTATTTATACTTTATTAAAGGATTAAACAACCTAAATAGAGGTATGGTACTTGGTAGTTTAGCTGCCACAGTACGTCTACA
++
+AAAAAEEEEEEE/EEAEEEEAEEEEEAEEEEEEEEAEEEEEEEEEEEEEAE/EEEEEEAEE/EEEEEEEEEEEEEEEEEAAAEA/EEEEEEEAAEEEEE/EEEEAEEEEEAAEEEE/AAAE<A<EEEE6AEEAAA<<<<AA<AE/EEAEEA
+@ERR5069949.870926 NS500628:121:HK3MMAFX2:4:11510:4100:9432/2
+AGCTAACACTGAAAAAGTCTGTCCTGGTTGAATGCGAACAAACTTATACTTAGGTGTCTTAGGATTGGCTGTATCAACCTTAAGCTTAAGTACACAATTTTGCATAGAATGTCCAATAACCCTGAGTTGAACATTACCAGCCTGTACCAAG
++
+AAAA6/EEEE6EEEEEEEEEEEE/EEEEEEEE/A/E/EEEEEEEEEEAEEEAEEEEEEE/AAE/AE<AE/<AEEE//AEEEEEEE<EEEE<<E<EEAEEE/E/<EEE/</</</<AAE/EA<AAAAEEEE/EE/</6AE/6<////EAAAA
+@ERR5069949.2257580 NS500628:121:HK3MMAFX2:3:11512:13762:8305/2
+TCTATAGGTCATGTCCTTAGGTATGCCAGGTATGTCAACACATAAACCTTCAGTTTTGAATTTAGTGTCAACACTGAGGTGTGTAGGTGCCTGTGTAGGATGTAACCCAGTGATTACCTTACTACAATCTTTAAAGAGTCCTGTTACATT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEAEEEAEAEEEEAEEEEEEAAAAEEEEAAEEAEEEEEE
+@ERR5069949.3249622 NS500628:121:HK3MMAFX2:2:21101:6979:17698/2
+ATCATGACGTTCGTGTTGTTTTAGATTTCATCGAAACGAACAAACAAAAATGTCTGATAATGGACCCCAAAATCATCGAAATGCACCCCGCATTACGGTTGGTGGACCCTCCGATT
++
+AAA/AE//EEE/EE6AE/A</EE//6AE6EE//EE/AE//A/EE//EEEE<EAA/EE//<E/A/E/EE//E/A/E/E//EE/<A/A<EE/A//</EE//E/E//A/EEE///A//6
+@ERR5069949.611123 NS500628:121:HK3MMAFX2:4:21402:10868:14846/2
+ATTATACTTAAACCAGCAAATAATAGTTTAAAAATTACAGAAGAGGTTGGCCACACAGATCTAATGGCTGCTTATGTAGACAATTCTAGTCTTACTATTAAGAAACCTAATGAATTATCTAGAGT
++
+AAAAAEEEEEA6EEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEE<EEEEEEEAEEEEEEEEEEEEAEEEEEEEEEEE/EEEEEEEA/AAEAAEAAEAE
+@ERR5069949.651338 NS500628:121:HK3MMAFX2:2:11307:9572:4021/2
+ATTCTGTGCTGGTAGTACATTTATTAGTGATGAAGATGCGAGAGACTTGTCACTACAGTTTAAAAGACCAATAAATCCTACTGACCAGTCTTCTTACATCGTTGATAGTGTTACAGTGAAGAATGGTTCCATCCATCTTTACTTTGATA
++
+AAAAAE/EAEEE/AEAEEE/EEEAAEEEEAEEEEE/EEEEAEEEEEEAEE/EEEE/EEE</EE/AEAE/<E/EEAEE<EEEE//AEEEEEE<EEAEE/EEE//E/<EE<A<A/EAA<AA/AEEA//A<A/A<A<6A6/AEE/AEEA<AE
+@ERR5069949.169513 NS500628:121:HK3MMAFX2:2:21303:6012:2171/2
+CATTTGGTGACGCAACTGGATAGACAGATCGAATTCTACCCATAAAGCCATCAAGCTTTTTCTTTTCAACCCTTGGTTGAATAGTCTTGATT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEE/EEEEEEEEEEEEEEEEEEEEEEE
+@ERR5069949.155944 NS500628:121:HK3MMAFX2:2:21111:24075:19469/2
+GTACACGGAACGTTCTGAAAAGAGCTATGAATTGCAGACACCTTTTGAAATTAAATTGGCAAAGAAATTTGACACCTTCAATGGGGAATGTCCAAATTTTGTATTTCCCTTAAATTCCATAATCAAGACTATTCAACCAAGGGTTGAAAA
++
+AAAA/EEEEEEAEEEEEEEEEEEE/EEEEEEEEEE/EAEEEEEEEEEEEEEEEAEEEAEE/AEEEEEEAAEEEEEEAEAEEEE/AEE/<EAE/E<EEA<<<AAEEAEEE<AA<EE/EAAEEEE<<<EEEA/AEAEE6</EEA<AEEE<<E
+@ERR5069949.2033605 NS500628:121:HK3MMAFX2:4:21508:23814:9228/2
+GTACTGGTAAGAGTCATTTTGCTATTGGCCTAGCTCTCTACTACCCTTCTGCTCGCATAGTGTATACAGCTTGCTCTCATGCCGCTGTTGATGCACTATGTGAGAAGGCATTAAAATATTTGCCTATAGATAAATGTAGTAGAATTATA
++
+AAAAAEAEEEEEEEEEEEEAEEEEEEEEEEEEE<EEEEEEEEEEE<EEEEEEEEAAEAEEEEEEEEEEEE/EAAEEEA/EEEEEEAE<EEEEEEEEE<AEEEEAAAEAE<EAEEEEEE//</A/AEAAAEA/<E<AEEEAEE<EEEEEE
+@ERR5069949.2730382 NS500628:121:HK3MMAFX2:4:21605:26351:3903/2
+ACTCATATGAGTGTGACATACCCATTGGTGCAGGTATATGCGCTAGTTATCAGACTCAGACTAATTCTCCTCGGCGGGCACGTAGTGTAGCTAGTCAATCCATCATTGCCTACACTATGTCACTTGGTGCAGAAAATTCAGT
++
+AAAAAEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE<EE/EEEEEEE/EEAEEEEEEEEEEEEEEEEEEAEEEEA<AEA<<EA<A<AEEEEEA<EAE<66A/AEEEEEEEAE<AAEA
+@ERR5069949.2125592 NS500628:121:HK3MMAFX2:1:11307:14726:9430/2
+TGTACACATAGTGCTTAGCACGTAATCTGGCATTGACAACACTCAAATCATAATTTGTGGCCATTGAAATTTCATCAAAGACAACTATATCTGCTGTCGTCTCAGGCAATGCATTTACAGTACAAAAGACATACTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAEE/EEE<EEEEEEEEEEEE//E<EE//A/<//6AA//</</</<EE//E//<//A
+@ERR5069949.1062611 NS500628:121:HK3MMAFX2:1:11304:13814:5808/2
+TGGTAATGCTTTAGATCAAGCCATTTCCATGTGGGCTCTTATAATCTCTGTTACTTCTAACTACTCAGGTGTAGTTACAACTGTCATGTTTTTGGCCAGAGGTATTGTTTTTATGTGTGTTGAGTATTGCCCTATTTTCTTCATAACTGGT
++
+AAAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAAEEE<AAEEEEEEEEAAEEAAA6AEEEEEEEEEEEEEEEEEA
+@ERR5069949.1778133 NS500628:121:HK3MMAFX2:1:11203:21125:7007/2
+TGCCACAACTGCTTATGCTAATAGTGTTTTTAACATTTGTCAAGCTGTCACGGCCAATGTTAATGCACTTTTATCTACTGATGGTAACAAAATTGCCGATAAGTATGTCCGCAATTTACAACACAGACTTTATGAGTGTCTCTATAGAAA
++
+AAAAAEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAAEEEEEE<AEEEEEEE/AAAE<AAEEAAEEEA<EAAEEEA<AAEEEEEE/EEAAAEE/EAAAAEEEEEAEAEE
+@ERR5069949.3057020 NS500628:121:HK3MMAFX2:2:21204:6342:14557/2
+GCGTACAACCAAGCTAAAGTTACTGGCCATAACAGCCAGAGGAAAATTAACTTAATTATATACAAAAACCTATTCCTGTTGGCATAGGCAAATTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@ERR5069949.2972968 NS500628:121:HK3MMAFX2:4:21405:4877:14897/2
+CTTTGAAGCAAGGTGAAATCAAGGATGCTACTCCTTCAGATTTTCTTCGCGCTACTGCAACGATACCGATACAAGCCTCACTCCCTTACGGATGGCTTATTGTAGGCGTTGCACTTCTAGCTGTTTTTCAGAGCGCTTCCAAAAACA
++
+AAAAAEEEEEEA/E/EEEEEEE/EEEEEEEAEEEEE/EEE/E/E/<66<6/EEEE///EEE<<E/<AEAEEEE/EEEEE6A/EEA/E//AE<EEE</6E/E/E/A<AAEA/A//AE<A/E<EE//6A/A<AAA6<E<<A66</E/AA
+@ERR5069949.2734474 NS500628:121:HK3MMAFX2:4:11507:10751:13427/2
+CCCATTGGTGCAGGTATATGCGCTAGTTATCAGACTCAGACTAATTCTCCTCGGCGGGCACGTAGTGTAGCTAGTCAATCCATCATTGCCTACACTATGTCACTTGGTGCAGAAAATTCAGTTGCTTACTCTAATAACTCTATTGCCA
++
+AAAA/EEEEEEEEE/E/EE6EEEEAEEEEEEAEEEEE/EEEEEEEEEEEAE/EAEE/EEEEEAE/EE<EAEEEEEEA/E<EEEEAE/EA<EEEEAEE/E/EE<EEEEE</EE/E//<<<AA6A<A<A/<AE/AE/EEEA6<A6A/</A
+@ERR5069949.856527 NS500628:121:HK3MMAFX2:2:11301:22893:2784/2
+ATACTTAGGTGTCTTAGGATTGGCTGTATCAACCTTAAGCTTAAGTACACAATTTTGCATAGAATGTCCAATAACACTTAGTTGAACATTACCAGCCTGTACCAAGAAATTATGATTAGACTTACGAATGAGTAATCTTCATAATTAGGGT
++
+AAAAAAAE/EEEEAEEEEEEEEAEEEEEE<E/<<EEEEEEAEEAAAEEEAE</EA<<EE/EAEEE6E//EE/EEE/6A/EEEEEEEEEAE<E/E<AAA/EEEEA/A<EEEEE<6EE<A/A<A<</6A<A<E<EEEE/AE//<6A/A6///A
+@ERR5069949.2098070 NS500628:121:HK3MMAFX2:3:11406:20552:1667/2
+CATAGTGCTTAGCACGTAATCTGGCATTGACAACACTCAAATCATAATTTGTGGCCATTGAAATTTCATCAAAGACAACTATATCTGCTGTCGTCTCAGGCAATGCATTTACAGTACAAAAGACATACTGTTCTAATGTTGAATTCACT
++
+AAAAA/EEEEEEEEEEEEEEEEEEEEE/EEEEEEEAAEEEEAAEEEEA/E/E<<EEEEEEAEEEE<AE/AEEEEAA<EEEE/EE6A<AEAEEAEA/<E/AEEA/AAA6A/E<A/AEAAEEE/EE/AE<<//EAAEE6/</6EEE6/A6A
+@ERR5069949.1552198 NS500628:121:HK3MMAFX2:1:11108:11797:17621/2
+GGTCAAGGTTAATATAGGCATTAACAATGAATAATAAGAATCTACAACAGGAACTCCACTACCTGGCGTGGTTTGTATGAAATCACCGAAATCATACCAGTTACCATTGAGATCTTGATTATCTAATGTCAGTACACCAACAATACCAGCA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEA<AA/<EE//AEA/EEE</E/AEA/EE/EEAEEEEE/</EEA<AE/E/E<EAAA6AA<///<A/AAAEA<EAA/EA
+@ERR5069949.2385514 NS500628:121:HK3MMAFX2:3:21611:13560:16508/2
+GTGCATTACCATGGACTTGACAATACAGATCATGGTTGCTTTGTAGGTTACCTGTAAAACCCCATTGTTGAACATCAATCATAAACGGATTATAGACGTAATCAAATCCAATAGAATGATGCCAACAGGCATAAGTGTCTGAAGCAGTGG
++
+A/AAAEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEEEAAEEEEAAEEEEEEEEEEAAEEEEAEAE/EEEEEAAEEEEAEEEEAEEEAEAEAEEEAA<<EEAAEEAA/EEAAEE<AA<EAAE<AA
+@ERR5069949.2270078 NS500628:121:HK3MMAFX2:4:21510:16801:19554/2
+CTTTATGACAAGTTGCAATTTACAAGTCTTGAAATTCCACGTAGGAATGTGGCAACTTTACAAGCTGAAAATGTAACAGGACTCTTTAAAGATTGTAGTAAGGTAATCACTGGGTTACATCCTACACAGGCACCTACACACCTCAGTGTTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEAEEEEEEEEEEEAEEAAEEEEAEEEEEAEEAEEEAAEAEEEEEAEAEEAEA<EAA
+@ERR5069949.114870 NS500628:121:HK3MMAFX2:4:11410:23953:14505/2
+GCATGAAGCTTTACCAGCACGTGCTAGAAGGTCTTTAATGCACTCAAGAGGGTAGCCATCAGGGCCACAGAAGTTGTTATCGACATAGCGAGTGTATGCCCCTCCGTTAAGCTCACGCATGAGTTCACGGGTAACACCACTGCTATGTTT
++
+AAAAAAEEEEEEEAEAEE/EEEEEEEEEEEEEEEEAEEEAAEEEEEEAAEAAEAEEEEEEEAEE<EEEEEEEEAE/EEEEAEAEEAEE/A/AEEEE<A<E/EEEA<AEE<EEAA<E/A<EA<EA/<<<<<<<<AA/E/EAAA<<<AAEEA
+@ERR5069949.2668880 NS500628:121:HK3MMAFX2:1:21104:5191:15623/2
+CTGTGGATCACGGACAGCATCAGTAGTGTCAGCAATGTCTCTGCCAAATTGTTGGAAAGGCAGAAACTTTTTGTTAGACTCAGTAAGAACACCTGTGCCTGTTAAACCATTGAAGTTGAAATTGACACATTTGTTTTTAACCAAATTAGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEAEEAEEEEEEEEEEEEAEEEEEAEEEEEEEAEEEAAE<EEEEEE6A<AEEEEEEAEAEAE/EA
+@ERR5069949.257821 NS500628:121:HK3MMAFX2:2:21307:14790:19018/2
+GCCTATACAGTTGAACTCGGTACAGAAGTAAATGAGTTCGCCTGTGTTGTGGCAGATGCTGTCATAAAAACTTTGCAACCAGTATCTGAATTACTTACACCACTGGGCATTGATTTAGATGAGTGGAGTATGGCTACATA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEAEEEEEE
+@ERR5069949.2243023 NS500628:121:HK3MMAFX2:2:21301:7521:3094/2
+TGGTAAGAGAATTCCTTACACGTAACCCTGCTTGGAGAAAAGCTGTCTTTATTTCACCTTATAATTCACAGAATGCTGTAGCCTCAAAGATTTTGGGACTACCAACTCAAACTGTTGATTCATCACAGGGCTCAGAATATGACTATGTCAT
++
+AAAAAEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEEEEE/EEEEAAEEAEEA<EEEEEEEEEEEEEEAEEEEE<<AA6AAEEEAEE
+@ERR5069949.2605155 NS500628:121:HK3MMAFX2:1:21311:18812:15592/2
+CCAAAAATCCAGCCTCTTATTATGTTAGACTTCTCAGTGGAAGCAAAATAAACACCATCATTAAATGGTAGGACAGGGTTATCAAACCTCTTAGTACCATTGGTCCCAGAGACATGTATAGCATGGAACCAAGTAACATTGGAAAAGAAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEAEEEEEAEEAEEEEEEAEAEAEEE/EEAEEEEEEAEEEAAAEEEAEEAEEEEEAEAEE/EEEE/A
+@ERR5069949.1340552 NS500628:121:HK3MMAFX2:2:11309:7320:9079/2
+GGTAATGCAACAGAAGTGCCTGCCAATTCAACTGTATTATCTTTCTGTGCTTTTGCTGTAGATGCTGCTAAAGCTTACAAAGATTATCTAGCTAGTGGGGGACAACCAATCACTAATTGTGTTAAGATGTTGTGTACACACACTGGTA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEEEEAEAEEEEEEEEEE<EEEEEEEEEEEEEEAAA<AEEEEEEEEEEEEEE
+@ERR5069949.1561137 NS500628:121:HK3MMAFX2:2:11303:23159:8428/2
+GTTAAAAACAGTACAATTCTGTGATGCCATGCGAAATGCTGGTATTGTTGGTGTACTGACATTAGATAATCAAGATCTCAATGGTAACTGGTATGATTTCGGTGATTTCATACAAACCACGCCAGGTAGTGGAGTTCCTGTTGTAGATT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEE<EAEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEE<EEEEAEEEE/EEEEAAAEEEA<AEEEAEEEEEAEAEEA/AA<A
+@ERR5069949.2361683 NS500628:121:HK3MMAFX2:2:21103:8659:7427/2
+ATAGACAACAGGTGCGCTCAGGTCCTATTTTCACAAAATACTTCATAGATGTCAACTCAAAGCCATGTGCCCATAAGACAAATACGACTCTGTCAGAGAGATTTTTAAGTGTGTCACTTAACATTTGTACAATCTTTATACGCACTAC
++
+AAAAAAEEAEE/EEAEEEEEE/AEEEEEEAE<EEEEEEEEEEE/E/<EEEE<EEE/EEAA/EA/E<E/<<6AE<EE/EEE6/EAE/EEEEA/EEEEEEAE/</E/AA/AAE6EEAE/A<E//AAAAEEEE<EEA</AAEEA/<AE6<A
+@ERR5069949.2521353 NS500628:121:HK3MMAFX2:3:11508:26680:5334/2
+GTCCACACCCAAATTATTGAGTATTTTCACCTCTGGTACTGGTTTAATGTTGCGCTTAGCCCAAAGCACAAAAGCTACATTAACAGGTAATGTTGTTTTATTTTCAAACAATTCTACATCAACACCATCAACTTTTGTGAAAACAGTGAT
++
+AAAAAEEEEEE/EEEE6EEEEEAE/EE<EEE<EEA/EEE/AAE/E//EEEE/A/AEEE/EAEEA/<//EEA//AEE/<EAEE6AA/<//AEE/AE<<//A/<AAA/EEEEEAE6<///A<</A<E/EEE<E//AEEA<A/EE/</<A//A
+@ERR5069949.1261808 NS500628:121:HK3MMAFX2:3:11601:6846:14627/2
+TTGTCATCAGTGCAAGCAGTTTGTGTAGTACCGGCAGCACAAGACATCTGTCGTAGTGCAACAGGACTAAGCTCATTATTCTGTAATTTGACAGCAGAATTGGCC
++
+AA//AEEEA6EAE/EAEE/EAEEEEEEEA//E/EE/EEEEEE/EAEE/A/AEEEEAEA//EEEE/E<EAE/E/EE////6E</EAEA<A/E<EE<</EEE/<//<
+@ERR5069949.2734873 NS500628:121:HK3MMAFX2:4:21611:11682:9300/2
+CATTGGTGCAGGTATATGCGCTAGTTATCAGACTCAGACTAATTCTCCTCGGCGGGCACGTAGTGTAGCTAGTCAATCCATCATTGCCTACACTATGT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@ERR5069949.3017828 NS500628:121:HK3MMAFX2:1:21111:18284:17510/2
+ATTAACTATTAACGTACCTGTCTCTTCCGAAACGAATGAGTACATAAGTTCGTACTCATCAGCTTGTGCTTACAAAGGCACGCTAGTAGTCGTCGTCGGTTCATCA
++
+AAAAAEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEE<EEEEEEEEEEEE<E/EE<EAEEEEEE<EAEEEEEEEEE/<E<EAEEAEEEAAE/A
+@ERR5069949.573706 NS500628:121:HK3MMAFX2:3:21407:23793:9689/2
+TTTCTGTACAAACAACACCATCCAATTTATAAGTAACTGGTTTTATGGTTGTTGTGTAACTGTTTTCTTTGTAGAAAACATCCGTAATAGGACCTTTGTATTCTGAGGACTTTGTAAGTAAAGCACCGTCTATGCAATACAAAGTTTCT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEAEEEEEEAEAEAAAAEAEEEEEEA<EEE<<AAEEAEEEA<AEEEEEEEE<E<AA
+@ERR5069949.1980512 NS500628:121:HK3MMAFX2:1:21303:8321:13259/2
+GTAAAGTACAAATAGGAGAGTACACCTTTGAAAAAGGTGACTATGGTGATGCTGTTGTTTACCGAGGTACAACAACTTACAAATTAAATGTTGGTGATTATTTTGTGCTGACATCACATACAGTAATGCCATTAAGTGCACCTACACTAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEAEEEEEEEEEA<EEEEEEAEEEEAAAEEAAE<EEEAAAA<AA<EEEE/AE
+@ERR5069949.1014693 NS500628:121:HK3MMAFX2:1:11109:13494:20222/2
+ATCATACACAGTTCTTGCTGTCATAAGGATTAGTAACACTACAGCTGATGCATACATAACACAGTCTTTTAGCTTAAAACCAGACAAACTAGTATCAACCATATCCAACCATGTCATAATACGCATCACCCAACTAGCAGGCATATAGCC
++
+AAAAAEEEEAE/EEAEEEEEEEE/EEEEEEEEE/EEEEEEEEEEAAEEEAEAEAAEEEEEEAEEE<EEEEEAA<6EEEEEEEA<EEEAEAEA/EE/E/EEEEE/6AE/AAEAA<EEEAE///<AE<AA<A<6<EAEA<AA<<AAEEE</A
+@ERR5069949.3184655 NS500628:121:HK3MMAFX2:3:21606:7576:13093/2
+TTTATCTAAGTCACTAACTGAGAATAAATATTCTCAATTAGATGAAGAGCAACCAATGGAGATTGATTAAACGAACATGAAAATTATTCTTTTCTTGGCACTGATAACACTCGCTACTTGTGAGCTTTATCACTACCAAGAGTGTGTTAG
++
+AAAAAEEEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEE<EEEEE/EEEEEAEAEEEEE/EEEAEEE<EEEEEE<EEAAEEAEEEEEAAAEEE/E<AAEEAAAE6A/A<<A<AAAEE/AA6AE/A
+@ERR5069949.29668 NS500628:121:HK3MMAFX2:2:21202:17958:8710/2
+TATAATTAATAACTAATTACTGTCGTTGACAGGACACGAGTAACTCGTCTATCTTCTGCAGGCTGCTTACGGTTTCTTCCGTGTTGCAGCCGATCATCAGCACATCTAGGTTTTGTCCGGGTGTGACCGAAAGGTAAGATGGAGAGCCTT
++
+AAA/E/EEEEEEEEEAEEEEEEEEE/</E/E/EE<E/EEAAEA/E/EE//EA/EEEEEA/AEEE/EEEEE/E/EA/EE/EEE<E/E///E<AEE<<EEE/<EEEAA///AE/6A///A/AE/EAEE</EAEAE///AA/EEAEE/AAEAA
+@ERR5069949.3258358 NS500628:121:HK3MMAFX2:3:21508:15641:1056/2
+ATACTGCGTCAAGGTTCACCGCACTCACTCAACATGCCAAGGAAGACATAAAA
++
+AAAAAA/A/////EE/EEE/A</EEEEEE/EAEEE//AEEEAEEEEA///6E6
+@ERR5069949.1476386 NS500628:121:HK3MMAFX2:2:21303:26270:9400/2
+GTATACGACATCAGTACTAGTGCCTGTGCCGCACGGTGTAAGACGGGCTGCACTTACACCGCAAACCCGTTTAAAAACGATTGTGCATCAGCTGACTGAAGCATGGGTTCGCGGAGTTGATCACAACTACAGCCATAACCTTTCCACA
++
+AA6AAEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEAE//AAEEEEE/AEEEEEEEEEEEEEEEEEEEEE<<EEAAEEEEEEEEAEEE<EEE/EAE/EA6E/EEA/E/EEEAAA/AEA<AEEEAA
+@ERR5069949.2415814 NS500628:121:HK3MMAFX2:1:11107:11754:7444/2
+CAATAGTCCAGTCAACACGCTTAACAAAGCACTCGTGGACAGCTAGACACCTAGTCATGATTGCATCACAACTAGCTACATGTGCATTACCATGGACTTGACAATACAGATCATGGTTGCTTTGTAGGTTACCTGTAAAACCCCATTGTTG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAAEAEAEEEEEEEEEEEEAEEEE/EEEEE<AEEEA

tests/software/gatk4/splitncigarreads/main.nf

@@ -0,0 +1,18 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { GATK4_SPLITNCIGARREADS } from '../../../../software/gatk4/splitncigarreads/main.nf' addParams( options: [:] )
+
+workflow test_gatk4_splitncigarreads {
+
+    def input = []
+    input = [ [ id:'test' ], // meta map
+              [ file("${launchDir}/tests/data/bam/sarscov2_aln.bam", checkIfExists: true)] ]
+
+    fasta = [ file("tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.fna", checkIfExists: true),
+              file("tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.fna.fai", checkIfExists: true),
+              file("tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.dict", checkIfExists: true)]
+
+    GATK4_SPLITNCIGARREADS ( input, fasta )
+}

tests/software/gatk4/splitncigarreads/test.yml

@@ -0,0 +1,8 @@
+- name: gatk4 splitncigarreads
+  command: nextflow run ./tests/software/gatk4/splitncigarreads -entry test_gatk4_splitncigarreads -c tests/config/nextflow.config
+  tags:
+    - gatk4
+    - gatk4_splitncigarreads
+  files:
+    - path: output/gatk4/test.split_cigar.bam
+      md5sum: 9228ac79882a4929c797c891b1bd9f6d