Mercurial > repos > jjohnson > bamutil_diff
changeset 0:2cafa8420c04 draft default tip
"planemo upload for repository https://github.com/jj-umn/galaxytools/tree/master/bamutil/ commit c1945909ca200610f128577b68a82d9228905f3d-dirty"
| author | jjohnson |
|---|---|
| date | Fri, 26 Mar 2021 13:16:53 +0000 |
| parents | |
| children | |
| files | bamutil_diff.xml macros.xml test-data/diff.sam test-data/diff.txt test-data/diff_only1_in1.sam test-data/diff_only2_in2.sam test-data/in1.sam test-data/in2.sam test-data/in3.sam |
| diffstat | 9 files changed, 866 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /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="<\t1a3\t74M74N1M" /> + <has_text_matching expression=">\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>
--- /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>
--- /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:-
--- /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
--- /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
--- /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
--- /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:-
--- /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
--- /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:-