# HG changeset patch
# User jjohnson
# Date 1616764613 0
# Node ID 2cafa8420c046f99fb60cf9e723fc72e8637feaa

"planemo upload for repository https://github.com/jj-umn/galaxytools/tree/master/bamutil/ commit c1945909ca200610f128577b68a82d9228905f3d-dirty"

diff -r 000000000000 -r 2cafa8420c04 bamutil_diff.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/bamutil_diff.xml	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,265 @@
+<tool id="bamutil_diff" name="BamUtil diff" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" python_template_version="3.5">
+    <description>two coordinate sorted SAM/BAM files</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements" />
+    <command detect_errors="exit_code"><![CDATA[
+    bam diff 
+    --in1 '$in1'
+    --in2 '$in2'
+    #if $fields.choice == 'all':
+        --all
+    #elif $fields.choice == 'select':
+        $fields.flag
+        $fields.mapQual
+        $fields.mate
+        $fields.isize
+        $fields.seq
+        $fields.baseQual
+        $fields.noCigar
+        $fields.noPos
+        #if $fields.samtags.tagchoice == 'everyTag':
+            --everyTag
+        #elif $fields.samtags.tagchoice == 'specify':
+            --tags '$fields.samtags.tags'
+        #end if
+    #end if
+    --posDiff $posDiff
+    --recPoolSize -1
+    $onlyDiffs
+    --params
+    --noPhoneHome
+    --out $output_as
+    ]]></command>
+    <inputs>
+        <param argument="--in1" type="data" format="sam,bam" label="Input BAM 1"/>
+        <param argument="--in2" type="data" format="sam,bam" label="Input BAM 2"/>
+        <param argument="--posDiff" type="integer" value="100000" min="0" label="max base pair difference between possibly matching records"/>
+        <param argument="--onlyDiffs" type="boolean" truevalue="--onlyDiffs" falsevalue="" checked="false" label="only print the fields that differ"/>
+        <conditional name="fields">
+            <param name="choice" type="select" label="BAM fields to diff">
+                <option value="default" selected="true">Read Name, Flag Fragment bit, Position, Cigar</option>
+                <option value="all">Diff all the SAM/BAM fields</option>
+                <option value="select">Select SAM/BAM fields to diff</option>
+            </param>
+            <when value="default"/>
+            <when value="all"/>
+            <when value="select">
+                <param argument="--flag" type="boolean" truevalue="--flag" falsevalue="" checked="false" label="diff the flags."/>
+                <param argument="--mapQual" type="boolean" truevalue="--mapQual" falsevalue="" checked="false" label="diff the mapping qualities."/>
+                <param argument="--mate" type="boolean" truevalue="--mate" falsevalue="" checked="false" label="diff the mate chrom/pos."/>
+                <param argument="--isize" type="boolean" truevalue="--isize" falsevalue="" checked="false" label="diff the insert sizes."/>
+                <param argument="--seq" type="boolean" truevalue="--seq" falsevalue="" checked="false" label="diff the sequence bases."/>
+                <param argument="--baseQual" type="boolean" truevalue="--baseQual" falsevalue="" checked="false" label="diff the base qualities."/>
+                <param argument="--noCigar" type="boolean" truevalue="--noCigar" falsevalue="" checked="false" label="do not diff the the cigars."/>
+                <param argument="--noPos" type="boolean" truevalue="--noPos" falsevalue="" checked="false" label="do not diff the positions."/>
+                <conditional name="samtags">
+                    <param name="tagchoice" type="select" label="Tags to diff">
+                        <option value="none">Do not diff tags</option>
+                        <option value="everyTag">Diff every tag</option>
+                        <option value="specify">Specify tags to diff</option>
+                    </param>
+                    <when value="none"/>
+                    <when value="everyTag"/>
+                    <when value="specify">
+                        <param argument="--tags" type="text" label="diff the specified Tags formatted as Tag:Type,Tag:Type,Tag:Type...">
+                            <validator type="regex" message="SAM 2-char Tag:type">^([A-Za-z][A-Za-z0-9]:[AifZHB])(,[A-Za-z][A-Za-z0-9]:[AifZHB])*$</validator>
+                        </param>
+                    </when>
+                </conditional>
+            </when>
+        </conditional>
+        <param name="output_as" type="select" label="Output format">
+            <option value="diff.txt">ASCII text diff file</option>
+            <option value="diff.bam">BAM files: diff, only_in_file1, only_in_file2</option>
+            <option value="diff.sam">SAM files: diff, only_in_file1, only_in_file2</option>
+        </param>
+    </inputs>
+    <outputs>
+        <data name="diff_bam" format="bam" from_work_dir="diff.bam" label="${tool.name} on ${on_string}: diff.bam">
+            <filter>output_as == 'diff.bam'</filter>
+        </data>
+        <data name="diff_only1_bam" format="bam" from_work_dir="diff_only1_*.bam" label="${tool.name} on ${on_string} only in: ${in1.element_identifier}">
+            <filter>output_as == 'diff.bam'</filter>
+        </data>
+        <data name="diff_only2_bam" format="bam" from_work_dir="diff_only2_*.bam" label="${tool.name} on ${on_string} only in: ${in2.element_identifier}">
+            <filter>output_as == 'diff.bam'</filter>
+        </data>
+        <data name="diff_sam" format="sam" from_work_dir="diff.sam" label="${tool.name} on ${on_string}: diff.sam">
+            <filter>output_as == 'diff.sam'</filter>
+        </data>
+        <data name="diff_only1_sam" format="sam" from_work_dir="diff_only1_*.sam" label="${tool.name} on ${on_string} only in: ${in1.element_identifier}">
+            <filter>output_as == 'diff.sam'</filter>
+        </data>
+        <data name="diff_only2_sam" format="sam" from_work_dir="diff_only2_*.sam" label="${tool.name} on ${on_string} only in: ${in2.element_identifier}">
+            <filter>output_as == 'diff.sam'</filter>
+        </data>
+        <data name="diff_txt" format="txt" from_work_dir="diff.txt" label="${tool.name} on ${on_string}: diff.txt">
+            <filter>output_as == 'diff.txt'</filter>
+        </data>
+    </outputs>
+    <tests>
+        <!-- Test-1 --> 
+        <test>
+            <param name="in1" ftype="sam" value="in1.sam"/>
+            <param name="in2" ftype="sam" value="in2.sam"/>
+            <param name="posDiff" value="100000"/>
+            <param name="onlyDiffs" value="true"/>
+            <conditional name="fields">
+                <param name="choice" value="default"/>
+            </conditional>
+            <param name="output_as" value="diff.txt"/>
+            <output name="diff_txt" file="diff.txt"/>
+            <output name="diff_txt">
+                <assert_contents>
+                    <has_text text="NB500964:249:HHLFNBGX7:3:21407:1974:9687" />
+                    <has_text_matching expression="&lt;\t1a3\t74M74N1M" />
+                    <has_text_matching expression="&gt;\ta3\t74M66N1M" />
+                </assert_contents>
+            </output>
+        </test>
+
+        <!-- Test-2 --> 
+        <test>
+            <param name="in1" ftype="sam" value="in1.sam"/>
+            <param name="in2" ftype="sam" value="in2.sam"/>
+            <param name="posDiff" value="100000"/>
+            <param name="onlyDiffs" value="true"/>
+            <conditional name="fields">
+                <param name="choice" value="select"/>
+                <param name="flag" value="true"/>
+                <param name="seq" value="true"/>
+                <conditional name="samtags">
+                    <param name="tagchoice" value="specify"/>
+                    <param name="tags" value="AS:i,MD:Z"/>
+                </conditional>
+            </conditional>
+            <param name="output_as" value="diff.sam"/>
+            <output name="diff_sam">
+                <assert_contents>
+                    <has_text text="NB500964:249:HHLFNBGX7:4:12608:21020:10228" />
+                    <not_has_text text="NB500964:249:HHLFNBGX7:4:11510:10074:3541" />
+                    <not_has_text text="NB500964:249:HHLFNBGX7:1:12312:5087:3846" />
+                </assert_contents>
+            </output>
+            <output name="diff_only1_sam">
+                <assert_contents>
+                    <has_text text="NB500964:249:HHLFNBGX7:1:12312:5087:3846" />
+                    <not_has_text text="NB500964:249:HHLFNBGX7:4:11510:10074:3541" />
+                    <has_text text="TGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCG" />
+                    <has_text text="AS:i:0" />
+                    <has_text text="MD:Z:75" />
+                </assert_contents>
+            </output>
+            <output name="diff_only2_sam">
+                <assert_contents>
+                    <has_text text="NB500964:249:HHLFNBGX7:4:11510:10074:3541" />
+                    <not_has_text text="NB500964:249:HHLFNBGX7:1:12312:5087:3846" />
+                    <has_text text="ATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCA" />
+                    <has_text text="AS:i:0" />
+                    <has_text text="MD:Z:75" />
+                </assert_contents>
+            </output>
+        </test>
+        <!-- Test-3 --> 
+        <test>
+            <param name="in1" ftype="sam" value="in1.sam"/>
+            <param name="in2" ftype="sam" value="in3.sam"/>
+            <param name="posDiff" value="100000"/>
+            <param name="onlyDiffs" value="true"/>
+            <conditional name="fields">
+                <param name="choice" value="default"/>
+            </conditional>
+            <param name="output_as" value="diff.txt"/>
+            <output name="diff_txt">
+                <assert_contents>
+                    <not_has_text text="NB500964:249:HHLFNBGX7:3:21407:1974:9687" />
+                </assert_contents>
+            </output>
+        </test>
+    </tests>
+    <help><![CDATA[
+**bamUtil diff**
+
+The diff option on the bamUtil executable prints the difference between two coordinate sorted SAM/BAM files. This can be used to compare the outputs of running a SAM/BAM through different tools/versions of tools.
+The diff tool compares records that have the same Read Name and Fragment (from the flag). If a matching ReadName & Fragment is not found, the record is considered to be different.
+diff assumes the files are coordinate sorted and uses this assumption for determining how long to store a record before determining that the other file does not contain a matching ReadName/Fragment. If the files are not coordinate sorted, this logic does not work.
+By default, just the chromosome/position and cigar are compared for each record.
+Note: The headers are not compared.
+
+Options are available to compare::
+
+ - all fields
+ - flags
+ - mapping quality
+ - mate chromosome/position
+ - insert size
+ - sequence
+ - base quality
+ - specified tags
+ - all tags
+ - turn off position comparison
+ - turn off cigar comparison
+
+**Inputs**
+  Two BAM or SAM alignment files
+
+**Outputs**
+  Choice of 2 Output Formats:
+
+  ::
+
+    **Diff Format**
+    There are 2 types of differences.
+    ReadName/Fragment combo is in one file, but not in the other file within the window set by recPoolSize & posDiff
+    ReadName/Fragment combo is in both files, but at least one of the specified fields to diff is different
+    Each difference output consists of 2 or 3 lines. If the record only appears in one of the files, the diff is 2 lines, if it appears in both files, the diff is 3 lines.
+    The first line of the difference output is just the read name.
+    The 2nd and 3rd line (if present) begin with either a '<' or a '>'. If the record is from the first file (--in1), it begins with a '<'. If the record is from the 2nd file (--in2), it begins with a '>'.
+    The 2nd line is the flag followed by the diff'd fields from one of the records.
+    The 3rd line (if a matching record was found) is the flag followed by the diff'd fields from the matching record.
+    The diff'd record lines are tab separated, and are in the following order if --onlyDiffs is not specified::
+
+      - '<' or '>'
+      - flag
+      - chrom:pos (chromosome name ':' 1 based position) - if --noPos is not specified
+      - cigar - if --noCigar is not specified
+      - mapping quality - if --mapq or --all is specified
+      - mate chrom:pos (chromosome name ':' 1 based position) - if --mate or --all is specified
+      - insert size - if --isize or --all is specified
+      - sequence - if --seq or --all is specified
+      - base quality - if --baseQual or --all is specified
+      - tag:type:value - for each tag:type specified in --tags or for every tag if --all or --everyTag specified
+
+
+    **BAM Format**
+    In SAM/BAM format there will be 3 output files::
+
+      1. the specified name with record diffs
+      2. specified name with _only_<in1>.sam/bam with records only in the in1 file
+      3. specified name with _only_<in2>.sam/bam with records only in the in2 file
+
+    Records that are identical in the two files are not written in any of these output files.
+    When a record is found in both input files, but a difference is found, the record from the first file is written with additional tags to indicate the values from the second file, using the following tags::
+
+      - ZF - Flag
+      - ZP - Chromosome:1-based Position
+      - ZC - Cigar
+      - ZM - Mapping Quality
+      - ZN - Chromosome:1-based Mate Position
+      - ZI - Insert Size
+      - ZS - Sequence
+      - ZQ - Base Quality
+      - ZT - Tags
+
+    If --onlyDiffs is not specified, all fields that were compared will be printed in the tags. If --onlyDiffs is specified, then only the differing compared fields will be printed in the tags.
+
+
+
+
+https://genome.sph.umich.edu/wiki/BamUtil:_diff
+
+    ]]></help>
+    <expand macro="citations" />
+</tool>
diff -r 000000000000 -r 2cafa8420c04 macros.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,22 @@
+<macros>
+    <token name="@TOOL_VERSION@">1.0.15</token>
+    <token name="@VERSION_SUFFIX@">0</token>
+    <xml name="requirements">
+        <requirements>
+        <requirement type="package" version="@TOOL_VERSION@">bamutil</requirement>
+            <yield/>
+        </requirements>
+    </xml>
+    <xml name="citations">
+        <citations>
+            <citation type="bibtex">@online{bamUtil,
+              author = {iStatistical Genetics},
+              title = {bamUtil},
+              year = 2015,
+              url = {https://github.com/statgen/bamUtil},
+              urldate = {2021-03-01}
+            }</citation>
+            <yield />
+        </citations>
+    </xml>
+</macros>
diff -r 000000000000 -r 2cafa8420c04 test-data/diff.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/diff.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,59 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075444.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZF:i:163
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZF:i:163
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZF:i:163
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	ZF:i:163	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:74M66N1M
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	ZF:i:163	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:74M74N1M
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	ZF:i:163	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZC:Z:74M66N1M
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	ZF:i:163	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZC:Z:74M66N1M
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	ZF:i:163	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:74M74N1M
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	ZF:i:163	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:74M66N1M
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZF:i:83	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:68M74N2M5S
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	ZF:i:83	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:68M74N2M5S
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	ZF:i:83	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZC:Z:68M74N1M6S
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	ZF:i:83	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:68M66N2M5S
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZF:i:83	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZC:Z:68M74N2M5S
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	ZF:i:163	AS:i:-5	ZT:Z:AS:i:-5;MD:Z:70A4	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP	ZC:Z:75M
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZF:i:83
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZF:i:83
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZF:i:83	AS:i:-10	ZT:Z:AS:i:-10;MD:Z:7A62	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZC:Z:70M5S
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZF:i:83
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZF:i:83
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZF:i:83
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	ZF:i:99	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
diff -r 000000000000 -r 2cafa8420c04 test-data/diff.txt
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/diff.txt	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,109 @@
+NB500964:249:HHLFNBGX7:3:21407:1974:9687
+<	1a3	74M74N1M
+>	a3	74M66N1M
+NB500964:249:HHLFNBGX7:3:22501:2652:1890
+<	1a3	74M66N1M
+>	a3	74M74N1M
+NB500964:249:HHLFNBGX7:3:22503:25677:1737
+<	1a3	74M74N1M
+>	a3	74M66N1M
+NB500964:249:HHLFNBGX7:3:22611:13268:8039
+<	1a3	74M74N1M
+>	a3	74M66N1M
+NB500964:249:HHLFNBGX7:4:13607:25752:7238
+<	1a3	74M66N1M
+>	a3	74M74N1M
+NB500964:249:HHLFNBGX7:4:21503:16243:7076
+<	1a3	74M74N1M
+>	a3	74M66N1M
+NB500964:249:HHLFNBGX7:3:13604:7137:3437
+<	153	68M66N2M5S
+>	53	68M74N2M5S
+NB500964:249:HHLFNBGX7:3:23501:15949:4416
+<	153	68M66N2M5S
+>	53	68M74N2M5S
+NB500964:249:HHLFNBGX7:4:12404:22663:7685
+<	153	68M66N1M6S
+>	53	68M74N1M6S
+NB500964:249:HHLFNBGX7:4:12608:21020:10228
+<	153	68M74N2M5S
+>	53	68M66N2M5S
+NB500964:249:HHLFNBGX7:4:23603:22568:16607
+<	153	68M66N2M5S
+>	53	68M74N2M5S
+NB500964:249:HHLFNBGX7:3:23508:17777:5737
+<	1a3	70M5S
+>	a3	75M
+NB500964:249:HHLFNBGX7:3:23508:17777:5737
+<	153	75M
+>	53	70M5S
+NB500964:249:HHLFNBGX7:4:11510:10074:3541
+>	a3	6:52995646	75M
+NB500964:249:HHLFNBGX7:2:13209:7691:13203
+>	a3	6:52995675	75M
+NB500964:249:HHLFNBGX7:4:11510:10074:3541
+>	53	6:52995715	75M
+NB500964:249:HHLFNBGX7:3:21512:8039:1524
+>	a3	6:52995730	75M
+NB500964:249:HHLFNBGX7:2:13209:7691:13203
+>	53	6:52995789	75M
+NB500964:249:HHLFNBGX7:4:21507:2226:9226
+>	a3	6:52995797	75M
+NB500964:249:HHLFNBGX7:4:21602:10644:5676
+>	a3	6:52995802	75M
+NB500964:249:HHLFNBGX7:4:21507:2226:9226
+>	53	6:52995834	75M
+NB500964:249:HHLFNBGX7:4:21602:10644:5676
+>	53	6:52995834	75M
+NB500964:249:HHLFNBGX7:3:21512:8039:1524
+>	53	6:52995836	75M
+NB500964:249:HHLFNBGX7:3:13604:19654:15542
+>	63	6:53004731	75M
+NB500964:249:HHLFNBGX7:3:13604:19654:15542
+>	93	6:53004936	75M
+NB500964:249:HHLFNBGX7:1:12312:5087:3846
+<	a3	6:52995649	75M
+NB500964:249:HHLFNBGX7:2:23304:14989:7965
+<	81	6:52995665	75M
+NB500964:249:HHLFNBGX7:4:13611:22454:15447
+<	1a3	6:52995667	75M
+NB500964:249:HHLFNBGX7:3:12409:4568:1161
+<	a3	6:52995706	75M
+NB500964:249:HHLFNBGX7:4:21406:23062:17974
+<	a3	6:52995706	75M
+NB500964:249:HHLFNBGX7:4:13611:22454:15447
+<	153	6:52995712	75M
+NB500964:249:HHLFNBGX7:3:23507:21978:1662
+<	a3	6:52995758	75M
+NB500964:249:HHLFNBGX7:2:11209:14749:12659
+<	a3	6:52995766	75M
+NB500964:249:HHLFNBGX7:3:12409:4568:1161
+<	53	6:52995785	75M
+NB500964:249:HHLFNBGX7:4:21406:23062:17974
+<	53	6:52995825	75M
+NB500964:249:HHLFNBGX7:2:13103:3678:13730
+<	59	6:52995829	75M
+NB500964:249:HHLFNBGX7:2:11209:14749:12659
+<	53	6:52995834	75M
+NB500964:249:HHLFNBGX7:3:23507:21978:1662
+<	53	6:52995835	1S73M1S
+NB500964:249:HHLFNBGX7:4:12611:5420:13168
+<	51	6:52995838	74M1S
+NB500964:249:HHLFNBGX7:1:12312:5087:3846
+<	53	6:52995848	73M2S
+NB500964:249:HHLFNBGX7:3:12402:8344:2452
+>	63	6:73519463	54M366N21M
+NB500964:249:HHLFNBGX7:3:12402:8344:2452
+>	93	6:73519916	75M
+NB500964:249:HHLFNBGX7:3:21607:1536:15142
+<	63	6:73518202	63M89N12M
+NB500964:249:HHLFNBGX7:3:21607:1536:15142
+<	93	6:73518250	2S15M89N58M
+NB500964:249:HHLFNBGX7:2:21106:19199:5423
+<	63	6:73519036	75M
+NB500964:249:HHLFNBGX7:2:21106:19199:5423
+<	93	6:73519096	75M
+NB500964:249:HHLFNBGX7:1:23112:23126:11688
+<	163	6:73519370	75M
+NB500964:249:HHLFNBGX7:1:23112:23126:11688
+<	193	6:73519434	75M
diff -r 000000000000 -r 2cafa8420c04 test-data/diff_only1_in1.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/diff_only1_in1.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,58 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075444.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	163	6	52995649	60	75M	=	52995848	274	TGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCG	AAAA/EEEE6EEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEE6EEEEEEEE/EEEEEEEEEEEEEEEEEE/EE/	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:-1	XO:i:0	XG:i:0	YS:i:-2	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:23304:14989:7965	129	6	52995665	1	75M	2	32916255	0	GCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCT	AAAAAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEE<AE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YT:Z:UP
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	419	6	52995667	1	75M	=	52995712	120	CAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCC	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	163	6	52995706	60	75M	=	52995785	154	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AA/AAEEAEEAAAAEEAEEEEEEEEEEE/EEAEEEEEEEAEEEE<EEAEEEEEEEAA<EEAE<AA/AEAEAEEAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	163	6	52995706	60	75M	=	52995825	194	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AAAAAEEEEEAEEEEAEEEAEEAAEEEEE/AEEEEE<EAEAEEEEEEEEEA/EEEEEEAEEEEAE<EE6E<E6EE	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	339	6	52995712	1	75M	=	52995667	-120	CCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	163	6	52995758	60	75M	=	52995835	151	GTCGAAGAGGACGACCATCCCCGATAGAGGAGGACCTGTCTTCGGTCAAGGGTATACTAGTAGCTGCGCTCCCCT	AA/AAEEE6EAEEEAEEEEEEAEEE6EEE/EEA/EA/E/EEEEEEEEEEAEEEEEEE/E/66/<E<A<A6A////	AS:i:-6	MD:Z:36G20G17	NM:i:2	NH:i:1	XN:i:0	XM:i:2	ZS:i:-6	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	163	6	52995766	60	75M	=	52995834	143	GGACGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAAC	AAAAA/EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEAE<EEAEEEEEEE<6E	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	83	6	52995785	60	75M	=	52995706	-154	AGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAG	</EA<AEEE/EAAEEAAEAE/EEEEEEEEEAEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	83	6	52995825	60	75M	=	52995706	-194	GCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAG	EAEEA<<A/EEEEE/EEEEEEEEEEEEEE/EEAE/EEEEEE/EAEE<AEAEEEEEEEEEEEAAA/EEEEEA6AAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:13103:3678:13730	89	6	52995829	60	75M	=	52995829	0	CCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTC	EEEEEEEAAEEEEAAEEEEAEAEEEEEEEEEE/EEEEEEEEEEAEEEAEEEEAEEEEEEEE6EAEEEAAEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YT:Z:UP
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	83	6	52995834	60	75M	=	52995766	-143	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	83	6	52995835	60	1S73M1S	=	52995758	-151	GTAGAACCTCCAAACAATCTCTCAAGGTCCATTTGTAGGAGAACGTAGTGTAGTCAAGCTTCCAAGACTCCAGAA	//E/EE/A/6/E/E/EE//E/E/E/<6EE/EAEE//E<A/EE///EEA/EEEEEEEAE6EEEEAEAE6EEAAA6A	AS:i:-8	MD:Z:16G30G25	NM:i:2	NH:i:1	XN:i:0	XM:i:2	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12611:5420:13168	81	6	52995838	60	74M1S	2	32916257	0	AACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACAC	EEEEAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-1	MD:Z:74	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YT:Z:UP
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	83	6	52995848	60	73M2S	=	52995649	-274	CAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAAG	<EEEEEEEEEEEEAEEEEEEEEAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-2	MD:Z:73	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	99	6	73518202	60	63M89N12M	=	73518250	118	GTGGCAATCCAATACAGGGGCATAGCCGGCGCTTATTTGGCCTGGATGGTTCAGGATAATCACCTGAGCAGTGAA	AAAAAEEEEEEEEEEEAEEA/EEEEAEAEAEEEEEEEAEEE<EAA/EEEEE<EEAE/EA</A/AAAAAEAA/E/<	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:-2	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	147	6	73518250	60	2S15M89N58M	=	73518202	-118	TTGTTCAGGATAATCACCTGAGCAGTGAAGCCAGCTGCTTCCATTGGTGGGTCATTTTTGCTGTCACCAGCAACG	A/<EAEEA//<//EA/AEAEE/EEA//EE<EEE/EEA/AEEEEEAAAEE/E/EEEEEEEEEEE/EEEEEEAAAAA	AS:i:-2	MD:Z:73	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:-2	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	99	6	73519036	1	75M	=	73519096	135	CCTTAACAATTTCCTCATATCTCTTCTGGCTGTAGGGTGGCTCAGTGGAATCCATTTTGTTAACACCGACAATTA	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	147	6	73519096	1	75M	=	73519036	-135	TAACACCGACAATTAGTTGTTTCACACCCAGTGTGTAAGCCAGAAGGGCATGCTCTCGGGTCTGCCCATTCTTGG	EEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	355	6	73519370	1	75M	=	73519434	139	GTCTCTGTGTCCTGGGGCATCAATGATAGTCACATAGTACTTGCTGGTCTCAAATTTCCACAAGGAGATATCAAT	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	403	6	73519434	1	75M	=	73519370	-139	GAGATATCAATGGTGATACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
diff -r 000000000000 -r 2cafa8420c04 test-data/diff_only2_in2.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/diff_only2_in2.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,51 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075445.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075449.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075449.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075453.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075522.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075453.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075522.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:4:11510:10074:3541	163	6	52995646	60	75M	=	52995715	144	ATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCA	AAAAAEE/EEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEE/	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:13209:7691:13203	163	6	52995675	60	75M	=	52995789	189	ATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTG	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:11510:10074:3541	83	6	52995715	60	75M	=	52995646	-144	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAGGAG	EEEAEEE/EAEAEEEAEEE<EAEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21512:8039:1524	163	6	52995730	60	75M	=	52995836	181	GTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTC	AAAAAEEAEEEEEEEEEEEEEAE/EEE/EEEEEEEAEEEEE/EEEEEEEEEE/EAEEEEAEEEEEAEAEAA/EEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:2:13209:7691:13203	83	6	52995789	60	75M	=	52995675	-189	GGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCC	EAAAA<<AEEAEEEEEAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21507:2226:9226	163	6	52995797	60	75M	=	52995834	112	CTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAG	AAAAAEEEEEEEEEEE/EEAEEEEEEEE/EEEEE/EEEEEAA//EEEE/E/EE/EEEE/E<E/AEE/6EEEEE/E	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21602:10644:5676	163	6	52995802	60	75M	=	52995834	107	GTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAA	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21507:2226:9226	83	6	52995834	60	75M	=	52995797	-112	CTAGAACCTCCAAACAACCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	AEEEEE<EEEEEEEAEA/<A//EEEEEEEEEEEAEE/EEAAE/EAE/EE</AE/EEEEAEEEEEEEEEEEAAAAA	AS:i:-3	MD:Z:17G57	NM:i:1	NH:i:1	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21602:10644:5676	83	6	52995834	60	75M	=	52995802	-107	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6AAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21512:8039:1524	83	6	52995836	60	75M	=	52995730	-181	AGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACAC	/EA/EEEEAE/E/A</E/E/6/EEEEEEEE/EEEEE/EEEEEAEEEEEEEE6EEEEAE6EEEEEEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:19654:15542	99	6	53004731	60	75M	=	53004936	280	CTTAAATTGCACCCCTGTTGCAATACTGCATTAATATTATTACATTTTTCAAAATCTGTAGAAAATATTATCTGG	AAAAAEEEEEEEA//EEEEEEEEEEEEEEAEEEEEEEEEEEEAEE/EEEEEEEEEEEEEEEEE/EEEEEEE/EEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:19654:15542	147	6	53004936	60	75M	=	53004731	-280	AATAATTTTATACAAAAAAGCCTACTGCATTCAAATCAATATTTTAAAAAAAAACTTTAAAAAAGAATATTGACA	EEEEEAEEEEEEEEEEEEEAEE/EAEEEEEEEEEEEEEEEEEEEEEEEEEEAEAAEEAEEEEAAEEEEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12402:8344:2452	99	6	73519463	60	54M366N21M	=	73519916	162	CTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCGCATCTCAGCAGCCTCCTTCTCAAA	A6AAAEE6EEE6/EEEE/<EEEE/<EEEEEEEEEEEE//EEE/E<E/////<E/A/AA///AAEEE<//E</</A	AS:i:-3	MD:Z:50C24	NM:i:1	NH:i:1	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-9	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:3:12402:8344:2452	147	6	73519916	60	75M	=	73519463	-162	TCTTTTGGCGATGCCACCGCATTTATAGATCAGATGGCCAGTAGTGGTGGACTTGCCCGAAGCTACGGGTCCAAT	/EEAE6</EE//EEE//AEEEA/EE/<//AE/EE/E//AEAA/EAAAA/AE6EE///E//A/EAE6E/EE/AAAA	AS:i:-9	MD:Z:7T53T5T7	NM:i:3	NH:i:1	XN:i:0	XM:i:3	ZS:i:-9	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
diff -r 000000000000 -r 2cafa8420c04 test-data/in1.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/in1.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,103 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075444.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:142	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995648	ZS:Z:CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	ZQ:Z:AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	163	6	52995649	60	75M	=	52995848	274	TGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCG	AAAA/EEEE6EEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEE6EEEEEEEE/EEEEEEEEEEEEEEEEEE/EE/	AS:i:0	ZS:i:-1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:-2	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:23304:14989:7965	129	6	52995665	1	75M	2	32916255	0	GCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCT	AAAAAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEE<AE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YT:Z:UP	NH:i:2
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	419	6	52995667	1	75M	=	52995712	120	CAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCC	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	163	6	52995706	60	75M	=	52995785	154	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AA/AAEEAEEAAAAEEAEEEEEEEEEEE/EEAEEEEEEEAEEEE<EEAEEEEEEEAA<EEAE<AA/AEAEAEEAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	163	6	52995706	60	75M	=	52995825	194	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AAAAAEEEEEAEEEEAEEEAEEAAEEEEE/AEEEEE<EAEAEEEEEEEEEA/EEEEEEAEEEEAE<EE6E<E6EE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995832	ZI:i:132	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	ZT:Z:AS:i:-6;MD:Z:23G44C6;NM:i:2;NH:i:2;XN:i:0;XM:i:2;ZS:i:-6;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995788	ZI:i:80	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	ZQ:Z:AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995785	ZI:i:85	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995848	ZI:i:148	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	ZT:Z:AS:i:-3;MD:Z:68C6;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995789	ZI:i:81	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	ZQ:Z:AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	ZT:Z:AS:i:-3;MD:Z:71G3;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995838	ZI:i:138	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	339	6	52995712	1	75M	=	52995667	-120	CCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995637	ZI:i:-153	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995636	ZI:i:-154	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N1M6S	ZF:i:83	ZM:i:1	ZN:Z:6:52995648	ZI:i:-142	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	ZQ:Z:</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	ZT:Z:AS:i:-6;MD:Z:69;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M66N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	163	6	52995758	60	75M	=	52995835	151	GTCGAAGAGGACGACCATCCCCGATAGAGGAGGACCTGTCTTCGGTCAAGGGTATACTAGTAGCTGCGCTCCCCT	AA/AAEEE6EAEEEAEEEEEEAEEE6EEE/EEA/EA/E/EEEEEEEEEEAEEEEEEE/E/66/<E<A<A6A////	AS:i:-6	ZS:i:-6	XN:i:0	XM:i:2	XO:i:0	XG:i:0	NM:i:2	MD:Z:36G20G17	YS:i:-8	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	163	6	52995766	60	75M	=	52995834	143	GGACGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAAC	AAAAA/EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEAE<EEAEEEEEEE<6E	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:0	ZN:Z:6:52995832	ZI:i:138	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP	ZP:Z:6:52995769	ZS:Z:CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	ZT:Z:AS:i:-5;MD:Z:70A4;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-10;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M1S	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-85	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995785	ZS:Z:AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	ZQ:Z:EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	ZT:Z:AS:i:-5;MD:Z:11T62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	83	6	52995785	60	75M	=	52995706	-154	AGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAG	</EA<AEEE/EAAEEAAEAE/EEEEEEEEEAEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	83	6	52995825	60	75M	=	52995706	-194	GCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAG	EAEEA<<A/EEEEE/EEEEEEEEEEEEEE/EEAE/EEEEEE/EAEE<AEAEEEEEEEEEEEAAA/EEEEEA6AAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:13103:3678:13730	89	6	52995829	60	75M	=	52995829	0	CCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTC	EEEEEEEAAEEEEAAEEEEAEAEEEEEEEEEE/EEEEEEEEEEAEEEAEEEEAEEEEEEEE6EAEEEAAEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YT:Z:UP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-132	ZS:Z:TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995832	ZC:Z:75M	ZQ:Z:EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:70M5S	ZF:i:83	ZM:i:0	ZN:Z:6:52995769	ZI:i:-138	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995832	ZS:Z:TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	ZQ:Z:EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-10;MD:Z:7A62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-10;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995719	ZI:i:-116	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:0	ZN:Z:6:52995717	ZI:i:-118	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	ZT:Z:AS:i:-5;MD:Z:47G27;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	83	6	52995834	60	75M	=	52995766	-143	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	83	6	52995835	60	1S73M1S	=	52995758	-151	GTAGAACCTCCAAACAATCTCTCAAGGTCCATTTGTAGGAGAACGTAGTGTAGTCAAGCTTCCAAGACTCCAGAA	//E/EE/A/6/E/E/EE//E/E/E/<6EE/EAEE//E<A/EE///EEA/EEEEEEEAE6EEEEAEAE6EEAAA6A	AS:i:-8	XN:i:0	XM:i:2	XO:i:0	XG:i:0	NM:i:2	MD:Z:16G30G25	YS:i:-6	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:12611:5420:13168	81	6	52995838	60	74M1S	2	32916257	0	AACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACAC	EEEEAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:74	YT:Z:UP	NH:i:1
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-148	ZS:Z:CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995848	ZC:Z:75M	ZQ:Z:AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	ZT:Z:AS:i:-3;MD:Z:26G48;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	83	6	52995848	60	73M2S	=	52995649	-274	CAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAAG	<EEEEEEEEEEEEAEEEEEEEEAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:73	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	99	6	73518202	60	63M89N12M	=	73518250	118	GTGGCAATCCAATACAGGGGCATAGCCGGCGCTTATTTGGCCTGGATGGTTCAGGATAATCACCTGAGCAGTGAA	AAAAAEEEEEEEEEEEAEEA/EEEEAEAEAEEEEEEEAEEE<EAA/EEEEE<EEAE/EA</A/AAAAAEAA/E/<	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:-2	YT:Z:CP	XS:A:-	NH:i:1
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	147	6	73518250	60	2S15M89N58M	=	73518202	-118	TTGTTCAGGATAATCACCTGAGCAGTGAAGCCAGCTGCTTCCATTGGTGGGTCATTTTTGCTGTCACCAGCAACG	A/<EAEEA//<//EA/AEAEE/EEA//EE<EEE/EEA/AEEEEEAAAEE/E/EEEEEEEEEEE/EEEEEEAAAAA	AS:i:-2	ZS:i:-2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:73	YS:i:0	YT:Z:CP	XS:A:-	NH:i:1
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	99	6	73519036	1	75M	=	73519096	135	CCTTAACAATTTCCTCATATCTCTTCTGGCTGTAGGGTGGCTCAGTGGAATCCATTTTGTTAACACCGACAATTA	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:3
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	147	6	73519096	1	75M	=	73519036	-135	TAACACCGACAATTAGTTGTTTCACACCCAGTGTGTAAGCCAGAAGGGCATGCTCTCGGGTCTGCCCATTCTTGG	EEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:3
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	355	6	73519370	1	75M	=	73519434	139	GTCTCTGTGTCCTGGGGCATCAATGATAGTCACATAGTACTTGCTGGTCTCAAATTTCCACAAGGAGATATCAAT	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	403	6	73519434	1	75M	=	73519370	-139	GAGATATCAATGGTGATACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZP:Z:6:73519451	ZC:Z:66M366N9M	ZF:i:99	ZM:i:1	ZN:Z:6:73519983	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-	ZI:i:241	ZS:Z:ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:3;XN:i:0;XM:i:0;ZS:i:0;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP;XS:A:-
diff -r 000000000000 -r 2cafa8420c04 test-data/in2.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/in2.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,96 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075445.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075449.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075449.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075453.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075522.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075453.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075522.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	163	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	163	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:11510:10074:3541	163	6	52995646	60	75M	=	52995715	144	ATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCA	AAAAAEE/EEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEE/	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	163	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:142	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995648	ZS:Z:CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	ZQ:Z:AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:2:13209:7691:13203	163	6	52995675	60	75M	=	52995789	189	ATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTG	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEAEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	163	6	52995709	1	74M66N1M	=	52995832	132	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	163	6	52995709	1	74M74N1M	=	52995788	80	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	163	6	52995709	1	74M66N1M	=	52995785	85	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	163	6	52995709	1	74M66N1M	=	52995848	148	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	163	6	52995709	1	74M74N1M	=	52995789	81	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	163	6	52995709	1	74M66N1M	=	52995838	138	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995832	ZI:i:132	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	ZT:Z:AS:i:-6;MD:Z:23G44C6;NM:i:2;NH:i:2;XN:i:0;XM:i:2;ZS:i:-6;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995788	ZI:i:80	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	ZQ:Z:AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995785	ZI:i:85	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995848	ZI:i:148	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	ZT:Z:AS:i:-3;MD:Z:68C6;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995789	ZI:i:81	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	ZQ:Z:AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	ZT:Z:AS:i:-3;MD:Z:71G3;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995838	ZI:i:138	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	83	6	52995715	1	68M74N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	83	6	52995715	1	68M74N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	83	6	52995715	1	68M74N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	83	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	83	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995637	ZI:i:-153	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995636	ZI:i:-154	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N1M6S	ZF:i:83	ZM:i:1	ZN:Z:6:52995648	ZI:i:-142	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	ZQ:Z:</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	ZT:Z:AS:i:-6;MD:Z:69;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M66N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:11510:10074:3541	83	6	52995715	60	75M	=	52995646	-144	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAGGAG	EEEAEEE/EAEAEEEAEEE<EAEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21512:8039:1524	163	6	52995730	60	75M	=	52995836	181	GTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTC	AAAAAEEAEEEEEEEEEEEEEAE/EEE/EEEEEEEAEEEEE/EEEEEEEEEE/EAEEEEAEEEEEAEAEAA/EEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	163	6	52995769	0	75M	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70A4	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:0	ZN:Z:6:52995832	ZI:i:138	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP	ZP:Z:6:52995769	ZS:Z:CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	ZT:Z:AS:i:-5;MD:Z:70A4;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-10;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	83	6	52995785	1	74M1S	=	52995709	-85	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M1S	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-85	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995785	ZS:Z:AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	ZQ:Z:EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	ZT:Z:AS:i:-5;MD:Z:11T62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:2:13209:7691:13203	83	6	52995789	60	75M	=	52995675	-189	GGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCC	EAAAA<<AEEAEEEEEAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21507:2226:9226	163	6	52995797	60	75M	=	52995834	112	CTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAG	AAAAAEEEEEEEEEEE/EEAEEEEEEEE/EEEEE/EEEEEAA//EEEE/E/EE/EEEE/E<E/AEE/6EEEEE/E	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:-3	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21602:10644:5676	163	6	52995802	60	75M	=	52995834	107	GTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAA	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	83	6	52995832	1	75M	=	52995709	-132	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	83	6	52995832	0	70M5S	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-132	ZS:Z:TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995832	ZC:Z:75M	ZQ:Z:EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:70M5S	ZF:i:83	ZM:i:0	ZN:Z:6:52995769	ZI:i:-138	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995832	ZS:Z:TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	ZQ:Z:EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-10;MD:Z:7A62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-10;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	83	6	52995834	1	75M	=	52995719	-116	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	83	6	52995834	0	75M	=	52995717	-118	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995719	ZI:i:-116	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:0	ZN:Z:6:52995717	ZI:i:-118	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	ZT:Z:AS:i:-5;MD:Z:47G27;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21507:2226:9226	83	6	52995834	60	75M	=	52995797	-112	CTAGAACCTCCAAACAACCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	AEEEEE<EEEEEEEAEA/<A//EEEEEEEEEEEAEE/EEAAE/EAE/EE</AE/EEEEAEEEEEEEEEEEAAAAA	AS:i:-3	XN:i:0	XM:i:1	XO:i:0	XG:i:0	NM:i:1	MD:Z:17G57	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21602:10644:5676	83	6	52995834	60	75M	=	52995802	-107	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6AAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21512:8039:1524	83	6	52995836	60	75M	=	52995730	-181	AGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACAC	/EA/EEEEAE/E/A</E/E/6/EEEEEEEE/EEEEE/EEEEEAEEEEEEEE6EEEEAE6EEEEEEEEEEEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	83	6	52995848	1	75M	=	52995709	-148	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-148	ZS:Z:CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995848	ZC:Z:75M	ZQ:Z:AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	ZT:Z:AS:i:-3;MD:Z:26G48;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:19654:15542	99	6	53004731	60	75M	=	53004936	280	CTTAAATTGCACCCCTGTTGCAATACTGCATTAATATTATTACATTTTTCAAAATCTGTAGAAAATATTATCTGG	AAAAAEEEEEEEA//EEEEEEEEEEEEEEAEEEEEEEEEEEEAEE/EEEEEEEEEEEEEEEEE/EEEEEEE/EEE	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:13604:19654:15542	147	6	53004936	60	75M	=	53004731	-280	AATAATTTTATACAAAAAAGCCTACTGCATTCAAATCAATATTTTAAAAAAAAACTTTAAAAAAGAATATTGACA	EEEEEAEEEEEEEEEEEEEAEE/EAEEEEEEEEEEEEEEEEEEEEEEEEEEAEAAEEAEEEEAAEEEEEEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	99	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZP:Z:6:73519451	ZC:Z:66M366N9M	ZF:i:99	ZM:i:1	ZN:Z:6:73519983	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-	ZI:i:241	ZS:Z:ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:3;XN:i:0;XM:i:0;ZS:i:0;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP;XS:A:-
+NB500964:249:HHLFNBGX7:3:12402:8344:2452	99	6	73519463	60	54M366N21M	=	73519916	162	CTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCGCATCTCAGCAGCCTCCTTCTCAAA	A6AAAEE6EEE6/EEEE/<EEEE/<EEEEEEEEEEEE//EEE/E<E/////<E/A/AA///AAEEE<//E</</A	AS:i:-3	ZS:i:-3	XN:i:0	XM:i:1	XO:i:0	XG:i:0	NM:i:1	MD:Z:50C24	YS:i:-9	YT:Z:CP	XS:A:-	NH:i:1
+NB500964:249:HHLFNBGX7:3:12402:8344:2452	147	6	73519916	60	75M	=	73519463	-162	TCTTTTGGCGATGCCACCGCATTTATAGATCAGATGGCCAGTAGTGGTGGACTTGCCCGAAGCTACGGGTCCAAT	/EEAE6</EE//EEE//AEEEA/EE/<//AE/EE/E//AEAA/EAAAA/AE6EE///E//A/EAE6E/EE/AAAA	AS:i:-9	ZS:i:-9	XN:i:0	XM:i:3	XO:i:0	XG:i:0	NM:i:3	MD:Z:7T53T5T7	YS:i:-3	YT:Z:CP	NH:i:1
diff -r 000000000000 -r 2cafa8420c04 test-data/in3.sam
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/in3.sam	Fri Mar 26 13:16:53 2021 +0000
@@ -0,0 +1,103 @@
+@HD	VN:1.0	SO:coordinate
+@SQ	SN:1	LN:248956422
+@SQ	SN:2	LN:242193529
+@SQ	SN:3	LN:198295559
+@SQ	SN:4	LN:190214555
+@SQ	SN:5	LN:181538259
+@SQ	SN:6	LN:170805979
+@SQ	SN:7	LN:159345973
+@SQ	SN:8	LN:145138636
+@SQ	SN:9	LN:138394717
+@SQ	SN:10	LN:133797422
+@SQ	SN:11	LN:135086622
+@SQ	SN:12	LN:133275309
+@SQ	SN:13	LN:114364328
+@SQ	SN:14	LN:107043718
+@SQ	SN:15	LN:101991189
+@SQ	SN:16	LN:90338345
+@SQ	SN:17	LN:83257441
+@SQ	SN:18	LN:80373285
+@SQ	SN:19	LN:58617616
+@SQ	SN:20	LN:64444167
+@SQ	SN:21	LN:46709983
+@SQ	SN:22	LN:50818468
+@SQ	SN:X	LN:156040895
+@SQ	SN:Y	LN:57227415
+@SQ	SN:MT	LN:16569
+@PG	ID:hisat2	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:hisat2-55424A4-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools	PN:samtools	PP:hisat2-55424A4-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 0.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075444.dat
+@PG	ID:hisat2-3A2CCEF5	PN:hisat2	VN:2.1.0	CL:"/home/galaxy/galaxy/miniconda2/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/rissdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools-6ADB4A65	PN:samtools	PP:hisat2-3A2CCEF5	VN:1.11	CL:samtools sort -@ 1 -m 10240M -T sorttemp -O sam -o 1.sam /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat
+@PG	ID:samtools.1	PN:samtools	PP:samtools	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:samtools.2	PN:samtools	PP:samtools-6ADB4A65	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075446.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075450.dat 0.sam 1.sam
+@PG	ID:hisat2-6ADB4A65	PN:hisat2	VN:2.1.0	CL:"/panfs/roc/galaxy/shared/miniconda3/envs/mulled-v1-e7321ba46fa5ea4c6b9a06b78e6cd5182b33c0a47c2c86b5d610e1361f8b1686/bin/hisat2-align-s --wrapper basic-0 -p 4 -x /panfs/roc/risdb/galaxy/genomes/GRCh38_canon/hisat2_index/GRCh38_canon/GRCh38_canon --pen-cansplice 0 --pen-noncansplice 12 --pen-canintronlen G,-8.0,1.0 --pen-noncanintronlen G,-8.0,1.0 --known-splicesite-infile splice_sites.txt --min-intronlen 20 --max-intronlen 500000 --dta -1 input_f.fastq -2 input_r.fastq"
+@PG	ID:samtools.3	PN:samtools	PP:samtools.1	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+@PG	ID:samtools.4	PN:samtools	PP:samtools.2	VN:1.11	CL:samtools merge -@ 1 -s 1 -f -h /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075454.dat /panfs/roc/galaxy/PRODUCTION/database/files/001/075/dataset_1075523.dat 0.sam 1.sam
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEAEEEE<EEE//<AAEAAEEEE<AAE/EEEEAE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	419	6	52995635	1	75M	=	52995715	155	CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:155	ZS:i:-11	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995635	ZS:Z:CTGGCTGCGACATCTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCC	ZQ:Z:AAAAAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA/EAEEAAAEEEEE	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-11;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	419	6	52995648	1	75M	=	52995715	142	CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:75M	ZF:i:163	ZM:i:1	ZN:Z:6:52995715	ZI:i:142	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995648	ZS:Z:CTGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCATCACC	ZQ:Z:AAAAAEEE//EEEEEE/EAEEEEEEEEEE/AEE/EEEEEEEEEEEEAEAEEA/AEEAEEE<<EE/EEEE/EEAEE	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	163	6	52995649	60	75M	=	52995848	274	TGTCACCCCATTGATCGCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCG	AAAA/EEEE6EEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEE6EEEEEEEE/EEEEEEEEEEEEEEEEEE/EE/	AS:i:0	ZS:i:-1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:-2	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:23304:14989:7965	129	6	52995665	1	75M	2	32916255	0	GCCAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCT	AAAAAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEE<AE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YT:Z:UP	NH:i:2
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	419	6	52995667	1	75M	=	52995712	120	CAGGGTTGATTCGGCTGATCTGGCTGGCTAGGCGGGTGTCCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCC	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	163	6	52995706	60	75M	=	52995785	154	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AA/AAEEAEEAAAAEEAEEEEEEEEEEE/EEAEEEEEEEAEEEE<EEAEEEEEEEAA<EEAE<AA/AEAEAEEAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	163	6	52995706	60	75M	=	52995825	194	CCCCTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCG	AAAAAEEEEEAEEEEAEEEAEEAAEEEEE/AEEEEE<EAEAEEEEEEEEEA/EEEEEEAEEEEAE<EE6E<E6EE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	419	6	52995709	1	74M74N1M	=	52995832	124	CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	AS:i:-6	MD:Z:23G44C6	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995832	ZI:i:132	ZS:i:-6	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTCCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:A//AAAE/EE//EE/EEEAAEA//EEEEE/AEEEEEEEEEEEE/E/E6/A/E//A/E/EEEA/AAEAE//EEA/A	ZT:Z:AS:i:-6;MD:Z:23G44C6;NM:i:2;NH:i:2;XN:i:0;XM:i:2;ZS:i:-6;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22501:2652:1890	419	6	52995709	1	74M66N1M	=	52995788	88	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	AS:i:-3	MD:Z:69C5	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995788	ZI:i:80	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCACGATA	ZQ:Z:AAAAAE/EEEA/EAAEEEEEEEEEEEE/<EEEEAA<EEEE/6EEEEEEEAEEEAEAAE/EEEAEE<EAE/EEE<A	ZT:Z:AS:i:-3;MD:Z:69C5;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	419	6	52995709	1	74M74N1M	=	52995785	77	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995785	ZI:i:85	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/E/AAEAEE/EE/AAAEEE/EEEEEEEEAEAEEEEEAEEEEEAEE/EE//EEE/AEE/AE//E/E/EA/E/	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	419	6	52995709	1	74M74N1M	=	52995848	140	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	AS:i:-3	MD:Z:68C6	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995848	ZI:i:148	ZS:i:-3	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCACCGATA	ZQ:Z:AAAAAEEAEAEEEA/EEEEAAE/AEE/<AEEEE/EEEEAEEE<AAEEAEEEAEEEEEEEEE/AEAEAE/EEAEE/	ZT:Z:AS:i:-3;MD:Z:68C6;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13607:25752:7238	419	6	52995709	1	74M66N1M	=	52995789	89	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	AS:i:-3	MD:Z:71G3	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M74N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995789	ZI:i:81	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCCATA	ZQ:Z:AAAAA/EEE/EEEEEEEEEAE/A6AEEEEEEEEEEE/<EEEEAEAEEEE/EEE/EEE/EAEE/AAE/EEEE///A	ZT:Z:AS:i:-3;MD:Z:71G3;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:21503:16243:7076	419	6	52995709	1	74M74N1M	=	52995838	130	CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	AS:i:-3	MD:Z:67C7	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M66N1M	ZF:i:163	ZM:i:1	ZN:Z:6:52995838	ZI:i:138	ZS:i:-3	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995709	ZS:Z:CTTCCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATACCCGATA	ZQ:Z:AAAA/EEE/EE/AEAEA/EEAA/AAE/EEE/EAEEA/AEEEE<EEEE<<//E/AE/E6EE/</EAAE/<EE/EEE	ZT:Z:AS:i:-3;MD:Z:67C7;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-3;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:13611:22454:15447	339	6	52995712	1	75M	=	52995667	-120	CCTCCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:13604:7137:3437	339	6	52995715	1	68M66N2M5S	=	52995637	-153	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995637	ZI:i:-153	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:<EEEEEEEEEEEEEE<EEEE6EEEE<E/EEEEEAEEEEEEAEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23501:15949:4416	339	6	52995715	1	68M66N2M5S	=	52995636	-154	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995636	ZI:i:-154	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEAEEAEEEEEEEEEEEEEAEEEEEEEEAAEEEAEEEEAEEEEEEEEAEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12404:22663:7685	339	6	52995715	1	68M66N1M6S	=	52995648	-142	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	AS:i:-6	MD:Z:69	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N1M6S	ZF:i:83	ZM:i:1	ZN:Z:6:52995648	ZI:i:-142	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATACAGGAG	ZQ:Z:</EEA/EEA/EEA/EEEEAEEEEE/EEAEEEEEEEAE<//EEAEE/EEEEEEAE<AE/6EEE/EEEEEAEAAAAA	ZT:Z:AS:i:-6;MD:Z:69;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:12608:21020:10228	339	6	52995715	1	68M74N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M66N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:AEEEEEEA<AEEAEEEEE<<EEEEEEEAAEAEEEEEAEEEEEE/EEEEEAEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:4:23603:22568:16607	339	6	52995715	1	68M66N2M5S	=	52995635	-155	CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:68M74N2M5S	ZF:i:83	ZM:i:1	ZN:Z:6:52995635	ZI:i:-155	ZS:i:-5	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	ZP:Z:6:52995715	ZS:Z:CCCTCACCGCTCCATGTGCGTCCCTCCCGAAGCTGCGCGCTCGGTCGAAGAGGACGACCATCCCCGATAAAGGAG	ZQ:Z:EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-5;MD:Z:70;NM:i:0;NH:i:2;XN:i:0;XM:i:0;ZS:i:-5;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	163	6	52995758	60	75M	=	52995835	151	GTCGAAGAGGACGACCATCCCCGATAGAGGAGGACCTGTCTTCGGTCAAGGGTATACTAGTAGCTGCGCTCCCCT	AA/AAEEE6EAEEEAEEEEEEAEEE6EEE/EEA/EA/E/EEEEEEEEEEAEEEEEEE/E/66/<E<A<A6A////	AS:i:-6	ZS:i:-6	XN:i:0	XM:i:2	XO:i:0	XG:i:0	NM:i:2	MD:Z:36G20G17	YS:i:-8	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	163	6	52995766	60	75M	=	52995834	143	GGACGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAAC	AAAAA/EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEAEAE<EEAEEEEEEE<6E	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	419	6	52995769	0	70M5S	=	52995832	138	CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	AS:i:-5	MD:Z:70	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZC:Z:75M	ZF:i:163	ZM:i:0	ZN:Z:6:52995832	ZI:i:138	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-10	YT:Z:CP	ZP:Z:6:52995769	ZS:Z:CGACCATCCCCGATAGAGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGATCCTC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAE	ZT:Z:AS:i:-5;MD:Z:70A4;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-10;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22503:25677:1737	339	6	52995785	1	74M1S	=	52995709	-77	AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	AS:i:-5	MD:Z:11T62	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZC:Z:74M1S	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-85	ZS:i:-5	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995785	ZS:Z:AGGAGGACCGGGCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAC	ZQ:Z:EA<E//EEA</6EAEAEAEEAEE/EEE6EEEEEEEAEAEE/EEEEEEA/EEEAEEEEEEEEEE6EE/AEE/AAA/	ZT:Z:AS:i:-5;MD:Z:11T62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-5;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12409:4568:1161	83	6	52995785	60	75M	=	52995706	-154	AGGAGGACCGGTCTTCGGTCAAGGGTATACGAGTAGCTGCGCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAG	</EA<AEEE/EAAEEAAEAE/EEEEEEEEEAEEEEEEEEAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:21406:23062:17974	83	6	52995825	60	75M	=	52995706	-194	GCTCCCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAG	EAEEA<<A/EEEEE/EEEEEEEEEEEEEE/EEAE/EEEEEE/EAEE<AEAEEEEEEEEEEEAAA/EEEEEA6AAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:2:13103:3678:13730	89	6	52995829	60	75M	=	52995829	0	CCCTGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTC	EEEEEEEAAEEEEAAEEEEAEAEEEEEEEEEE/EEEEEEEEEEAEEEAEEEEAEEEEEEEE6EAEEEAAEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YT:Z:UP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21407:1974:9687	339	6	52995832	1	75M	=	52995709	-124	TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-132	ZS:Z:TGCTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAG	XO:i:0	XG:i:0	YS:i:-6	YT:Z:CP	ZP:Z:6:52995832	ZC:Z:75M	ZQ:Z:EEEEEEE/EEE<EEE/EEAEEEE/A</EEEE6EEAE/EEEAEEEEEEEEEEEEEA6EEEEEEE/EAEEEEAAAA<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-6;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:23508:17777:5737	339	6	52995832	0	75M	=	52995769	-138	TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-10	MD:Z:7A62T4	NM:i:2	NH:i:2	XN:i:0	XM:i:2	ZC:Z:70M5S	ZF:i:83	ZM:i:0	ZN:Z:6:52995769	ZI:i:-138	ZS:i:-10	XO:i:0	XG:i:0	YS:i:-5	YT:Z:CP	ZP:Z:6:52995832	ZS:Z:TGCTAGATCCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACCCCAG	ZQ:Z:EEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	ZT:Z:AS:i:-10;MD:Z:7A62;NM:i:1;NH:i:2;XN:i:0;XM:i:1;ZS:i:-10;XO:i:0;XG:i:0;YS:i:-5;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:12605:20423:18895	339	6	52995834	1	75M	=	52995719	-124	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	AS:i:0	MD:Z:75	NM:i:0	NH:i:2	XN:i:0	XM:i:0	ZF:i:83	ZM:i:1	ZN:Z:6:52995719	ZI:i:-116	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:/EEEEAEAEAEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEAEEEAAAAA	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:2;XN:i:0;XM:i:0;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:3:21607:9911:11221	339	6	52995834	0	75M	=	52995717	-126	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	AS:i:-5	MD:Z:47G27	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:0	ZN:Z:6:52995717	ZI:i:-118	ZS:Z:CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAAGGTAGTCAAGCTTCCAAGACTCCAGAC	XO:i:0	XG:i:0	YS:i:-8	YT:Z:CP	ZP:Z:6:52995834	ZC:Z:75M	ZQ:Z:EEEEEAEAEEEEEEEE</E/EEEE/E<A//EE6EEE/EEEEEEEEEEEEEEEE/A//EEEEEEEEEA/EEAAAAA	ZT:Z:AS:i:-5;MD:Z:47G27;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-8;YT:Z:CP
+NB500964:249:HHLFNBGX7:2:11209:14749:12659	83	6	52995834	60	75M	=	52995766	-143	CTAGAACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGAC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:23507:21978:1662	83	6	52995835	60	1S73M1S	=	52995758	-151	GTAGAACCTCCAAACAATCTCTCAAGGTCCATTTGTAGGAGAACGTAGTGTAGTCAAGCTTCCAAGACTCCAGAA	//E/EE/A/6/E/E/EE//E/E/E/<6EE/EAEE//E<A/EE///EEA/EEEEEEEAE6EEEEAEAE6EEAAA6A	AS:i:-8	XN:i:0	XM:i:2	XO:i:0	XG:i:0	NM:i:2	MD:Z:16G30G25	YS:i:-6	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:4:12611:5420:13168	81	6	52995838	60	74M1S	2	32916257	0	AACCTCCAAACAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACAC	EEEEAEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-1	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:74	YT:Z:UP	NH:i:1
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP
+NB500964:249:HHLFNBGX7:3:22611:13268:8039	339	6	52995848	1	75M	=	52995709	-140	CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	AS:i:-3	MD:Z:26G48	NM:i:1	NH:i:2	XN:i:0	XM:i:1	ZF:i:83	ZM:i:1	ZN:Z:6:52995709	ZI:i:-148	ZS:Z:CAAGCTCTCAAGGTCCATTTGTAGGATAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAGG	XO:i:0	XG:i:0	YS:i:-3	YT:Z:CP	ZP:Z:6:52995848	ZC:Z:75M	ZQ:Z:AE/EE/<EE//AEA<AEEEEEEEE6E/EEAEEEEAE6EE/EE///EE//EE<EEEAEAAEAAE/EEAEEEAAA6A	ZT:Z:AS:i:-3;MD:Z:26G48;NM:i:1;NH:i:2;XN:i:0;XM:i:1;XO:i:0;XG:i:0;YS:i:-3;YT:Z:CP
+NB500964:249:HHLFNBGX7:1:12312:5087:3846	83	6	52995848	60	73M2S	=	52995649	-274	CAAGCTCTCAAGGTCCATTTGTAGGAGAACGTAGGGTAGTCAAGCTTCCAAGACTCCAGACACATCCAAATGAAG	<EEEEEEEEEEEEAEEEEEEEEAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA	AS:i:-2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:73	YS:i:0	YT:Z:CP	NH:i:1
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	99	6	73518202	60	63M89N12M	=	73518250	118	GTGGCAATCCAATACAGGGGCATAGCCGGCGCTTATTTGGCCTGGATGGTTCAGGATAATCACCTGAGCAGTGAA	AAAAAEEEEEEEEEEEAEEA/EEEEAEAEAEEEEEEEAEEE<EAA/EEEEE<EEAE/EA</A/AAAAAEAA/E/<	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:-2	YT:Z:CP	XS:A:-	NH:i:1
+NB500964:249:HHLFNBGX7:3:21607:1536:15142	147	6	73518250	60	2S15M89N58M	=	73518202	-118	TTGTTCAGGATAATCACCTGAGCAGTGAAGCCAGCTGCTTCCATTGGTGGGTCATTTTTGCTGTCACCAGCAACG	A/<EAEEA//<//EA/AEAEE/EEA//EE<EEE/EEA/AEEEEEAAAEE/E/EEEEEEEEEEE/EEEEEEAAAAA	AS:i:-2	ZS:i:-2	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:73	YS:i:0	YT:Z:CP	XS:A:-	NH:i:1
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	99	6	73519036	1	75M	=	73519096	135	CCTTAACAATTTCCTCATATCTCTTCTGGCTGTAGGGTGGCTCAGTGGAATCCATTTTGTTAACACCGACAATTA	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:3
+NB500964:249:HHLFNBGX7:2:21106:19199:5423	147	6	73519096	1	75M	=	73519036	-135	TAACACCGACAATTAGTTGTTTCACACCCAGTGTGTAAGCCAGAAGGGCATGCTCTCGGGTCTGCCCATTCTTGG	EEEEEEEEEEEEE<EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:3
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	355	6	73519370	1	75M	=	73519434	139	GTCTCTGTGTCCTGGGGCATCAATGATAGTCACATAGTACTTGCTGGTCTCAAATTTCCACAAGGAGATATCAAT	AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEE	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:1:23112:23126:11688	403	6	73519434	1	75M	=	73519370	-139	GAGATATCAATGGTGATACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCC	EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEAEEEEEEEEEEEEAAAAA	AS:i:0	ZS:i:0	XN:i:0	XM:i:0	XO:i:0	XG:i:0	NM:i:0	MD:Z:75	YS:i:0	YT:Z:CP	NH:i:2
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-
+NB500964:249:HHLFNBGX7:4:13509:13430:12762	355	6	73519451	1	66M366N9M	=	73519983	241	ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	AS:i:0	MD:Z:75	NM:i:0	NH:i:3	XN:i:0	XM:i:0	ZP:Z:6:73519451	ZC:Z:66M366N9M	ZF:i:99	ZM:i:1	ZN:Z:6:73519983	ZS:i:0	XO:i:0	XG:i:0	YS:i:0	YT:Z:CP	XS:A:-	ZI:i:241	ZS:Z:ACCACGTTCACGCTCAGCTTTCAGTTTATCCAAGACCCAGGCATACTTGAAGGAGCCCTTTCCCATCTCAGCAGC	ZQ:Z:AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEE/EE/EEEEEAEEAEEEAAEEEEEEEE<E<	ZT:Z:AS:i:0;MD:Z:75;NM:i:0;NH:i:3;XN:i:0;XM:i:0;ZS:i:0;XO:i:0;XG:i:0;YS:i:0;YT:Z:CP;XS:A:-