Mercurial > repos > iuc > bamutil_diff
changeset 0:c28c4c132459 draft default tip
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bamutil commit 9f971aa7d32d756aa4763ab29d1770178a5ece6b"
| author | iuc |
|---|---|
| date | Fri, 02 Apr 2021 20:20:54 +0000 |
| parents | |
| children | |
| files | diff.xml macros.xml test-data/in1.sam test-data/in2.sam test-data/in3.sam test-data/input.bam test-data/input.sam test-data/input_qname_sorted.bam test-data/output.sam test-data/output_qname_sorted.bam test-data/output_stats.txt |
| diffstat | 11 files changed, 956 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/diff.xml Fri Apr 02 20:20:54 2021 +0000 @@ -0,0 +1,265 @@ +<tool id="bamutil_diff" name="BamUtil diff" version="@WRAPPER_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@"> + <description>two coordinate sorted SAM/BAM files</description> + <macros> + <import>macros.xml</import> + </macros> + <expand macro="requirements" /> + <expand macro="edam"/> + <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 expect_num_outputs="1"> + <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"> + <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 expect_num_outputs="3"> + <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 expect_num_outputs="1"> + <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 Apr 02 20:20:54 2021 +0000 @@ -0,0 +1,25 @@ +<macros> + <token name="@WRAPPER_VERSION@">1.0.15</token> + <token name="@VERSION_SUFFIX@">1</token> + <token name="@PROFILE@">20.09</token> + <xml name="requirements"> + <requirements> + <requirement type="package" version="@WRAPPER_VERSION@">bamutil</requirement> + </requirements> + </xml> + <xml name="edam"> + <edam_topics> + <edam_topic>topic_0091</edam_topic> + <edam_topic>topic_0080</edam_topic> + </edam_topics> + <edam_operations> + <edam_operation>operation_0369</edam_operation> + </edam_operations> + </xml> + <xml name="citations"> + <citations> + <citation type="doi">10.1101/gr.176552.114</citation> + </citations> + </xml> +</macros> +
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/in1.sam Fri Apr 02 20:20:54 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 Apr 02 20:20:54 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 Apr 02 20:20:54 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/input.sam Fri Apr 02 20:20:54 2021 +0000 @@ -0,0 +1,178 @@ +@SQ SN:>chrM LN:0 +@SQ SN:GTATTCTTACTCCATAAACACATAGGCTTGGTCCTAGCCTTTTTATTAGT LN:0 +@SQ SN:CTAAATCACGTCTCTACGATTAAAAGGAGCAGGTATCAAGCACACTAGAA LN:0 +@SQ SN:ATAAAAATTAAGCTATGAACGAAAGTTCGACTAAGTCATATTAAATAAGG LN:0 +@SQ SN:ATAAATCTCCGGCGTAAAGCGTGTCAAAGACTAATACCAAAATAAAGTTA LN:0 +@SQ SN:AAGTGACTTTAATACCTCTGACTACACGATAGCTAAGACCCAAACTGGGA LN:0 +@SQ SN:AAGCTATTCGCCAGAGTACTACTAGCAACAGCCTAAAACTCAAAGGACTT LN:0 +@SQ SN:CCGATAAACCCCACCATCCCTTGCTAATTCAGCCTATATACCGCCATCTT LN:0 +@SQ SN:AACGTTAGGTCAAGGTGTAGCCCATGGGATGGAGAGAAATGGGCTACATT LN:0 +@SQ SN:TGGAGACTAAAGGAGGATTTAGCAGTAAATTAAGAATAGAGAGCTTAATT LN:0 +@SQ SN:CACAAATCATAACATAACATAAAACCGTGACCCAAACATATGAAAGGAGA LN:0 +@SQ SN:TGTAGCTTAAACAAAGCATCCAGCTTACACCTAGAAGATTTCACTCAAAA LN:0 +@SQ SN:TTAGTCACTTAACTAAAACATTCACCAAACCATTAAAGTATAGGAGATAG LN:0 +@SQ SN:ATGCATTAAAAGTACTAAACAGCAAAGCTTACCCCTTTTACCTTTTGCAT LN:0 +@SQ SN:CGAAACCAGACGAGCTACCTATGAACAGTTACAAATGAACCAACTCATCT LN:0 +@SQ SN:AGCCTGGTGATAGCTGGTTGTCCAGAAACAGAATTTCAGTTCAAATTTAA LN:0 +@SQ SN:AAAGGTACAGCTTTTTAGATACAGGTTACAACCTTCATTAGAGAGTAAGA LN:0 +@SQ SN:TTCAAGCTCAACGACACATCTATCTTAATCCCAACAATCAACCCAAACTA LN:0 +@SQ SN:TTAATATGAGTAACAAGAATTATTTCTCCTTGCATAAGCTTATATCAGAA LN:0 +@SQ SN:TCATCTATTTAAACCATTGTTAACCCAACACAGGCATGCATCTATAAGGA LN:0 +@SQ SN:ACCAAAAACATCACCTCTAGCATTTCCAGTATTAGAGGCACTGCCTGCCC LN:0 +@SQ SN:taatcacttgttccctaaatagggacttgtatgaatggccacacgagggt LN:0 +@SQ SN:cgggaatgactaaataagacgagaagaccctatggagcttTAATTAACTG LN:0 +@SQ SN:TTGATTGAATCAGCAATTTCGGTTGGGGTGACCTCGGAGAACAAAACAAC LN:0 +@SQ SN:TTGATCCAAACCATTGATCAACGGAACAAGTTACCCTAGGGATAACAGCG LN:0 +@SQ SN:TTGGATCAAGACATCCTAATGGTGCAACCGCTATTAAGGGTTCGTTTGTT LN:0 +@SQ SN:CGGTTTCTATCTATTCTATACTTTTCCCAGTACGAAAGGACAAGAAAAGT LN:0 +@SQ SN:AATCTAACTAATTTATAACTTCTACCGCCCTAGAACAGGGCTCgttaggg LN:0 +@SQ SN:caactcctctccctaacaacaTGTTCATAATTAACGTCCTCCTCCTAATT LN:0 +@SQ SN:CTTAGGCTATATGCAACTTCGCAAAGGACCCAACATCGTAGGCCCCTATG LN:0 +@SQ SN:CTACAACCACTAACATCATCGACATCCATATTCATCATCGCACCAATCCT LN:0 +@SQ SN:CACTAATCAACATAAACCTAGGAATTCTATTCATACTAGCCATGTCCAGC LN:0 +@SQ SN:CGCCCTAATTGGAGCTCTACGAGCAGTAGCACAAACCATCTCATACGAAG LN:0 +@SQ SN:ACATTATCAACACTTATTATTACCCAAGAATACCTCTGATTAATCTTCCC LN:0 +@SQ SN:ACCGAGCTCCATTTGACCTAACAGAAGGAGAATCAGAACTCGTCTCTGGA LN:0 +@SQ SN:ATACGCAAACATCATCATGATAAACATCTTCACAACAACCCTATTTCTAG LN:0 +@SQ SN:ATTAAAGCTCTCCTTCTAACATGTTCCTTCCTATGAATCCGAGCATCCTA LN:0 +@SQ SN:TACCACTCACACTAGCCCTCTGCATATGACACGTCTCACTTCCAATCATA LN:0 +@SQ SN:ACTTTGATAGAGTAAAACATAGAGGCTCAAACCCTCTTATTTctagaact LN:0 +@SQ SN:ttacaccatgtcctaCAAGTAAGGTCAGCTAAATAAGCTATCGGGCCCAT LN:0 +@SQ SN:CTTCACAACTATTCTAATAACAGTTCTTCTAGGAACTATAATCGTTATAA LN:0 +@SQ SN:GCCATTATCCCTATCCTAATAAAAAAGTACAATCCCCGAACCATAGAAGC LN:0 +@SQ SN:TAGCGATCATCATTAACCTCATACACTCAGGCCAATGAACAATCACAAAA LN:0 +@SQ SN:ACTTGGACTCACACCATTCCACTTCTGAGTACCCGAAGTCACACAGGGCA LN:0 +@SQ SN:ATATCAATCCTATATCAAATCTCACCCTCAATTAACCTAAATATCTTATT LN:0 +@SQ SN:AAACCCAACTACGAAAAATCATAGCATACTCGTCAATCGCGCATATAGGA LN:0 +@SQ SN:ATTAATTTACATTATAATAACACTCACAATATTCATACTATTTATCCACA LN:0 +@SQ SN:CTAACCACTACACTAATCTTAATTACCTTACTATCCATAGGAGGCCTCCC LN:0 +@SQ SN:AAAATAGCAGCATCATCCTCCCCACACTAATAGCCATTATAGCACTACTC LN:0 +@SQ SN:CCCATCCACAAACAACATAAAAATAAAATGACAATTCGAAACCAAACGAA LN:0 +@SQ SN:ACCCCCATACTATCAATTTTGGACTAGGAATTTAGGTTAACATCCCAGAC LN:0 +@SQ SN:TAAGGACTGCGAGACTCTATCTCACATCAATTGAACGCAAATCAAACTCT LN:0 +@SQ SN:TTTAGTTAACAGCTAAATACCCTAATCAACTGGCTTCAATCTACTTCTCC LN:0 +@SQ SN:TCCTTTGAATTTGCAATTCAATGTGAAAATTCACCACGGGACTTGATAAG LN:0 +@SQ SN:CCATCTTACCTATGTTCATCAACCGCTGACTATTTTCAACTAACCACAAA LN:0 +@SQ SN:AACTGCCCTAAGCCTCCTAATCCGTGCTGAATTAGGCCAACCTGGGACCC LN:0 +@SQ SN:GTAATAATTTTCTTTATGGTCATACCCATTATAATCGGAGGATTCGGAAA LN:0 +@SQ SN:TAAACAACATAAGCTTCTGATTACTTCCCCCATCATTCCTACTTCTTCTC LN:0 +@SQ SN:TCCTCTAGCTGGAAATCTGGCGCATGCAGGAGCCTCTGTTGACTTAACCA LN:0 +@SQ SN:TTTATTACCACAATCATTAACATAAAACCACCAGCCCTATCCCAATATCA LN:0 +@SQ SN:TAGCCCTCCCGGTCCTAGCAGCAGGCATTACCATGCTTCTCACAGACCGT LN:0 +@SQ SN:TTATCAACACCTATTCTGATTCTTCGGACACCCCGAAGTCTATATTCTTA LN:0 +@SQ SN:AAAAAGGAACCTTTTGGCTACATGGGTATAGTGTGAGCTATAATATCCAT LN:0 +@SQ SN:TAGACGTTGACACACGAGCATACTTCACATCAGCTACCATAATCATCGCT LN:0 +@SQ SN:AAATATCAAATGATCTCCAGCTATACTCTGAGCTCTAGGCTTCATCTTCT LN:0 +@SQ SN:GATATTGTTCTCCACGATACTTATTATGTAGTAGCACATTTCCATTATGT LN:0 +@SQ SN:TCCCTCTATTCTCAGGATACACACTCAACCAAACCTGAGCAAAAATCCAC LN:0 +@SQ SN:CCTTGGCCTCTCAGGAATGCCACGACGCTATTCTGATTATCCAGACGCAT LN:0 +@SQ SN:GCAGTGATACTAATAATTTTCATAATTTGAGAAGCGTTCGCATCCAAACG LN:0 +@SQ SN:GATGCCCCCCACCATACCACACATTTGAAGAACCCACCTACGTAAACCTA LN:0 +@SQ SN:tcataaccactatgtctttctcCATCAATTGAGGTATTAGTAAAAATTAC LN:0 +@SQ SN:GCCTACCCCTTCCAACTAGGATTCCAAGACGCAACATCCCCTATTATAGA LN:0 +@SQ SN:GCTCTCTAGTATTATATATTATCTCATCAATACTAACAACTAAATTAACC LN:0 +@SQ SN:ACCAGCCATCATCCTTATTCTAATCGCCCTCCCATCCCTACGAATTCTAT LN:0 +@SQ SN:CACCAATGATACTGAAGCTACGAGTATACCGATTACGAAGACTTGACCTT LN:0 +@SQ SN:TTCTAGAAGTCGACAATCGAGTGGTTCTCCCCATAGAAATAACCATCCGA LN:0 +@SQ SN:AGGCCTAAAAACAGACGCTATCCCTGGGCGCCTAAATCAGACAACTCTCG LN:0 +@SQ SN:TCAAACCACAGCTTTATACCAATTGTCCTTGAACTAGTTCCACTGAAACA LN:0 +@SQ SN:TAGCATTAACCTTTTAAGTTAAAGATTGAGGGTTCAACCCCCTCCCTAGT LN:0 +@SQ SN:AATCCTAACTCTATTTATTGTATTTCAACTAAAAATCTCAAAGCACTCCT LN:0 +@SQ SN:CCTTGAGAATCAAAATGAACGAAAATCTATTCGCCTCTTTCGCTACCCCA LN:0 +@SQ SN:CCTATTCCCCTCACCCAACCGACTAATCAACAATCGCCTAATCTCAATTC LN:0 +@SQ SN:AGCAAAGGACAAACCTGAACTCTTATACTCATATCACTGATCCTATTCAT LN:0 +@SQ SN:CACAACTATCAATAAACCTAGGCATAGCTATTCCCCTATGGGCAGGGACA LN:0 +@SQ SN:ACCTCAAGGGACGCCCATTTTCCTCATCCCCATACTAGTAATTATCGAGA LN:0 +@SQ SN:AACATTACCGCCGGACACCTCCTAATACACCTCATCGGAGGGGCAACACT LN:0 +@SQ SN:TAATTCTACTAACTATCCTCGAATTCGCAGTAGCTATAATCCAAGCCTAC LN:0 +@SQ SN:CACCAAACCCACGCTTACCACATAGTAAACCCCAGCCCATGACCACTTAC LN:0 +@SQ SN:ACTTTAACTCAACCTTACTTCTAGCTATAGGGCTATTAACTAACATCCTT LN:0 +@SQ SN:CCATCACACATCAATCGTTCAAAAGGGACTCCGATATGGCATAATCCTTT LN:0 +@SQ SN:CACTCAAGCCTAGCCCCCACACCCGAACTAGGCGGCTGCTGACCACCCAC LN:0 +@SQ SN:TGCTCCTAGCATCTGGAGTCTCTATCACCTGAGCCCACCATAGCCTAATA LN:0 +@SQ SN:AGGCGTATACTTCACCCTTCTCCAAGCCTCAGAATACTATGAAGCCTCAT LN:0 +@SQ SN:TTCCACGGACTACACGTAATTATCGGATCTACCTTCCTCATTGTATGTTT LN:0 +@SQ SN:AAGCAGCCGCTTGATACTGACACTTCGTCGACGTAGTCTGACTATTCTTG LN:0 +@SQ SN:CAATTGACTTCCAATCAATCAGCTTCGGTATAACCCGAAAAAGAATAATA LN:0 +@SQ SN:ACTCATCGCATTCTGACTACCACAACTAAACATCTATGCAGAAAAAACCA LN:0 +@SQ SN:TCAATAAAATTTTTCTTAGTGGCCATTACATTTCTGCTATTCGACTTAGA LN:0 +@SQ SN:ACACTATACTTATCATAGCACTAGTCCTAATCTCTCTTCTAGCCATCAGC LN:0 +@SQ SN:TTAGTTTAAACCAAAACAAATGATTTCGACTCATTAAACTATGATTAACT LN:0 +@SQ SN:ACAGTATCCCTCGTAGGCCTACTAATGTACCGATCCCACCTAATATCCTC LN:0 +@SQ SN:TAATAGTCCTAAACACCCACTTCACACTAGCTAGTATAATACCTATCATC LN:0 +@SQ SN:CATAGTCTCCAATACTTATGGAGTAGACCACGTACAAAACCTTAACCTCC LN:0 +@SQ SN:GACTATCAAAAAAGAATATAATCTGAATCAACACTACAACCTATAGTCTA LN:0 +@SQ SN:CCTAAACTTCTCACTAATATTCTTCTCCGATCCCCTATCAGCCCCACTTC LN:0 +@SQ SN:CATCTATCTAAGGAACCACTAATCCGAAAAAAACTCTACATCACCATGCT LN:0 +@SQ SN:TCTCCTTCTACATCCTATTTGAAGCCACATTAGTTCCAACACTAATTATC LN:0 +@SQ SN:CCTATTCTACACACTAATAGGTTCCCTCCCACTCTTAGTTGCACTAATCT LN:0 +@SQ SN:AACCAAGCACTACCCGACTCTTGATCCAATATTTTCCTATGACTAGCATG LN:0 +@SQ SN:TCCCAAAAGCCCATGTAGAAGCCCCAATTGCCGGATCCATAGTGCTAGCA LN:0 +@SQ SN:ACTAAACCCCCAAACTAGCTTTATAGCCTACCCCTTCCTCATACTATCCC LN:0 +@SQ SN:AAATCACTTATTGCATACTCCTCTGTCAGCCACATAGCCCTAGTAATCGT LN:0 +@SQ SN:TAATCGCTCACGGCCTTACATCATCAATACTATTCTGCCTGGCAAACTCA LN:0 +@SQ SN:AACACTTCTTCCCCTTATAGCAGCCTGATGACTATTAGCCAGCCTAACCA LN:0 +@SQ SN:ATATCATCATTCTCATGATCAAATATTACCATTATCCTAATAGGAGCCAA LN:0 +@SQ SN:GAGGGAAATACACACACCATATCAACAGCATTAAACCTTCATTTACACGA LN:0 +@SQ SN:TAACCCTAAAATTATCCTAGGCTTTACGTACTGTAAATATAGTTTAACAA LN:0 +@SQ SN:CGAGAAAGTATGCAAGAACTGCTAATTCATGCCCCCATGTCCAACAAACA LN:0 +@SQ SN:CCAAAAAATTGGTGCAACTCCAAATAAAAGTAATCAACATGTTCTCCTCC LN:0 +@SQ SN:CTTCAATACCTACAAAAACAGCACGTTCCCGCATCATGTAAAAAACACTA LN:0 +@SQ SN:TCTGGACAAGAAACAATTATCTCAAACTGACACTGAATAACCATACAAAC LN:0 +@SQ SN:TACCAGTAGCCCTATTCGTAACATGATCTATTATGGAATTCTCCCTATGA LN:0 +@SQ SN:ATTCCTCATCACTATAATAATTCTAGTCACAGCTAACAACCTTTTCCAAC LN:0 +@SQ SN:TGATGATACGGCCGAACAGATGCCAACACCGCGGCCCTTCAAGCAATCCT LN:0 +@SQ SN:TATTCAACACCAACACATGAGACCTCCAACAAATCTTCATACTCGACCCC LN:0 +@SQ SN:ATCCGCTCAATTTGGACTCCACCCATGACTTCCTTCAGCCATAGAGGGCC LN:0 +@SQ SN:GTCTTCCTGCTAATCCGCTTCCATCCACTAATAGAAAACAACAAAACAAT LN:0 +@SQ SN:TCTGCGCACTCACTCAAAACGATATCAAAAAAATCATTGCTTTCTCCACC LN:0 +@SQ SN:CCTAGCATTCCTCCACATTTGCACTCACGCATTCTTCAAAGCTATACTAT LN:0 +@SQ SN:CGAAAAATAGGCGGACTATTTAATGCAATACCCTTCACCACCACATCTCT LN:0 +@SQ SN:ACTCCAAAGACCTCATCATCGAAACCGCCAACACATCGTACACCAACGCC LN:0 +@SQ SN:TACCCGAATCATCTTCTTTGCACTCCTAGGGCAACCCCGCTTCCTCCCTC LN:0 +@SQ SN:CTCTTAATTGGCAGCATTTTTGCCGGATTCTTCATCTCCAACAATATCTA LN:0 +@SQ SN:CCCTCGCAGTAACCATCCTAGGATTTACACTAGCCCTAGAACTAAGCTTG LN:0 +@SQ SN:CAACCTCCTAGGATACTACCCAACAATTATACACCGACTCCCACCGCTCG LN:0 +@SQ SN:TGACTAGAAAACATCCTGCCAAAATCTATCTCCCAGTTCCAAATAAAAAC LN:0 +@SQ SN:CATTCCTCATCACCCTTACCCTAAGCATACTACTTTTTAATCTCCACGAG LN:0 +@SQ SN:ATCACAACCCAAGCCCCATAACTATACAATGCAGCAGCCCCTATAATTTC LN:0 +@SQ SN:CACCACTAAACTTAAACACTACCCCCACTTCCTCACTCTTCAGAACATAT LN:0 +@SQ SN:AGTCGTATTAGACACCCATACCTCAGGATACTGCTCAGTAGCCATAGCCG LN:0 +@SQ SN:ATCAACCCCAAAAAGGACCCTCCAAAATTCATAATAATACCACAACCTAC LN:0 +@SQ SN:AAGAAAACCCCACAAAACTAACAACAAAAATAACACTCAAAATAAACACA LN:0 +@SQ SN:GAAAAATCATCGTTGTATTTCAACTATAAGAACACCAATGACAAACATCC LN:0 +@SQ SN:CCAGCCCCCTCAAACATTTCATCATGATGAAACTTCGGCTCCCTCCTAGG LN:0 +@SQ SN:ACACATCAGACACGACAACTGCCTTCTCATCCGTCACTCACATCTGCCGA LN:0 +@SQ SN:AATATTTTTTATCTGCCTCTTCATTCACGTAGGACGCGGCCTCTACTACG LN:0 +@SQ SN:ACAGTTATAGCTACAGCATTCATGGGCTATGTCCTACCATGAGGCCAAAT LN:0 +@SQ SN:ACATCGGTACTACCCTCGTCGAGTGAATCTGAGGTGGATTCTCAGTAGAC LN:0 +@SQ SN:CATCACAGCCCTGGTAGTCGTACATTTACTATTTCTTCACGAAACAGGAT LN:0 +@SQ SN:CCATATTATACAATTAAAGACATCCTAGGACTCCTCCTCCTGATCTTGCT LN:0 +@SQ SN:ACTACACCCCAGCTAACCCTCTCAGCACTCCCCCTCATATTAAACCAGAA LN:0 +@SQ SN:AGGCGGCGTATTAGCCCTAATCCTCTCCATCCTGATCCTAGCACTCATCC LN:0 +@SQ SN:CAATGCGTATTCTGACTCTTAGTGGCAGACTTACTGACACTAACATGAAT LN:0 +@SQ SN:CAATCCTCTACTTCTCCCTAATTCTCATTTTTATACCACTCGCAAGCACC LN:0 +@SQ SN:ACCCTGGTCTTGTAAACCAGAAAAGGGGGAAAACGTTTCCTCCCAAGGAC LN:0 +@SQ SN:TACTTAAACTATTCCTTGATTTCTTCCCCTAAACGACAACAATTTACCCT LN:0 +@SQ SN:CTGACATGCAATATCTTATGAATGGCCTATGTACGTCGTGCATTAAATTG LN:0 +@SQ SN:GTACATTATATTATTGATCGTGCATACCCCATCCAAGTCAAATCATTTCC LN:0 +@SQ SN:GCGGGAAATCAGCAACCCTCCCAACTACGTGTCCCAATCCTCGCTCCGGG LN:0 +@SQ SN:TCTTTCTTCAGGGCCATTCCCACCCAACCTCGCCCATTCTTTCCCCTTAA LN:0 +@SQ SN:CTGTGATTTCATGCATTTGGTATCTTTTTATATTTGGGGATGCTATGACT LN:0 +@SQ SN:CTTAAATTGAACGTTATTCCTCCGCATCAGCAACCATAAGGTGTTATTCA LN:0 +@SQ SN:ctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacct LN:0 +@SQ SN:gcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgc LN:0 +@SQ SN:ctgtgcacctACCCGCGCAGTAAGCAAGTAATATAGCTTTCTTAATCAAA LN:0 +@SQ SN:GCCAAACCCCAAAAACAAGACTAAACAATGCACAATACTTCATGAAGCTT LN:0 +@SQ SN:GAACTTTCCCCCCGCCATTAATACCAACATGCTACTTTAATCAATAAAAT LN:0 +@SQ SN:TTCTTCCCCC LN:0 +HWI-EAS91_1_30788AAXX:1:1:1698:516 16 * 2735 25 36M * 0 0 TTTACACTCAGAGGTTCAACTCCTCTCNNTAACAAC I9IIIII5IIIIIIIIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1398:854 16 * 3921 25 36M * 0 0 CACCCTTCCCGTACTAATAAATCCCCTNNTCTTCAC IIIII=AIIIIIIIIIIIIIIBIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1716:413 0 * 6040 25 36M * 0 0 GATCCAANNCTTTATCAACACCTATTCTGATTCTTC IIIIIII""IIIIIIIIIIIIIIIIIIIIIIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1634:211 0 * 9350 25 36M * 0 0 GAAGCAGNNGCTTGATACTGACACTTCGTCGACGTA IIIIIII""IIIIIIIIIIIIIIIIIIIIII9IIDF NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1513:715 16 * 9563 25 36M * 0 0 CTGACTACCACAACTAAACATCTATGCNNAAAAAAC I+-II?IDIIIIIIIIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1310:991 16 * 10002 25 36M * 0 0 CTCCTATGCCTAGAAGGAATAATACTANNACTATTC I:2IEI:IIDIIIIII4IIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1711:249 16 * 10617 25 36M * 0 0 ACCAAACAGAACGCCTGAACGCAGGCCNNTACTTCC IIIIIIIIIIIIIIIIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1491:637 16 * 10864 25 36M * 0 0 TGTAGAAGCCCCAATTGCCGGATCCATNNTGCTAGC DBAIIIIIIIIIIIFIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1630:59 16 * 12387 25 36M * 0 0 TCATACTCGACCCCAACCTTACCAACCNNCCGCTCC FIIHII;IIIIIIIIIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27 +HWI-EAS91_1_30788AAXX:1:1:1218:141 16 * 14062 25 36M * 0 0 ACAAAACTAACAACAAAAATAACACTCNNAATAAAC I+IIII1IIIIIIIIIIIIIIIIIIII""IIIIIII NM:i:1 X1:i:1 MD:Z:7N0N27
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.sam Fri Apr 02 20:20:54 2021 +0000 @@ -0,0 +1,178 @@ +@SQ SN:>chrM LN:0 +@SQ SN:GTATTCTTACTCCATAAACACATAGGCTTGGTCCTAGCCTTTTTATTAGT LN:0 +@SQ SN:CTAAATCACGTCTCTACGATTAAAAGGAGCAGGTATCAAGCACACTAGAA LN:0 +@SQ SN:ATAAAAATTAAGCTATGAACGAAAGTTCGACTAAGTCATATTAAATAAGG LN:0 +@SQ SN:ATAAATCTCCGGCGTAAAGCGTGTCAAAGACTAATACCAAAATAAAGTTA LN:0 +@SQ SN:AAGTGACTTTAATACCTCTGACTACACGATAGCTAAGACCCAAACTGGGA LN:0 +@SQ SN:AAGCTATTCGCCAGAGTACTACTAGCAACAGCCTAAAACTCAAAGGACTT LN:0 +@SQ SN:CCGATAAACCCCACCATCCCTTGCTAATTCAGCCTATATACCGCCATCTT LN:0 +@SQ SN:AACGTTAGGTCAAGGTGTAGCCCATGGGATGGAGAGAAATGGGCTACATT LN:0 +@SQ SN:TGGAGACTAAAGGAGGATTTAGCAGTAAATTAAGAATAGAGAGCTTAATT LN:0 +@SQ SN:CACAAATCATAACATAACATAAAACCGTGACCCAAACATATGAAAGGAGA LN:0 +@SQ SN:TGTAGCTTAAACAAAGCATCCAGCTTACACCTAGAAGATTTCACTCAAAA LN:0 +@SQ SN:TTAGTCACTTAACTAAAACATTCACCAAACCATTAAAGTATAGGAGATAG LN:0 +@SQ SN:ATGCATTAAAAGTACTAAACAGCAAAGCTTACCCCTTTTACCTTTTGCAT LN:0 +@SQ SN:CGAAACCAGACGAGCTACCTATGAACAGTTACAAATGAACCAACTCATCT LN:0 +@SQ SN:AGCCTGGTGATAGCTGGTTGTCCAGAAACAGAATTTCAGTTCAAATTTAA LN:0 +@SQ SN:AAAGGTACAGCTTTTTAGATACAGGTTACAACCTTCATTAGAGAGTAAGA LN:0 +@SQ SN:TTCAAGCTCAACGACACATCTATCTTAATCCCAACAATCAACCCAAACTA LN:0 +@SQ SN:TTAATATGAGTAACAAGAATTATTTCTCCTTGCATAAGCTTATATCAGAA LN:0 +@SQ SN:TCATCTATTTAAACCATTGTTAACCCAACACAGGCATGCATCTATAAGGA LN:0 +@SQ SN:ACCAAAAACATCACCTCTAGCATTTCCAGTATTAGAGGCACTGCCTGCCC LN:0 +@SQ SN:taatcacttgttccctaaatagggacttgtatgaatggccacacgagggt LN:0 +@SQ SN:cgggaatgactaaataagacgagaagaccctatggagcttTAATTAACTG LN:0 +@SQ SN:TTGATTGAATCAGCAATTTCGGTTGGGGTGACCTCGGAGAACAAAACAAC LN:0 +@SQ SN:TTGATCCAAACCATTGATCAACGGAACAAGTTACCCTAGGGATAACAGCG LN:0 +@SQ SN:TTGGATCAAGACATCCTAATGGTGCAACCGCTATTAAGGGTTCGTTTGTT LN:0 +@SQ SN:CGGTTTCTATCTATTCTATACTTTTCCCAGTACGAAAGGACAAGAAAAGT LN:0 +@SQ SN:AATCTAACTAATTTATAACTTCTACCGCCCTAGAACAGGGCTCgttaggg LN:0 +@SQ SN:caactcctctccctaacaacaTGTTCATAATTAACGTCCTCCTCCTAATT LN:0 +@SQ SN:CTTAGGCTATATGCAACTTCGCAAAGGACCCAACATCGTAGGCCCCTATG LN:0 +@SQ SN:CTACAACCACTAACATCATCGACATCCATATTCATCATCGCACCAATCCT LN:0 +@SQ SN:CACTAATCAACATAAACCTAGGAATTCTATTCATACTAGCCATGTCCAGC LN:0 +@SQ SN:CGCCCTAATTGGAGCTCTACGAGCAGTAGCACAAACCATCTCATACGAAG LN:0 +@SQ SN:ACATTATCAACACTTATTATTACCCAAGAATACCTCTGATTAATCTTCCC LN:0 +@SQ SN:ACCGAGCTCCATTTGACCTAACAGAAGGAGAATCAGAACTCGTCTCTGGA LN:0 +@SQ SN:ATACGCAAACATCATCATGATAAACATCTTCACAACAACCCTATTTCTAG LN:0 +@SQ SN:ATTAAAGCTCTCCTTCTAACATGTTCCTTCCTATGAATCCGAGCATCCTA LN:0 +@SQ SN:TACCACTCACACTAGCCCTCTGCATATGACACGTCTCACTTCCAATCATA LN:0 +@SQ SN:ACTTTGATAGAGTAAAACATAGAGGCTCAAACCCTCTTATTTctagaact LN:0 +@SQ SN:ttacaccatgtcctaCAAGTAAGGTCAGCTAAATAAGCTATCGGGCCCAT LN:0 +@SQ SN:CTTCACAACTATTCTAATAACAGTTCTTCTAGGAACTATAATCGTTATAA LN:0 +@SQ SN:GCCATTATCCCTATCCTAATAAAAAAGTACAATCCCCGAACCATAGAAGC LN:0 +@SQ SN:TAGCGATCATCATTAACCTCATACACTCAGGCCAATGAACAATCACAAAA LN:0 +@SQ SN:ACTTGGACTCACACCATTCCACTTCTGAGTACCCGAAGTCACACAGGGCA LN:0 +@SQ SN:ATATCAATCCTATATCAAATCTCACCCTCAATTAACCTAAATATCTTATT LN:0 +@SQ SN:AAACCCAACTACGAAAAATCATAGCATACTCGTCAATCGCGCATATAGGA LN:0 +@SQ SN:ATTAATTTACATTATAATAACACTCACAATATTCATACTATTTATCCACA LN:0 +@SQ SN:CTAACCACTACACTAATCTTAATTACCTTACTATCCATAGGAGGCCTCCC LN:0 +@SQ SN:AAAATAGCAGCATCATCCTCCCCACACTAATAGCCATTATAGCACTACTC LN:0 +@SQ SN:CCCATCCACAAACAACATAAAAATAAAATGACAATTCGAAACCAAACGAA LN:0 +@SQ SN:ACCCCCATACTATCAATTTTGGACTAGGAATTTAGGTTAACATCCCAGAC LN:0 +@SQ SN:TAAGGACTGCGAGACTCTATCTCACATCAATTGAACGCAAATCAAACTCT LN:0 +@SQ SN:TTTAGTTAACAGCTAAATACCCTAATCAACTGGCTTCAATCTACTTCTCC LN:0 +@SQ SN:TCCTTTGAATTTGCAATTCAATGTGAAAATTCACCACGGGACTTGATAAG LN:0 +@SQ SN:CCATCTTACCTATGTTCATCAACCGCTGACTATTTTCAACTAACCACAAA LN:0 +@SQ SN:AACTGCCCTAAGCCTCCTAATCCGTGCTGAATTAGGCCAACCTGGGACCC LN:0 +@SQ SN:GTAATAATTTTCTTTATGGTCATACCCATTATAATCGGAGGATTCGGAAA LN:0 +@SQ SN:TAAACAACATAAGCTTCTGATTACTTCCCCCATCATTCCTACTTCTTCTC LN:0 +@SQ SN:TCCTCTAGCTGGAAATCTGGCGCATGCAGGAGCCTCTGTTGACTTAACCA LN:0 +@SQ SN:TTTATTACCACAATCATTAACATAAAACCACCAGCCCTATCCCAATATCA LN:0 +@SQ SN:TAGCCCTCCCGGTCCTAGCAGCAGGCATTACCATGCTTCTCACAGACCGT LN:0 +@SQ SN:TTATCAACACCTATTCTGATTCTTCGGACACCCCGAAGTCTATATTCTTA LN:0 +@SQ SN:AAAAAGGAACCTTTTGGCTACATGGGTATAGTGTGAGCTATAATATCCAT LN:0 +@SQ SN:TAGACGTTGACACACGAGCATACTTCACATCAGCTACCATAATCATCGCT LN:0 +@SQ SN:AAATATCAAATGATCTCCAGCTATACTCTGAGCTCTAGGCTTCATCTTCT LN:0 +@SQ SN:GATATTGTTCTCCACGATACTTATTATGTAGTAGCACATTTCCATTATGT LN:0 +@SQ SN:TCCCTCTATTCTCAGGATACACACTCAACCAAACCTGAGCAAAAATCCAC LN:0 +@SQ SN:CCTTGGCCTCTCAGGAATGCCACGACGCTATTCTGATTATCCAGACGCAT LN:0 +@SQ SN:GCAGTGATACTAATAATTTTCATAATTTGAGAAGCGTTCGCATCCAAACG LN:0 +@SQ SN:GATGCCCCCCACCATACCACACATTTGAAGAACCCACCTACGTAAACCTA LN:0 +@SQ SN:tcataaccactatgtctttctcCATCAATTGAGGTATTAGTAAAAATTAC LN:0 +@SQ SN:GCCTACCCCTTCCAACTAGGATTCCAAGACGCAACATCCCCTATTATAGA LN:0 +@SQ SN:GCTCTCTAGTATTATATATTATCTCATCAATACTAACAACTAAATTAACC LN:0 +@SQ SN:ACCAGCCATCATCCTTATTCTAATCGCCCTCCCATCCCTACGAATTCTAT LN:0 +@SQ SN:CACCAATGATACTGAAGCTACGAGTATACCGATTACGAAGACTTGACCTT LN:0 +@SQ SN:TTCTAGAAGTCGACAATCGAGTGGTTCTCCCCATAGAAATAACCATCCGA LN:0 +@SQ SN:AGGCCTAAAAACAGACGCTATCCCTGGGCGCCTAAATCAGACAACTCTCG LN:0 +@SQ SN:TCAAACCACAGCTTTATACCAATTGTCCTTGAACTAGTTCCACTGAAACA LN:0 +@SQ SN:TAGCATTAACCTTTTAAGTTAAAGATTGAGGGTTCAACCCCCTCCCTAGT LN:0 +@SQ SN:AATCCTAACTCTATTTATTGTATTTCAACTAAAAATCTCAAAGCACTCCT LN:0 +@SQ SN:CCTTGAGAATCAAAATGAACGAAAATCTATTCGCCTCTTTCGCTACCCCA LN:0 +@SQ SN:CCTATTCCCCTCACCCAACCGACTAATCAACAATCGCCTAATCTCAATTC LN:0 +@SQ SN:AGCAAAGGACAAACCTGAACTCTTATACTCATATCACTGATCCTATTCAT LN:0 +@SQ SN:CACAACTATCAATAAACCTAGGCATAGCTATTCCCCTATGGGCAGGGACA LN:0 +@SQ SN:ACCTCAAGGGACGCCCATTTTCCTCATCCCCATACTAGTAATTATCGAGA LN:0 +@SQ SN:AACATTACCGCCGGACACCTCCTAATACACCTCATCGGAGGGGCAACACT LN:0 +@SQ SN:TAATTCTACTAACTATCCTCGAATTCGCAGTAGCTATAATCCAAGCCTAC LN:0 +@SQ SN:CACCAAACCCACGCTTACCACATAGTAAACCCCAGCCCATGACCACTTAC LN:0 +@SQ SN:ACTTTAACTCAACCTTACTTCTAGCTATAGGGCTATTAACTAACATCCTT LN:0 +@SQ SN:CCATCACACATCAATCGTTCAAAAGGGACTCCGATATGGCATAATCCTTT LN:0 +@SQ SN:CACTCAAGCCTAGCCCCCACACCCGAACTAGGCGGCTGCTGACCACCCAC LN:0 +@SQ SN:TGCTCCTAGCATCTGGAGTCTCTATCACCTGAGCCCACCATAGCCTAATA LN:0 +@SQ SN:AGGCGTATACTTCACCCTTCTCCAAGCCTCAGAATACTATGAAGCCTCAT LN:0 +@SQ SN:TTCCACGGACTACACGTAATTATCGGATCTACCTTCCTCATTGTATGTTT LN:0 +@SQ SN:AAGCAGCCGCTTGATACTGACACTTCGTCGACGTAGTCTGACTATTCTTG LN:0 +@SQ SN:CAATTGACTTCCAATCAATCAGCTTCGGTATAACCCGAAAAAGAATAATA LN:0 +@SQ SN:ACTCATCGCATTCTGACTACCACAACTAAACATCTATGCAGAAAAAACCA LN:0 +@SQ SN:TCAATAAAATTTTTCTTAGTGGCCATTACATTTCTGCTATTCGACTTAGA LN:0 +@SQ SN:ACACTATACTTATCATAGCACTAGTCCTAATCTCTCTTCTAGCCATCAGC LN:0 +@SQ SN:TTAGTTTAAACCAAAACAAATGATTTCGACTCATTAAACTATGATTAACT LN:0 +@SQ SN:ACAGTATCCCTCGTAGGCCTACTAATGTACCGATCCCACCTAATATCCTC LN:0 +@SQ SN:TAATAGTCCTAAACACCCACTTCACACTAGCTAGTATAATACCTATCATC LN:0 +@SQ SN:CATAGTCTCCAATACTTATGGAGTAGACCACGTACAAAACCTTAACCTCC LN:0 +@SQ SN:GACTATCAAAAAAGAATATAATCTGAATCAACACTACAACCTATAGTCTA LN:0 +@SQ SN:CCTAAACTTCTCACTAATATTCTTCTCCGATCCCCTATCAGCCCCACTTC LN:0 +@SQ SN:CATCTATCTAAGGAACCACTAATCCGAAAAAAACTCTACATCACCATGCT LN:0 +@SQ SN:TCTCCTTCTACATCCTATTTGAAGCCACATTAGTTCCAACACTAATTATC LN:0 +@SQ SN:CCTATTCTACACACTAATAGGTTCCCTCCCACTCTTAGTTGCACTAATCT LN:0 +@SQ SN:AACCAAGCACTACCCGACTCTTGATCCAATATTTTCCTATGACTAGCATG LN:0 +@SQ SN:TCCCAAAAGCCCATGTAGAAGCCCCAATTGCCGGATCCATAGTGCTAGCA LN:0 +@SQ SN:ACTAAACCCCCAAACTAGCTTTATAGCCTACCCCTTCCTCATACTATCCC LN:0 +@SQ SN:AAATCACTTATTGCATACTCCTCTGTCAGCCACATAGCCCTAGTAATCGT LN:0 +@SQ SN:TAATCGCTCACGGCCTTACATCATCAATACTATTCTGCCTGGCAAACTCA LN:0 +@SQ SN:AACACTTCTTCCCCTTATAGCAGCCTGATGACTATTAGCCAGCCTAACCA LN:0 +@SQ SN:ATATCATCATTCTCATGATCAAATATTACCATTATCCTAATAGGAGCCAA LN:0 +@SQ SN:GAGGGAAATACACACACCATATCAACAGCATTAAACCTTCATTTACACGA LN:0 +@SQ SN:TAACCCTAAAATTATCCTAGGCTTTACGTACTGTAAATATAGTTTAACAA LN:0 +@SQ SN:CGAGAAAGTATGCAAGAACTGCTAATTCATGCCCCCATGTCCAACAAACA LN:0 +@SQ SN:CCAAAAAATTGGTGCAACTCCAAATAAAAGTAATCAACATGTTCTCCTCC LN:0 +@SQ SN:CTTCAATACCTACAAAAACAGCACGTTCCCGCATCATGTAAAAAACACTA LN:0 +@SQ SN:TCTGGACAAGAAACAATTATCTCAAACTGACACTGAATAACCATACAAAC LN:0 +@SQ SN:TACCAGTAGCCCTATTCGTAACATGATCTATTATGGAATTCTCCCTATGA LN:0 +@SQ SN:ATTCCTCATCACTATAATAATTCTAGTCACAGCTAACAACCTTTTCCAAC LN:0 +@SQ SN:TGATGATACGGCCGAACAGATGCCAACACCGCGGCCCTTCAAGCAATCCT LN:0 +@SQ SN:TATTCAACACCAACACATGAGACCTCCAACAAATCTTCATACTCGACCCC LN:0 +@SQ SN:ATCCGCTCAATTTGGACTCCACCCATGACTTCCTTCAGCCATAGAGGGCC LN:0 +@SQ SN:GTCTTCCTGCTAATCCGCTTCCATCCACTAATAGAAAACAACAAAACAAT LN:0 +@SQ SN:TCTGCGCACTCACTCAAAACGATATCAAAAAAATCATTGCTTTCTCCACC LN:0 +@SQ SN:CCTAGCATTCCTCCACATTTGCACTCACGCATTCTTCAAAGCTATACTAT LN:0 +@SQ SN:CGAAAAATAGGCGGACTATTTAATGCAATACCCTTCACCACCACATCTCT LN:0 +@SQ SN:ACTCCAAAGACCTCATCATCGAAACCGCCAACACATCGTACACCAACGCC LN:0 +@SQ SN:TACCCGAATCATCTTCTTTGCACTCCTAGGGCAACCCCGCTTCCTCCCTC LN:0 +@SQ SN:CTCTTAATTGGCAGCATTTTTGCCGGATTCTTCATCTCCAACAATATCTA LN:0 +@SQ SN:CCCTCGCAGTAACCATCCTAGGATTTACACTAGCCCTAGAACTAAGCTTG LN:0 +@SQ SN:CAACCTCCTAGGATACTACCCAACAATTATACACCGACTCCCACCGCTCG LN:0 +@SQ SN:TGACTAGAAAACATCCTGCCAAAATCTATCTCCCAGTTCCAAATAAAAAC LN:0 +@SQ SN:CATTCCTCATCACCCTTACCCTAAGCATACTACTTTTTAATCTCCACGAG LN:0 +@SQ SN:ATCACAACCCAAGCCCCATAACTATACAATGCAGCAGCCCCTATAATTTC LN:0 +@SQ SN:CACCACTAAACTTAAACACTACCCCCACTTCCTCACTCTTCAGAACATAT LN:0 +@SQ SN:AGTCGTATTAGACACCCATACCTCAGGATACTGCTCAGTAGCCATAGCCG LN:0 +@SQ SN:ATCAACCCCAAAAAGGACCCTCCAAAATTCATAATAATACCACAACCTAC LN:0 +@SQ SN:AAGAAAACCCCACAAAACTAACAACAAAAATAACACTCAAAATAAACACA LN:0 +@SQ SN:GAAAAATCATCGTTGTATTTCAACTATAAGAACACCAATGACAAACATCC LN:0 +@SQ SN:CCAGCCCCCTCAAACATTTCATCATGATGAAACTTCGGCTCCCTCCTAGG LN:0 +@SQ SN:ACACATCAGACACGACAACTGCCTTCTCATCCGTCACTCACATCTGCCGA LN:0 +@SQ SN:AATATTTTTTATCTGCCTCTTCATTCACGTAGGACGCGGCCTCTACTACG LN:0 +@SQ SN:ACAGTTATAGCTACAGCATTCATGGGCTATGTCCTACCATGAGGCCAAAT LN:0 +@SQ SN:ACATCGGTACTACCCTCGTCGAGTGAATCTGAGGTGGATTCTCAGTAGAC LN:0 +@SQ SN:CATCACAGCCCTGGTAGTCGTACATTTACTATTTCTTCACGAAACAGGAT LN:0 +@SQ SN:CCATATTATACAATTAAAGACATCCTAGGACTCCTCCTCCTGATCTTGCT LN:0 +@SQ SN:ACTACACCCCAGCTAACCCTCTCAGCACTCCCCCTCATATTAAACCAGAA LN:0 +@SQ SN:AGGCGGCGTATTAGCCCTAATCCTCTCCATCCTGATCCTAGCACTCATCC LN:0 +@SQ SN:CAATGCGTATTCTGACTCTTAGTGGCAGACTTACTGACACTAACATGAAT LN:0 +@SQ SN:CAATCCTCTACTTCTCCCTAATTCTCATTTTTATACCACTCGCAAGCACC LN:0 +@SQ SN:ACCCTGGTCTTGTAAACCAGAAAAGGGGGAAAACGTTTCCTCCCAAGGAC LN:0 +@SQ SN:TACTTAAACTATTCCTTGATTTCTTCCCCTAAACGACAACAATTTACCCT LN:0 +@SQ SN:CTGACATGCAATATCTTATGAATGGCCTATGTACGTCGTGCATTAAATTG LN:0 +@SQ SN:GTACATTATATTATTGATCGTGCATACCCCATCCAAGTCAAATCATTTCC LN:0 +@SQ SN:GCGGGAAATCAGCAACCCTCCCAACTACGTGTCCCAATCCTCGCTCCGGG LN:0 +@SQ SN:TCTTTCTTCAGGGCCATTCCCACCCAACCTCGCCCATTCTTTCCCCTTAA LN:0 +@SQ SN:CTGTGATTTCATGCATTTGGTATCTTTTTATATTTGGGGATGCTATGACT LN:0 +@SQ SN:CTTAAATTGAACGTTATTCCTCCGCATCAGCAACCATAAGGTGTTATTCA LN:0 +@SQ SN:ctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacct LN:0 +@SQ SN:gcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgcacctgtgc LN:0 +@SQ SN:ctgtgcacctACCCGCGCAGTAAGCAAGTAATATAGCTTTCTTAATCAAA LN:0 +@SQ SN:GCCAAACCCCAAAAACAAGACTAAACAATGCACAATACTTCATGAAGCTT LN:0 +@SQ SN:GAACTTTCCCCCCGCCATTAATACCAACATGCTACTTTAATCAATAAAAT LN:0 +@SQ SN:TTCTTCCCCC LN:0 +HWI-EAS91_1_30788AAXX:1:1:1698:516 16 * 2735 25 36M * 0 0 TTTACACTCAGAGGTTCAACTCCTCTCNNTAACAAC I9IIIII5IIIIIIIIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1398:854 16 * 3921 25 36M * 0 0 CACCCTTCCCGTACTAATAAATCCCCTNNTCTTCAC IIIII=AIIIIIIIIIIIIIIBIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1716:413 0 * 6040 25 36M * 0 0 GATCCAANNCTTTATCAACACCTATTCTGATTCTTC IIIIIII""IIIIIIIIIIIIIIIIIIIIIIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1634:211 0 * 9350 25 36M * 0 0 GAAGCAGNNGCTTGATACTGACACTTCGTCGACGTA IIIIIII""IIIIIIIIIIIIIIIIIIIIII9IIDF MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1513:715 16 * 9563 25 36M * 0 0 CTGACTACCACAACTAAACATCTATGCNNAAAAAAC I+-II?IDIIIIIIIIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1310:991 16 * 10002 25 36M * 0 0 CTCCTATGCCTAGAAGGAATAATACTANNACTATTC I:2IEI:IIDIIIIII4IIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1711:249 16 * 10617 25 36M * 0 0 ACCAAACAGAACGCCTGAACGCAGGCCNNTACTTCC IIIIIIIIIIIIIIIIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1491:637 16 * 10864 25 36M * 0 0 TGTAGAAGCCCCAATTGCCGGATCCATNNTGCTAGC DBAIIIIIIIIIIIFIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1630:59 16 * 12387 25 36M * 0 0 TCATACTCGACCCCAACCTTACCAACCNNCCGCTCC FIIHII;IIIIIIIIIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1 +HWI-EAS91_1_30788AAXX:1:1:1218:141 16 * 14062 25 36M * 0 0 ACAAAACTAACAACAAAAATAACACTCNNAATAAAC I+IIII1IIIIIIIIIIIIIIIIIIII""IIIIIII MD:Z:7N0N27 NM:i:1 X1:i:1
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output_stats.txt Fri Apr 02 20:20:54 2021 +0000 @@ -0,0 +1,8 @@ +Overlap Statistics: +Number of overlapping pairs: 0 +Average # Reference Bases Overlapped: 0 +Variance of Reference Bases overlapped: 0 +Number of times orientation causes additional clipping: 0 +Number of times the forward strand was clipped: 0 +Number of times the reverse strand was clipped: 0 +Completed ClipOverlap Successfully.