From cc4a7110282bb2cf1ff14925f7ceda1cea5e15c5 Mon Sep 17 00:00:00 2001 From: "Rahamim, Ben" Date: Wed, 18 Nov 2020 17:03:28 +0200 Subject: [PATCH 1/3] Added an optimisation to not spare not required shuffle of the reads when they are unpaired --- .../adam/cli/TransformFragments.scala | 2 +- .../org/bdgenomics/adam/rdd/ADAMContext.scala | 2 +- .../adam/rdd/read/AlignmentDataset.scala | 10 +- .../adam/rdd/read/SingleReadBucket.scala | 44 +- adam-core/src/test/resources/unpaired.sam | 592 ++++++++++++++++++ .../rdd/fragment/FragmentDatasetSuite.scala | 18 +- .../adam/rdd/read/MarkDuplicatesSuite.scala | 2 +- 7 files changed, 648 insertions(+), 22 deletions(-) create mode 100644 adam-core/src/test/resources/unpaired.sam diff --git a/adam-cli/src/main/scala/org/bdgenomics/adam/cli/TransformFragments.scala b/adam-cli/src/main/scala/org/bdgenomics/adam/cli/TransformFragments.scala index 125f62830a..7af05922c9 100644 --- a/adam-cli/src/main/scala/org/bdgenomics/adam/cli/TransformFragments.scala +++ b/adam-cli/src/main/scala/org/bdgenomics/adam/cli/TransformFragments.scala @@ -130,7 +130,7 @@ class TransformFragments(protected val args: TransformFragmentsArgs) extends BDG val rdd = if (args.loadAsAlignments) { sc.loadAlignments(args.inputPath) - .toFragments + .toFragments() } else { sc.loadFragments(args.inputPath) } diff --git a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/ADAMContext.scala b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/ADAMContext.scala index a793573143..ec4fb6812e 100644 --- a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/ADAMContext.scala +++ b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/ADAMContext.scala @@ -3680,7 +3680,7 @@ class ADAMContext(@transient val sc: SparkContext) extends Serializable with Log .querynameSortedToFragments } else { info(s"Loading $pathName as BAM/CRAM/SAM and converting to Fragments.") - loadBam(pathName, stringency).toFragments + loadBam(pathName, stringency).toFragments() } } else if (isInterleavedFastqExt(trimmedPathName)) { info(s"Loading $pathName as interleaved FASTQ and converting to Fragments.") diff --git a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/AlignmentDataset.scala b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/AlignmentDataset.scala index 08f5219bc6..0304ac9f76 100644 --- a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/AlignmentDataset.scala +++ b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/AlignmentDataset.scala @@ -463,11 +463,17 @@ sealed abstract class AlignmentDataset extends AvroReadGroupGenomicDataset[Align /** * Convert this set of reads into fragments. * + * @param mergePairedReads Merge paired reads into a fragment. Can be set to false when only unpaired + * reads are in the AlignmentDataset. Defaults to true. * @return Returns a FragmentDataset where all reads have been grouped together by * the original sequence fragment they come from. */ - def toFragments(): FragmentDataset = { - FragmentDataset(groupReadsByFragment().map(_.toFragment), + def toFragments(mergePairedReads: Boolean = true): FragmentDataset = { + val fragmentRDD = if (mergePairedReads) groupReadsByFragment().map(_.toFragment) + else + rdd.map(SingleReadBucket.buildFragmentFromSingleAlignment) + + FragmentDataset(fragmentRDD, sequences, readGroups, processingSteps) diff --git a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala index 6ddecdaa97..1ff459da50 100644 --- a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala +++ b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala @@ -105,6 +105,27 @@ private[read] object SingleReadBucket { def apply(fragment: Fragment): SingleReadBucket = { fromGroupedReads(fragment.getAlignments.toIterable) } + + def buildFragmentFromSingleAlignment(alignment: Alignment): Fragment = { + buildFragmentFromAlignments(List(alignment), alignment.getReadName, Option(alignment.getReadGroupId), Option(alignment.getInsertSize).map(_.toInt)) + } + + def buildFragmentFromAlignments(reads: Seq[Alignment], readName: String, readGroupId: Option[String], insertSize: Option[Int]): Fragment = { + + val builder = Fragment.newBuilder() + .setName(readName) + .setAlignments(seqAsJavaList(reads)) + + // is an insert size defined for this fragment? + insertSize.foreach(is => { + builder.setInsertSize(is) + }) + + // set read group name, if known + readGroupId.foreach(n => builder.setReadGroupId(n)) + + builder.build() + } } /** @@ -139,24 +160,17 @@ private[adam] case class SingleReadBucket( // want to pay the cost exactly once val unionReads = allReads - // start building fragment - val builder = Fragment.newBuilder() - .setName(unionReads.head.getReadName) - .setAlignments(seqAsJavaList(allReads.toSeq)) + val firstRead = unionReads.head - // is an insert size defined for this fragment? - primaryMapped.headOption - .foreach(r => { - Option(r.getInsertSize).foreach(is => { - builder.setInsertSize(is.toInt) - }) - }) + val insertSize = primaryMapped.headOption match { + case Some(alignment) => Option(alignment.getInsertSize).map(_.toInt) + case None => None + } - // set read group name, if known - Option(unionReads.head.getReadGroupId) - .foreach(n => builder.setReadGroupId(n)) + SingleReadBucket.buildFragmentFromAlignments(unionReads.toSeq, + firstRead.getReadName, Option(firstRead.getReadGroupId), + insertSize) - builder.build() } } diff --git a/adam-core/src/test/resources/unpaired.sam b/adam-core/src/test/resources/unpaired.sam new file mode 100644 index 0000000000..5db227d6de --- /dev/null +++ b/adam-core/src/test/resources/unpaired.sam @@ -0,0 +1,592 @@ +@PG ID:bwa PN:bwa VN:0.7.5a-r405 +@RG ID:SRR062634 LB:2845856850 SM:HG00096 PI:206 CN:WUGSC PL:ILLUMINA DS:SRP001294 +@SQ SN:1 LN:249250621 +@SQ SN:10 LN:135534747 +@SQ SN:11 LN:135006516 +@SQ SN:12 LN:133851895 +@SQ SN:13 LN:115169878 +@SQ SN:14 LN:107349540 +@SQ SN:15 LN:102531392 +@SQ SN:16 LN:90354753 +@SQ SN:17 LN:81195210 +@SQ SN:18 LN:78077248 +@SQ SN:19 LN:59128983 +@SQ SN:2 LN:243199373 +@SQ SN:20 LN:63025520 +@SQ SN:21 LN:48129895 +@SQ SN:22 LN:51304566 +@SQ SN:3 LN:198022430 +@SQ SN:4 LN:191154276 +@SQ SN:5 LN:180915260 +@SQ SN:6 LN:171115067 +@SQ SN:7 LN:159138663 +@SQ SN:8 LN:146364022 +@SQ SN:9 LN:141213431 +@SQ SN:GL000191.1 LN:106433 +@SQ SN:GL000192.1 LN:547496 +@SQ SN:GL000193.1 LN:189789 +@SQ SN:GL000194.1 LN:191469 +@SQ SN:GL000195.1 LN:182896 +@SQ SN:GL000196.1 LN:38914 +@SQ SN:GL000197.1 LN:37175 +@SQ SN:GL000198.1 LN:90085 +@SQ SN:GL000199.1 LN:169874 +@SQ SN:GL000200.1 LN:187035 +@SQ SN:GL000201.1 LN:36148 +@SQ SN:GL000202.1 LN:40103 +@SQ SN:GL000203.1 LN:37498 +@SQ SN:GL000204.1 LN:81310 +@SQ SN:GL000205.1 LN:174588 +@SQ SN:GL000206.1 LN:41001 +@SQ SN:GL000207.1 LN:4262 +@SQ SN:GL000208.1 LN:92689 +@SQ SN:GL000209.1 LN:159169 +@SQ SN:GL000210.1 LN:27682 +@SQ SN:GL000211.1 LN:166566 +@SQ SN:GL000212.1 LN:186858 +@SQ SN:GL000213.1 LN:164239 +@SQ SN:GL000214.1 LN:137718 +@SQ SN:GL000215.1 LN:172545 +@SQ SN:GL000216.1 LN:172294 +@SQ SN:GL000217.1 LN:172149 +@SQ SN:GL000218.1 LN:161147 +@SQ SN:GL000219.1 LN:179198 +@SQ SN:GL000220.1 LN:161802 +@SQ SN:GL000221.1 LN:155397 +@SQ SN:GL000222.1 LN:186861 +@SQ SN:GL000223.1 LN:180455 +@SQ SN:GL000224.1 LN:179693 +@SQ SN:GL000225.1 LN:211173 +@SQ SN:GL000226.1 LN:15008 +@SQ SN:GL000227.1 LN:128374 +@SQ SN:GL000228.1 LN:129120 +@SQ SN:GL000229.1 LN:19913 +@SQ SN:GL000230.1 LN:43691 +@SQ SN:GL000231.1 LN:27386 +@SQ SN:GL000232.1 LN:40652 +@SQ SN:GL000233.1 LN:45941 +@SQ SN:GL000234.1 LN:40531 +@SQ SN:GL000235.1 LN:34474 +@SQ SN:GL000236.1 LN:41934 +@SQ SN:GL000237.1 LN:45867 +@SQ SN:GL000238.1 LN:39939 +@SQ SN:GL000239.1 LN:33824 +@SQ SN:GL000240.1 LN:41933 +@SQ SN:GL000241.1 LN:42152 +@SQ SN:GL000242.1 LN:43523 +@SQ SN:GL000243.1 LN:43341 +@SQ SN:GL000244.1 LN:39929 +@SQ SN:GL000245.1 LN:36651 +@SQ SN:GL000246.1 LN:38154 +@SQ SN:GL000247.1 LN:36422 +@SQ SN:GL000248.1 LN:39786 +@SQ SN:GL000249.1 LN:38502 +@SQ SN:MT LN:16569 +@SQ SN:NC_007605 LN:171823 +@SQ SN:X LN:155270560 +@SQ SN:Y LN:59373566 +@SQ SN:hs37d5 LN:35477943 +SRR062634.10022079 163 22 16082635 0 100M = 16082719 179 CTGGAGTGCAGTGGCATGATTTCAGCTCACTGTCGTCTCTGCCTCCCTGACTCAAGTGATTCTCCTGCCTCAGCCTCCCACGTCGCTCGGACTCCACGCC A:=D5D5E?D?DDD:.@@@@=?EE=DADDB@D=DD??ED=:CCCC?D:E=EEB=-C>C=@=EEEEB5EC-?A>=C-C?DC+34+4A>-?5:=/-A=@>>: RG:Z:SRR062634 XT:A:R NM:i:5 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:5 XO:i:0 XG:i:0 MD:Z:80A2A3G5A5G0 XA:Z:14,-19760295,100M,5;GL000244.1,-31171,100M,5; +SRR062634.10072753 163 22 16090555 0 99M1S = 16090636 181 CCCAGGCCCCAGGTGAACACCAGGTTTTAGGTGAACACGAGGTCCTAGGTGGGTGTCTATGCTTCTGGTGAACCTCAGGCCCTAGTGGACACTCAGGCCC GGGGGGFGGGEEECEGGGGGGDGE?EGGEG?CF?FGGGGEEFEFGGEEBEEED9AACBBBEEEGEEED?@CCEE:CB=CE@CCA@>?AACAB5AA8???# RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:14,-19752379,100M,0;GL000244.1,-23245,100M,1;14,+19823345,100M,1; +SRR062634.10091189 147 22 16056515 15 100M = 16056440 -175 AGGGTTCAGGCAGCTGTCCAAACACCTTATCTTTTCATCTCTGACCATTTCTTCAGAGTCGACTATGAGCAGAAATCACCCAGACACCCCACAGGTAAGA EEEEEBE=EEEEGGEEDBGEFEFAAFEFEBGGFGEGEDGEEEFEGGGFGFGGFGGGGGGGGGGGGGGGGGGGGGGFFDFFGGGGGGGDGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:99T0 XA:Z:14,+19786394,100M,1; +SRR062634.10151245 147 22 16086391 0 10S90M = 16086297 -184 TAGATCTGTGTCCCACTGGGGCCTTGTGTTCACCGGGGACTGGTATCCAGCTGTGGCCTGATGATCTACCGCATCCTGTTGCTCACCTATGGGCTGGTGT ##########????7??=@?-B+@B?BCB5B=BBEBB=?BDFEEB=CEB??FEFFEEDFF?:CACEFBFFEEB?CCEAEAAEEEBDDDDADFDFDEFFFE RG:Z:SRR062634 XC:i:90 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 XA:Z:GL000244.1,+27418,100M,0;14,+19756552,100M,1; +SRR062634.1016484 147 22 16063682 23 100M = 16063611 -171 ACAAGTGGGTTTCTCTAAAGCACAGTTTGGCTTTTTACTAGAATATCAAGTCACCTTTTTTGGGTGAGGAGAGTACAAAGGGAGTGGGCAGCTGGTGAAA BEEEEEAEEE?EEEEDDCDCFEBEEFCEEAEEEEEGEFGGEGGGFFGFEEEFBFGFGGGGFGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19779190,100M,1; +SRR062634.10182517 147 22 16064048 0 100M = 16063977 -171 AAAAACTGCAATATGATGTGTAAAGTCAGTGTAATCTTATGGAAGTCATATAAAAGTCACAACATGAGGCATGAGATACCTCAGATGTTGTTCCCTCTAA GEEFEEEEEEEEBEFEGGEGEGFGGEEEGGGGEGFEGGGBGGFGFFGFGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:39C60 XA:Z:14,+19778824,100M,1; +SRR062634.10194383 163 22 16070721 0 42M58S = 16070835 199 TCAGCTGATACAAGAAATCCAGAGAACTCCCGTGAGACTGAAGGTAGGACGTGACGCATATCTTTGCCTTCTTGTGTGGGTCTGAGTATTCTCTGGCTAT --;55,676(-4,956+9;.?A5A5@C-5C>-CC>?=,=@C?########################################################## RG:Z:SRR062634 XC:i:42 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:0C41 XA:Z:14,-19772208,100M,1;2,-132538703,100M,1; +SRR062634.10210679 163 22 16056340 22 75M25S = 16056452 192 AGGAAGCAGTCAGCAAGTGCCATCCACCCCAGCTTTTCTCTCCCTGGGGACTGCTGGACTGAGGGTCCAGGTGGCCCAGCTCCAAAGGCTGTGGCGGGCC FFDFFFFFFEEFEEFDDBEEFFFFFBFFFF?FFEEFFBFEFEEEEE?FDAB:B:@=A@B@A@:AA2A??5?;?7?######################### RG:Z:SRR062634 XC:i:75 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:75 XA:Z:14,-19786594,100M,0; +SRR062634.10257696 163 22 16086456 18 100M = 16086559 203 CTGTTGCTCACCTATGGGCTGGTGTCTACCTGGGGCTTGGTGATCACCTGGGAGCCAGATATCAACTTGGGGCCTGGGTGTCCACTTAAGGCCTGATGTG GGGFGGGGGGGGGGGGGGDFGFGGGGGFGGGGGGGGGGGFEGDFEGFGGGDF?FDGDFFGGBEGGGGDEGFEEBEFEB7FEDFDDEGEEEED5BF?DFDE RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:18 AM:i:0 X0:i:1 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19756477,100M,1;14,+19817593,100M,1;14,-19367444,100M,1; +SRR062634.10262546 147 22 16086560 15 100M = 16086434 -226 TGGGGCCTGATTGCCCACCCGGGGACTGGGTGTCCACCTGGGGTCTGATGTCCACCTGAAGTTAGGTATCTACCTAAGGCTTGGTGTCTACCTGTGGCCT 9DCAAB?CCB=BB?CA=AAAECCE>@EECEDEEEEEEEDGDEEGEGGEGGDDDBDEEGEEGFFEGFGGGEGBFFGEGGGEGGGGGGGGGFFGFGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.10437037 163 22 16058717 23 99M1S = 16058762 145 AGAGATCTCCATGGCAGTCCATCACATCACAGACCTTGAGGCCTAAGAGGAAAGAATGCTTTTGAGGGCCAGGACTTTACTGCTCCACACAGCCCGGGGA GGGEGEGGGGGGGGGGGGGEGGGFG?GGAGGGGFGFFFEG=FBFEBGEFEEGDEAEEE?EEGGEEBDBCEACDCFEAC?BEDEDCD=ECD?CBDDA@A5) RG:Z:SRR062634 XC:i:99 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:14,-19784191,100M,1; +SRR062634.10448889 117 22 16079761 0 * = 16079761 0 TTTCTTTCTTTTATATATATATACACACACACACACACACACACACATATATGTATATATACACGTATATGTATGTATATATGTATATATACACGTATAT @DF>C;FDC=EGEGGEFDGEFDD?DFDEEGFGFGGGDGGGGGGGEGGGGFGGGFGGGGGGFGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 +SRR062634.10483312 147 22 16092813 0 100M = 16092712 -201 TGTCCCCAGGTTGACATCACTCTCAAGGTGTACATTAGGCAACAGATGTACACCCAGGTCCAAGGCAGACACGAGTCCCCAGTAAAACTCAAGGCCCCAG 57:AAE=CCEDCCCECCE?EAECDEAEEEEDBAEDDGBEFBEBCBFDFFDEAFBF?FFEGDGGFGGDGGFEEGGAAGFEEGGGGGGGGFGGGGGFGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.10583394 147 22 16074056 0 60S40M = 16073897 -199 CGCGCTTACCTATATTCCCTATGCTTCCTCCCACAAAAGAGCTTAGGGTGGCTGCCAACTTGTGGGCAGAGACCCTCACTTTCCAATCCCCACGAGGGGC ############################################################@@=(C@:=A=?&6DBEDADBDDD RG:Z:SRR062634 XC:i:40 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:5 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:40 +SRR062634.1064731 163 22 16096744 23 69M31S = 16096883 208 CCACCATGGGGTACCCCAAGAGGCCTGAAGTGCACCCTGGTCTTGAACCAGGGATGCCAGGGTCCCCTGGGCCCAGCTCCGGGGCTGGTGGGAAGACACT GGGGGGGGGFGGGGGGGDGBGGGGFF=DFBACEEEED=CCCEED=CDEEB=CBCAEEEEED?7EABBB?<>C?*>EB:BE?ABC5A<:A-=C@@>-:5@::DC=DE+EED=DDA:EDEEDEBEEBD?DCDCC?5C=AD=D RG:Z:SRR062634 XC:i:88 XT:A:U NM:i:1 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:58C29 +SRR062634.10780186 163 22 16070732 0 100M = 16070817 185 AAGAAATCCAGAGAACTCCCGTGAGACTGAAGGTAGGCAGTGAAGCAAATGTTTGCATTCTTGTGTGGCTCTGATTAGCATCAGGATATCTACATTTTAC GGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGGFGGGGGGGGGGGGGGGGGEFGGGBGGGGFGEFEFGGEEEEEEEEEDECEFBE?CDDDFDEFEEEEEF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19772139,100M,0;9,+45410516,100M,1;2,-132538634,100M,1; +SRR062634.10902474 163 22 16069603 55 84M16S = 16069708 194 TCTGAATTACACCACCAAACAAAGATACATCACCACACCAAGCCACCACAGGAGCCTATCTGCAACCAAGGAACTTGTACAGAGCCTTGGTCATCTGAAA EEEDAFDFFFFEFFFDDFFBFFCCEEEEAEFDFF?FFDFDEBBCEEE:E?F:B?DEEEEFD?DDBFEEFB=?@B=?B?B?@?@@################ RG:Z:SRR062634 XC:i:84 XT:A:U NM:i:1 SM:i:18 AM:i:18 X0:i:1 X1:i:3 XM:i:1 XO:i:0 XG:i:0 MD:Z:0C83 XA:Z:14,-19773285,100M,2;9,-42282437,100M,2;9,+45409384,100M,2; +SRR062634.10940591 147 22 16056291 0 100M = 16056178 -213 GTGCTGGAATCAGCCCATAGTCTTTGCGTGTCCCATGGGGCAGGGCCAAAGGAAGCAGTCAGCAACTGCCATCCACCCCAGCTTTTCTCTCCCTGGGGAC AA?=9?AA4:B:B?=@@,CC7BAB:CCC:C>>>>=+EEEEBDGFGCBEABFF:FFDFFBF=EEFEFDFDFECC>FE?GDEEDEDEDA:EEDEGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:65G34 XA:Z:14,+19786618,100M,1; +SRR062634.11087705 147 22 16069500 0 100M = 16069388 -212 GCTTTTGGCACTTCGTCACTTCTGCACCTGTCGAGGACAAGTACACAGAGGAGAAGACCCTGCAGCCAGCTGTTACTCTTAAGCACCGTCTACTGGACTG =@:A@=?5=6<=5<:=5,:>BCCBC?DC:DCBD=EEEEGEEECECBEFEGEFEGEFADBFFAFBGGEFGGGGFEGBGGGGGGGGEGEGDGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:0 AM:i:0 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:86T13 XA:Z:21,+10629843,100M,0; +SRR062634.11091461 147 22 16050505 0 100M = 16050471 -134 CAACTGAGATGCCCAAAGGGGGCCACTCTCCTTGCTTTTCCTCCTTTTTAGAGGATTTATTTCCCATTTTTCTTAAAAAGGAAGAACAAACTGTGCCCTA =@@?CCCD?@DCDCC@CEFEEGEEEFFFFDEFEGAGGFGGGGGGGGGGGGGFGGGGGFGGDGEGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19792397,100M,0; +SRR062634.11146114 147 22 16068990 23 10S90M = 16068896 -184 GAATTAGGCTCTCTCCTCTTACAGGACTGTAGTGGGAAGAGAGTTGCTAAAGCAGAGGTTTCTCAGGGATGGTGAGACTTGCAGCCAAGAACAGTTTTGC ##########C,@>=5B??5ABA:A>@(C5?D=C:5BBA:5EEDEEE5DEEEE?BE?BBA-BE?EEECC:5DAADBEEDEEDDE=EEEEEEEE?AEEE:E RG:Z:SRR062634 XC:i:90 XT:A:U NM:i:1 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:10G79 XA:Z:14,+19773892,100M,2; +SRR062634.11150957 147 22 16089407 0 8S92M = 16089311 -188 ACCCATGCTTCAGGTGCACACTAGGCCCAAAGTGTACACCCAAGCCCAGGTGGGCATCAGGCCCAAGGTGTACACCAGACCCCAAGTGGACATCAGGTTC ########??:731.5?-?=B?5:AEBDDCC?:?B>A?A5>>,9@.@?)(7:4<:-:=D5A?EED:EEDEEBEBDEEEED:EEFFACE?AEEE?EDDDDE RG:Z:SRR062634 XC:i:92 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:34T57 +SRR062634.11155401 147 22 16072370 0 4S96M = 16072277 -189 TCTGTTAACAAAATGAAATCGATCTTTTTGTACTGTGGTGTAAAGATTTTCTGGACATAGTTCAGGTAAAAAAAAAAAAAAAGACCACTGTTGAATAATG ####?03;,;C??--?CB5?@330<@DBC:AA=A?@@B@?:CDD-CDEEC?BC@??A:A-EC=CCEAGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:96 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:96 XA:Z:14,+19770501,100M,0; +SRR062634.11163767 163 22 16073705 0 99M1S = 16073791 186 TGGGACCTCCCAGGAGGTGTAGAATATTCCACTAGGGACACCTCATACCCTTCCAGGATTAGACCTTGGGGCCTGGAGATCCCCAGGCAATTAGTATTGA GGGGGGGGGGEFGGDEECECGGGGBGGFGGGGEGE?GEFGEGEFFGEFEGD:FFDC=CAACCCE@:CC:C?EBAA?C?AACBCC?:CC>?7>16=..0.0AA87?5BCAC:?CEFECE:EDBEBAE?BAF=AGECFADFFDGFGG=EGFGGGGFGGEDGGFGFFGEEGGGGGGFGGGGGGG RG:Z:SRR062634 XC:i:84 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:84 XA:Z:14,+19790367,100M,0; +SRR062634.11307291 147 22 16070204 0 2S98M = 16070127 -175 GACCAGACCTTTGAATCAACCCAGTCAGACAAAAATAAAAAAGGAATTTAAAAATAATTAACCACACGTTTGAGAAATATAGGATTATGTAAAGGAACAA ##@?>-5=C:55:5-A:C:=@C=;.E-:?CEBEEE?AA:ADD=DD:EACCEE6DB@C@6>,CCCCC?-C??DD=DBEE:BB?DD?DB:AD?AD:AEEAE: RG:Z:SRR062634 XC:i:98 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:56G8T32 XA:Z:14,+19772669,100M,2; +SRR062634.11329415 163 22 16063043 36 100M = 16063112 169 TGGGCAACAAAGTGAGACCCAATCTCTACAAAAACTAGAGAAAAAATAGCTGGGTGTGGTGGCACTCACCTGTACTAAGGAGGCTGGGCTGGGAGGATCT GGGFGGGEGGGGFGGGGGGGGGGGBGGGGDGGGGGGDEFGEGFDGGFGGGGGGBEGFEGDFGBEFGGDEBDGEEEGEEFECFBBCEEEACEDD=DC?C?C RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:2,-132546315,100M,0;14,-19779829,100M,1; +SRR062634.11340799 147 22 16068541 0 100M = 16068458 -183 GGAAAGCAGCCTCTGCGATAGCATTTTGTTCATAAAAGGGAGTGAAACTGCAGTACGAGAGAAGAACATCCAGTCTCCCTCTGTGACTCACCTTTCCACT ?0?=,?CC@CEEEEEEECEEEE=EEDDCD@?E=GGEEEEEGEEEGDEFEGEGGGEDGGGGGGGGGEGGFGGGGFGEGEGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19774331,100M,0; +SRR062634.11413841 147 22 16072769 0 100M = 16072689 -180 GAACCCCGAATCTGTCAGGACATATGTGTGATCTTGGATAACGAATTCAACCTACTTGGCACTGAATTTTCTTGACTGTAAACTGTTTTAGAGATGTTCC CB?@@########################### RG:Z:SRR062634 XC:i:73 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:68T4 XA:Z:14,-19786774,100M,1;18,-14988210,100M,1;2,-132576152,100M,2; +SRR062634.11502577 147 22 16091746 37 100M = 16091651 -195 GTGGATACTAGACCCCAGGGGTACACCAGACTCCTGGTGGGCATAAGGCCCCACGAGGACACTAGAATCCAGGTGTATGTAAAGCCACTGGTTGACACCA DCCCBCCD5EEEEEE?DFDDEAEGEGEEGEEEEGEEGDEGGEGGGFGFFFGGEGGGGGGGEGGGFEGDGGGGGGGGGGGBGGGGGGGGGGFGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.11519149 163 22 16082932 0 96M4S = 16083020 182 AAGAAACTCAGGAATATGCATTGAAAATGTCTCCCTCAGGTGATTCTGATGTGATATGTGGTCTGAGTTTAAAATGTAAGGAAAATTACCTTCCCTGCCC FGGGGFFFEEGDGGFGFGGGGGGFDFFFGDBDEEFDECACCECDADEEC@CA:ACBDDCDBC=?DC=:CC@@A@:@AC?CA:ABB?5???>CAAAA#### RG:Z:SRR062634 XC:i:96 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:96 XA:Z:14,-19760002,100M,0;GL000244.1,-30878,100M,0; +SRR062634.11577003 147 22 16058706 0 100M = 16058583 -223 ACATGCATTCCAGAGATCTCCATGGCAGTCCATCACATCACAGACCTTGAGGTCTAAGAGGAAAGAATGCTTTTGAGGGCCAGGACTTTACTGCTCCACA BDD@DBCCBCBDDDC=CDEEEEEEGEBDDB?>>->GFEGFGGEFGGE=GGGGEEDGGGGGGGGGGGGGCGGGGFFGGGGGGGGGGGGGGGGGGGGGGGEG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:52C47 XA:Z:14,+19784201,100M,1; +SRR062634.11605012 147 22 16055272 22 18S82M = 16055159 -195 CTGAGCCAAAGTGAGCCCCAGCTTAAGGGCCTGGGCAAAGGCCTCTTTGCCCTCACCTTGTCCTTCAGGCTCAGCAGCAGGCAACCCCAATTGCAGAAGT ##################@7>@=>=<<?><4=D5DDCB>>5?DD=ABE=EBBDD:DDDA?DA-FFFEFFFEEBEEA RG:Z:SRR062634 XC:i:82 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:82 XA:Z:14,+19787655,100M,0; +SRR062634.11637216 147 22 16061051 60 100M = 16060966 -185 AAACTTCAAAAGTCTATTGACAAGGAAGCTATGCATGTTGCAAGTTTCACTAAAGAATATGATGATGCATGTGGGGACTTTTTTGGCTGGTTGTGTTGTT ??:??A+CCC???:=@C=CC=CCAEEEEEEBAED=EEGFGEEEFGECEC:DEGFGGEDEGGFGGGFFFFEGBGGGGGEGGGGGFGGEEGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.11714131 163 22 16051397 23 100M = 16051489 192 AAGGTCAGAGTTTAATCCTCATTCCTAAGTTATGTAAACCAAAAATAAAATTCTGAAGATGTCCTGATCATCTGAATGGACCCTTCCTCTGGACCAGGGC EBDDEEEEEEDEEBEEADEEDEDEEEBBEDEEEEEDE=E=ADECEBACCC??AE?5AA?:EDEEE?EEEEEBEEBEEEADDEAD=BAB+C==?AA?@=BC RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19791505,100M,1; +SRR062634.11751606 147 22 16070184 0 100M = 16070108 -176 CTAATTTAAGAGCTCAAAGACCAGACCTTTGAATCAACCCAGTCAGACAAAAATAAAAAAGGAATTTAAAAATAATGAACCACACTTTTGAGAAATATAG @DFDEB@DFEEFB:EDDEFDCEEEA:GEEEFFBFBECFFDEFF:FFBAGFGGGDGGGGGEGGGGGE>GGGGGGFGFFFEFGEGFGGGFGGDGGGGEEGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19772687,100M,0; +SRR062634.11782831 147 22 16074220 0 100M = 16074118 -202 GAAGGTCACACCTGTCAGGGACCAGTGCTGGAGGGCAAACCTCACCTTTAAAAGCTCCCACTGTTTATTTTAAATTTATTTTTATTTTTAATTTCTATGA ED?DBCECE?ADEDEGEGGGEGGEGGFFGGGFGGGGGGDGGGGCGBGGGEGEFGGGFGGEDGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19768647,100M,0;GL000244.1,+39587,100M,0; +SRR062634.11783511 163 22 16083698 0 100M = 16083751 153 CCCAGATTAGCTAAAAGCACCAAAATACCAACTTGTACACCAAATTTTGGAAATGGATGCTAAAACTTTGATGGCAAAAGGTGACTGATCTGCCTGAGAA GGGGGGGGGGGGGGGFGGFGGGFGGGEGDGGGGGGGDGGGFGGGGGGGGFGBGGGFFGFGGGGGGEGGGGBGGGGFFBGDGCGBGCFEEEGGEFFBED:D RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:16A83 +SRR062634.11793884 147 22 16079816 29 100M = 16079715 -201 TGGCTCATGACTGTAATCGCAGCACTTTGGGAGGCTGAGACAGGTGGACCATCTGAGGTCAGGATTGAGAGCAGCCTGGCCAACATGGTGAAACCCTGTC 57B?AC:CD?CCC?CC5BBC@A@==D=BEFE:DEGEGGGEFFGGGGG=FGGGGEGGGGGGFGGGFGGGGGAGG=GGGGGGGGGGGFGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+33985,100M,0;14,+19763051,100M,0; +SRR062634.11800397 163 22 16059480 0 100M = 16059548 168 AATGCCTGTACCCCAATTGTGTCTTGAAGTAAATAACTTGTTTTTGATTTTACAGGCTCATAGGTGGAAAAGACTTGCCTTGCCTCAGATGAGAACTTGG GGGGGGGGEGGGGGEFGGFGFDGGGGGGGEGGGGGGGGGGEBGFGGBGGGGFGEGGFEGFGEEEEGECDFFEFFECDDDCADEECEA=CDDCECCACBDA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783427,100M,0; +SRR062634.11810243 147 22 16067017 0 100M = 16066952 -165 GTGAGTTCTCATGACATCTGGCTATTTAAAAGTGTGTGGCACCTCTCCCTGCTTGCTCCTGCTCTGGCCATGTGACATACCTTCTTTTGCTTCACCTTCT A=A8>>8A?CA@A=A2AAC5ACA=C?C??@?AA@###################################### RG:Z:SRR062634 XC:i:62 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:62 XA:Z:14,-19776943,100M,0;2,-132543432,100M,1; +SRR062634.12009427 163 22 16063939 0 100M = 16064013 174 CTAGCAAGACTTACAAGCAGAGCTCAGGGAAAAGTCTGCATGGGAAAACAGAGCCAGACTTAGAAGGCATTGCTGATGAATCCCATTTTCTTATTACAGG GGGGGGGFGGGGGGGGGGGGGGGGGGGEGGGGGGGGGGGGGGFGGGGGGGFGGGGGGGGGGGGGGGGFGGGGGGGGGGFGGGGGGEGGGGGFFGFBFFGE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19778933,100M,0;2,-132545419,100M,0; +SRR062634.12017865 163 22 16073760 0 100M = 16073883 208 GGATTAGACCTTGGGGCCTGGAGATCCCCAGGCAATTAGTATTGAAGGTCGAATGGCCAATGACAGGAATGGGAAGGCCCACTGTGTCATTCACAAAGCA FFFFFFGGGGDGGGGDGGEGGGGFDGGGGEGGGGGGGDFEEFG?GGDGCGGEEFGFDGFBDDEEDGFBEEBFA?EEEEAFDDCBAECDD:BDDDDADBDC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19769107,100M,0; +SRR062634.12080663 163 22 16090872 0 100M = 16091029 220 ACCCAGAGTCTACCTGAAATCTTAGTGTCAACCTGGGGCCTATGTGTCCACTTGGAGTCTGACGTGCACCTGAAGCCTGAGTTTTCACCTAGGATCTGAT GGGGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGEGGGFGGGGGGGFGFGEGEGGGEGGEGFCEFGFEGGAG?FEECDEFDEDDDEEEE:DED5D RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-22927,100M,0;14,-19752061,100M,0;14,+19823661,100M,1; +SRR062634.12099057 161 22 16095533 0 85M15S GL000211.1 45004 0 CCAGTCCCCAGCTGAACTCTGGGCCCCAGCTGAACGTCATAACCCAGATGGTCACCCAGGCTCCAGGGGAACACATAGTCTTAGGTGGACATCAGGCCCC AACA>=DDDDAGDEGGDCGG:EDE:=DDDD5?EEE:BE=EFFF=EC5=5?>=?C?*;43;?<>A:2-%)0>:;<=@5A@?C?5:?############### RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:67T17 XA:Z:14,+19827864,100M,1;14,-19747870,100M,1;GL000244.1,-18281,100M,1; +SRR062634.12130956 147 22 16065121 15 100M = 16065049 -172 AGAGACCTAGAGGTCCAGGATGAAAGAATGGTTTCCTGTGTAGTGGGATTTTAAAAATTAATTAATTTATTTACTTATTTATTCTGAGATGGAGTCTCAC DFEFABDAEGEGFAEAFGGGGFGGGGGBBEGGGGBFGGGGGFFGGFGGGGGGGGGFGGGGGGGGGGGFGGGGBGGGGGFGFFGGGGDGGGGGFGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19777751,100M,0; +SRR062634.12148752 147 22 16068783 0 100M = 16068678 -205 CTGGGACTGAGAGAAGGATGCCATTTTTAATCTGGGCTCGTACAAAGTCAATCATTCTTTGGTGCCTGGCAATGGTGGCCACTGCAGACACTTTAATCTT =@EEECC:?CECEC=DDDEBB5DFFFBEGE@F?FEFDEDEDDDGGEF:GGDGEEGBGGFFGGGGGDDGGFGGFGGGGGGGGGGGFGGGGFGGGGGGEGGF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19774089,100M,0; +SRR062634.12153358 163 22 16052684 60 100M = 16052752 168 AGTTTGTTCTGCCTAGTTTAAGGACATGCCCATGACACTGCCTCAGGAGGTCCTGACAGCATGTGCCCAAGGTGGTCAGGATACAGCTTGCTTCTATATA EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEDEEEEDEEEEEDEEEDEEEEDD:D@ RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.12268287 163 22 16065775 0 100M = 16065860 185 GTCTAGTGACCTTGCAGCTGCACAGCTAGAGAAATAGAGTCTTCACAATGCCTGGGAAAGGAAGAGAGATGAGGCTCACTAGCCACAGAAAAACAGGCAG GGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGFGEGGGGGGFFGFGGGGBGEFGGGEGGBGFGEEGGDFGEGGFEGEEGFEFEEDBB@C77?A?A RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19777097,100M,0; +SRR062634.1227834 147 22 16094053 0 100M = 16093960 -193 CATTCAGGCCCCCAGTAGTCATTTGGCCCCAAGTGGACACTCAGGCACCAGGTTCACATGATGTCTTAACTGGACACCAAGGGGCCAGTTTGATACCCAA ?B@AABEEEEEBEEE?EEEEACDCDEB?EEECEEEEEGEEEEEGGGEEEEEGGFEEGEGGGGGGFGGGGGFGGEFEEGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.12325816 147 22 16086461 41 100M = 16086455 -106 GCTCACCTATGGGCTGGTGTCTACCTGGGGCTTGGTGATCACCTGGGAGCCAGATATCAACTTGGGGCCTGGGTGTCCACTTAAGGCCTGATGTTGGCCT =C1C=CB?C=F@BBECDDDCDD@ADEFFEEGEEEEEFDBECFBBEEGGGGEEGFGGGGGEGGGGGGGGEGGGFGGGFGFGGGGGGBFGGGFGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:23 AM:i:18 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:94G5 XA:Z:14,+19756472,100M,2; +SRR062634.12330584 147 22 16052336 0 100M = 16052270 -166 TAGGGAGATATGAGACTTCAGTCAAATACATTTTTAAAAAATACATTGGTTTTGTCCAGAAAGCCAGAACCACTCAAAGCAGGGGTTTCCAGGTTATAAG EEGEEGGGGEEGGEEECDEAGFFGEGFGFFGGGGEGGGGFGFEGFEGGGGGGGGEGGGGGGFEGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19790574,100M,0; +SRR062634.12334727 163 22 16063459 23 100M = 16063537 178 TCAAGACAAGACGAGAATATGTTAAGATACCACAAAAACAAAAAATATACATAGTTAAGAGCTTTCATGAGTAAGCAAATTTTATACATATAAAAATAAT GGFGGFGFGFGFEEGGGGGDGGGGGFGGGGEGGGGFFGDFFGGAGGFFG?EGGFGGGGGGEDGFGFAFDAEBCEEEFEFFAEEEEECGG:EEEEEGGCCA RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19779413,100M,1; +SRR062634.12347760 163 22 16072614 0 100M = 16072713 199 AAAATCAATAACATACAAGTAAAGAAGATCAAGAAGCCACAAAAAGCAATTCCAAGTATACAATGCATGGTGAAAAGTGATGGAAGAAAACTGCAGGAAT FDFDBFBDFDDEDFDFFFDDEEDEEFFDBDDFFDFEFFFDFFFD>AFFDD=EFFAFAF?CFEFFFFFDBD@EBEDFF@EE@DB@DEFE?>EEDEEEE?=E RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19770253,100M,0; +SRR062634.12349883 147 22 16081667 12 100M = 16081546 -221 GCAAAAGGAGGAAGCTCTCTGCAGCAGACAAGGGCGTCCAAGTGGATTGCCGTTTTTACAGTTGAATCAAAAAGCTTTTATAAGAAACTCCTCTCAGCTA @EDFGEEEEFEEGEEEGEEEGFEAEGFEGGGFEGFGFDGGGEFGGG?FGGGGDFGGGGGGGGDGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:8 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.12389979 163 22 16083814 0 100M = 16083907 192 CCTCTTTCCGCAAAGCTGGAGGGAGAGTCAACAAAAGCAAAAATAGGTCAAAGTCTTCTAAAAGCCATACCTGAAGGGTTTTCAATATCACTTGATCAGA GGGGGGGGGGGGGGGGGGGEGGGFGGGEGGGGGFGFGGGGGGGGGEGEGGFFF@BGGGGGGGGGFGGEGGGFEEGEDFA?EEEEAGGFGFGEBG@BFEDA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-29992,100M,0;14,-19759116,100M,0; +SRR062634.12397906 163 22 16054169 15 100M = 16054263 194 CCCTCCAGCTTGGCCCTCACATCTGTGAAACCCACAGTAATGGGGCTGACATCCTCTGCCCTATGCAAGAGAGGTTTCCCAAGCACTTGCAGCAAGTGAG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGGEEFGGGGEEGGGGGGFGGGGEEEGGEEEEFDFEBDECEC@CCDBGEEEEADCBCCBCAD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19788741,100M,0; +SRR062634.12447434 147 22 16087881 0 6S94M = 16087774 -201 GCCTGATGTCAACCTTGAGCTTAGGTATTCACAAGAGGACTAGTGTCCAGCTGGGGCCAGATGTTCACTTGGGGCCTGGTGTCAACTTGAAGCATGGTTG ######AAA=;5888:;;>3:C?E?A@AA:CDEADCCCABCAE=CEDEEGEEEGEBFDEGGGGGGGCGGFGGED?GGGGGGGGGGGFGGFGGGGFGGGGG RG:Z:SRR062634 XC:i:94 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:94 XA:Z:GL000244.1,+25924,100M,0;14,+19755058,100M,0; +SRR062634.12482352 147 22 16063086 37 100M = 16062991 -195 AAATAGCTGGGTGTGGTGGCACTCACCTGTACTAAGGAGGCTGGGCTGGGAGGATCTCCTGAGCCTAGGTGGTTGAGGCTGCAGTGAGCCAAGATCATGC @@@4:E?EEGBBBC??ADAA5<;;,><:>>???=?CC:EB=CGEGE?GGGFDBDGGBDDEEEGEGFDGGFEGGGEGGGFDGGGFFBDFDEFBEFD?EEEE RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:68T31 +SRR062634.12505449 163 22 16062163 22 100M = 16062251 188 CCCAGGGCTACAAAGCCTTCACTGTATGAGCTGCAAATCAAGCCAAAATGTTATCAGAGATGAAGAGATGCGGATAGCCATCAGGCCTGGGGATGGACCT GFGEGFGDDGGGGGGEGGGGGGGGGFGGGGGGGFGFGEFFFGBGEEEEFGGGGGGBGEGEGDEGGGFFFEEEAECEEDFDEEEA=CDDEDDD?9CCACAA RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:70G29 XA:Z:14,-19780709,100M,1; +SRR062634.12572740 147 22 16064137 23 100M = 16064046 -191 GTTCCCTCTAAATCTCATGTTGCATTGTGTTATCCCCAATGTTGGTGGTAGGTCTGGAGGCAGGTGATTGGATCATGGGGATGGATTTCTCAGGAATGGT FBEGEEEEEEEEDGEFEEEEEGGEEGGFGFEGEGDGGGGBGFGGGGGGGGGFGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19778735,100M,1; +SRR062634.12584179 163 22 16078457 0 100M = 16078597 183 TATTATTTCCCCAAGAAAACTCTGCACTTGGCCTTGTGAGGGAAGCACTATAGTTACTAAGTGGTGGGGGTGGGAAGGTGAATTTGGAAACAACATTTGG DDDDDDDDDDDDDDDDDDCADDDDDCDDDDDDDDDDBDBDDDBCDDBBBBBBCBCBBDBBCACCACCD>>6@?>3>BBBB?ABB@B?@B?@ RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:0C56C42 XA:Z:14,-19764410,100M,2;GL000244.1,-35350,100M,2; +SRR062634.12606172 65 22 16081155 0 100M GL000244.1 32652 0 TATACAGATGCTCACATTGTTTTGTTTATATGAAAATGTTTGTAACTACTACGCGCATTTTTGCTACTTAAGCCTTTTGGTCTTGCTGCCATAGCAAATA GGGGGEGGGFGGGFGFGGGGGGGDGGGAGFGDDGGGGGGGGGFGFGGGGGGFEGGGFGGGFGEFGGEGGGDEGEFEGED?CCFC:EEEBEEEEEEDDCB@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-32651,100M,0;14,-19761775,100M,0; +SRR062634.12658957 163 22 16055031 37 76M24S = 16055171 185 TGATTAGGCCTGAGAGGTTGTTTGTCTCCCATCAGTGGAGCTGGCAATCTGACGATGCCTCTCTTTGGACCCTGGCCTTCAGCCTCTTCAGGCCATGGGG GGFFEGDBGGGGGEGDEAEEBEEEEGGGGGGGGGBDGDGAGEGGGABFFDDFEE?BAEECEBEGECBECEFCA=D@######################## RG:Z:SRR062634 XC:i:76 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:76 +SRR062634.12675072 147 22 16080023 0 100M = 16079936 -187 GCAGTGGGCCAAGATAACACCACAGCACTACAGCCTGGGTGACAGGGTGAGACTGTCTCAAAAAACAAAAACAAAAACAAAAAAACACACATTTAATTGA GBFECDDF:EEDEFEDEFDFEBGAEFCF?BGFEGGGGGGGGGDGDFGGGGGGGFGGEGCGGGGGGGGGGGGGGGGGGGGGGGGGEGGGGGGGGGGGEGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+33778,100M,0;14,+19762844,100M,0; +SRR062634.12733781 163 22 16058903 29 100M = 16058999 196 CAAGCTATAAGCCTCGGAAACTTCCATGTGGTGTTAAGCCTATGGGCATGCAGAGTGCAAGAGTGATGGCTTGGCATTCTCTGCCTATATTTCAGAGGAC GBGGGGGGGGFGGFGGEG=EGGGEGFG=EFFD?FEGGFGGFEGG?GGDBGDGDBDEGGEEEBE@FEEBGBDDDEC?AC=CAB@?C:3<(4?>5,-4>;:> RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19784004,100M,0; +SRR062634.12758522 163 22 16064266 0 100M = 16064360 194 TGATGGTGAGTGAATTCTTGCAAGATCTATTTGTTTAAAAGTGTGTGGTGTTTCCCCTGCCTTGCCAGGTGATGTGCCTGCCTGCTCCTTCTTCATCTTC GGFGFGDDGBEGGGGEGGGGGGGFFGGGDGGGFGGGGGGG?CEEEFFFDFFFGGGGG=DGGEB?EFFEEABDDB?EBDC=CB5=C:BB:=CC?CC=CC:B RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19778606,100M,0; +SRR062634.12768799 163 22 16055306 23 100M = 16055410 204 TCACCTTGTCCTTCAGGCTCAGCAGCAGGCAACCCCAATTGCAGAAGTCTCTATCCATCTTGGGAGGTTGTCATCCCCTGAACCTTCCTAGAGCACATTG GGGGGGGGGGGGGGGGGGGGGGGEGGGGGGGFGGGGGGGGGGFGGGGEFGGGGGGGGGGGGGGGEGGEFFEFBFEEGEGEEEGBEEDE?EEEEE=ECCD@ RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19787603,100M,1; +SRR062634.12775255 163 22 16051427 23 100M = 16051498 171 TATGTAAACCAAAAATAAAATTCTGAAGATGTCCTGATCATCTGAATGGACCCTTCCTCTGGACCAGGGCATTCCAAAGTTAACCTGAAAATTGGTTTGG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGDGGGDEFGGGGGGGGGGEGGEGGGGEGEFGGGGEGGGB RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19791475,100M,1; +SRR062634.12777147 163 22 16059766 0 100M = 16059838 172 GACCTAATGGGAGGCATTTGGGCCAGGGGGACAGATACCTCATGGCTTGGTGCTATCCTTGCAATAGTGAGTGACTACTCATGAGATCTTGCTGTTTAAG DDDDDAEEEEEFFFFGFFGFGFGGGGDGGFGGGGGGDGGBGGGGGFGGEFEGFFEFFGGGFGGGGDEEEFFBGAFEEEGGEEEBGEEEFFEEECCCFDDD RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:74G25 XA:Z:14,-19783141,100M,1; +SRR062634.12857943 147 22 16081406 0 100M = 16081339 -167 CTAATGTGGCTGCTGAAGATGTTACCTAGAGCAGAGGACTTTGGGTTCAACCTAAGCACTTTATATCCTTTATTTTCAATTGGGTAGGAGATAAAATAAT ?=C?@GEEEEFGEBFEEEDECECACD:EEEEECECFDAEFEDFGGGGGEDFFFFEFF@FFGDAGGFGGGGGGGFGFGGFGGGGGGFGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:6 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.12975192 147 22 16086550 15 7S93M = 16086452 -191 AAGGCCTGATGTGGGCCTGGGGCCTGATTGCCCACCCGGGGACTGGGTGTCCACCTGGGGTCTGATGTCCACCTGAAGTTAGGTATCTACCTAAGGCTTG #######A?AAAAB:=@BDDEE=EEBABB8?CE:?==>=:@>(=:.=>=+@?6)? RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:87T12 XA:Z:14,-19787241,100M,1; +SRR062634.13334738 147 22 16065703 34 100M = 16065643 -160 ATTAATTGAGACGTGCATTATATCATTTCTTACTTTTCAAGGAAAAACATGTTTTACGACTTGAGTTTATCAGTCTAGTGACCTTGCAGCTGCACAGCTA FEDFGDEFEED=FDGEFEGCFGGEGGGGGGFCFFEEDGGFGDGGGDGFFGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGEGGFGG RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:16 AM:i:16 X0:i:1 X1:i:5 XM:i:1 XO:i:0 XG:i:0 MD:Z:71T28 +SRR062634.13381370 147 22 16053893 15 100M = 16053804 -189 CATGGACTGGCTCAGGGGTCTGTCAGATGATTCTAGAGCTGAGTTGGAGGTAAGGGCAGGGGGTTTGTTCCTGGGTTCAAGACCATGGAAGGAAGGGGTA AA=DDC=DECEEB@GEGEEEBFGEEGDFGGEEEEEEEGFFGEGEFGGDGGFEDFGGGGGFGEGGGGGGGFGGGFGGGGGGGDGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19789017,100M,0; +SRR062634.13567718 163 22 16052362 15 100M = 16052473 211 TACATTTTTAAAAAATACATTGGTTTTGTCCAGAAAGCCAGAACCACTCAAAGCAGGGGTTTCCAGGTTATAAGTAGATTTAAAATTTTTCTGATTGACA GGGDGGGGGGGGGGGEGEGGGGGEGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGCF?FDDEGGGGGFCGDGFBFDGGGFDFEFFGEGGGGGGGFEEDEGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19790548,100M,0; +SRR062634.13630324 147 22 16070963 15 100M = 16070895 -168 TCATCTCAGGACAGGCAGAGGGATGATTTCTGGTCCTGACCTTGTCCCAGAACTATGAAAATAACCTGGTAATTGAGACCGTCAGGGTGAAATTGAGCCC DDDDBBFBEGEEEGGEEGEGGGGFFEBGGGGFEFFGGDEEFGGFFGFGGGGGGGGGGGGGGGGEGGGGGGGGGGGEGEGGGGGGGGGGGGGGGGGGGEGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19771908,100M,0; +SRR062634.13661630 147 22 16075418 0 100M = 16075295 -223 CCTGAGGCCTTCCTCCCTTCCATCAGGTGACCCTACCATGTGGCCTCAGCTCTAGGGAGATGGGCCCTAGCTGGAGGCACTGCACAGCAGCATCCTGGGT A<5A::6;:9'CC)CCCA:=A::>>>=>;55A?>:-BEC@@5@>B=<:.?6?=::3@DB?CCC@>49@@5D?DDEDDFEFEGFE RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:3 XM:i:1 XO:i:0 XG:i:0 MD:Z:77A7 +SRR062634.13878407 147 22 16051373 0 100M = 16051315 -158 GAAATCTGGGCAAGAGAAGAAAAAAAGGTCAGAGTTTAATCCTCATTCCTAAGTTATGTAAACCAAAAATAAAATTCTGAAGATGTCCTGATCATCTGAA =@=>5@@A=CDDDBF@DD;CEEEGEEEEEGEEEDBEEGFEGFFGGGEFFGGGGFFGGFGGGGGGGGGGGEGGGGGFGFGGGGGGGFEFGGEGGFGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19791529,100M,0; +SRR062634.1391342 163 22 16089957 0 72M28S = 16090062 205 TCTGTGCCCCACTTAAACACCGGGTCTTAGGTAGACATCAGGCCTCAAATGGATGCCCAGGCCCCAGGCGGATATAAGGCCTCAGGCAAACAGCAGGCCC CBC???4>@4<7=BBDDD??<,5C?C-?BBC85CA0;16=ABB5BCDCD:==@@-5>4>=:-=A############################ RG:Z:SRR062634 XC:i:72 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:4 XM:i:1 XO:i:0 XG:i:0 MD:Z:68T3 +SRR062634.13925390 147 22 16063124 60 100M = 16063047 -177 GGCTGGGCTGGGAGGATCTCCTGAGCCTAGGTGGTTGAGGCTGCAGTGAGCCAAGATCATGCCTCTGCACTTGAGCTTGGGTGACAGAGTAAGACCCTAT BA-5BDGEBFFFDBDB=EDD@EEEEE:?=FEFGGGFAEGF?BEE=GFFDF?EFFGGGGGGGDFGGGGGFGGDGGGEGGGGFGGGGEFFFFFFFFFDGGGG RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:30T69 +SRR062634.14050328 147 22 16064082 23 100M = 16064025 -157 TCTTATGGAAGTCATATAAGAGTCACAACATGAGGCATGAGATACCTCAGATGTTGTTCCCTCTAAATCTCATGTTGCATTGTGTTATCCCCAATGTTGG EEEFDGFEEEEEDCEABCB5FGEGFGFEEGGGEFGEFGGFGGGEEGEGGGGGGGGGGEGFGGGGGGDFGGGGGFGGGGFFGGGGGGGEGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:2 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:2 XO:i:0 XG:i:0 MD:Z:5C13A80 XA:Z:14,+19778790,100M,3; +SRR062634.14065678 163 22 16092966 0 100M = 16092979 113 CAGGTGGATATCAGATTCCAGATGAACATTAGGCCCCAAATGGATACCTAGGCACCAGGTAGACATCAGACCCCAGGTGCATCCCCCAGTCTCAGGTGCA FFFFDGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGD:DDBACCCCFGGEEEFDFED?EEBDDCDD?CDCDEDDC:@C=BCAADBCCC4A>B=B?C>-B=A?@???@8+??)?#################### RG:Z:SRR062634 XC:i:80 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:80 XA:Z:14,-19776571,100M,0; +SRR062634.14233499 163 22 16076047 0 100M = 16076123 176 CCTACCCTGCAGGGAGTGCCTCTCCTCGGCCATCAGCTGATCCAGTGCCCAGAAGGTGTCTTCCTCTGGCAGATACAGGAGGAGGATGGCAGTTAGGCAG FFFFFEEFFEGGGGGGEFGGGEGGGEFEFEE:BEEDD:AECCCB??@@?C=????>,@>>CCCC:CCCCCCB@:C>->>3>@/==2???5>>=>8AB:AC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.14238941 163 22 16073434 17 55M45S = 16073521 187 CAGGGCTCACACTCCTGACACCCAGACTCCACGATACTCCTCTGTCCCCACCCAGGGCAGATCCCCGCCCCAAAAATTTTCCCCCTTTTTTCCCGCAAAA GGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGFFCCCA############################################# RG:Z:SRR062634 XC:i:55 XT:A:U NM:i:0 SM:i:17 AM:i:0 X0:i:1 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:55 +SRR062634.14242393 163 22 16064133 23 100M = 16064250 216 TGTTGTTCCCTCTAAATCTCATGTTGCATTGTGTTATCCCCAATGTTGGTGGTAGGTCTGGAGGCAGGTGATTGGATCATGGGGATGGATTTCTCAGGAA GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGFGGFGGEGGGGGGGEGGEEGGEGGGDEGGFEEE4E?=AAC=CDCDEEEE*=BA@ABDCCCCCC:C RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19778739,100M,1; +SRR062634.14250468 147 22 16058499 15 100M = 16058465 -134 AATGACTTAAAGCTGGAATTTATGTTTAAACAGGAAGCAGACCAGAAAAATTTGGAAAATTTGCAGTCTGGCCATGTGACAAAGAAAGAAAAAAGCTTTT DFDECEDBGEGEEEEGEGGDGGGFEFEGGDEFEGGEBGGGGF>FGGGGGGGFGGGGGGGGFGFFGGGGGGGGGGGFGGGGGGGGGGGGGGGGFGGGGGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19784408,100M,0; +SRR062634.14321084 163 22 16073360 23 97M3S = 16073465 204 CTTCCCTCTAAAGGAAGCTTGTGCAGGAAACATCATGCATGGAGCAGCATGGGGACAGGGGCTGGCCAGCTGGGCAGGGCTCACACTCCTGACACCCAGA GFGGGGGGGFGGBGGGGGGGAAGGGFFFFGGGEGGGGGGCGGGBGGFDFCADEB?BDDCCC:@CCAB?CC3>>64=??@?C### RG:Z:SRR062634 XC:i:97 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:97 XA:Z:14,-19769510,100M,1; +SRR062634.14347911 163 22 16088289 0 100M = 16088410 221 AACTTGGGACCTAAGTATTTACCTAGGGCCTGGGTGTCCACCTGGGGCCTGACTTCCAACTGGATCTTGTGTCAACATGGGGCCTGATGTCCACTTTGGG GGGGGGGGGGGGGGGGGGGGGGGGGEEFGGDGGDEGGGGGGGGAGFFGFGDFFFDFFEAFBEBEEGFEGEGEFEGFFEFEEFFFEBFBEBEDEFDEBBED RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:0 AM:i:0 X0:i:1 X1:i:6 XM:i:1 XO:i:0 XG:i:0 MD:Z:7A92 +SRR062634.14482419 163 22 16081594 20 100M = 16081696 202 AGTTCTTCATCTCACCACCAAGATGATTAAGGAACGGGGATACACGGGTGAGGTGGGAGTGAAAGTTTAATAAGCAAAAGGAGGAAGCTCTCTGCAGCAG FGEGEGGGFGGGGGGGGGGGGGGGGFGGGGGDGGGGBGF:GEGGGEDDAECEEBEEEBEBECEEDDDDFFGFGGGAEDEDFDFDBFG=DEE?DDDCGEBF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:20 AM:i:0 X0:i:1 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-32212,100M,1;14,-19761336,100M,1; +SRR062634.14485327 163 22 16092077 0 100M = 16092163 186 TAGCAAGTCCCCAGGTGGACATCAGGTCAAAAGTGAATATAAGTCTCTAGGAAGACATCTGGCCCCAGGTGGATACTGAACTAGAGGTTTACATCGGGCC GGGGGGGGGGGGGGGEGGGGGGGGGFFGGGGGGEGGGGGGGGGGGGGGGGGGGFGGGGGGFGGGGGBGFEGFEEFGDEEEGDEBEFGEEBCCCDEEGBEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.14568549 147 22 16054726 29 100M = 16054606 -220 CAGGAACAAGGCCCGGGGATGCTGCTTCTGGGTCTGAACCCCAGGTATGGGTCAAGAGGGCTCCAACTTCAAGGGTGAGGTTATTCATGGTAGAAGAAAC C;BAA=C?CDECEBEEEEEBGEEFEEEGBFEEEEEFEEFEEGEGGFGGGEGGEGFGGGGFGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:14A85 XA:Z:14,+19788183,100M,1; +SRR062634.14575835 163 22 16078586 37 99M1S = 16078704 211 AGGAATTCCATACACATCAGATCTGAAAATATGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTACATGTGTCTGTCTGTGTGTGTTTGTGTGTGTGT GGGGGGGGGGGGGGGGGGGGGGGGGGEGGGGGGEGEFEFFFEFDFDFEFDECECEDFEFCECECB?:-??BBBB5@?@-@?ABACCAC5CAAAAA?@?@- RG:Z:SRR062634 XC:i:99 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 +SRR062634.14579669 147 22 16086312 0 100M = 16086211 -201 GCCGAATGTCCCCCTCAGGGTGAATTCCACCTCAGGCCTGTATGTCCACCTGGGGCCTGATGTCTGCCTTAGATCTGTGTCCCACTGGGGCCTTGTGTTC CC?C:-AB:@ACCAC=?DCBEDCC@AEEEEDDDB?GEEGEEEDFEDABEDDGFFFFEGGFFFGGFFEFFGGEFDGGGEGEGFFEGGGGGDGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:6 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.14627816 147 22 16097683 15 19S81M = 16097564 -200 CGCGAGCCAGGAACGCCGGGGTCCCCTGGGCCCAGCGCAAGGGCTAATGGGAAGACACTTTCGTCCGTTGGGGACCCAGGCTCCGCTTCTCCGTGGTGCG ###################@=+==?:-BAA@?=-AB?BB=CCCCBA?:EE:CAAAD=ADA=A?AB:E?:EDDD5DD=:D?DCD=B:EDD?EBE=EDBGBG RG:Z:SRR062634 XC:i:81 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:81 XA:Z:GL000244.1,+16131,100M,0;14,+19745569,100M,0;14,-19829759,100M,0; +SRR062634.14683012 163 22 16059412 0 100M = 16059505 193 CATGTGGGGCATGTAACCTCTTTCTTTTGGCCAAGTTACCCCTTTGGGGATGGGAATGCCTGTACCACAATGCCTGTACCCCAATTGTGTCTTGAAGTAA FDFFEGGGGGGGGGGGGGBGDGGGGGGG:GGGGFEEGGGFFGFDGADFF=DFFDDEGFFEBDCDEE?EEE@EEEEB@DDBEECCCC?AB?AC?5:C@=:C RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783495,100M,0; +SRR062634.14697022 147 22 16059621 0 100M = 16059544 -177 TGGAGGACTATTGGGAAGGCATGATGGTATTTTGCAATGTGAGAAAAACATGAGATTTGGGGCTTGTCCTGGGGTGGAATGATATAGTTTGGATGTTTGT DEGEEDCDEDEEFEGFFEDGDGFGGEGGGFGEGGGEGGGGGGEGGGGFGGGGGGGEEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19783286,100M,0; +SRR062634.14706235 163 22 16096260 0 74M26S = 16096428 214 TTTCCTCTGCCCCAGGTGGGTCACCTTCCCCTCATGGGCCTTCTGCCTGCTTTGGGGTACCCCTAGCGGCCCCAGGCAGGGGCTTGGCTCCAATAAGGGG GGAGBGGFFGGFGDGACEDDFFFFDE?EEEE=EBBCDCEEEFF?EBC=BC:BAA-:=??@BECAAAB?AB@???########################## RG:Z:SRR062634 XC:i:74 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:47C26 XA:Z:GL000244.1,-17565,100M,1;hs37d5,+18470739,100M,1;hs37d5,-5437444,100M,1; +SRR062634.14729112 147 22 16067672 41 100M = 16067587 -185 AAAGGGAGACAGAGAGAAAAACGGTATGTATTTCTCTTTGAGAAAGTCCACTGGCACCAGTAAGGTTTTTGAGAGCTGGCAGGCTCTGATAGGTGAGGGA CC@CCCC@A:EEEEECEEECBB=CDEEEEEEFEEEEEEEEGFGEGGEGGEFEGGGEFEGGGGGGGGGGGFGGGGGFGGGGGGFEGGGFGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:18 AM:i:18 X0:i:1 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19775200,100M,1;21,+10631656,100M,1;9,+68332996,100M,1; +SRR062634.14785697 147 22 16065718 17 100M = 16065674 -144 CATTATATCATTTCTTACTTTTCAAGGAAAAACATGTTTTACGACTTGAGTTTATCAGTCTAGTGACCTTGCAGCTGCACAGCTAGAGAAATAGAGTCTA DEEEGF?GGGGGEGGECFDEE-FGGEEEFGGGFFFGGFEGEGGEGDGGGGGFGFGGGGGGFGGGGBGGGGGGGGEGGGGGGGGEGGFGEGGDGEGGDGGG RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:56T42T0 XA:Z:14,+19777154,100M,2; +SRR062634.14801948 163 22 16062504 0 100M = 16062575 171 ACAAAATTAGTTTAAGAACACATATAATTAAACAAAATAAAAATATTATCTCATCTTTTAAAGAATATTTACTAATTCACAGATATTACCCGAAGTTTAG DGGGGFGGGGDGGFGGGDGGGGGGFFGGGGGGGEGGGGFGEGGGDEGEGGEGGGGGEAEGGGGGGEFGFFGGGGGGBGGGDGGGAGEGGEFEDEDBBBFD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19780368,100M,0; +SRR062634.14825518 163 22 16055745 60 90M10S = 16055841 196 GGGCCCTGGGTAGAAGACCCCTCCCTGCAATTCCCACACCCCTGGAAGCTGTAAAGCCCAGCTCAGACCCTTGGAGGGCAAGTGGCCAACTCAGCACACG GGFFGGGGGGFFGDGGGGGGGGGGGEDG=FGDGFEDGGGGEGE=AAB=DCAFEAFEBFB?CECEBDDEEF?5?5?C99BB?5>>5>5:=@########## RG:Z:SRR062634 XC:i:90 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 XA:Z:14,-19787174,100M,1; +SRR062634.1484613 163 22 16051435 23 100M = 16051518 183 CCAAAAATAAAATTCTGAAGATGTCCTGATCATCTGAATGGACCCTTCCTCTGGACCAGGGCATTCCAAAGTTAACCTGAAAATTGGTTTGGGCCATGAT GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGCGGFGEGEGGGGFEGGGGEGGFEAEEEGGGGEEEEEEEEEEEB@EEEEED=CCC RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19791467,100M,1; +SRR062634.14985224 163 22 16098700 0 100M = 16098805 205 TCAGGAGAATTGCTTGAGTCTGGGAAACCGAGGTTGCAGTGAGTCAAGAGATCACAGCACTGCACTCCAGCCTGGACAAAAGAGTAAGAGTCCATCTTAA GGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGEFGGDGEGGGEGGGGGBGBGGGGGGFGGGGGGEGGGGFEEGGEEBFEEBCBCDFEFE@DDE?FEBCD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19830794,100M,0;GL000244.1,-15100,100M,0;14,-19744522,100M,0; +SRR062634.15045508 163 22 16059789 0 93M7S = 16059866 177 CAGGGGGACAGATACCTCATGGCTTGGTGCTATCCTTGCAATAGTGAGTGACTACTCATGAGATCTTGCTGTTTAAGGATGTTGCACCTCCCCCCCCACT FDGGDGGGGFEGGGGGGGGFDDGAGEGBGAFDFFFEEGGDFEAEBCADCDEEGEFEGBAEEG=ABEGEE?F@CEDFEED=A=AA?EDE@EDFC####### RG:Z:SRR062634 XC:i:93 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:51G41 XA:Z:14,-19783125,100M,1; +SRR062634.15136938 147 22 16056549 15 100M = 16056476 -173 TCATCTCTGACCATTTCTTCAGAGTCGACTATGAGCAGAAATCACCCAGACACCCCACAGGTAAGTCCACTTGCCTCACATGGCTCTTCCTCTTGTCAAG BFEDFDFDFCEEEBCBFBCFGFEDEFFGEEDFGGGGGGGGFGGEGGGGFFGCGGFGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19786360,100M,0; +SRR062634.15188632 147 22 16060807 22 100M = 16060773 -134 CAGGCTAGTCTTGAACTCCTGACCTCATGATTCACCGGCTTTGGCCTCCCAAAGTGCTGGGATTACAGCCTCTTTTCTTATACATTACCCAATCTCATGG @=C?CDCBDD????0C5@AAD2EDDDCDDDECFABEBFEDFGFGGEEEGGEEGGFEEGGGFFGFGDGGFGGGGGGGGEGGGGGDGGDGGGFGGEGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19782065,100M,0; +SRR062634.15382701 147 22 16074523 0 100M = 16074444 -179 GTTACCTCTCCCTGGGTCCCTCCCGTGACACATGAGGCTTATGGAAATTACAGTTCAAGATGACATTTGGGTGGCAACACAGCCAAACCATATCACCTTA 5@5(C>7>;:5--C>=0:B4+66>>:??47D5:CDD5FBFD?EEEE?DAC=@=?@?-C@CD5CBDE5A6DDD=DB:A=C=CAC=A:C::@DE?ABFFFED RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:10A89 XA:Z:14,+19768344,100M,1;GL000244.1,+39284,100M,1; +SRR062634.15395619 147 22 16055993 15 5S95M = 16055902 -186 CTGCTCCGCAGCAAAGGACTCTCTCCAACAGGCCCTATGATACAGCCACTCAGCCAGTTCACCCCACCCACATTCCACAGGCTGCTTTAGGCTTTAGGAC #####?:,AAABCE@?=>@=?5@,@EEDDB??B?@:5?BCCCEDDDFEAEGEG?GGAGFG=CEEFCE?ECEFF5FGGFGDDDDFFGGFFGEGGGGGGGFG RG:Z:SRR062634 XC:i:95 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:2T92 XA:Z:14,+19786921,100M,1; +SRR062634.15413978 163 22 16052970 0 81M19S = 16053078 202 TGGCTGACTTGAGATACCTGCCAGATGAGTCACTGGCAAAAGGTGCTGCTCCCTGGTGAGGGAGAAACACCAGGGGCTGGGAGAGGCCCAGAAGGCCCTG GGGGGGGGGGGGEGGGGGGGGGGGGGFDGGGFGGGEGGGGGEEEGEGGGFFGGEGDCFCFDD?DCBEFEDEEEEEB@?EE@################### RG:Z:SRR062634 XC:i:81 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:0G80 XA:Z:14,-19789959,100M,1; +SRR062634.15491942 163 22 16056252 23 100M = 16056385 223 CTCAGGGGGACTGGCAGAGTGGGCGCATCCTGAAGGGCTGTGCTGGAATCAGCCCATAGTCTTTGCGTGTCCCATGGGGCAGGGCCAAAGGAAGCAGTCA EDBFFFEF=DFA=FBEDEEBFFD=FBBDFFEFFDEFA=DF>=C@CEEEBE=5DBC:?@@BFDAEECD9C@:C?@?EE@ED-???@B=@B@EEE:=B:?BB RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19786657,100M,1; +SRR062634.15566617 147 22 16094844 0 100M = 16094766 -178 AAGTACCTAAGTGAAACAAAGCCTTAGGTGATTACCAGGCCATAGGTAGTCATTAGTCTCCAGCTGGACAATAGTCCCTAGGTGGATACCTAGGCCCCAG DCDD=EBBFEFEEEEDEECC:CBBB@CCDC:FF@FFEFGFDEFDDCGGFFGEGGGGGGGGGGFGGGFGFGGGGGGEGFGGGGGGGFGEGGGGGFFGGDEG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:6 X1:i:6 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.15670949 147 22 16070245 0 100M = 16070139 -206 GAATTTAAAAATAATGAACCACACTTTTGAGAAATATAGGATTATGTAAAGGAACAAAATCTATGACTTACTATTGGCATTTTTTTTTGAGACAGAATCT EEG?DEEGGGGEGGFEGDFGDFAFDFEGGGGGGGEGFGGGGDFGGGGGGGGGGGGGGGGGGGGGGEGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19772626,100M,0; +SRR062634.15736242 147 22 16086293 0 100M = 16086204 -189 GGCCTGATGCGCACCTGGGGCCGAATGTCCCCCTCAGGGTGAATTCCACCTCAGGCCTGTATGTCCACCTGGGGCCTGATGTCTGCCTTAGATCTGTGTC @CC?CCAC?AAA;A?@FEDEFEGEFDEEEEFEEEEEGEFGEEEBE?GBEDDGFGGGGGGGEGGGGGEFGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:5 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.15761612 163 22 16069444 23 100M = 16069544 200 TAGGTTTCTCCAGCAGCTCTGCCCCTGCTGAAGACTAGTACACATGGGGACAGGGTGCTTTTGGCACTTCGTCACTTCTGCACCTGTCGAGGACAAGTAC GGGFGGGGGGGGGGGGGGGGFGGGGDGGGGGGFFGGGFEGFFGCGGFAGEG=EFFEGGEGGFGBEEGEGFBEABEBEE?DDEDFCDCDBCCB@C@BBC@BBEAE:GGGGGGGGEGFGGGGGGEG RG:Z:SRR062634 XC:i:63 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:6 XM:i:0 XO:i:0 XG:i:0 MD:Z:63 +SRR062634.15865584 147 22 16091388 0 3S97M = 16091277 -208 CATCATGATTGTAAGTGTCCCGAGCATGCCCAGTCATGCTCCCTATTAAGCCTGAAGAACTGTGGGTCAGTTAAACCTCTTTTCCTCATAGGTTGCCCAA ###B=???5:B@E?E:B@ABBF@EEB?EDEADADEDEF?DBEEGGFFGEGAFAFCGGBFGDGGGGDDEGDFFFBFFF?BFBDGFGBFFGFFGGGGGEGDD RG:Z:SRR062634 XC:i:97 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:97 XA:Z:14,+19751548,100M,0;14,-19824177,100M,0;GL000244.1,+22414,100M,0; +SRR062634.1592431 147 22 16093950 0 100M = 16093857 -193 ACTAGAGGTTTACACCAGGCCCCAGGTTGACACCCAGGCTCAGGTGGACATTGGGCCCCAGGTGGATACCTAGGCCCCTAGTAAACCTCAGATTCTAGGT BCCBC=EEEE@EEEED?CDC=CDCDDC@DAGECEEEEEEEEEGEFEFEFGGGEGGEEDGGGGDFFGGEGGGGGGGGGGGGGGGGEGFGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:5 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.15969153 163 22 16059583 0 100M = 16059617 123 TGGGACTTTTGAGTTAATGCTGAAATGAGTTAAGGCTTTGGAGGACTATTGGGAAGGCATGATGGTATTTTGCAATGTGAGAAAAACATGAGATTTGGGG GGGGAGGGGGGBGGG=EEEEGGGBDGGGGEGGGGGGGGGGFCGDEGGEFGGGG?EFEBGFEEGGGCBEDFGGE:=GGEG?FAAAC???>ACB=?B:-?-4,:<85:@5@############################################# RG:Z:SRR062634 XC:i:55 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:6 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:55 +SRR062634.16287193 147 22 16052110 60 100M = 16052017 -193 ACCTGAGATTTTGCCATTGCACTCCAGCCTGGGCAACAGAGCAAGACTCTATCTCAAAAAACAAACAAACAAACAAACAAACAAACAAACTGTCAAAATC ?=:<@??5@?=:2?:;;5-<3:1>;,@B@@?;>>:?5A5ABEBEFAAA==:EFECC?>:FG>EGGGGGG@GGGGGGGGGGGFEGGDGGGFGFEGGGGGFG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.16316709 163 22 16062868 0 94M6S = 16062949 178 GGGAGAGGAGGCTCCTCATCCTGCAGACAAGAAGCAGCACCCACTGGTTCCTGCTCCAAAAGCCATTATGAACATTATAAACTGGCCAGGGGCGGTGGCT GGGFFGGGGGGEGGGGGGGGFGGFGGGGGGDDGGGGGGFGGGGFEGDEEGFCFE=EFDDDFEFEDEEEAAEEFEEEEEEEFEE@B@B?AB0???###### RG:Z:SRR062634 XC:i:94 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:20T69T3 XA:Z:14,-19780010,100M,2; +SRR062634.16336737 147 22 16072277 0 100M = 16072186 -191 TGCATGTTTATTACTGCATTGTTTGAACTTAAATGTCCACAACAGGATAACAGGATAATTAAACAGTTAATTGCACTATGGCACATCCATCTGTTAACAA C?EDDAA@CECCC@?C@AC? RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.16358075 163 22 16080903 0 100M = 16080957 154 TTTTTTAATCTACCGCACTCTAAGAACACATATTTTAAATCAATTTGTACATTCAGTGCCTAGAACAAAACCAGCATTTTGTAGATTCCAAAGAATTATT GGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGFGGGGGEGGGEGEEFFGEEGGGGEGFEEDEEEEEGEEGEGEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.16389448 163 22 16097253 20 86M14S = 16097354 196 GAGGCGCACCCTGGTTTCGAGCCAGGGACGCTAGGGTCTCTGGGGCCCAGTGTAGGGCTGATGGGTAGGGACGTTGGTCCGTGGGGGACCCAGGGGGCAC GGGGFGEGGGGGGGDFGGGFGGGGGGGAFGGGB?AFCFEEEDE:E5@=5?ACA=BCCA?CACC=C6?@??:?@?35>>3<;=@# RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:GL000244.1,-28238,100M,0;14,-19757372,100M,1; +SRR062634.16769670 117 22 16081625 0 * = 16081625 0 GGTGAGATGATTGCTGGGATTACAGGCGTGAGCCACCGCGCCTGGCCGTATGTTTATTCTTATGATAGTACCATACTGTTTTGTAGTATGTTTTATAGCT CCCCBCC:CA:CDCBDCBCCA?C@ABBCCBADFECFGGGGGGGGEGDGGGGGGGGFDGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 +SRR062634.16950120 163 22 16080662 0 60M40S = 16080799 203 TCGGCCCACTGCCGCCTCCACCTCCCGGGCTCAAGCGATTGTCCTGCCTCTGTCTCGTGAGTAGCTGGGACTACAGGCATGTGCCACCACCCCCGGCCAA C??A:A?A?<5=>==C:DDD?DDAA5@CC'+;8C?:A;:?@=C:::A?>C=>-@######################################## RG:Z:SRR062634 XC:i:60 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:29T26C3 XA:Z:GL000244.1,-33184,100M,2;14,-19762308,100M,2; +SRR062634.17022019 147 22 16079228 0 100M = 16079138 -190 ATAAATTAAGACGAAGTCCGAGGAGTTATGGTTCAATTTATAATTTAACAATATGGATAGCACTTATGTTTTGATACGAGCAGCTTCGTTGGAGTGAGTA ?6@@?>5@=5,>(<;==ECAA??,?@,@@@A:<:@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-38665,100M,0;14,-19767725,100M,0; +SRR062634.17459387 147 22 16093291 0 100M = 16093198 -193 CCAGCCCAGCTGAACATCAGGCTGCAGGTGGATGCCTAGGTTCCAGGTGCACAACAGGTCACAGTTGGACATTCAGCCCCAGGTGAACATCAGGCCATGG A,,:?C?=BC:B:AA;8CAE?ACDCAEEDDE:CEACC8BEGEEEEEGGGGFEGBGFFDF5EFGGFGGGGFGD?GGGFBFEGGEFFGG=GGEGGFFGFGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:5 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.1749702 163 22 16081352 0 80M20S = 16081446 194 CACCCTGTCAGAAGCAAAGAGACATCTACTCCACCATTACCACCCATCCCTCTACTAATGTGGCTGCTGAAGATGTTACCTAGAGCAGAGGACTTTGGGG GGGGGGGGGGGGGGGGEGGGGGGGFGGGGGGGGGGFGGGGGGGGGGEGFGEGEDEEEEEEEEEGEB?4:4=;-5==>>:=:5=;5;*83:=??@7:21+88748:@CBBC:CCC<<>17,>C>>@@C??>>A???A:GGGGGGGEGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 +SRR062634.17705940 147 22 16095914 0 100M = 16095814 -200 CCGACTTGGGGTATCCCTAGCCACCCTAGGCCCACACTGGACTCGAAACAGGGACGCCAGGGTCCCTGGGGCTCAACGCAAGGGTTCATGGGAATACACT @?@:@=+>>>A?:)BBBB:?BB5CB?=-?=?BC?EBC=?BEBC:EDEA?EGGGFGGGGDGGGFFGGGFGGGFEFGFGGEGGGGGGFGGFGGGGGFGGFGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:31G68 XA:Z:GL000244.1,+17885,100M,1;14,+19747474,100M,2; +SRR062634.17707542 147 22 16070980 15 100M = 16070912 -168 GAGGGATGATTTCTGGTCCTGACCTTGTCCCAGAACTATGAAAATAACCTGGTAATTGAGACCGTCAGGGTGAAATTGAGCCCACTCACTTTTAGGGCTA EEEGFEGGEE?EFGFDE?FGGEDGEFFEGBGGGGEGGGEEGFGFFFEGGGGGGGGGGGGGEGGGGGEGGGGGGGGGGGGFGGGDGGGEGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19771891,100M,0; +SRR062634.17753450 147 22 16061525 0 100M = 16061440 -185 GTTGGCCTTCACATACATTCCATAAAATGATCTTTGGCCAATAGTGGGGACTGAGGCCTCAACTCTGACCTTTTATTCTCTTGCCAAAAATCCTGTCTAA DA@B@>=:5@=B=@8?@A@@EDDFFB?BB@@EBDEEBD=@B@DEEEEDF@DB?CBBBEFFDEFBFFF@AAEFFEFFFDEDFFDEFFFFFFF=DFFFFFEF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19781347,100M,0;2,+132547828,100M,0; +SRR062634.17763088 163 22 16086911 0 100M = 16086968 156 CCGGGTGTTTACATAGGACCTGATGTCCAGCTGGTGCCTAGGTGCCCACCGGGGGCCTTGTGTTAACCTGGGGACTGGTATCCAGCTTGGTCCAAATGAA GFGGFDGGGGGGGGGGGGGGGGGGGGGGGGGGGBEGDGGGGFEGGGFGGFEFBDEDFEEGDFEDACADDC?CC<=?@@7@BABABC=DCC7@A,@AACCC RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:3 XM:i:2 XO:i:0 XG:i:0 MD:Z:87C5T6 +SRR062634.17792264 163 22 16067723 0 100M = 16067829 206 TGGCACCAGTAAGGTTTTTGAGAGCTGGCAGGCTCTGATAGGTGAGGGACAGTGGTGGGTAAAACTAGTCTTAGTGTTGCAAGAGGTTGTTTCAGTAACC GGGGGGFGFGEGGGG?EGGGGGFGGGF=G?GDFGBGGFGGGGEEFGGFEGFFEFFDFFF?EEEECEEC:GEEAA?CBCBCCCB=DE=BCBACBCBEEEEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:2,-132541641,100M,0;14,-19775149,100M,0; +SRR062634.17821762 163 22 16065989 20 83M17S = 16066088 199 TTGGACCAGGGCCAGGGCACCACTGCCCTGTACCACCCCCAAGAGGGTGTTCCTTGCATCCTGGCCCTTCTAGCTCCAGCCTTGGTTCAAAGGGCTCCAG EEB?=D=DDC=EEE???AEEEEEA=CDDB=AABBBADBA;1>A===>=C=EBEBBD:@EEF:=DE??5?=??==?@@?,=@??################# RG:Z:SRR062634 XC:i:83 XT:A:U NM:i:1 SM:i:20 AM:i:0 X0:i:1 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:41G41 XA:Z:2,-132543389,100M,2;14,-19776900,100M,2; +SRR062634.17902563 163 22 16096119 20 60M40S = 16096233 193 GGGGTACCCCTAGCGGCCCGAGGCACACCCTGGGCTCGAAACAGGGATGCCAGGGTCCCCGGGGCACAGCGCAAGGGCTGGTGGGGAGGCAGTTTCTTCC GGGGEGGGGGGGFFD=CEE=EA:?:FAFFF@BFFBCDDCDDCC@C:@@ACEC@=C?CC?B######################################## RG:Z:SRR062634 XC:i:60 XT:A:U NM:i:0 SM:i:20 AM:i:0 X0:i:1 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:60 XA:Z:14,-19747309,100M,1;GL000244.1,-17720,100M,1; +SRR062634.17958337 147 22 16079279 0 100M = 16079167 -212 TATGGATAGCACTTATGTTTTGATACGAGCAGCTTCGTTGGAGTGAGTAGGTGAGAAACCATATTGGAGTGTTTGCTGAAATGAAAGTCTGATCTTTGGA 9A:CCAAAA=>5>=4??B?=CEC@>:D@CCABCCC:EABGEEGGDBAFEBF?GGGGGEGGGGGGGGFGGGGGG=FGGGGGGGGGGGGGFGGGGFGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+34518,100M,0;14,+19763582,100M,0; +SRR062634.17969132 133 22 16096818 0 35M65S = 16096818 0 TGGGCCCCGTGGACCCTGGCGACCCCCGGGGGAGGCCCGGGGGGGCCCCCTGGCCCCCAAGGGGGGGCCCCAACGGGGAGAAGGGTCCCTAGGGGGGGGG GGGGGGGGGFGGGGGGEEGEAA<@A########################################################################### RG:Z:SRR062634 XC:i:35 +SRR062634.18017266 147 22 16066333 0 49S51M = 16066186 -198 TCTTCTAGGACAGTGCCAAGTGGAAAGGGGGGGTTGGGGGCCCCACACGGGGTCCTGACTGGGGGGGCGGGGGCTGCCTAGTGGAGCTGTGGGAGGGGGG #################################################@BB:@::BAB6GGGC@GEBGGGGGDGGGGGEGGGGFGGGGEGGGFGGGGGG RG:Z:SRR062634 XC:i:51 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:51 XA:Z:14,+19776588,100M,0; +SRR062634.18118350 147 22 16080800 0 100M = 16080709 -191 GTCAGGCTGGTCTCGAACTCCTGACCTCAAGTGATCCACCTGCCTCGGCCTCCCAAACTGCTGGGACTACAGGCATGAGCCACTGTGTTCAGTCTATTTT >3=9?=5=?1(:<>>>2:-=:CC8CC=CDC@??>+EC@@>FFFEFDGGGAEE5B?@EBEEEA:EBEEEABB:@EE5==@@?5A@?E=BEB5@=BB############### RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:85 XA:Z:GL000244.1,-23878,100M,0;14,-19753012,100M,0;14,+19822726,100M,1; +SRR062634.18478179 163 22 16055377 23 100M = 16055483 184 CATCCCCTGAACCTTCCTAGAGCACATTGTCCAGGCCCCCAAAGCCTCCAAACTCATGCAGTTTGATCAGGCAGCCAGCAGATGCAGGTGACAGGCACTG GGGGGGGFGGFGGGGGGGGFGGGDGGGGDEGGFGGGAGGGGFGDAGEEFEBFFEE5ECEED?BBCBBEDDBD:A=BBFF:BB-?: RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.18842008 163 22 16065289 43 100M = 16065387 198 GCCATTCTCCTGCCTCAGCCCCCCAAGTAGTTGGGACTACAGGCACCTGCCACCATGCCTGGCTAATTTTTTGTATTTTTAGTAGAGACGGGTTTTCACT GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGEGGGGGFGGGGGGGGGGGDEEGEGEEGEE?FGGGGFGGEEGGGGEEBEEEEEEEEDD>5>6?5??:)=BBB=AAA@=@+:BEEAABA???C=EC?=BA/@B@:E?ADCEEFBFFBDCFBABEBDGGEDGGGEFDEAAFDFFFEDFFFF RG:Z:SRR062634 XC:i:96 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:96 XA:Z:14,+19769173,100M,0; +SRR062634.18923396 163 22 16078528 0 100M = 16078617 181 GGGAAGGTGAATTTGGAAACAACATTTGGAATTACTACACTTTAAAGGTAAAGCTGGCAGGAATTCCATACACATCAGATCTGAAAATATGTGTGTGTGT EEEEEEEAECGGGGGGGGFGGDGGGGEEDGGEEFGEFFEGEEEEEDBDBCDD?ECEBEEEC:DA?CCAAA?CC?:A:ACAB@B9C-=5C-:B5??>:5:42><,(,6:'11)-=DDBCC52>-?A?,?A5;6CA??CA?DD@CC5ACC==:C?5DD?6DDDB>@C@=65C:CDD5=? RG:Z:SRR062634 XC:i:94 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:94 XA:Z:14,+19786286,100M,1; +SRR062634.19083820 147 22 16055317 46 100M = 16055213 -204 TTCAGGCTCAGCAGCAGGCAACCCCAATTGCAGAAGTCTCTATCCATCTTGGGAGGTTGTCATCCCCTGAACCTTCCTAGAGCACATTGTCCAGGCCCCC =>;,?:CBB=CA5B:BDE??DFECADFEFDFFGDDEDGGEEDFFGEFEDGGEGGFFGDGDDGGEABGFGFGFGDDGFGEGGDEGFGGGGGFFGGGGBFFG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19787592,100M,1; +SRR062634.19136817 163 22 16062908 0 80M20S = 16063007 159 CCACTGGTTCCTGCTCCAAAAGCCATTATGAACATTATAAACTGGCCAGGTGCGGTGGCTCAAACCTGTAATCCCAGCACCCTTTGGGAGGTTGAGGCAG DDDDDDDDDDDDD?DDDDDDCDDCDDDDDDDDDDDDCDCDDDCCDBCBBCBC?BB?BBBBB=@B@BA?B:?-??@A?@@B#################### RG:Z:SRR062634 XC:i:80 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:80 XA:Z:14,-19779984,100M,0; +SRR062634.19138156 147 22 16088456 0 100M = 16088347 -209 TTAGAGACTGGTATCCACCTGGGGTCCAGGTATCCACTTGGGACCTGATGTTCACCTGGAGTGTAGGAATTCACGTGGGGCCTGGTGTCCACCTTGAGTG ?=3>/=9=6:>6:2=;5::*CCCCBCC4CCBA?CEAC?EEEFADEEECFDFFDADEDEGDGAEEEEGDEEAFAFBFFGGBFGGGGGGGGGEGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.19144807 163 22 16096725 23 100M = 16096883 194 CCCCTCATAGGCTTTCTGCCCACCATGGGGTACCCCAAGAGGCCTGAAGTGCACCCTGGTCTTGAACCAGGGATGCCAGGGTCCCCTGGGCCCAGCTCAG GGGGGGGGFGGGGGGDGGEEGFGGDFDDFGGGGFGFDFDFBDGEGEEBEAEEDEFE@?E:?@@:?ECC@@CC=,C:C??=B?,@?A=?-A-><1;???3: RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19746666,100M,0;GL000244.1,-17071,100M,1;14,+19828645,100M,1; +SRR062634.19197346 147 22 16055045 60 100M = 16054958 -187 AGGTTGTTTGTCTCCCATCAGTGGAGCTGGCAATCTGACGATGCCTCTCTTTGGACCCTGGCCTTCAGCCTCTTCAGGCCATGGTGATCGGGGTCCACCA =B7BBC9C?C==:@-CCABCFDEECGE=EEDEDFBFG@EEEFDFAEEE:EGEGGDEDGGGGGFGGGGGGFFGGGDFGGGGGGGGGGGGFGGGGGDFFFFF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.19224030 147 22 16055071 60 100M = 16054985 -186 CTGGCAATCTGACGATGCCTCTCTTTGGACCCTGGCCTTCAGCCTCTTCAGGCCATGGTGATCGGGGTCCACCAGGAGAAAGATCTCGGCTGTCATCCAA CC@EECBCEFECFED=BEDFDDEFECEC1ADEBF?DEEGGEGAGGFGGGDEFGGEGGGGGGGGGGGGGGGEGFGGGGGGFGGGFGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.19263473 163 22 16065783 0 100M = 16065871 188 ACCTTGCAGCTGCACAGCTAGAGAAATAGAGTCTTCACAATGCCTGGGAAAGGAAGAGAGATGAGGCTCACTAGCCACAGAAAAACAGGCAGTTAATTTT GGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGEGGGGGGGGEGGGEFFGEEEFGEGGGFGGEFBFFDFEGEGEGEGGDGEEFEEEEEEE@EFEEEECCDCC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19777089,100M,0; +SRR062634.19284303 163 22 16082944 0 100M = 16083048 204 TATATGCATTGAAAATGTCTCCCTCAGGTGATTCTGATGTGATATGTGGTCTGAGTTTAAAATGTAAGGAAAATTACCTTCCCTGCCCTCCAGTTGGCCC GGGGGGGGGGGGGGGGGGGGGGGGGCFFDFEGGGGGEGGGGGGFGGGGGGGGGEDGGGGFGGGGEGGGGGGFGGGGFEEGGFGFGFEEGEGEFEEAAAA: RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:1 XO:i:0 XG:i:0 MD:Z:0A99 +SRR062634.19348749 163 22 16096984 0 100M = 16097082 198 TTCCCCTCATGGGCGTTCTGTCAGACTTGGGGTACCCCTAGTGGCCAGACGCACACCCTGGGTTCGAAACTGGGACACCAGGTTCCCTGGGGCCCAGCGC FGGGFGGGGGBEFFEFGGGGGGGGFGGGGGGG=FGGGGE:FFD:FEEBEEEGEEBDGB=A=A@@CDBDBA@AD?@BBBB@5@4-:=@:==?????????@ RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:31T68 XA:Z:14,-19746243,100M,0;14,+19829065,100M,0;GL000244.1,-16811,100M,1; +SRR062634.19371464 147 22 16083729 0 100M = 16083665 -164 CTTGTACACCAAATTTTGGAAATGGATGCTAAAACTTTGATGGCAAAAGGTGACTGATCTGCCTGAGAAATGAACCTACATGATCCCTCTTTCCGCAAAG FDDEECFCECDFDEDEDGGGGGGGGDGGCGGGFEFGFFGGAGGGFGDGGGGGGGGFGGGGGFGGGGGGGGGGGEGGGFGGGFGGGGGGFGGGGDGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+30077,100M,0;14,+19759201,100M,0; +SRR062634.19383963 147 22 16053784 60 100M = 16053667 -217 CACCAACCCTACCTGGCCCCTAACCTGCTGCCTGGAGTAGGCAGGTACCAGAGGCATGGGGTGAGGCATGTTGCAGGTCGAGGACCAGGGCCATCTCACT <,>==(=??:,@A?@ABB@B?=8==5C=EEACAFGEGEEFFEEGGDCBFFFEEGEEBGGGBFGDFDFFGGGGGFAGGGEFGGGEG?GGGDGGGGGGGBAG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.19429719 147 22 16054354 15 100M = 16054261 -193 GAGATGAGACCAGTGAAAACCAGGGAACCCTGCTCATTAGGGGAACCCTTGGTGGTCTCTGCTTGGAAGCCGGGCTGCAGAGGACTGTAAGATCAACCGC DCDE?EEBEEEEEEEEE@BEECEEE@7EAEEEEBGFDF?FGGGGBBGFGFGGGGFGGEFFFGGGGGGGFGGGFGGGGGGGGGGGGGGGGGGGGGGGFGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19788556,100M,0; +SRR062634.19449912 147 22 16072983 23 100M = 16072925 -158 TATGTTGTGTCTTTTTTTAATATAAAAGGCTTTGAAAATGTCCCACAGTTTCTTCAGTAACTGACTGTCAAAAGGGGCAGCCTTCAAAAGTAGAATTCTG ?DDDDDDE=C4;B9AADDEEEEDEBE@EEABEEEEEDEDEEBDEEDEEEDCBAEAEBEEBC?EBBEE=EEEEEECDECEEBEEDEEEEEEEEEEEDEEEE RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19769884,100M,1; +SRR062634.19460776 163 22 16088345 0 100M = 16088457 212 CAACTGGATCTTGTGTCAACATGGGGCCTGATGTCCACTTTGGGCCTAGGTAACTTCCTGACGACTAATGCCCACATGGCTCCTAAGGACCATCTGAGGC GGGGGGGGGGGGGFGGGFFGGGGGGGGGGGGGGGGGGFGGGGGFGGGEGGGDGGGGGGGGGGGGFGGEFFGEGEFGEGGGFGGG=DDECFEDCFFE@EGE RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:6 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:83G16 +SRR062634.19467803 147 22 16088892 18 100M = 16088790 -202 AGGGGTCTACCTAGTCCTCAAGTGAACACCAGGCACTAGATTGACACACAGGTACCAGGTGGATATCAGGCCACAGGTGAACACCAGGCCCCAGGTGGTT ?CACA?A@4AAA@AA?C4AAC??BBCC??DEEEEEEFEEEEEFEEEEEEEEEFEEFEFEEEEFFEFEEEGEEEEEGEGGGEFEGGEGFFEEGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:18 AM:i:0 X0:i:1 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19821682,100M,1;14,+19754041,100M,1;GL000244.1,+24907,100M,1; +SRR062634.19542791 163 22 16074324 0 100M = 16074412 188 ATAGTAGGTGTATATATTTATGGGGTACATGAGATAGTTTGACACAGGCATGCAATGTGTAATAATCACATCAGGGTATAAATGGGGCATTCATCACCTC CCC=5C;;6DDEACECEEAEEEEE:B?CC?DC?BCBCBCCB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-27263,100M,0;14,-19756397,100M,0;14,+19819331,100M,1; +SRR062634.19800690 163 22 16089498 0 100M = 16089596 198 CCAGGTTGACACCAGTCTCTAGGTAGATCCTTAAGCCCCAATTGGTCATCAGGCCCAAGGTGGATACCTTGACCCCAGGGGGTCACTAGGTCCCAGGCAG FFFFFDGEFFGGGGGFGGGGGGGGFGGGGFGGFGGGGGGGDGGGDEGFGGDGCGGGFAGEEGDEDGDGEGFEEEGE@BADEC@C?CE=?CCCED5CEC@@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-24301,100M,0;14,-19753435,100M,1; +SRR062634.19904246 163 22 16073440 0 46M54S = 16073514 170 TCACACTCCTGACACCCAGACTCCACGATACTCCTCTGTCCCCACCCAGGGGCGGTCCCCGCCCCCCAAAGTTTTCCCCTTTTTTTCCCCCAAAAGCGCC GFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG-CECC@###################################################### RG:Z:SRR062634 XC:i:46 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:46 +SRR062634.19916147 163 22 16094211 0 100M = 16094329 220 TAGGTGAATAACACAACCCACGTGGTCATTAGGCCCCAGAATGACACCAGTCCCCAGGTTAACAGGAAGCCCCCAGTGGGCACCTATGCCCCAGCTGGAC GGGGEGGGGGGGGGGGGGGGGFEGFEGFGGGGEEGGGGGGGDGGFGGGGGEGGGGGFGFGGGGGEEGGDEEGGEED@DCDBCAE@BDDEEFG@EBCBCAC RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:5 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:80A19 +SRR062634.1999937 163 22 16081437 19 100M = 16081528 191 CAGAGGACTTTGGGTTCAACCTAAGCACTTTATATCCTTTATTTTCAATTGGGTAGGAGATAAAATAATTCAGCAGCAATAAAAGTCACACTTCTTAAAG GGGFGGGGGGGBGEGGGGGGGGGDGGGGGGGGGGGGGGGGDGGGGFBEGGGGFEFGGDGFEGBBFFGEGFEEDEG=ECCEGBGGGEGGDFGGGDEGDEDG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.20049351 147 22 16063667 0 100M = 16063575 -192 GGATTTGAAGAAGTAACAAGTGGGTTTCTCTAAAGCACAGTTTGGCTTTTTACTAGAATATCAAGTCACCTTTTTTGGGTGAGGAGAGTACAAAGGGAGT CA=A==2?=0@A1?AA::AAAC*DDC?E:EEEGFGFEEAFEFGGDGGGGGGGGGGGGGGEEGFGGGGGDFGGGDGGGGGGGEEFFF RG:Z:SRR062634 XC:i:93 XT:A:U NM:i:2 SM:i:0 AM:i:0 X0:i:1 X1:i:1 XM:i:2 XO:i:0 XG:i:0 MD:Z:29A41G21 XA:Z:GL000244.1,+16292,100M,1; +SRR062634.20482498 147 22 16052786 60 100M = 16052750 -136 TTAGGGAGAAAATACATCAGCCTGTAAACAAAAAATTAAATTCTAAGGTCCCTGAACCATCTGAATGGGCTTTCTTCTAGGCCAGGGCACTCTAAAATTG DEFFFBE9FF?EE=EEDGDGDGGFGGFEGGGGFFGGGGGFEDGGGDGGFGFGGGGEGFGGGFGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGEGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19790124,100M,1; +SRR062634.20563591 99 22 16050128 0 100M = 16050189 161 GGACAACATTCACCTTTAAAAGTTTATTGATCTTTTGTGACATGCACGTGGGTTCCCAGTAGCAAGAAACTAAAGGGTCGCAGGCCGGTTTCTGCTAATT GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGFGGGGGGGFGGGGGEGGGGGGEFGFGGFGGGGFBGGGDE>EEECBCCEEC>E:A@=ADBCAB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19792774,100M,0; +SRR062634.20624314 163 22 16067888 0 100M = 16067954 166 TGCTTGGAGTGGTGTTTTTGCCCTGAATGCCTCCTCCCCCTGGCTTCTTGATGGTTTTATTTGGTGTGACAAGAATGACCCAATTCATATGATCAATTTT GEGGGFDAFEEECEEEEEEEGGEGFGGGFBGGGGGGGGGEEBFEFBEFG?=CEEAEEFCECG=BCFBF:AC:CACED?C=DC=BCDADCDE?EEEDD?EE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19774984,100M,0; +SRR062634.20671502 163 22 16064710 20 61M39S = 16064800 190 AAGGATAGGCTGATGCGGTCTCAGATGGAAATGAGGAACTTACTGGGACTTCGAAAAAAGGGCATGCCTGTTATGCCTTTGCAAAGAACTTGGCTGCCGG GGGGGGGGGGGGGGGGGGGGFGDGGGGGFFFFGGGFGFGGGGGGGGGFGGGGFGGGDGG:?####################################### RG:Z:SRR062634 XC:i:61 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:61 XA:Z:14,-19778201,100M,0; +SRR062634.20745047 147 22 16092785 0 100M = 16092670 -215 GTACACACTCAGACCCCAAGAGGACAACTGTCCCCAGGTTGACATGACTCTCAAGGTGTACATTAGGCAACAGATGTACACCCAGGTCCAAGGCAGACAC AA?A?@75=BCE=AEA=CFF=EE5EEDFDF5FFEDGGGGFB:=<=5@EDDGGGFEGGGGGGFGGGFEFGEDGFGGGGGADGDGEGGEGDGGGGGGGGGGF RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:45C54 XA:Z:14,+19750148,100M,1;GL000244.1,+21014,100M,1;14,-19825572,100M,2; +SRR062634.20796691 147 22 16076022 0 6S94M = 16075936 -180 GGGTGCGTGAGGCCCTGGGGGAAGCTGTCTGCCTACCCTGCAGGGAGTGCCTCTCCTCGGCCATCAGCTGATCCAGTGCCCAGAAGGTGTCTTCCTCTGG ######B@B4AA@@?+AEEEC?@EEBE?DC?BCE?DEEEGEAGGGEGGGDEEEBEEFFGGDGGDFGGDGGDFFGGGEGGGFEGGEGGGGFGGFGGFGGGG RG:Z:SRR062634 XC:i:94 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:94 XA:Z:GL000244.1,+37791,100M,0;14,+19766851,100M,1; +SRR062634.20799768 147 22 16062369 60 100M = 16062277 -192 AATAAAAATGCATGCACTTTCAAATGTACTTAGTGTAACATATTGAACTTAAATTCCAGTTTTCCTGGAATTACTTGTGTCTTGAGCTAAAGACTGTATT GGFGGGGGGGGGGFGAFGGFGFGGGGGEGGGGGDGGGGGGGGGGGFGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19780503,100M,1; +SRR062634.20838869 147 22 16070801 0 100M = 16070708 -193 TCTGATTAGCATCAGGATATCTACATTTTACATGGGAAAGGAGGGAGCAGAGGAAACAGTTAATTATGCGTTCATCTCACATTCAGTAAATCTTCATTTA :DEAEDFGEFE=EGDFFDF?EDEEEFE=AFDFDFGGGFFGGGGGFEGFGFGGGGGGGGEGGGGGGGDGGGDFGGFGFGGGDGFFGGGGDGFFGGFGGGFG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:99T0 XA:Z:14,+19772070,100M,1; +SRR062634.20839193 147 22 16065234 23 100M = 16065140 -194 CTGGAGTGCAGTGGCACAATCTCAGCTCACTGCAAGCTCTGCCTCCCAGGTTCATGCCATTCTCCTGCCTCAGCCCCCCAAGTAGTTGGGACTACAGGCA DCFCCDBFECFEDEFE=FBFFDFE?F:GDGDGGGGEGEFFBGGAGDGGGGGCGGDGGGGGAFGFGGGGGEGFGGGEGGGGGGGGGGGGGGGEFGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19777638,100M,1; +SRR062634.2087100 177 22 16087194 0 100M GL000244.1 26601 0 CCGATCTGGTGTCAAATTCACTGGGAGCCTAGGTATCCACCTGGGGCCTGATGTCCACCTGGGACTAGGTGTCAACGTGTGGCCTGATGTAAACCTCGAG @EEEFEEEEFEEEEEEEEGEGEEFEEGEEEDFEEEEEGEDGEGGEGGGGGGGGBGGEGGFGGGGGFFFFAGGGGGGGGGGGGFGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:4 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:4 XO:i:0 XG:i:0 MD:Z:0T1A1C92T2 XA:Z:GL000244.1,+26605,100M,4;14,+19755739,100M,4; +SRR062634.20876170 163 22 16052205 0 100M = 16052299 194 AAATCTGTACAGTATGTGAAGAGATTTGTTCTGAGCCAAATATGAATGACCATGGTCCATGACACAGCCCTCAGAAGACCCTGAGAACATGTGCCCAAGG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFEGGGGGGGGGGGGGGEGDFF:BBEEGFFEDEEGEEDADFD@@-@@ RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:34A65 XA:Z:14,-19790705,100M,1; +SRR062634.20905629 147 22 16075062 0 100M = 16074951 -211 GGAAAAGAAAACCCAACCCCAAACCCATGGAGAATCCTCACCTTGGGCGAGGCCCTGGCTGGGACTCAGCCTCTTGTCAGACCCTCCAGTAGCTCCATCT @-;4?25:<:/@@?5:??:?A=)=BAA?A??=>3:@A>A3DCB=CFEFDDFCFGGFFFAEEEE:DFF?AFGDEFGEADEEEEEECGFFGGGGGDAGDDFF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19767805,100M,0;GL000244.1,+38745,100M,0; +SRR062634.20911784 133 22 16060584 0 35M65S = 16060584 0 TGTAGTGGCAGGGGCCCGTTATCCCAAACTACCTGGGGGGGGGGGGGGGGGGGAACACCTAAAACCCGGGGGGGGGGGGGTTGGTGGGGGCTTTATCGCA GGGGGGGDGG@######################################################################################### RG:Z:SRR062634 XC:i:35 +SRR062634.20933829 163 22 16089055 0 100M = 16089132 177 CAGGATCCAGGTGGATACCCATGTCCTAGGTGAACACCAGGTTCCAAATGAACATCAGGCTCCAAGTGAACACACAGGCCCCAGTTCAATACCAGCCTCA GGFBFGFEDE?CEEDDFEFEFGFEFEFFGF@EEEAEEBE?@>>>CDCDAA?ECEBCEC?E=FEDDCABC=CCBECEEE=ECCCC9CCCBCEE@EEAC:C? RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:6 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.20936213 147 22 16075609 0 5S95M = 16075515 -189 CCGAGGGGCTACATGGCCATCTCCCAGACATGGGGCCCGGCTGGGGGACATGGGATGGGCAGACACTGCCATCTTGACTTCATTGGCCCATCTGTGGGCT #####AC=C9A?A?@@A@63C::CC5?=A=>-*DCDE@EGDFGGEGEFGGGGFGGGCGGEGEGGGGGGGGGGGEGEFGGGGGDGGGGFGGGGGGGGG RG:Z:SRR062634 XC:i:95 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:20G9T64 XA:Z:14,+19767263,100M,2;GL000244.1,+38203,100M,2; +SRR062634.20951918 147 22 16073419 17 100M = 16073317 -202 GGCTGGCCAGCTGGGCAGGGCTCACACTCCTGACACCCAGACTCCACGATACTCCTCTGTCCCCACCCAGGGCAGATCCCTGCCCTAAAAGTTTTCCCTC ?=AAABBCBCE?D@ECBCCDA@A=A8A@?BCECFADEECECECBF?FDDDEFEDGDEDF@CGDFBFFEGGFGGGGEEGEGGEGDGGFGGGGGGGGGGGGF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:17 AM:i:0 X0:i:1 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.21037423 147 22 16066924 22 100M = 16066834 -190 TTGGAGTTGGGGTCTGGTGGGAGGTGATGGGGTCATGAGGGCTGATTTCTCATGTTTAGTACCATCCTCTTGGCACTGTCCTTGCAAGAGTGAGTGAGTT C@ECADF=EEDEDC?GGGGGFCFFBEEEFGEGGFEFFEGGFGFFGGGEGGGGGGDGGGGGEFGGGFGGGGGGGGGFFGGGGGGGGGFGGGGGGGGGFGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19775948,100M,0; +SRR062634.21053410 147 22 16082624 0 25S75M = 16082494 -205 TTTTTTTTATGAGACGGAGTCTCCCTGTTGCCCAGGCTGGAGTGCAGTGGCATGATTTCAGCTCACTGTCGTCTCTGCCTCCCTGACTCAAGTGATTCTC #########################C==@=5??:?5BCAD:?DAADB?-C:CB?:?C55>@A-@CA:A::C?,8C=>-.= RG:Z:SRR062634 XC:i:75 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:75 XA:Z:GL000244.1,+31207,100M,0;14,+19760331,100M,0; +SRR062634.21129408 147 22 16062666 0 65S35M = 16062543 -158 TAGAAAACCACCTAAGAACAATTGTTTAAAAGTTATTTAGGAAAAATGAAGCAAAATTGTTGTTACAATCTGAGATTTTAACAGCAAGTGCACTCCTGTT ###############################################################################################@>.?C RG:Z:SRR062634 XC:i:35 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:9 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:20C14 +SRR062634.21175705 147 22 16053357 60 100M = 16053268 -189 ACCTGGTCCTACACCTGGCTCCTGGGTTGTCAGCAGGTGTTTGTTGGGCCAACGAGTGCATGGATGGAAACACAGACAGAAGGACAGATGGAGAGATGGT CBEA?CCBAEDEBBFEE:FEG:GGFEBEEBGEGDEEGGGGGFFEDGGFDEGGFGGGFGGGGGGGGGFGGGGDGGGGFGGGGGGGFGGGDGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.21215759 163 22 16065809 0 100M = 16065879 170 TAGAGTCTTCACAATGCCTGGGAAAGGAAGAGAGATGAGGCTCACTAGCCACAGAAAAACAGGCAGTTAATTTTTAAAGGGCTCCAGCTCTTTCTCTTTC GGGDFGGGGGGGGGGGGGEGGGDGGGGGGGGGGGFGGGGGGGGFGGGGGGGGGGFFBGGGGGFGGGEFGBGFDGGGFAEEDGEGFGGGEGDEFGEFEEGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19777063,100M,0; +SRR062634.21220378 147 22 16097882 0 100M = 16097777 -205 CTCGAACCAGGGTCCACAGGGCCCAGCGCAGGGGTTGATGGGAAGGCATTTTCCTCCGTGGGGGACCCAGGCCCAGCTTCTCCTAGGCGCGGCTTTTTTT =4AA?(;6;@:1<5=*A4A,?:@C::;=::;:B3<;:-=ADCC=5C>;C:A;C?665D:=DD:DDAA=ABBB=:==.@;== RG:Z:SRR062634 XC:i:90 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:3 XM:i:1 XO:i:0 XG:i:0 MD:Z:13T76 +SRR062634.21434459 147 22 16057183 60 10S90M = 16057075 -198 GAGCGGAGCAGCATGGGGGGAAGCAGCCCTGGGCCCCCTGAAGTCCTAGGTGGGCAAACACAGACTACAGAGGGAGCTCCCTTTTCAGATTAGTGGTGAA ##########BA:7CCCDDE:EDEEDECD=GEGFEEDEFFB5EEGGEAFEGGGGGFGGGGGGGGEGGGDGGGGGGGGFGFGGGGGEGGGGFGGGGGGGGG RG:Z:SRR062634 XC:i:90 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 XA:Z:14,+19785734,100M,1; +SRR062634.21445731 147 22 16070372 0 100M = 16070263 -209 TGGCATGATCTTGGCTCACTGAAACCTCCACCTCCTGGGTTCAAGCAATTCTTGTGCCTCAGCCTCCCAAGTATCTGAGATTACAGGCGTATGCCACCAT 5CFDDED=DDEEFGEDFCFEECACEFDDFADFDFFDGGGGFGGGGGGGGFGGGGGGGGGGFGGGGGGGGGGGGGGGGGGFGGEGGGFGGGGGEGGEGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19772499,100M,0; +SRR062634.2147406 147 22 16089112 0 100M = 16088992 -220 GGCTCCAAGTGAACACACAGGCCCCAGTTCAATACCAGCCTCAGGTAGACATCAGGACCCAGGTGGACCCCAGGCCCAATGTGCATGCCTAGTCTCTTGG ?A;A<:A<5A5A@CCCEEEEE@CABB?A?DDDB:;6@;CE=EEBE?B?EA?EEEEE?DEEEEE=DD55-DD:CD:--AA5AAB;@0ADD4,.5:@7,6?>=9>><-AA->@6*>==-C= RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:71T28 XA:Z:14,-19788025,100M,1; +SRR062634.2196891 163 22 16059675 23 100M = 16059791 189 ATTTGGGGCTTGTCCTGGGGTGGAATGATATAGTTTGGATGTTTGTCCCTGCCAAATCTCATGCAATTTTAATTCCCGATGCTGGTGGTGGGACCTAATG EEEEE-=E:EBEEEB::D:=@A@8?CC>C@A5D@DDCAC-:>?2?@=CD=:EDBACBAEBDFA?A=C-ECEEBEE-C=>=:@-5?@B@B:+'1>;>C=:: RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783232,100M,1; +SRR062634.21996593 163 22 16055446 22 100M = 16055531 185 GGCAGCCAGCAGATGCAGGTGACAGGCACTGTTGCCACCAGCCAGGATGGCCAGTAGTCAGTAGCCCAGTGCCCTTCCCAGGATGCCCTGCATCCAGGAG EGDFGGGGGGGDGGGGEGGFFFFFFGGGGGFGGGGGGGGEDGGFGFGGGGGGFGGGGCGEDFFFFFDEBCFDFFE?EEEEC:5?CCEB=5A>:AA::A== RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19787463,100M,0; +SRR062634.22021967 147 22 16072554 0 100M = 16072466 -188 CTAGGGTAATTATTAAAACATTAAGAAATACTCTTTAATCATATATGTTTAAACTTGTAGAAAATCAATAACATACAAGTAAAGAAGATCAAGAAGCCAC DBBFF5EECEB=FBFDD@=?DBDC=AC:@:C?CEE=C@:5DA?DBDADDCDD?BD::AAAEEEE=DEE:E?AAAA:DAA;EEDDECCCA,G=DGGGFGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19770313,100M,0; +SRR062634.22055833 147 22 16072644 0 1S99M = 16072556 -187 CAAAAAGCCACAAAAAGCAATTCCAAGTATACAATGCATGGTGAAAAGTGATGGAAGAAAACTGCAGGAATTTAAAATGCAAAATTAAAATACAGCCCTA #A:5CEBC88@GEGGGEEEGBEEFGGFGGEEGGFEEGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:2G96 XA:Z:14,+19770224,100M,1; +SRR062634.22171331 163 22 16079996 0 92M8S = 16080065 169 ATCACTTGAATCTGGAAGGCAGAGACTGCAGTGGGCCAAGATAACACCACAGCACTACAGCCTGGGTGACACGGTGAGACTGTCTCAAAAAACAAAAACA GGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGEGGGGGGGGGGGGGGFFGFGDGGGFFFAGFFFEEACECE(+)=6111/5CD8CAA?A<@?B?DCBBEEECCCAAEGEEEDDDADEC@EEFE?EDEBGGEEEGGFGGEGGGGFGGFGGFG RG:Z:SRR062634 XC:i:75 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:75 XA:Z:14,+19745767,100M,0;14,-19829567,100M,0;GL000244.1,+16328,100M,0; +SRR062634.22390291 163 22 16077676 0 100M = 16077737 161 AGTAGGAAATCTATGGTTAACCAAACAAATTAAATAATGTTTATGAGGAAATTATACTGCTAATGTTTTATTATTTGGCAAGCTATTCAGTGGAACACAA FFFFFEGGGGGGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGGGGGDGGGGGGGGFEDGGGG=DFFFGFGGDGEGGGGEEDABEEDGGGBEDA=CAFFEFF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19765191,100M,0;GL000244.1,-36131,100M,0;2,-132531723,100M,0; +SRR062634.2244701 147 22 16069673 57 100M = 16069591 -182 GAACTTGTACAGAGCCTTGGTCCTCTGAAAGTACCCAGAAATGAAGCAAACTGATCCTATACAACATACACCACAGTCATACCCTTAAGAGAAAGACAGA GEEEEEG?EGEFGEEEEEEGEEEFDEEGGGFGGGFGGGGGGGGGGGGGGGGEGGFGGGGFGFGGGGFGGFGGGGGGGFGGEGFGGGGGGGGGGGGDDGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:20 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.22473110 163 22 16050790 15 100M = 16050896 206 TTGTTAAGGGATTTCCTCTAGGGCCACTGCACGTCATGGGGAGTCACCCCCAGACACTCCCAATTGGCCCCTTGTCACCCAGGGGCACATTTCAGCTATT GGGGGFGGGGGGGGGGGGGGGGGGGGGGFGGGGEGGGGGGGFG?GGGGGGGGGBGGGEGGGFBEBDCFDFFDCB?EBEEEEEE@FD5@=AD@??CDAB=C RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19792112,100M,0; +SRR062634.22479007 147 22 16058484 22 100M = 16058447 -137 ATCCAGGAGCAAAGAAATGACTTAAAGCTGGAATTTATGTTTAAACAGGAAGCAGACCAGAAAAATTTGGAAAATTTGCAGTCTGGCCATGTGACAAAGA DCCEEFFGGE@EGGGEFEEEGEBGFGGGFFGGGGGGGGGFGGGGGGGGGGFGGGFGFEGGGGGGGGGGGGGGGGEGGGGGGGGGFFEGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19784423,100M,0; +SRR062634.2253669 147 22 16095712 0 10S90M = 16095596 -206 CTCCAGATGGCCATAGGTGGATAACTAAGCCTCTCCTGGATATCAGGCCCCAGGTAGACACCAGGCTCCAGGCGAACATCTAGCCCCAGGGGGACATCCA ##########@@9?6=A4@@>>A8?=A?==A2CA@CCAC:?AEEEEEEEEEEEFDECCE5DFDFFFEFGEGGGGAGGGFGGGGGGGGGGGGGGGGFGGGG RG:Z:SRR062634 XC:i:90 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:6 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 +SRR062634.22573022 147 22 16095255 0 100M = 16095150 -205 GACATCACAAACCAGGTGGACACTCAGGCCCCTACTGAATATCTGTCCCCAGGTGGACATCCATCCCAAGGGGAACATCAGGCCACAGATGTACACTTAA EEEBEECEDBABFBDDEEGDGEFEGEFEGEEFEEGGGGGGEGGFDFGDGCGGGGGGFGGDGGBGFGGGEGGFFGGGGGGGGGGFGGGGGGGGGGFGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+18544,100M,0;14,+19748133,100M,0; +SRR062634.22582258 147 22 16064268 22 15S85M = 16064158 -195 GGTGCTGTCCTGTTGATGGTGAGTGAATTCTTGCAAGACCTATTTTTTTAAAAGTGTGTGGTGTTTCCCCTGCCTTGCCAGGTGATGTGCCTGCCTGCTC ###############@@ABA:???=?B?BA>A>:5A>;>93&4@@,>CC@AAA>5:<;>:E=BAEF=GGGFGFDGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:23T6G54 XA:Z:14,+19778619,100M,2; +SRR062634.22601302 147 22 16063217 29 100M = 16063132 -185 ACCCTATCTCAAAAAATAAATTACAAACTGATTCACAATAACTTTGGTTTTGTCACTAATATGCTGAATATTTTTGTTACAACTAATATGCAAAATGCAA EGBGFDEFGAAEEGBGFFDGEBEEEEEFDFFFDFFGGDGGFGEGGDEGGGGGGEGGFGGCGGGGBGEFGGGGGGGGGGGGDGGFDGGGGGEFGGGGGDGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19779655,100M,0;21,+10636123,100M,1; +SRR062634.22665881 147 22 16054450 0 21S79M = 16054322 -207 GGCAGGGGACTGTAAGATCAACCGCCGGGGCTGTAGGTGATGGCTCAGGCACTGGATTCTCCATGGAGGACCCCATGAATGCACTTTGAAAGGTCAGCTC #####################EBEEEGGGFFEGFGGFGGGGGFFGGGGGGEAGFGFGGDEGDGGGGGGGEGBBGGGGGGGGGGGGGGGGGGGFGGGGGGG RG:Z:SRR062634 XC:i:79 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:79 XA:Z:14,+19788481,100M,0; +SRR062634.22717174 163 22 16053387 60 100M = 16053474 187 CAGCAGGTGTTTGTTGGGCCAACGAGTGCATGGATGGAAACACAGACAGAAGGACAGATGGAGAGATGGTGGGTGGCCAGACAAAGGAGTAACTTGGTGA GGGGGGGCGEFFFFFGGGGGFGGGGGGGGGGGGFFGGGGGGGGGGEGEGEGGGEGEGEEGFCFCGEDFEBGEGEGEGE=EBEDEEEDBD9?CCCADC>=5?:?5:CA>CCEEEEEEAEEEEAEEEED?GCEAEEE=CCEEFGGEGFGFGGEECB?@>CAC RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:85 XA:Z:14,+19767270,100M,0;GL000244.1,+38210,100M,0; +SRR062634.2296540 147 22 16060462 29 56M2I42M = 16060351 -209 TGAGCCAATTAACTCTCTTTTTTCTTTTTCTTTTTCTTTTTCTTTTTCTTTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTT CE?B5BECACECGFG>FEEECGFF?G;FGGGFFGFGGFEGFGGGFGGGGGGGGFGEGGCGGGFGFGGGGGGFGGGGFGGEGGGGGGGGEGGGGGGDGGGG RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:1 XG:i:2 MD:Z:98 XA:Z:14,+19782445,100M,5; +SRR062634.22999200 163 22 16054881 22 100M = 16054984 203 CCAAGACAAAAAGGTATGAGGCTGGTTGTCCAGCACCGACCCACCCTGGGGCTGGAAAGGCTGGGAGGGCAGGTTGGCACATTCCCAGTTTGAGGGAGCC GGGGFEGEGGGGFGGGBFGGGFGGGGGGGFGGGGGGGGGGGGGGGGGGGFGFGFGEFGGGGFGGFEGFE;DFFCFFDD@GB:FFDFBFEEEECFDBAEBD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19788028,100M,0; +SRR062634.23011367 163 22 16061776 15 95M5S = 16061834 158 TCTACTGACTTCATGTCTTTACACAATGGTTTAACTCCCTTGACCAACTGTCAACTAAAGAATTCCTAAAACTCACCTATGACTTGTAAGACCCCTCGTT DDDDDDD?DBDADDDBDDDDDDDDDADDD??ACAABDDD?BD??=?C??A;A?A?CCBD?A?>AAA??237?:?C@:B##### RG:Z:SRR062634 XC:i:95 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:95 XA:Z:14,-19781101,100M,0; +SRR062634.23055198 163 22 16080130 0 100M = 16080200 170 TACAAAGACTGATACTAGAGGTTGTTGGAGAGGACATTGTTCCACATAATATCTTATGAGTTGCTGATGGACAAGTAAGATGGTAAAATGCCTTTGAAAA GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGEGGGGGGGGGGGEFGGGGFGGGGGGFGEGGGEFGFFGDDGEEFEEGGGEEFG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:35G64 XA:Z:14,-19762737,100M,1;GL000244.1,-33671,100M,1; +SRR062634.23138785 147 22 16088988 10 100M = 16088891 -197 GGTTGGGTTACTTATAGCATAGGTGGCCATCAGGTCCCAGGTCTATAGCCACTCCCCACCTGAAAATCAGGATCCAGGTGGATACCCATGTCCTAGGTGA DDBFGGFFECFEFGGBGFGFGFG?EDDEDFBEFGEFGGGGGDGGGGGBFGGDGFGGGEGGGGGGGGFGGFGEGEGGGGGGGGGGGGFGGGDGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.23184540 147 22 16080500 17 37S63M = 16080398 -165 TATATGTATATATACACATACATATATATTATATATGTATGTGTATATATATTAGACAGATATATATGTATATAATACATATACATATATACATATATGT DGGGDFGGFGGGGFGGGGDFGGEGGGGGGGGGGGGGGGGGGGGGGGGFDGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:M NM:i:0 SM:i:17 AM:i:17 XM:i:0 XO:i:0 XG:i:0 MD:Z:63 +SRR062634.23228418 147 22 16075217 0 100M = 16075121 -196 TCTTCCACTCAGCAGGTACCGCATGATCTTTGGGGAGGATTTGTGTTATACCTGCTCCTGGTGGGATAGGAGCCTCCGGAGCTGGGCAGTATTTGGGCTG ?A?AAC@D9CEEEAEEEEEEGEEDEEEEEEEGEGGEGEEEFEGGGGEGGEGGEEGEEEEGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19767650,100M,0;GL000244.1,+38590,100M,0; +SRR062634.23243072 147 22 16070403 0 100M = 16070294 -209 CTCCTGGGTTCAAGCAATTCTTGTGCCTCAGCCTCCCAAGTATCTGAGATTACAGGCGTATGCCACCATGCCTGGCTAATTTTTGTATTTTTAGTAGAGA B=BGFGGEDGGGGEAGEFGGGGFGGFGGGFGGGGGGGGFGGGGGGGGGGGFGGGGGGGGGGGGGEGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19772468,100M,0; +SRR062634.23278600 147 22 16063002 0 100M = 16062947 -155 AGGCAGAAGGATTGCTTGAGCCCAGGAGTTTGAGAACAGCCTGGGCAACAAAGTGAGACCCAATCTCTACAAAAACTAGAGAAAAAATAGCTGGGTGTGG ??:,?BBBDE@B:B@=BB?BE-DADE?GFEBCA?=?ABE?BCECD,:FFFFFDGGFGEGGGG?EFGGFCFGG?EEGDDGGGCGGGGFGGGFFGFBEGFGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.2332992 147 22 16063157 60 100M = 16063064 -193 GTTGAGGCTGCAGTGAGCCAAGATCATGCCTCTGCACTTGAGCTTGGGTGACAGAGTAAGACCCTATCTCAAAAAATAAATTACAAACTGATTCACAATA FEEEEEEEAEEEEEEEE?GEDFBEEEE,FFEGEAGDFEFGG?GGGGFGGGGGGGGGGGGGFGGGGGFGGGGGGGGGGGGGGGGGGGGFGGGEGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19779715,100M,1; +SRR062634.23345893 163 22 16053256 23 88M12S = 16053361 205 TGGGTCTTCTGCTGCAACTCAGCCCCTCTGTACCAGCCCTGGCCTCATTCCCTTGGTTCCAGGACACCCAGCTGACAAAAGGGACTTGCCTGTACCCCTG GGBDGGGGDDAFFFFDGGGEFGG?DB:EEDEGGDEGEBEBBEADDE?DEA:@?AAC@CC=AB@;@CCC@/CC,AACA*?@?,??,3??############ RG:Z:SRR062634 XC:i:88 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:88 XA:Z:14,-19789666,100M,1; +SRR062634.2335459 163 22 16083297 0 100M = 16083360 163 GCGAATACTTATAGGTGATTTTGGATTTCCTCTGATTGCATAAGGCTCTGGGTAACCTGATGCTATAGACCAAATGTTATTGTCCCCCTCAAATTCAAGG FFGGGGFGGGGGGFGCGGGGGFGGFGDGGGGGGGGGGGGGGGFGEFGGGGGFCEGGFDFDEGGEEEGDFEFBEEE?AEEEFDEGEEEECD@D:ACED=CA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19759633,100M,0;GL000244.1,-30509,100M,0; +SRR062634.23393747 147 22 16096579 18 10S90M = 16096488 -181 AGCCAGGGATGACAGGGTCCCCGGGGCCCAGCGCAGGGGCTGATGAGAAGGCACTTTCGTGCGTGGGGGACCCTGGCCCTGCTTCTCTGTGGCGCGGTTT ##########?=:9)A<:>@?AA>?1A7<=C,A?EEECE@DBDC=E@DE;88=('1:;8GGGGEGFGGGGGGGGGGGFGGGFGGGFGGGGGGGE RG:Z:SRR062634 XC:i:90 XT:A:U NM:i:1 SM:i:18 AM:i:0 X0:i:1 X1:i:3 XM:i:1 XO:i:0 XG:i:0 MD:Z:50C39 XA:Z:GL000244.1,+17230,100M,2;14,+19746823,100M,2;14,-19828499,100M,2; +SRR062634.23531594 147 22 16058911 36 1S99M = 16058824 -186 TAAGCCTCGGAAACTTCCATGTGGTGTTAAGCCTATGGGCATGCAGAGTGCAAGAGTGATGGCTTGGCATTCTCTGCCTATATTTCAGAGGACGTATGAG #?BAD6>@.>AB5AD:>:5:?:,?A5.B9.B6;@.>+55>?-5DDD5EE5A?5AC::=5D=5=FF:D:FFDDDDD5=-EDDEEED>@>@;5B?:?C>ACC RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:14,+19783997,100M,0; +SRR062634.2355614 163 22 16052165 25 88M12S = 16052260 195 AACAAACAAACAAACAAACAAACAAACAAACAAACTGTCAAAATCTGTACAGTATGTGAAGAGATTTGTACTGAGCCAAATATGAGTGACCATGGTCCAT GGGGGGGGGGGGGGGGGFGGGGGGGGGGBGGDECG;55??=??>5/3286:A=<:/256>B?A@BB-=5+;;9<43=@==?@@@@?0?############ RG:Z:SRR062634 XC:i:88 XT:A:U NM:i:4 SM:i:25 AM:i:0 X0:i:1 X1:i:0 XM:i:4 XO:i:0 XG:i:0 MD:Z:2A66T4A10A2 +SRR062634.23628948 163 22 16075170 0 100M = 16075251 181 TTCCCTGCCCCATGGAACCCTGGGCCTCTAGGAAAGAGTTGAGGGTGTCTTCCACTCAGCAGGTACCGCATGATCTTTGGGGAGGATTTGTGTTATACCT GGFGGGGGGGGGGGFGFGGGGGGGGGEGGGGGEGGGDGEFFFFFF?EEEDGGGGGGGGDGGGGEGEDEEECDDAEGEFEGFGEGE;?CBE@BCDCEEBEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-38637,100M,0;14,-19767697,100M,0; +SRR062634.23661209 147 22 16090777 0 100M = 16090704 -173 TTGGGCCTAGGAATCTACCTGGGGCCAGATGTCAATCTGGGGCCTGAAGTCTACTCAGGTTCAGCTGTCCACCTAGGGCGGGGTGTTTTTCTGGGACCCA B?DCDDAC?CA@?@?C9DEDGEEFFEEEDEAEEBFBGEGGFGEFGGGEGGGGEFEGGGGGGGFGGGGGGGCGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19752156,100M,0;GL000244.1,+23022,100M,0; +SRR062634.23735128 163 22 16075046 0 100M = 16075129 183 CTCCTCCTCCAGGGCAGGAAAAGAAAACCCAACCCCAAACCCATGGAGAATCCTCACCTTGGGCGAGGCCCTGGCTGGGACTCAGCCTCTTGTCAGACCC GGGFGGGGGFGGGGGGGGGGGGGGGGGGGGFEGGFGGGGGGGGDGGGGFDBGGGGEGGEEEDEDGEEEDEGEFAEDEDEBDDDADEA@E@CCA RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19786051,100M,1; +SRR062634.23811334 147 22 16068763 0 100M = 16068686 -177 CTGCAGGCATTTCCGCAGACCTGGGACTGAGAGAAGGATGCCATTTTTAATCTGGGCTCGTACAAAGTCAATCATTCTTTGGTGCCTGGCAATGGTGGCC F?GGFBFFE=ACEEGGEFGFEGGGGFGGGGGGGGGFFFFEEGGGFGGGGGGGGGGGGGGGGGGFGGFGGGFGFGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19774109,100M,0; +SRR062634.2381660 163 22 16052452 23 100M = 16052555 203 CTGATTGACAATTGGTTGAAAGAGTTGTCAATAGAAAGGAATGTCTGCATTGTGACAAGAGGTTGTGGAGACCAAGTTTCTGTCATGCAGATGAAGCCTT FFFFFFDBFFFFFFFFFFEEEDCEEEEEEEDEEE=BDDDDEFFFFFFDFFFDAF5:9>>7CC=CCCEE@C5ABBAF@DDFFE@@?B=B5C@?B@@BB@?B RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19790458,100M,1; +SRR062634.23864416 147 22 16060745 0 100M = 16060660 -185 CTGCCACTGCACCCAGCTAATTTTTTGTATTTTTAATAGAGATGGGGTTTCACCATGTTAGCCAGGCTAGTCTTGAACTCCTGACCTCATGATTCACCGG B-A5B==>=@(-?-6:BAAA,AF5FFF:FAFEBBEEEF-BFFGD=GFGGDBGE?GGGFBD?FEDE:EFDDFDEDE=?E=E?ECACDD RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:52A47 XA:Z:14,+19782127,100M,1; +SRR062634.23874540 147 22 16086410 0 100M = 16086339 -171 TCACCGGGGACTGGTATCCAGCTGTGGCCTGATGATCTACCGCATCCTGTTGCTCACCTATGGGCTGGTGTCTACCTGGGGCTTGGTGATCACCTGGGAG AC>=CCDCBBCCCCE?=@BBBDEEDD=EEEEEEEEBEECEEEFEEFEEEFDGEGGEDEGEFGFGGGGGFFGGGEGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+27389,100M,0;14,+19756523,100M,1; +SRR062634.23933318 163 22 16098033 0 100M = 16098160 208 CCTCACCTCTCCTCTCTCAAACGTTAACTTCCCATCATGGGCTTTCTGTCCGACTTGGGGTATTCCTAGCGGCCCAAGGCGCTCCCTTGACTCGAACCAT GGGGGGGGEFGGGGEGGGFGDFEGGFBGGGGFGGGGGFFGFGAGGGGGGGEBDEDFFFGGBDFFDACCCEEEGEED:DC=CDBB@CC@C,=>=>,:>>3> RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-15766,100M,0;14,+19830112,100M,0;14,-19745202,100M,0; +SRR062634.24104349 147 22 16064095 23 100M = 16064004 -191 ATATAAAAGTCACAACATGAGGCATGAGATACCTCAGATGTTGTTCCCTCTAAATCTCATGTTGCATTGTGTTATCCCCAATGTTGGTGGTAGGTCTGGA BEDEGBEBCBDBDDBDEEGFGGGFEEE?ECBEEADGGGGGEFGGGEEEAEFGGB=GFFGGGGGGGEGGGGGFGGDGDEGFGBGFGGDGGGGFFGGGGDFF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19778777,100M,1; +SRR062634.24188342 147 22 16086506 15 100M = 16086420 -186 GGAGCCAGATATCAACTTGGGGCCTGGGTGTCCACTTAAGGCCTGATGTGGGCCTGGGGCCTGATTGCCCACCCGGGGACTTGGTGTCCACCTGGGGTCT ??3@7//8/8999397:5?ADC?A?ACCAA=AC?C@@?ACC5?CE?DEBCGEEEEGGGFEEGGGGEAA?EAEFDFGGDEEFFFFFEEEEAFEGGGGGGFG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:81G18 XA:Z:14,+19756427,100M,1; +SRR062634.24208951 163 22 16061542 0 100M = 16061647 160 TTCCATAAAATGATCTTTGGCCAATAGTGGGGACTGAGGCCTCAACTCTGACCTTTTATTCTCTTGCCAAAAATCCTGTCTAAGGGGCCTGGAGATTCAC GGGBGDBFGDGGFGGGGFDGFGEAGBEEECF?=FDE=CBEGDGGFFFFAGEEGEDFDGEG>EDFEFGEEB?BABBFDD?EB==B==??B@D=ABA@B@5@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:2,-132547811,100M,0;14,-19781330,100M,0; +SRR062634.24280820 163 22 16071145 0 100M = 16071226 181 CTGAGGAGACATGCAGCCTTCTGTCACTGTGGGAACCTGGCTTATGGATGAGGCTATGACACAGGGTTGTGAAGTTATAGCTATGTGTTTGAGAACAAGG GFGGGGEGFGDGGGGGGG?GGGFEGFGGEEDGGBDGGGFGFGFEFFGEGFFGGGGFGGEFEGDGGECDECEC?ECDC?CDEEEDEFEEE3BCBB=FE=F@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19771726,100M,0;9,-68329546,100M,1; +SRR062634.24304065 147 22 16093338 0 100M = 16093253 -185 TGCACAACAGGTCACAGTTGGACATTCAGCCCCAGGTGAACATCAGGCCATGGGTGGATAAACAGTCCACAGGTGGACATCAGGTCAAAGGTGAACATCA C?CCC=FEEEDEFFDFEEGFEDF?:DGGGEEFDGEGEGFGFEFFGGGGEEGGGGGGGGGGGFGFDEGGEGGGGGGGGFDGGGGGGGGGGGGGGGGGGDGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:6 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.2431756 163 22 16074096 0 100M = 16074200 204 TGTGCAGTGGGGAGGACAAGGCCCCCTCCTCTGTCTCCTCCAAGACCCTGTTTTCTGAGGAAGGTCACTCTGGGAACTGTTGGCCTCTGCAGATGGGGGC GGGGGGGGGGGGGGGGGGGGEDGGGGFGGEFEEEGEGGBGGGEEEGFFEEEEGEEEEEEEEEFE?CCECEEEEE?BCCCCEEECBEE?EEAA?AAA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-39711,100M,0;14,-19768771,100M,0; +SRR062634.24369770 147 22 16063469 60 100M = 16063366 -203 ACGAGAATATGTTAAGATACCACAAAAACAAAAAATATACATAGTTAAGAGCTTTCATGAGTAAGCAAATTTTATACATATAAAAATAATTTACATGCAA E?GFGGGGGFGBFGGGGGGEGGGGGGGGGGGGGGGFGGGGGFGEGGGGGGGGFGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19779403,100M,1; +SRR062634.24381897 147 22 16063097 37 100M = 16063028 -169 TGTGGTGGCACTCACCTGTACTAAGGAGGCTGGGCTGGGAGGATCTCCTGAGCCTAGTTGGTTGAGGCTGCAGTGAGCCAAGATCATGCCTCTGCACTTG ;>5?<-5;A5>A=>:>5?<<,7:)4D?B=@C>A:=,=DCD@><4>:?;3=:==?5DD?=DAAAA?D:DBBDD?DDD=D=DDAADDC>C>A6>=?>>C?C@ RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.24396953 147 22 16073897 0 100M = 16073786 -211 CCCTGTGAGACCTGCAGAGCAGGGGCTCTCACAAAGGAGGAGTCGGGAATGTCAAGATTTTAACACCTTCTCCAAGTCAGGATCAGGAAATGCTGTCCCA EE=DBCCCC@DEBE=EEDGDEGEGFEFDDGFGGGGGGGGGGGGGGGGFGGGFGGGGGGGGFEFFCFGFFGGGGFGEDGGGGGDGGGGGGGGGGGGGEGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19768970,100M,0; +SRR062634.24399064 163 22 16097666 0 94M6S = 16097778 212 CGAACCAGGAACGCCAGGGTCCCCTGGGCCCAGCGCAAGGGCTAATGGGAAGACACTTTCGTCCGTTGGGGACCCAGGCTCCGCTTCTCCGTGGTGCGGG GGGGGGGFGGGGGGGGGFGFGGGFGD?GGGGFGGFFEDGFGBGG?GFGFDEGAGFFEEGEFEEGEEEEEDEDDDEBEC5CD:BDCBE=EEA@@C###### RG:Z:SRR062634 XC:i:94 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:94 XA:Z:GL000244.1,-16135,100M,0;14,-19745573,100M,0;14,+19829742,100M,0; +SRR062634.24456714 163 22 16052044 60 100M = 16052165 180 AATCTTAGCTACTGCAGAGACTGAGGCAGGAGAATCGCTTGAACCCAGGAGGTGGAGGGTGCAGTGACCTGAGATTTTGCCATTGCACTCCAGCCTGGGC DDDDD=BA?D?DD?DDDADDD?DA=DBDAD=BAA=A:B';?8?5:CC=::?5?--B@:=?CDB:BBAC?C5?B:B5:::: RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:58T41 +SRR062634.24472247 147 22 16089024 0 15S85M = 16088926 -183 GGTGGCCATCAGGGCCCAGGTCTATAGCCTCTCCCCACCTGAAAATCAGGATCCAGGTGGATACCCATGTCCTAGGTGAACGCCAGGTTCCAAATGAACA ###############@@@5:?@B=C=?AC&=4/=:5BEAE@BBB?-DAAAE5EEEDAC=B@?=CC?E?EE?EAEE?DEEED-DDAEEEA?DDEE?DDCDD RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:2 XO:i:0 XG:i:0 MD:Z:14A51A18 +SRR062634.2504971 147 22 16071810 0 100M = 16071729 -181 CCAACCAATTCTGCGTTTTGAGTGTAAATTACCATAGCCTTTGTAACAGAAATTAGGCAGCAACTCTCGTCATTCAAAATGCATAAGCAATTTGGCTAAG CADAD=EEEEEEFEEGGEEGDEEEDEEEEECFEEEEDBEGEFFFA@DGGGGGDGGFGFGGGGEGGGFGGGGEGGEGGGDGFGGGGGGGGGGEGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19771061,100M,0; +SRR062634.2621871 147 22 16063267 22 100M = 16063166 -201 TGTCACTAATATGCTGAATATTTTTGTTACAACTAATATGCAAAATGCAAACTGGCTATCACTACCATTCTGATAATGGTATTAGCTGTCCTACGTACCT ADAD?C,BEBABBC?EB?;A;AEEEEDEADEBEB?EEE=E:EEEECEEEEDEEEEDEEEE@DCB??EE?EEEEEEEEEEEEDEEEEEEEEDEEEEEEEEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19779605,100M,0; +SRR062634.2634503 147 22 16052430 0 13S87M = 16052319 -198 GGGGTTTCCAGGTTATAAGTAGATTTAAAATTTTTCTGATTGACAATTGGTTGAAAGAGTTGTCAATAGAAAGGAATGTCTGCATTGTGACAAGAGGTTG #############@AAB=5AA::A:,93A5A*>8175=??B@:98*>5>7C>692C?8CA?5C:@ACC:@@C>CC7997A3=D5?D:DDDD RG:Z:SRR062634 XC:i:87 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:87 XA:Z:14,+19790493,100M,0; +SRR062634.2650102 163 22 16053623 0 65M35S = 16053728 182 CATCCTTGTGGGAGGCACTCACCTTTTGCTCAGCCTCTTGTGGCTACAGCCCAGCAGGTCCCAGGGGGCACCAGCCAAGATGAAGGGGGCGTTGAGGGGT GGGGGGGGGGGGGFGDFGGGGEGGGDGGGGGGEGEGEGGFGGEFGE,C@:@?:<,=;>69>>45?################################### RG:Z:SRR062634 XC:i:65 XT:A:U NM:i:2 SM:i:0 AM:i:0 X0:i:1 X1:i:1 XM:i:2 XO:i:0 XG:i:0 MD:Z:1G34A28 XA:Z:14,-19789322,100M,1; +SRR062634.265500 163 22 16068206 0 100M = 16068285 179 CAATAATCCCAGGAGAGAAACATACCTTGGCAAATTCGAAGATATTTCTAATTTTATTTAACAAATTTAAACTTAATTAATGTATCAAAGATTCGCCTAA GGGGGGGGGGGGGGGB-CDDGFGEGGGEGDGGGGGFGGGEGEEFGGGFGGBGGGEEGGEEEGEEEGGEGEEEFEEEFEEEGGBEEEEEEECFEEEGEEEC RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:2G97 XA:Z:14,-19774666,100M,1; +SRR062634.265672 163 22 16060380 48 100M = 16060481 201 CTTCTGCCAGAATTGTGAACTTTCTGAGACCTCCCCAGAAATGGATGCCAGCATTATGCTTCCTATACAGCCTGCAGAACCATGAGCCAATTAACTCTCT GGGGGGGGFGGGGGGEGEGGFFGGFGGGGGGGGGGGEEGGGEDDG-ACDBDEDDDFFDGEFGEEB?EFDDEFCCFEEEEEEEEEBCDAA@A@C9CADCCA RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19782527,100M,1; +SRR062634.2718504 163 22 16070237 0 100M = 16070342 205 TAAAAAAGGAATTTAAAAATAATGAACCACACTTTTGAGAAATATAGGATTATGTAAAGGAACAAAATCTATGACTTACTATTGGCATTTTTTTTTGAGA GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGEGGGDGEEGGGEGGGGFEEEGEFEEGEEEGEGGGGEEEEEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19772634,100M,0; +SRR062634.2729515 147 22 16094634 0 100M = 16094524 -210 TGGATATCAGGCCTCAGGTGAACCCTAGGTCCAACATAGACATCAGGCCTTAGGTCGACACTCAAGCACCAGATGGACTGCTGCACCTAAGTAGAAAACA ACECCEEEEEEBEEECEDDD?EGGEEEGFEEGEEEEEGFEFEEGEGFGGGGFGGGFGEGEGGFGGGGEGGGGGGFGGGGGGGGGEFGGGGGGGGGGGGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19826963,100M,0;GL000244.1,+19165,100M,0;14,+19748756,100M,1; +SRR062634.2876822 147 22 16054124 0 100M = 16054052 -172 AGCCTGCCTACCCAGCACATCTCTCCTCATGGACATGAGAGAAAACCCTCCAGCTTGGCCCTCACATCTGTGAAACCCACAGTAATGGGGCTGACATCCT =?BDCDCCC:CEEEEEAEAAC?C9BFDED@FEEBFEEEFGBFDFCBFFCFGGEGBEFDGGGGGFGGFGEGGGGGFGFGGFGGGGGDGGGFGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19788786,100M,0; +SRR062634.2889923 163 22 16090821 0 100M = 16090927 206 TGAAGTCTACTCAGGTTCAGCTGTCCACCTAGGGCGGGGTGTTTTTCTGGGACCCAGAGTCTACCTGAAATCTTAGTGTCAACCTGGGGCCTATGTGTCC FFFFFFFFFFFFFFFEFFFFFFFFFFFFFFFFFFFFFFFEFEFFFEFEEBEEFFFFEEEEFEEEFEEEEEEEEFFEEEEE?EEEBEED@EEEEEEEEE@B RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19823610,100M,0;GL000244.1,-22978,100M,0;14,-19752112,100M,0; +SRR062634.2918482 163 22 16055807 22 100M = 16055899 190 TCAGACCCTTGGAGGGCAAGTGGCCAACTCAGGGCCTGGTGTCTTTGTCCACTGCTACTCTGCCTGGTCCTCAGATCAGGCACCTGGATTTACCTCCTCA GGGGGGGGGGGGGGGGFGGGBGGGGGGGGGFFGBGGDGGAFDDFFGGAFGEGEFFAFFDE:CEECDD?FECEECCCFBCDFBCFBCB::CC?DDADEACA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19787102,100M,0; +SRR062634.2947227 147 22 16086468 18 22S78M = 16086361 -185 CTACCGCATCCTGTTGCTCACCTATGGGCTGGTGTCTACCTGGGGCTTGGTGATCACCTGGGAGCCAGATATCAACTTGGGGCCTGGGTGTCCACTTAAG ######################A>9AA:A@CAC@ACC>CBEECEEEFEEEECDCD@AAEEGAEGEFGGGEBFGAEGGFGFGCBGFGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:78 XT:A:U NM:i:0 SM:i:18 AM:i:0 X0:i:1 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:78 XA:Z:14,-19817605,100M,1;14,+19367454,100M,1;14,+19756487,100M,1; +SRR062634.2987383 163 22 16069302 20 85M15S = 16069421 194 AAGGGGCACACCAAAGCCTGGGAACAGATTTAGAGAGGTGGTCATCTCTTGCCCTCTGTCACCTTCCTAGTGTACTGCTGCAGACCCAGCAGTGGACATC EEEEAEFGGGGDFGFGDGGFEGE?FDFCEEGEGBGEEEEE@ABBBCAEACEDCDAC?A@AEEE>?AC<4>:?3?;>=A>?@:<:@############### RG:Z:SRR062634 XC:i:85 XT:A:U NM:i:1 SM:i:20 AM:i:0 X0:i:1 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:72C12 XA:Z:14,-19773585,100M,2;2,-132540079,100M,2; +SRR062634.3138488 163 22 16079362 0 100M = 16079422 160 AAAGTCTGATCTTTGGATTCTTAAGGAAAAATTTTATTGGACTGCAGAGGCATAAATTGGATTTCTAATTCTGGAGTCAGGATGTGCTCAGTAAACATGG GFGGEGGGGGGGGGGEGGGGGGGGFFGGGGGGGGGGGGGFGFGDGGGGGGGFGEEGGGGGEGGGGGFFGGFGGEEEEGEEGEEGEEDEGEDDDDDAECCC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-34435,100M,0;14,-19763499,100M,0; +SRR062634.3149997 163 22 16058933 29 100M = 16058997 158 GTGTTAAGCCTATGGGCATGCAGAGTGCAAGAGTGATGGCTTGGCATTCTCTGCCTATATTTCAGAGGACGTATGAGAAAGCCTGGGTGTTTAGGCAGAA GGGGGGGGGGGGGGFGGGGGFGGEGDGG?GFEGDGGGGGGEEFG=BEDEEFEEEDEEGEEGEEEGEEEBEEACDFCC?CCECCCECB1B/???@?A@CAA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783974,100M,0; +SRR062634.3203184 163 22 16098661 0 99M1S = 16098752 191 TCTGCCCCAGGTGCCTCACCAGCTACTTGGGAGGCTGAGTCAGGAGAATTGCTTGAGTCTGGGAAACCGAGGTTGCAGTGAGTCAAGAGATCACAGCACT GGGGFGGGGGGFGGGGGGGGFGEEGGFGGGGDGFFFFDGEGEGDEGBB=EGEGGEDE?EDEEFCBED?DCDCBCDBBBAB?A8AA?CCC::CCAA@,AA) RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:GL000244.1,-15140,100M,0;14,+19830755,100M,1;14,-19744562,100M,1; +SRR062634.3351841 147 22 16055025 60 100M = 16054912 -213 GGCCTATGATTAGGCCTGAGAGGTTGTTTGTCTCCCATCAGTGGAACTGGCAATCTGACGATGCCTCTCTTTGGACCCTGGCCTTCAGCCTCTTCAGGCC C:C:=C=C?BCCEEED=GBGDF>FDFDGGFDGAGGFEDGGGGGGGEFGGGGGEGGGGGGGFFFGGFFFFFGDGGEDGGFGGEGGGDGGGEFGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:37 AM:i:23 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:45G54 +SRR062634.3373760 147 22 16059989 15 100M = 16059893 -196 ATGATAGCACCATGCTTGTACAGCTTGCATTAAAGCTTGTGAACCCATTAAACATATTTTCTTTATATTTACCCAGTTGAAGGTATTTCTTTGTAGCAAG EEGGDEGGEEEEGGGFGFFEFFGGBGGFGGGGGGFGGGGFGGGGGGGGGGFGGGGGEEGGGGGGGGGGGFEGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19782918,100M,0; +SRR062634.3402989 163 22 16051061 46 100M = 16051136 173 AAAGGCCTAGAAAAGGAAGGCCTGACTGCATTAATGGAGATTCTCTCCATGTGCAAAATTTCCTCCACAAAAGAAATCCTTGCAGGGCCATTTTAATGTG GGGEGEGFGEABFFFGGGGGBGGD?FFDFFE?FFFDGGAEGGG:GGGGGDBB?EEEAEECFGDFGEBDEBB:EBA;EEEBEEF?BD@EDA:@A=+AB=BB RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19791841,100M,1; +SRR062634.3413999 163 22 16050792 22 100M = 16050904 212 GTTAAGGGATTTCCTCTAGGGCCACTGCACGTCATGGGGAGTCACCCCCAGACACTCCCAATTGGCCCCTTGTCACCCAGGGGCACATTTCAGCTATTTG GGGGGGGGGEGGGGGGGGFGFGGGGGGGGGGEEEFGGGFBGBEEGFGGEEEEEEEEEEEEC@DADBDCD?:?>>:ACC(A??A:,?67;:A?9==?CEEE=@F?FF RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783415,100M,0; +SRR062634.3767224 147 22 16085247 0 100M = 16085154 -193 GTGGGGATGCAGCACAAGAACATCTTGCTCTCTTGAGCGTCTCCCACCAAGTGAGCTGGCTACGGGGCTAACTCTAGGATGTGGGTGCCTGGTTATCGGA A:CCCCA?CCA@EEEEEEEEE?EEEEDEEEEFEEEEEEEEDCDEDAEEEEEEEEGFEEFEEEGGGGGGGGFGGGGGGGGGGGGGGGGGFGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+28558,100M,0;14,+19757684,100M,1; +SRR062634.3836617 163 22 16063353 60 100M = 16063453 200 TGTCCTACGTACCTGACGACCTAATGCTTAACCTAATGTTCCTCTTTCCCACTTTGATTTGGAAACTTGTTGCAATAGAGTTCTTTACCCCAAAGTCTTA GGGGDGGGGGGGGGGEGGGGGGFGGGFGGGGGGGGEGGFFGEDBEGFDFD@DECFD?DEEDB=CCEEE@FCE@ABEE@EEEEEEED?;DA?6,?3398>: RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.383903 163 22 16076533 0 100M = 16076661 198 TGGGCCACTGATCACCACACAATGGGTCCTCTGATGGCCCAGGAGCTGCCTGCCAGGCACAGGAGGGCGGCTGGGTCCAGACCCCATGTGGGCAGCCCCT GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGGEGGGFGFDGGEDE5EEF:GGE??BE?BB=A@CCBC5AB6@BC(5C?>?(9 RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:98A1 XA:Z:14,-19766334,100M,1;GL000244.1,-37274,100M,1; +SRR062634.40183 163 22 16065713 0 100M = 16065802 182 ACGTGCATTATATCATTTCTTACTTTTCAAGGAAAAACATGTTTTACGACTTGAGTTTATCTGTCTAGTGACCTTGCAGCTGCACAGCTAGAGAAACAGG EEEBEFFFFDGEDGGGFGGGGGGFGAGFEFGDDFGGGGGGGEDGFEGGDGFGEFFEEGGGEFGBGFFB@EDFEEBEGGEFEEFBFFEDBEBDEDBEEEE= RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:6 X1:i:4 XM:i:2 XO:i:0 XG:i:0 MD:Z:96T2A0 +SRR062634.404230 147 22 16088646 0 56S44M = 16088481 -209 GACCAGGGTATACCGGGGACCTGATGTCAATATGGGACCTGGGCATCCATCTAGGACCTGATGTTCACATAAGGGCTGGCATTCTCCTGGCCTGGTTCCC ########################################################@B?@AB@C5@@AA?CEBE@EBECEA?EBEFGECFFDAGGFF=FF RG:Z:SRR062634 XC:i:44 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:44 XA:Z:14,-19821436,100M,0;14,+19754343,100M,0;GL000244.1,+25209,100M,0; +SRR062634.4055787 163 22 16094250 0 82M18S = 16094267 106 AATGACACCAGTCCCCAGGTTAACAGGAAGCCCCCAGTGGGAACCTATGCCCCAGCTGGACATCAGGCCTAAAGTGGACACCCAGGATCAAGATGGACAT GGGGGGGGGGGGGGGGGGGEEFFFFGGFGEFGGGGAEECEEEEFEEEEEEDDD>>><)A################## RG:Z:SRR062634 XC:i:82 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:82 +SRR062634.4088212 163 22 16090238 0 98M2S = 16090345 207 CCTCAAATGGACATGATGTTCCAGATGGATATGAGGCCCCAAGTGGATACTAGGCCCAGGTGGACCCCAGGTCTGAGGGGTACACCAGGTCCCAGAGGAA GGGEFGGGGGDGGGGGGGGFGGGGGGGGGGGGGGGEGGGGGGGEDGGGFGEFFEFAEBCCA@??A@==BGGGGDBGDF@FEAEEE:FEFDF?FFFFAFGEGFDDCFGGFGEAEEDEEDEEDDFEDFEDGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19778658,100M,0; +SRR062634.4549361 163 22 16070185 0 100M = 16070281 196 TAATTTAAGAGCTCAAAGACCAGACCTTTGAATCAACCCAGTCAGACAAAAATAAAAAAGGAATTTAAAAATAATGAACCACACTTTTGAGAAATATAGG EEEEEEEEEEDDEEEEEEEEEEDDE?EDDDDDDEAD=EBEDEEBEBEEEEDCAEAAACC?EE=DEEEDAEDEACA??DAC=D?BBEEEEEDEDBEE:EA: RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19772686,100M,0; +SRR062634.4608717 147 22 16065037 23 3S97M = 16064962 -172 TGAAAAAAAGCTAGGTGCTGATAGCCAAGACAATTAAAAAAAGGTCTTGAAGGCATTTCAGAAATCTTTGCAGCAGCTCATCCCAGCAGAGACCTAGAGG ###BBEFEEF=@CDDEE@EEDADDD:EFEEBFEEEEGGGGGEGGGGGGGEGFFGGGFFGGGGGGEGFGGGGGGGGGGGGFGGGGGGGGGGGFGGGGGGGG RG:Z:SRR062634 XC:i:97 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:97 XA:Z:14,+19777838,100M,1; +SRR062634.4650450 147 22 16054765 60 100M = 16054680 -185 CCCAGGTATGGGTCAAGAGGGCTCCAACTTCAAGGGTGAGGTTATTCATGGTAGAAGAAACTGTAAGTCTCCCAGGACATTCCCTCCATGTAAGCAGGAA ?CA;@CC<@C>?;A?AEE@CEAAEECBDCBCCCDFGEGEBFAE=EFEEEEGGEFGFEGGEGEGFGEGGGEFGGFGGGGGGEFGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.4656511 163 22 16055811 22 100M = 16055865 154 ACCCTTGGAGGGCAAGTGGCCAACTCAGGGCCTGGTGTCTTTGTCCACTGCTACTCTGCCTGGTCCTCAGATCAGGCACCTGGATTTACCTCCTCATGCT GGGGGGGGGGGGEGGGEGGGGGGGGGGGGGDGGGGEGEGGGGGGGGGGGEGEGFEFEEEECBEAEEEEEEEEEBEEFECEAED,AACAA?AEC=C9:9A@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19787098,100M,0; +SRR062634.4761788 163 22 16072413 0 84M16S = 16072513 200 TTTCTGGACATAGTTCAGGTAAAAAAAAAAAAAAAAACCACTGTTGAATAATGCTTCCGGTATAACACAATGTGTTTGAAAATGTAAACAATACACTGTA GGGGGGGGGGGGGGGGFGGGFGGGGGGGGGGGGGGGEED:??:?@:B@-B?3@?@BE@CC85?@@.3.<3:?5>>CA2=BBAA?################ RG:Z:SRR062634 XC:i:84 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:35G48 XA:Z:14,-19770470,100M,1; +SRR062634.4770057 147 22 16054425 0 64S36M = 16054296 -165 GACAAGGAAAACCCAGGGACCCCTGCTCTTTGGGGGACCCCTTGGTGGTCTCTGCTTGGAAGCCGGGCTGCAGAGGACTGTAAGATCAACCGCCGGGGCT ################################################################A?><5=-:AAC><:;>(856-5:DCFGGGGGGGFGG RG:Z:SRR062634 XC:i:36 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36 XA:Z:14,+19788549,100M,0; +SRR062634.4789722 133 22 16071485 0 35M65S = 16071485 0 ATAATATAATGAAAAATAGAGACTGAGAGAGAAAAAACAAGACCCTTTATCTTTATATTTTTTCATATGTGTTTTTTTTTCTGTCTGGTTTTTTTGTTTT #################################################################################################### RG:Z:SRR062634 XC:i:35 +SRR062634.4851733 147 22 16084401 0 41S59M = 16084276 -184 CCCCCCAGCCCCCCCGGCAGGCACGCAAAGCCCCCCAGCAGCAGTAGCAGCAGGCTCTTCAGCAGCAGGGCTGCTGCTCTGCTGAATGAGAGAAGTCCCT #########################################A4A;:==,6?,C@C?=A;C=?CCCACEC?F?EEEGFEDFGGGGGGGGGGFGGGGGGGGG RG:Z:SRR062634 XC:i:59 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:5 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:59 +SRR062634.4860515 163 22 16076256 18 100M = 16076373 192 AGAGGTCACCTGGGAGAGCTCCCTCTAGCCCTAGGGGATTTAGGGTGCAGACTCTGCACCCCTCTCCTGACCCTGGGCGTGAGGACTAAGCAAGTCCCCC FGGGGGGGGGGFFGGGGGEFGGGFGBEGGFGEECGGGFEBEGEEECGGEFEFEEEGGEDEDFCEBCC588>?@>=:CCC?CAC5CAA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19766611,100M,0;GL000244.1,-37551,100M,0; +SRR062634.4930603 147 22 16066899 46 100M = 16066822 -177 GCATGTTGAATTGTAATTCCCAGTGTTGGAGTTGGGGTCTGGTGGGAGGTGATGGGGTCATGAGGGCTGATTTCTCATGTTTAGTACCATCCTCTTGGCA @=ACC@C=9CBCCCBEEB=EBCCAC?C@E:DFDDEDE?ADEGEEGGCEEDEEEEGGGGEEDGGFGGGFFGGGGFFGGGFFGGGGGDGGEGGGGEGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19775973,100M,1; +SRR062634.4971548 163 22 16081162 0 100M = 16081211 149 ATGCTCACATTGTTTTGTTTATATGAAAATGTTTGTAACTACTACACGCATTTTTGCTACTTAAGCCTTTTGGTCTTGCTGCCATAGCAAATACTGTGCC GGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGFGGFGFEEGEGGEFFEEEDEFEF@FEGEECFE= RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:45G54 XA:Z:GL000244.1,-32644,100M,1;14,-19761768,100M,1; +SRR062634.4997602 163 22 16091471 0 100M = 16091580 209 CATAGGTTGCCCAATATCAGGTAGTTTTTTATAACAGTGTGAAAATGGACTAATACAAGGTCTTAGGATAACAACCATGCTTCAGGCCATAGGTGGACAT GGGGGGGGGGGGGGGGGGGGGEGGGGGGGGGGFGGGGGGDGFGGGGGGGGGGGFGGGGGFCGGGGGFFEGGGGFEGGGGGBEEDEEGEDEEEE?EDCDDB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-22328,100M,0;14,+19824260,100M,0;14,-19751462,100M,1; +SRR062634.5006015 163 22 16065688 15 100M = 16065768 175 AAAGACAGTTATTCAATTAATTGAGACGTGCATTATATCATTTCTTACTTTTCAAGGAAAAACATGTTTTACGACTTGAGTTTATCTGTCTAGTGACCTT FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFEFFFFFFCFFEEEFEFFFFFFFEFEDFBECFEEB RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:15 AM:i:0 X0:i:1 X1:i:6 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.5029545 147 22 16055853 60 100M = 16055771 -182 GTCCACTGCTACTCTGCCTGGTCCTCAGATCAGGCACCTGGATTTACCTCCTCATGCTCCGCCCCTGGCTGCCCCACTGCTCCACCCACCTGATCCAGCA 88>83?::?:6=4=>-B??AD??C@C@BC?BAB5ECEFFFEFDFFFDFE=BD?AAAAA:;;@C RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:56T43 XA:Z:14,+19781678,100M,2; +SRR062634.5613583 147 22 16050264 0 100M = 16050206 -158 TAACTTCCTGGAGGGAGGCTTATCATTCTCTCTTTTGGATGATTCTAAGTACCAGCTAAAATACAGCTATCATTCATTTTCCTTGATTTGGGAGCCTAAT @E=EAEEEEGF?EEGEGFFEDBGEGFGGGGFGFFGFGGGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19792638,100M,0; +SRR062634.5713832 147 22 16081079 0 100M = 16080983 -196 TACAGTGATAACTTTTGTTATGTATAAAATGATTTCTGTCTCAGTTAAAAAATATTTTTTGCATGAGTTATTGTCATATACAGATGCTCACATTGTTTTG EEEEGEEEGEBEEEGDGEBFGEFGGGFGFGEGGGEGGGEGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGFGGGGGGFGGGGGGGFGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19761851,100M,0;GL000244.1,+32727,100M,0; +SRR062634.5888676 147 22 16081039 0 100M = 16080973 -166 ATTTACAACATCTATATACAAACAATATGTTACCTCAGAATACAGTGATAACTTTTGTTATGTATAAAATGATTTCTGTCTCAGTTAAAAAATATTTTTT @ACC@E:EE?BEEFEEBED?EFFBGBEBEEF@DFDDEEFGEBEBBEEBDEGE;GGGFGGGFFGDGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19761891,100M,0;GL000244.1,+32767,100M,0; +SRR062634.5907528 163 22 16096331 0 95M5S = 16096458 201 CCAGGCAGGTGTTGGGCTCGAATAAGGGTCGCCGGGGACCCCGGGGCCCAGCGCAAGGGTCGATGGGAAGACACTTTCACTGGTGGGGGACCAAGGTCCC GGGGGF?FFDGGGFGDBGGFDBFFFGFGAGFFFEE:EBEEEF?AD:BEEEDFEB==A=:(>A==@BB5?B5AAA?7?@@BA?A=@7@-<+?????##### RG:Z:SRR062634 XC:i:95 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:95 XA:Z:GL000244.1,-17473,100M,0;14,-19747059,100M,0; +SRR062634.5927323 163 22 16072960 46 100M = 16072999 139 CCCTGCACACAGACCCACAGCAGTATGTTGTGTCTTTTTTTAATATAAAAGGCTTTGAAAGTGTCCCACAGTTTCTTCAGTAACTGACTGTCAAAAGGGG GGGGFGGGGGGGGGGGGGGGGGGEEGGEGGGGEGGGGGGFGEFGEGDEGGEEGEEGGEGF-?@>@FEDEDEEEEEEEEDEEEEEEAEEEEEEEEEEEAEE RG:Z:SRR062634 XT:A:U NM:i:1 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:60A39 XA:Z:14,-19769907,100M,2; +SRR062634.5992619 147 22 16072939 23 100M = 16072880 -159 GGCAGCTGTGCCCTTCCTGGTCCCTGCACACAGACCCACAGCAGTATGTTGTGTCTTTTTTTAATATAAAAGGCTTTGAAAATGTCCCACAGTTTCTTCA CFEDEEEGBEEEEEEEGGEFEFGGGGGEGEGGGDEGGEGGGGGGGGGGGGGEGGGGGGGGFGGGFGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19769928,100M,1; +SRR062634.6088536 163 22 16063860 0 100M = 16063971 211 TTTTCTAATCCCAACAATGCAAACCTCAGCCACTGTCAAGGGCAATCTCAGACTGGGGGTTGCTGCCTAGGGTGAGGTTCTAGCAAGACTTACAAGCAGA GGGGGGFGGGGDFGGGGGGGEGGFFGGGFFGGDGDFGGGGBDGEEFEGDEGDDF?EEBEA=ECBACDDCDDEAE@CBAACEEEEBCDAA@ACCCA@D??? RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:74G25 XA:Z:14,-19779012,100M,1;2,-132545498,100M,1;21,+9845102,100M,2; +SRR062634.6172826 163 22 16087370 12 100M = 16087481 211 TGTATTCACCCATGGCCTGATGGTCACCTGGGGTTGAATGTCCAACTGTGGCCAGATGTGCACCCGGAACCTGGGCATCCACCTGGGGCCTGATGTTTAG GGGGGGGGGGGGGGGGGGFGEGGFGGGGGGGGGGFGFGGAEGGFGFGGEGEEEGEDCD?DDEGDFGEEEEEEB@EDDCEDEEEEEBC@BC?C=CB>A# RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:0A16C81 XA:Z:GL000244.1,-15132,100M,2; +SRR062634.6274777 147 22 16079689 60 1S99M = 16079688 -100 GTATATATACATATATATATACATATACGTGTATATATACATATGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTATATATATATATAAAAGA -C5=CAAA?-C?CEEB?-?5D:D:D>?,:>AE5EDA?AA-D:D::D?GGGGGGDFFGGGGGGGGFGGGGDGFGGGEGGGBGFEGGGEGGGAGDFGFGGGG RG:Z:SRR062634 XC:i:99 XT:A:U NM:i:1 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:16C82 +SRR062634.6315188 163 22 16089135 0 100M = 16089229 194 CCAGTTCAATACCAGCCTCAGGTAGACATCAGGACCCAGGTGGACCCCAGGCCCAATGTGCATGCCTAGTCTCTTGGAATACATCATTTTCAAGGTGGAC GGEEFGGFGGGGFGDGGGGGFGEGGGGGGGGFGFFGGGGGEGGGGFGGGEEFFGEEEFEEGEEEEECDFCEEEEEEEBAEEEEEECE?EEE=5CAACCAC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:8 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.6454224 147 22 16068105 0 100M = 16068007 -198 TTGCTCCGTCCCACATCTGGGAAAGGTATGGAGGTCTACATCAGGTCTATGTCAAGATTTATGGCTTAAAACATCTATCCGAGGAAAAAAATAAAACTGA B?EB=@BA9A=>(@=@?B@EFEFFBF?EEFFEAFDCCBCA-FEFEEFFEDCF:FFFFFDFFFDDBFF?DAEEEDEF5FFFFFFFFCFEFFFFFDFEDBFB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19774767,100M,0; +SRR062634.6485861 147 22 16052967 0 100M = 16052923 -144 CAAGGGCTGACTTGAGATCCCTGCCAGATGAGTCACTGGCAAAAGGTGCTGCTCCCTGGTGAGGGAGAAACACCAGGGGCTGGGAGAGGCCCAGAAGGCT :>:CB@CACC?@:???@>/=CAC=AEDDBDCCDD9D@CDDGEBEED?A:FED@GFEDGGDGBGFGFFGGFE@DEBEGGGGFGGGGFGFAFGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:18A81 XA:Z:14,+19789943,100M,1; +SRR062634.6497971 147 22 16074511 0 15S85M = 16074402 -194 CGTGGGGAAACCGGCCGCCATGATTCAGTTACCTCTCACTGGGTCCCTCCCGTGACACATGAGGCTTATGGAAATTACAGTTCAAGATGACATTTGGGTG ###############D=?DBDB:CCCCCBC5BBAB5C@C?EBG=EEEEF:EFEEECCDEADEBED?AEEDAFGFGEEEEDAEEEDFGGFGGFDDFGFGGG RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:85 XA:Z:14,+19768371,100M,0;GL000244.1,+39311,100M,0; +SRR062634.6506287 147 22 16051121 46 100M = 16051051 -170 TCCTCCACAAAAGAAATCCTTGCAGGGCCATTTTAATGTGTTGGCCCTGTGACAGCCATTTCAAAATATGTCAAAAAATATATTTTGGAGTAAAATACTT BEEBDDEEEEGGEGFBBEGEEGGEEFGEDEGGGGGGEGGGGGGGFFGGGEGGGGGGGGGGDGGGGGGGGFFGGFGGGGGGGFFGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19791781,100M,1; +SRR062634.6513061 163 22 16059528 0 100M = 16059599 171 TTTACAGGCTCATAGGTGGAAAAGACTTGCCTTGCCTCAGATGAGAACTTGGACTTGGGACTTTTGAGTTAATGCTGAAATGAGTTAAGGCTTTGGAGGA GGGFGDGGEFGGGGGE5EEEEEEEDGEFGGEGE=EGFGBGEBFBGB=FEGGD?BFEF?F@EDEGD?CECEEEFFAEBBED-?=?>EBC?BCCEDBC>CC? RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19783379,100M,0; +SRR062634.658637 163 22 16064446 0 100M = 16064535 169 CAAACTTCTTTTCAAACCTCTTTTCTTTATAAATTACCCAGTCTCAGGTATTTCTTTATAGCAATGCAAGAAAATTGGTACTGAGACATGGGGCATTGCT GGGGGGFGFGGGGGGGGGGGGGGGGGGGGBGGGGGGGGGEDCEEGEEC;EE@EFEEEEGGDGGFCDFCDCD?F??CCC>ECDECCDECB?A1CCA?>CCA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19778426,100M,0; +SRR062634.6589014 163 22 16055723 23 78M22S = 16055810 185 CAGAGGGCCTGGCCATGGGTGAGGGCCCTGGGTAGAAGACCCCTCCCTGCAATTCCCACACCCCTGGAAGCTGTAAAGCCCAGCTCAGACCCTTGGAGGG FFFFFBBDFFGBDGGGGEA?FDFFADEEBEEEB-E?B?@BDEFDEEE5DAB?EE?AA:?BBB@C*@A@???,?3?@@?###################### RG:Z:SRR062634 XC:i:78 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:78 XA:Z:14,-19787208,100M,1; +SRR062634.6650508 147 22 16065858 0 4S96M = 16065764 -190 TATCCACAGAGAAACAGGCAGTTAATTTTTAAAGGGCTCCAGCTCTTTCTCTTTCTCAGGGGGGATTTGGGTTTTCTTACATACAACTGAGTTTCTGCTT ####?5:61>'=>5=?:4@?8A?C?@@A?A@??8?=ACEEEEDEECE=ECDCEBCCFEEEEEFGGGGBFEEGGADGGFDGGGFGEGGGFGEGGGGFGGGG RG:Z:SRR062634 XC:i:96 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:6A89 XA:Z:14,+19777018,100M,1; +SRR062634.6683814 163 22 16094487 0 99M1S = 16094592 205 TTGAACATTAGACCCCAGGTGAACACCAGGTCTTAGGTGGACATTAGGCCCCTGGTGGACATAAAGTACCAGTGGACATCCATGTTGCAGGTGGACACCC GGGGGGGGGGGDGFGGGGGEFGGFAGGGGEEGGDGFF@FFGGGGGFAAGEDEEFDAEBDEDBED?ECECCEEEEEC?DCE?CCA@C?A==<9>75=CCA# RG:Z:SRR062634 XC:i:99 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:0C98 XA:Z:14,-19748904,100M,1;GL000244.1,-19313,100M,1;14,+19826816,100M,1; +SRR062634.6711108 163 22 16073809 0 100M = 16073888 179 CGAATGGCCAATGACAGGAATGGGAAGGCCCACTGTGTCATTCACAAAGCACTTCCAAACCCATCACCACAGGTGACCCTCACAACAACCCTGTGAGACC GGGGGGGGGGGGGGGGGGGFEGGFGGEGEGGEEGGEGFEFEGEGGEGFEGEGEEFEEEFEEEGEEEEEEF@CA>CAACC@8B>DEEEEBB888';BB>(7.*7:??#################################################### RG:Z:SRR062634 XC:i:48 XT:A:M NM:i:0 SM:i:12 AM:i:12 XM:i:0 XO:i:0 XG:i:0 MD:Z:32 +SRR062634.7268093 147 22 16069163 0 100M = 16069059 -204 AAGAAACAGGAGGGAGGTGAACCCGTTTCCCCAGGAGACCTAAACTTTGGTGAGGACCAGTCCTAAGGGAAGGAGGAATACAGCTTGCCGAACACAGCCC ?29?@>=0<=*2BA:BCCCC:BDDEEEEFEBEBDE=EEEEEGEEDGDGGFDFEFG=DFEDGGEEGCGGGGGGGGGGGGGBGGFEGGGFGGGGFGFGGGFG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19773709,100M,0; +SRR062634.7274906 163 22 16090332 0 100M = 16090441 209 GAGGAACACCAGACCCTAGGTAAGCATGCAGTCCCAGGTGGACATCAGGTTCCAGGAGGACACCAGGACCCAGTTGGTCATCAATCCACAGTTGAACACC GGFGGGGGGGGGGGGGGDEGFGGGGGGGGGGGGGGGGGDGGGGGGGGFGFGGFGGGGGGFGGGGGDGGGEGEEEEEECFEEGEGFEGEGFEBDDDDEFDC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-23467,100M,0;14,-19752601,100M,1;14,+19823122,100M,1; +SRR062634.7301099 113 22 16070263 0 100M 14 19772605 0 TCTCACTTTTGAGAAATATAGGATTATGTAAAGGAACAAAATCTATGACTTACTATTGGCATTTTTTTTTGAGACAGAATCTCACACTGTTACCCAGTCT BC@C=DDDDEEGEEFFDD5DDGEEEGDGEFEGGGGEGGGGEEGGFGGEGGFEGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:0C1A97 XA:Z:14,+19772608,100M,2; +SRR062634.7386055 147 22 16080275 0 100M = 16080232 -143 AGAAATTCCATCTTGTCATTTATTTCCACAAGAGAAAACAAGAGACCTGTATTAAAATATCCACTTCTATTGAAAATAATCACAGTTTCCTGAAATTTTT EEGBG@?GFGEGAEGGGEEFGGGGEGEEGGGGGGGGGFGGGGGFEGGGGGGGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33A66 XA:Z:GL000244.1,+33526,100M,1;14,+19762592,100M,1; +SRR062634.7392094 163 22 16064899 0 100M = 16064968 169 AAGAGCAGACCATACAAGTATAGAAAATTTGCAGCCTGCCCATTTGGCAAAGGAAGGTAAAGCATTTTCAGGGGAAGAATCTAAACAGGCTGTGGCACTA GGGGGGGGGGGGGGGGGFEGGGGGGGEGGGGGGGGGEFGGGEEFEEBFEEEEFDBED;EDECCCBEEEEBEEEE=DD?C@E?CCDC@ACD@C?D5:>>6: RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19777973,100M,0; +SRR062634.7502452 147 22 16095910 0 100M = 16095815 -195 CTGTCCGACTTGGGGTATCCCTAGCCACCCTAGGCCCACACTGGACTCGAAACAGGGATGCCAGGGTCCCTGGGGCTCAACGCAAGGGTTCATGGGAATA A5>,5>A,?56:?@?A,5BBBB-ABC=CCC@@@@CCCAC=DCBEDFAEEEBEFEDEGGGGGFFGFFEFGGDGDDGGGGBGGEGGFGGGGFDGGGFGGGGF RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:2 XO:i:0 XG:i:0 MD:Z:35G22C41 XA:Z:GL000244.1,+17889,100M,2;14,+19747478,100M,3; +SRR062634.7532215 147 22 16056947 37 50S50M = 16056784 -213 CCTTGGCTGCGGGGCCCCTGCGGGGAACAGTAATGGTCTCTGCCCACCAGCCTGTGAGCAGCAAAAAGGCAGGACTGGGTCTTGAACCCAGTCTCCCTCC ##################################################@=?@@@=?-@:=BABB5@BAADF=EBDGEEBEDEEEBFFEC?B4;EEEEE RG:Z:SRR062634 XC:i:50 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:50 +SRR062634.754379 163 22 16094137 0 100M = 16094215 178 CCAGTTTGATACCCAAGTCCTATGTGGGTGCCAGGTCCAGGGTTACACTCAAGCCCCAAGTGGACACCAGGCCCTAGGTGAATAACACAACCCACGTGGT GGGEGGGGEGGGGDGGBGGGFGGDGGGFEDFFGDFGGGFAFDEFFGFGGGGFFGGGGEDAAEECEEEGEBEGBEEBCA>=BCCCCEEE=EEEEEED??CA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19749253,100M,0;GL000244.1,-19662,100M,0;GL000244.1,-20119,100M,1; +SRR062634.7566852 147 22 16063130 37 4S96M = 16063032 -194 CTGGGCTGGGAGGATCTCCTGAGCCTAGGTGGTTGAGGCTGCAGTGAGCCAAGATCATGCCTCTGCACTTGAGCTTGGGTGACAGAGTAAGACCCTATCT ####???C5BCDAA@CCCC@DCECCBAECBADDAE@BDEBEEEGEDDEEBEEF?EFE?EGGDGGGECE?BFGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:96 XT:A:U NM:i:1 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:24T71 +SRR062634.7623885 147 22 16072083 0 2S98M = 16071980 -201 GGGGCTACAAGTGCGCATCACCATACCCAGCTAATTTTTGTATTTTTTGTAGAAACGGGGTTTTGCCATGTTGGCCAGGCTGGTCTCAAACTTCTGCCCA ##?-A5:5<:>,>C=.=*?3:>@9.8?:B=7;4A.-6A::0=7>C)5+7/=:B98:::CC;>995--3A-6>6B6:CACCB.5B56C=>C;55;55-;;; RG:Z:SRR062634 XC:i:98 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:94A2T0 XA:Z:14,+19770786,100M,2;2,+132537314,100M,2; +SRR062634.7745521 147 22 16096853 0 51S49M = 16096720 -182 CGGGTCCCCGGGGCCCGCCTCAGGGGCTGGTGGGAGGACACTTTTGTCCGGGGGGGACACAGGCCCCGCTTCTCCGCGGCAGGGTTTTTTTTTTTTTTTT ###################################################@>A@>',*9@=:ABC@-:?5:@A?::ACA-DE@@FGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:49 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:49 XA:Z:14,-19828773,100M,0;GL000244.1,+16994,100M,0;14,+19746589,100M,1; +SRR062634.7769519 163 22 16081691 0 85M15S = 16081808 187 CAGACAAGGGCGTCCAAGTGGATTGCCGTTTTTACAGTTGAATCAAAAAGCTTTTATAAGAAACTCCTCTCAGCTATGTATAAAAAACTGTCTGCACAAT GGGDGGFFAFGGGGGDEEEGFEGFEFFFEEEED3ECECADACBCCDEDD?A@DDCC?C?CCCDA:CCACCC5=A=,4472:)???############### RG:Z:SRR062634 XC:i:85 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:85 +SRR062634.7774141 163 22 16067564 39 100M = 16067684 169 AATTAAAATGTTTTATTTGAGATACAAGTATTACAATTTGACACATACATGCAGAATGGGTGTCTTTAGTATGTCCGAAACATAAAGAGAAAGTTAGAGG FFFFF?EAEEFFFFFFFECFFFFFFDFFFFFEFEBDFFDFEB5EEEEEEBDEAFBEFFEDEFFEFFFDF@CC?B?EE?@E=BC@@@C?B@7?@:?@B@A@ RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:16 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19775308,100M,1; +SRR062634.7830631 163 22 16091015 0 100M = 16091120 205 GTTTTACATCTGGGCCCCAGGTATACACTAGACACCAAAAGAACTCCAGCCCCTATCTTAACATGAGGTCCTAGGTGGATGCCCAGGCCTCATGTCTACA GGGGGGGGGGGGGGGGGGGGGEGGGGGGGGGGGGGGGGGGGGGGGGGGEGEGFFEGGGGEFGGEFEEEFEFEE@FEEEEEDEFEBBCC?CEEEEEFEEEB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19751918,100M,0;GL000244.1,-22784,100M,1;14,+19823804,100M,1; +SRR062634.7884703 163 22 16098210 0 53M47S = 16098346 211 CTCTAAGGTGCGTTTTTTTTTTTTTTCTCTGCCCCAGGTCCCTCACCTTCCCCTCATTGGCCTTCTGCCCCCCTTGGGGTCCCCCTATCAGGCCCGGGGC GGGGFGGGGGGGGGGGGGGGGGGGG@6<0:*;AB?,:=,?AA5>)/:386?A?############################################### RG:Z:SRR062634 XC:i:53 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:24C28 XA:Z:14,+19830289,100M,1; +SRR062634.7937929 163 22 16069767 29 14M1D84M2S = 16069874 207 TACAGAATGAAAAATTAAAAAAAAAAAACTCCTTCCAAATGATATCAAATTCAACAAAAGAAGCATCAAGTCTCTCAGATGAGAAGAAACCAGCACAAGC FGGGFGGGGGGGGGGGGGGGGGGGGGGGEDGFGDGEGGG?DBFBEEGGBFGFBGEGGGD5EA;?ECBCBCBDDECB:CA@BC@:EDCEDAC?@DDDCB## RG:Z:SRR062634 XC:i:98 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:1 XG:i:1 MD:Z:0G13^T84 XA:Z:14,-19773106,70M1D28M,2; +SRR062634.7998607 163 22 16056562 12 100M = 16056647 185 TTTCTTCAGAGTCGACTATGAGCAGAAATCACCCAGACACCCCACAGGTAAGTCCACTTGCCTCACATGGCTCTTCCTCTTGTCAAGAACCTACAAGGGC GGGGGGFGGGGEGGGGGGGGGGGGFGGGGGGGGGGFGFGGGGGFGGGDEGGGEFDBGEEFGEEEDFD@CCFBCCCDECCCEEEEECCBC?DAAA2==9=@ RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19786347,100M,0;18,-14987783,100M,1; +SRR062634.8003424 163 22 16079678 57 100M = 16079775 197 ACATATACGTGTATATATACATATATACATACATATACGTGTATATATACATATGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTATATATAT GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGFGGGGGEGGGGFGFGEEEEEEFFFFFFFFFFFFFFEFEFEFEFFFFBEDFFGGGEF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:20 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.8004962 147 22 16054219 23 100M = 16054183 -136 ATCCTCTGCCCTATGCAAGAGAGGTTTCCCAAGCACTTGCAGCAAGTGAGACTGCACAGGATGGCGAATCCACAAAGAACACGTTGTTCTCATGCTCTTT @EEEEEBEEEEEBEEEEEEEEEEEEEEEEEEEGEEGEEGFEGEEEGEGGGGGGGGEFGGGGGGGGGGGGGGEFGGGGGGGGGGGFGGGGGGGGGGGGGGF RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19788691,100M,1; +SRR062634.8007234 147 22 16066081 15 100M = 16065984 -197 GGCTCCAGATGAGCTCAGGCCATGGCCTCAGAGGGTGCAAACTATAAGCCTTGGTGGCTTCCATGTGGTGTTAAGCCTGCAGGTGAGCAAAGTGCAAGAG CCAA9AC=?CC=,:=@AAAA RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19773516,100M,0;21,+9850686,100M,1; +SRR062634.8021494 147 22 16087542 39 1S99M = 16087435 -206 ACCTGGCGCCTGATGTCCACCTGGAGCCTGAGGACCCTTCTCAGGCCTGATGTCCACAATTGGCCTGGTATTCATCTGGGACCTTTGTGTTAATGTGGCA #CCCCD?DB@CCCDCCCDBDEEEGEBCCEFEFDBEEFEGGGFFGGEEGGFGGGGGFGGGGFGGGGGGGGGGGGGGGGGGGEGGGGGGFGGGGGGGGGGGG RG:Z:SRR062634 XC:i:99 XT:A:U NM:i:1 SM:i:23 AM:i:16 X0:i:1 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:98C0 XA:Z:GL000244.1,+26258,100M,2; +SRR062634.8061040 147 22 16051199 60 100M = 16051146 -153 TATATTTTGGAGTAAAATACTTTCATTTTCCTTCAGAGTCTGCTGTCGTATGATGCCATACCAGAGTCAGGTTGGAAAGTAAGCCACATTATACAGCGTT DDEDEFBBEBEDDFDFDF?EFEDE?DFBFGGE?EEGGGFGGGBGGGGGGGGGFDGGGGGEGGGGGGGEEGGGGGGGGGGGGGFGGGFFGGGGGGGGGGFG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.8108915 163 22 16072285 0 100M = 16072325 140 TATTACTGCATTGTTTGAACTTAAATGTCCACAACAGGATAACAGGATAATTAAACAGTTAATTGCACTATGGCACATCCATCTGTTAACAAAATGAAAT EGGGGGFGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGFGGGGGDGFGGGGGGGDFEGGGEGGGGGGFGGEGFGGGEEEGEEGEEEEEEAEEEGGEEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19770582,100M,0; +SRR062634.8204793 163 22 16072492 0 100M = 16072604 182 AAATGTACACAATACACTGTAATATCCAACACATACACCTATCAAAGTAAATATATTTTATTCTAGGGTAATTATTAAAACATTAAGAAATACTCTTTAA GGGGGDGGGGGFGGGGGGGGGGGGGGGBGGGGGGGGGGGGGFGGGGGEGGGGGGFGGGGGGGGGGE?GEEEGFGGGFDFGEFGEFEEEEGFEEBFECEEE RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19770375,100M,0; +SRR062634.8239925 163 22 16077961 0 100M = 16078030 169 TATTTAAATTCCATGCCTTTCATTTTAAAGTTTAGGAAATGCTTATAGCTTTCCAAACTGACATTTGTACTAACACAGTTAATCTAATCGAGTGATAACA GGGGGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGEGGGGGGGGGGGGGGGGGFGGGGGFGGGGEEGGF=FGEGGGEEGGGGGFEFEB=EAEEEEEGE RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:66A33 XA:Z:14,-19764906,100M,1;GL000244.1,-35846,100M,1; +SRR062634.8272386 147 22 16066698 0 100M = 16066605 -193 TGTTTTGATTTTACAGGCGCATCAGTGGAAGGGACTTGCCTTGTCTCAGAAGAAACTGCACTTTCGAGTGATGCCAAAACAGGTTGAGACTTTTGGGGGA EEEEEDDCDCC:CBCEDDDCBD:DEEEEFEGEEEGGEGGGGEGEGEEGGGGGGGEFFFFDFFGEGGFGGGFFFFBGGGFGGGGGEGGGFGGFGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:2 XM:i:1 XO:i:0 XG:i:0 MD:Z:75G24 +SRR062634.8303717 163 22 16066683 15 100M = 16066806 223 TAGAAGTAAATAACTTGTTTTGATTTTACAGGCGCATCAGTGGAAGGGACTTGCCTTGTCTCAGAAGAAACTGCACTTTCGAGTGATGCCGAAACAGGTT GGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGEGGEEGGGEGGGGGGEEGEEFGEBDEGFEDEEEEEGFEEGE5FCDCCDCDFE;EECCB?? RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.8331942 163 22 16077084 0 100M = 16077147 163 TACAGCCAGATCAGGATCAATGTCCTCTCTCTGGAATAAATGCAACGACCACGGTTCTTTGAGCATTATTTATCTTAAATTACTGTTTTAATTAGAAATG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGDGGGGGGGGGEGGGGGFGFDFEGFGGFEEGEEGFFGGEEEGGEEEEGEFFGEEGG?EEEEFBEEEG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19765783,100M,0;GL000244.1,-36723,100M,0; +SRR062634.8396408 163 22 16063308 60 100M = 16063415 207 AAAATGCAAACTGGCTATCACTACCATTCTGATAATGGTATTAGCTGTCCTACGTACCTGACGACCTAATGCTTAACCTAATGTTCCTCTTTCCCACTTT DADEBDED:EDEDDEEABD?EEEBAEEDEEDBE?CEDDBADD?BEE?DEBEEEEDEBB@=:5?:==7??=???=:??>>::BFFDDBABB=EE?EBBEDFF@CC-:DEBBEECCBC>B?B?DEEDDDDEFF RG:Z:SRR062634 XC:i:90 XT:A:U NM:i:0 SM:i:17 AM:i:17 X0:i:1 X1:i:4 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 +SRR062634.84766 147 22 16091435 0 100M = 16091352 -183 CCTGAAGAACTGTGGGTCAGTTAAACCTCTTTTCCTCATAGGTTGCCCAATATCAGGTAGTTTTTTATAACAGTGTGAAAATGGACTAATACAAGGTCTT CBCC?EEEEEEEFEEEEEEGEGFGEGGEFGGGEGGFGGGGGGGGGGGGGGGGGGGFGGDGGGGFGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19751498,100M,0;14,-19824224,100M,0;GL000244.1,+22364,100M,0; +SRR062634.8504458 147 22 16089656 0 100M = 16089543 -213 ATGTGGGCTTTAGTCATAAGGAGCTTACCTAGGCCCTAAGTGGACATCAGGTCCCAGGTTGACACAATGAACCATGTAGAAGTCAGGCTCTAAGTAGACA <6?BA?C@4AC RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19788064,100M,0; +SRR062634.8829645 147 22 16055333 46 100M = 16055217 -216 GGCAACCCCAATTGCAGAAGTCTCTATCCATCTTGGGAGGTTGTCATCCCCTGAACCTTCCTAGAGCACATTGTCCAGGCCCCCAAAGCCTCCAAACTCA AAC??DDCEDDCDCD@EDEEEFBEEBDEE?DEDEGAGBEGGGEFEBEGFGGFGGEGGGGGFDGGGGGEGGGGGFGGGGDGFGGGGGGGGGFGGFGFGDGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:23 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19787576,100M,1; +SRR062634.8881472 147 22 16058330 37 100M = 16058255 -175 AATTTCTTAGAGACTGGTTAAATGGCTGCATTCTGTTCATGCCCTAGGAATTTGTGGAAGTTCAAACTTAAGAATGATGATTTAGAGTGTCTGATGGAAG ?EEEEDECCEBDCD:DDDCDED@DDFEGAGFGGGEFEGGGGGGEGEGEGGGGGGGFGGGGFFGFGGGFGGGGGGGGGGGGGEGGGAFFGFGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.88820 163 22 16053839 46 100M = 16053928 189 ATGGGGTGAGGCATGTTGCAGGTCGAGGACCAGGGCCATCTCACTGCCTGAGCCCATGGACTGGCTCAGGGGTCTGTCAGATGATTCTAGAGCTGAGTTG GGGGGGFGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGFGGEEFEGGGGAGGDDEFGFGEEEEGGG?EEDDCCE@DDA@ECDDD=C=CCAE*,=,:>5=331049@?)27:@############# RG:Z:SRR062634 XC:i:87 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:83A3 XA:Z:14,-19788775,100M,1; +SRR062634.8909118 147 22 16073892 0 24S76M = 16073760 -208 CCCACCCCCACAGGTGACCCTCACAACAACCCTGTGAGACCTGCAGAGTAGCGGCTCTCACAAAGGAGGAGTCGGGAATGTCAAGATTTTAACACCTTCA ########################??A@5@DDADCCABABCCC?=CBC=BE:EDGEGEGEGGGGGGGGGGGGFGFGGGFGGFGGGGGGGGGGGEGGGGGG RG:Z:SRR062634 XC:i:76 XT:A:R NM:i:3 SM:i:0 AM:i:0 X0:i:2 X1:i:2 XM:i:3 XO:i:0 XG:i:0 MD:Z:24C2G47T0 XA:Z:14,+19768999,100M,3;9,+68326811,100M,4;2,-132578235,100M,4; +SRR062634.8911912 163 22 16060491 48 100M = 16060556 165 TTTTTTCTTTTTCTTTTTCTTTTTCTTTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTCTTTTTGTTT GGGGGGEGGGGFDGGGFGDGGGGGEGGGGGAGGGAGGGDDGGEGGEAFFEEEE8EFE8FFE>EEE:DD@2CD>3DEC>DEB:DDC8ECBA@0>@A RG:Z:SRR062634 XT:A:U NM:i:5 SM:i:25 AM:i:23 X0:i:1 X1:i:1 XM:i:5 XO:i:0 XG:i:0 MD:Z:0C89T1G3T0C2 XA:Z:22,+16060485,31M2D69M,3; +SRR062634.9015131 147 22 16083581 0 100M = 16083536 -145 CAAGCCCACACCAGACTTGGAATCTGCTAACCACTTGATCTTGAATTTTGCCATCTCCAGAACTGAGAGAAAGAAATTTGTTTACAAGCTTCCTCATCTA BCCCDEEEGEFFEECFEEGEFFFGEEGFGFGGEGEDGGGGGGGGGEEDGGGGGEGEGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,+30225,100M,0;14,+19759349,100M,1; +SRR062634.9071348 163 22 16053649 22 84M16S = 16053708 146 TGCTCAGCCTCTTGTGGCTACAGCCCAGCAGGTCCCAGGTGGCACCAGCCAAGATGAAGGTGGCACTGAGGGCTGAAGTCTCCCTCCCCATGAAGGGATG GGGGGGGGGGFGGGGGGGGGEF:EFEE=5:?C:@CC=5@,?@@???:ABB-<:;:><>>?9>>==;27/5>=3,?6585905@?################ RG:Z:SRR062634 XC:i:84 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:10A54T18 XA:Z:14,-19789277,100M,2; +SRR062634.911914 147 22 16097980 60 7S93M = 16097970 -103 GCTTTTTTTTGTTGTTTTGTTTTGTTTTGTTTTGTTTTGTTTTGTTTTCTGCCACAGATGCCTCACCTCTCCTCTCTCAAACGTTAACTTCCCATCATGG #######B@@/;:?68=E?C8>?/86AE;?A/<;A<:<;@(21;;?:*B-:DD?DB>@>;CC5C:DDCDCB@,@4C:?CC?CC=C??EE=EEEEEEE=E? RG:Z:SRR062634 XC:i:93 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:93 +SRR062634.9119161 133 22 16062369 0 * = 16062369 0 TTATATGTGTTTTTAAACTAAACTAATTTTATAGGAAAAATAATTTCTTTCCTTCCCTGTTATATCAAATACAGCCTTTAGCTCAAGACACAAGTAATTC CCAC?BD?DCCD:CCACD5=DCDAB=AA=5DBC=BC:??BB85:CBBB<7;'25????? RG:Z:SRR062634 +SRR062634.9146674 163 22 16085888 0 100M = 16085960 172 AAAGTGCTGGGATTACAGGTGTGAGCCAATGTGCCTGGCCTGGTTACCAGAATTCTTAGTTCTGTTAGGGTCTGTTGGCAAGGAAGTGAGGTCGCTTCTT GGGEEGFGDGGGGDGGAGGGGGGGGFGEGGGGEGGFFEEFGGGGFEGFGGEEGGFGFEGEFGGGGBGGFGEFGGEGGGEBDEFCDGCFBFCACECDDCDD RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-27911,100M,0;14,-19757045,100M,1; +SRR062634.9263376 163 22 16078724 0 94M6S = 16078857 183 AGAGAGGAGAGAGAAGGAGAGAATCAGACAGATACAGAAAGAGAGTGACAGGGCTACAGAGACAGAATTAAAAGGGGGCACCAAGAATACTGAGCTGAGA >>;@@==-6BA-ACBA=CC?C5DDDD?-D5>-@C@:CC,C?C>5A??:>AC-A5@;@>C-BCD@:DBDABDB5=,>*?:??5A###### RG:Z:SRR062634 XC:i:94 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:94 XA:Z:GL000244.1,-35079,100M,0;14,-19764143,100M,0; +SRR062634.9308239 147 22 16059040 37 35S65M = 16058928 -177 ATGGGAAGGGCTGGGTGTTTAGGGAGAAGCCTGCTGCAGGGGTGGAGCCATCACAGAGAACCTCTCCTAGGGCAGTACCAAGTGGAAATGTGGGGTTGGA ###################################A=%9CCCCDC4CCC@=BCBECCCC5FFCEDFFEDGFEFFBFBF?FGGFGGGFGGGGGEGGEGDGG RG:Z:SRR062634 XC:i:65 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:65 +SRR062634.9333081 163 22 16077138 0 100M = 16077197 159 TTCTTTGAGCATTATTTATCTTAAATTACTGTTTTAATTAGAAATGTATCTAACTTACTTATATTAGCCAAATTTCCTTTCAATGTAACCAAATTTCTTT GGGGGGFGGGGGGGGGGFGGGGGGGGGGGGGGGGGGFGGGGGDGGGGGGGGGGGGGGGGGGGGGGGGGFEEGGGGGGGGGGFEGGEGGFGEEAGEGGGEG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19765729,100M,0;GL000244.1,-36669,100M,0; +SRR062634.9341717 163 22 16065959 20 99M1S = 16066039 180 TTCTTTAATTTCTTTTAATTCCTGTTCCGCTTGGACCAGGGCCAGGGCACCACTGCCCTGTACCACCCCCAGGAGGGTGTTCCTTGCATCCTGGCCCTTC GGGGGGGEGFGGGGGGGGGGGGGGGGGFGGGDGEGGEGEEAEEEEBEFEGFEGEGGGGEDABCEEEDEEE2CBAA??0B?C@CC5=BAAAAA95=AA??# RG:Z:SRR062634 XC:i:99 XT:A:U NM:i:0 SM:i:20 AM:i:0 X0:i:1 X1:i:2 XM:i:0 XO:i:0 XG:i:0 MD:Z:99 XA:Z:14,-19776914,100M,1;2,-132543403,100M,1; +SRR062634.9369051 163 22 16081516 0 100M = 16081643 215 TAAAAGTCACACTTCTTAAAGTTGCAGTCTCGCCAAGGCAACACAATGTAGCAGTCTCTCTTGTGAGGTATCACCTGGAGTTCTTCATCTCACCACCAAG GGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGEEGGGGGGEGEFEGEGEEEEEEEEEEEE?CCC>CCCCDDCDCBDBCCC@CBCCBEEEEEEEEEB RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-32290,100M,0;14,-19761414,100M,1; +SRR062634.9370249 163 22 16055029 37 100M = 16055164 200 TATGATTAGGCCTGAGAGGTTGTTTGTCTCCCATCAGTGGAGCTGGCAATCTGACGATGCCTCTCTTTGGACCCTGGCCTTCAGCCTCTTCAGGCCATGG GGGGGGGGGEGGGGGGGGGEGGFGGGGGGGGGGGFGGEEGEGGFFGGFFGGFGGGGBEGGGGGGFGFGEACFEG?A=DEBEE=EEE?D?AC<@BCD=C?C RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:37 AM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 +SRR062634.938102 147 22 16073397 23 100M = 16073285 -212 CATGGAGCAGCATGGGGACAGGGGCTGGCCAGCTGGGCAGGGCTCACACTCCTGACACCCAGACTCCACGATACTCCTCTGTCCCCACCCAGGGCAGATC ;3.6;=59:?>;2ACCDCBCAEEEEDDBDD:FEGEGGF@EGEFEDEGEGEGGEFEGBDEFGGEGGGDEGGGDEGFGFGGEGDGFGFGGEFGGGGGGGGGG RG:Z:SRR062634 XT:A:U NM:i:0 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19769470,100M,1; +SRR062634.9399882 147 22 16052983 0 100M = 16052904 -179 ATACCTGCCAGATGAGTCACTGGCAAAAGGTGCTGCTCCCTGGTGAGGGAGAAACACCAGGGGCTGGGAGAGGCCCAGAAGGCTCTGAAGGAGTTTTGGT 9>=CC@CAA?A<>ABA=C?D?BC:DEEEECECBD?ECEEEDEG?GEGEEGFEGEFEGEEEGGFGGGGGGFGGGDGGGGGGGGGEGGGGGGFFFDFGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19789927,100M,0; +SRR062634.9401900 163 22 16063058 23 100M = 16063141 183 GACCCAATCTCTACAAAAACTAGAGAAAAAATAGCTGGGTGTGGTGGCACTCACCTGTACTAAGGAGGCTGGGCTGGGAGGATCTCCTGAGCCTAGGGGG GFGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGDEGEGGDGGGGGEGGGFFGEGGGGGEFEEEFCEAEECGGE=EEC5EADECABCBDAAC@5?B RG:Z:SRR062634 XT:A:R NM:i:2 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:2 XO:i:0 XG:i:0 MD:Z:96T0T2 XA:Z:2,-132546300,100M,2;14,-19779814,100M,3; +SRR062634.9476947 147 22 16096797 23 65S35M = 16096625 -207 TAGCCTTCCCCCCCCCCAGGGCGCCCCCCAGGGGGCTTAACCCCCCCCCTGTTCTCGAACCAGGGATGCCAGGGTCCCCTGGGCCCAGCTCAGGGGCTGA #################################################################@=>-@-@B?@A:?=EAEEFFEEE5C?AE:CEEFAF RG:Z:SRR062634 XC:i:35 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:10 X1:i:5 XM:i:0 XO:i:0 XG:i:0 MD:Z:35 +SRR062634.95155 147 22 16074504 0 100M = 16074417 -187 AAACTGCCGCCATGATTCAGTTACCTCTCACTGGGTCCCTCCCGTGACACATGAGGCTTATGGAAATTACAGTTCAAGATGACATTTGGGTGGCAACACA =>:?,ADDD?EEEEECACCCB@)BCCECECE?EGFEEGGEEEFFEGDGCFGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGFGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19768363,100M,0;GL000244.1,+39303,100M,0; +SRR062634.9524762 163 22 16098596 23 52M1I23M24S = 16098697 201 CGTCTGTGAGGGACCCAGGCCCCACTCCTCTGCACCCGAGGTTTTTGTTTTGTTTTTTTTTTTTTCTCTGCCCCAGGTGCCTCCCCAGCTCCTTGGGGGG GGEGGGGGGGGFGGGFGGGGGGGGGFGGGGFFGEGGGG?FFCEFFFD@CCCAFDGEGEFGEGEEG*;-=-=?@@(?######################## RG:Z:SRR062634 XC:i:76 XT:A:U NM:i:1 SM:i:23 AM:i:0 X0:i:1 X1:i:1 XM:i:0 XO:i:1 XG:i:1 MD:Z:75 XA:Z:GL000244.1,-15229,11M1I64M,2; +SRR062634.9559653 147 22 16097853 29 17S83M = 16097725 -211 CTGCCCGTCGTGGGGTACCCCTAGCGGGCCCAAGGCGCACCCGGGGCTCGAACCAGGGTCCACAGGGCCCAGCGCAGGGGTTGATGGGAAGGCATTTTCA #################?:@C?4DCAEEEEC@CDEDEE@BEBEFGEEFEGGEGGGGGGFGGFGGGGGFGGEGGGFGGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XC:i:83 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:4 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:82C0 XA:Z:14,+19745394,100M,1;14,-19829932,100M,1;GL000244.1,+15958,100M,1; +SRR062634.9625653 163 22 16068337 0 91M9S = 16068404 167 TAATATATATTTTATATGAGCATTCTTTTATTCAAACCATCCTTTTTTCCTAGAATGATGGATTAGGGGCTTTCGGTGTGCCTCAGCCACTTGGAAGTAG GGGGGGFGEGGGGDGGGGGGGGGGGGGGGFGGGEGGGFGGGGGGGGGGGGFGFGGGEDGDFEGDGGEEED5CC@:C@:?B:?>?:B@BAB?######### RG:Z:SRR062634 XC:i:91 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:91 XA:Z:14,-19774544,100M,0; +SRR062634.9676388 147 22 16087636 0 10S90M = 16087527 -199 CGTGTTAATGTGGCCTAATGTACTCCTGGGTTCTAGGGTCCTCTTGGGACCTGATGTCTACCAGGATCCTGGTATCCACCTGGGGCCTGGTATCCACCTA ##########@@?@?=;3<=8&318?-?A@=??:5AA:?5EC:ECBB@,?ABA@A=A@=5C@5C@=:F?FEEEEDDE?EEFGGGFEDDGFFGGGFFFDFF RG:Z:SRR062634 XC:i:90 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:4 X1:i:3 XM:i:0 XO:i:0 XG:i:0 MD:Z:90 +SRR062634.9683698 163 22 16071154 0 100M = 16071262 208 TATGCAGCCTTCTGTCACTGTGGGAACCTGGCTTATGGATGAGGCTATGACACAGGGTTGTGAAGTTATAGCTATGTGTTTGAGAACAAGGAAGGCAGTA GGGAGGGGGFGGGGGGGGGGGGGDGDGGGGGGEGGFGFBGFFGGFGGDGFGGFEDFFCEF:CBBCBEBEAEGGEGEEEBEEE=EEEBEEE@?CAB5,?:? RG:Z:SRR062634 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:1 XO:i:0 XG:i:0 MD:Z:0C99 XA:Z:14,-19771717,100M,1;9,-68329537,100M,2; +SRR062634.9716434 147 22 16062768 0 100M = 16062661 -207 CAGGGAAATATTCACTCCTGTTTTCTTCTCTCCCCAATCTCAGCCTTTCGGTGACTGTTTGGGCAAAGCCTCCCTTGTGGTAAAAGATGCCTCACTTCTG EEEEECCACC@AC>CEFDEE@ACBDEEF?EEEEGFCCAFCFEGEEEDDEDGGGEGEEGGGGGGGFGGFGGFFGFGGGGGGGGGFFGGGGEGGGEGGGFGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19780104,100M,0; +SRR062634.9741621 147 22 16064712 0 100M = 16064620 -192 GGATAGGCTGATGCGGTCTCAGATGGAAATGAGGAACTTACTGGGACTTCGAAAAAAGGTCATGCCTGTTATGCCTTTGCAAAGAACTTGGCTGCAGTCT EEBEEEEEEEEEEEEEFBEEEGGEEFEFEEGGGEGEGEGEGBGGGFGEGFGGGGGGGFGGGGGGGGGGFGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19778160,100M,0; +SRR062634.9744609 163 22 16096842 0 61M39S = 16096978 216 CTTTTGTCCGTGGGGGACACAGGCCCCGCTTCTCCGCGGCAGGGGTTTTTTTTTTTTTTTTTTTTTCCCCCTCCGCGGCGGGGGGTTTTTTTTTTTTTTT GGGGGGGGGGGGGGGGEGGGGGGGGGGFGGFGEGGFGF?ECEEE*ACABBDEEEGGGGEE@####################################### RG:Z:SRR062634 XC:i:61 XT:A:R NM:i:1 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:44T16 XA:Z:GL000244.1,-16993,100M,1; +SRR062634.9796441 147 22 16058918 22 15S85M = 16058782 -221 CAAGCTATAAGCCTCGGAAACTTCCATGTGGTGTTAAGCCTATGGGCATGCAGAGTGCAAGAGTGATGGCTTGGCATTCTCTGCCTATATTTCAGAGGAC ###############A??=;=48?A?AAAC>=>:A,AAB:?@?C@@CBB>?=@5F=8CFGGEGGGGGAGFGEEGF?FD?GFGGGDFCEBFGGGGGFGEDFBDDGFFEDFAGFGGFGEGG RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,+19757247,100M,0;GL000244.1,+28113,100M,1; +SRR062634.9961917 163 22 16065519 0 100M = 16065565 146 CGATGGGGTTGAGGAGGGTATTTATTATTTAGGTGCACCAGCCCAGTCAGATTAACATCCAAAGGACTGAGTCCTGAACAAAGACTTAAGTTACCTTTTA GGGGGGGFGDEBEEBEEDBEFFEAFGBGEGEEFEFGGDGFFDFEBEBEBFABD?EDFEBFEEDBABEDDDBA@ABBBAADD?AA?BB@BDAA?=BABAC? RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:2 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:14,-19777353,100M,0; +SRR062634.997014 163 22 16081633 0 100M = 16081734 201 ATACACGGGTGAGGTGGGAGTGAAAGTTTAATAAGCAAAAGGAGGAAGCTCTCTGCAGCAGACAAGGGCGTCCAAGTGGATTGCCGTTTTTACAGTTGAA GGGGGGGGGCEEEECFFFEFDGGGGGEGGGGGFGGGGGFEGFAFGEGGFEEEGGGGGGGGGGEG?DEEGECDDD:CACC?CDBBEBA?CCC;C7C?@?:; RG:Z:SRR062634 XT:A:R NM:i:0 SM:i:0 AM:i:0 X0:i:3 X1:i:1 XM:i:0 XO:i:0 XG:i:0 MD:Z:100 XA:Z:GL000244.1,-32173,100M,0;14,-19761297,100M,0;21,-10617533,100M,1; diff --git a/adam-core/src/test/scala/org/bdgenomics/adam/rdd/fragment/FragmentDatasetSuite.scala b/adam-core/src/test/scala/org/bdgenomics/adam/rdd/fragment/FragmentDatasetSuite.scala index 5fe63efc06..d97951d6f2 100644 --- a/adam-core/src/test/scala/org/bdgenomics/adam/rdd/fragment/FragmentDatasetSuite.scala +++ b/adam-core/src/test/scala/org/bdgenomics/adam/rdd/fragment/FragmentDatasetSuite.scala @@ -335,8 +335,8 @@ class FragmentDatasetSuite extends ADAMFunSuite { } sparkTest("union two genomic datasets of fragments together") { - val reads1 = sc.loadAlignments(testFile("bqsr1.sam")).toFragments - val reads2 = sc.loadAlignments(testFile("small.sam")).toFragments + val reads1 = sc.loadAlignments(testFile("bqsr1.sam")).toFragments() + val reads2 = sc.loadAlignments(testFile("small.sam")).toFragments() val union = reads1.union(reads2) assert(union.rdd.count === (reads1.rdd.count + reads2.rdd.count)) // all of the contigs small.sam has are in bqsr1.sam @@ -636,4 +636,18 @@ class FragmentDatasetSuite extends ADAMFunSuite { assert(convertedRdd.rdd.collect().toSet == convertedDataset.rdd.collect().toSet) } + + sparkTest("Compare unpaired FragmentDataset with optimisation to non-optimised") { + + val alignmentDataset = sc.loadAlignments(testFile("unpaired.sam")) + + val merged = alignmentDataset.toFragments() + + val notMerged = alignmentDataset.toFragments(false) + + assert(merged.rdd.count === notMerged.rdd.count) + + assert(merged.toDF().except(notMerged.toDF()).count() === 0) + + } } diff --git a/adam-core/src/test/scala/org/bdgenomics/adam/rdd/read/MarkDuplicatesSuite.scala b/adam-core/src/test/scala/org/bdgenomics/adam/rdd/read/MarkDuplicatesSuite.scala index daeaf2011d..273cf5e4e6 100644 --- a/adam-core/src/test/scala/org/bdgenomics/adam/rdd/read/MarkDuplicatesSuite.scala +++ b/adam-core/src/test/scala/org/bdgenomics/adam/rdd/read/MarkDuplicatesSuite.scala @@ -208,7 +208,7 @@ class MarkDuplicatesSuite extends ADAMFunSuite { private def markDuplicateFragments(reads: Alignment*): Array[Alignment] = { AlignmentDataset(sc.parallelize(reads), SequenceDictionary.empty, rgd, Seq.empty) - .toFragments + .toFragments() .markDuplicates() .toAlignments .rdd From d6a66783de76e23a184fbb2760a3d7d5ea79a7fb Mon Sep 17 00:00:00 2001 From: "Rahamim, Ben" Date: Wed, 18 Nov 2020 21:04:22 +0200 Subject: [PATCH 2/3] added mergePairedReads as a parameter to Python API --- adam-python/bdgenomics/adam/rdd.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/adam-python/bdgenomics/adam/rdd.py b/adam-python/bdgenomics/adam/rdd.py index 7c7ffa9111..6f8fc21dab 100644 --- a/adam-python/bdgenomics/adam/rdd.py +++ b/adam-python/bdgenomics/adam/rdd.py @@ -870,16 +870,18 @@ def _inferConversionFn(self, destClass): return "org.bdgenomics.adam.api.java.AlignmentsTo%s" % self._destClassSuffix(destClass) - def toFragments(self): + def toFragments(self, mergePairedReads=True): """ Convert this set of reads into fragments. + :param mergePairedReads: Merge paired reads into a fragment. Can be set to false when only unpaired + reads are in the AlignmentDataset. Defaults to true. :return: Returns a FragmentDataset where all reads have been grouped together by the original sequence fragment they come from. :rtype: bdgenomics.adam.rdd.FragmentDataset """ - return FragmentDataset(self._jvmRdd.toFragments(), self.sc) + return FragmentDataset(self._jvmRdd.toFragments(mergePairedReads), self.sc) def toCoverage(self, collapse = True): From e6332ef5bac0dc6ce37cb840631c62f3e67a8351 Mon Sep 17 00:00:00 2001 From: "Rahamim, Ben" Date: Thu, 19 Nov 2020 22:30:11 +0200 Subject: [PATCH 3/3] code format fix --- .../scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala index 1ff459da50..889fb302bd 100644 --- a/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala +++ b/adam-core/src/main/scala/org/bdgenomics/adam/rdd/read/SingleReadBucket.scala @@ -164,7 +164,7 @@ private[adam] case class SingleReadBucket( val insertSize = primaryMapped.headOption match { case Some(alignment) => Option(alignment.getInsertSize).map(_.toInt) - case None => None + case None => None } SingleReadBucket.buildFragmentFromAlignments(unionReads.toSeq,