BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,445,981	T→G	22.0%	G34G (GGT→GGG)	ydbL →	DUF1318 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,445,981	0	T	G	22.0%	184.7 / 7.6	98	G34G (GGT→GGG)	ydbL	DUF1318 family protein
Reads supporting (aligned to +/- strand): ref base T (17/55); new base G (23/3); total (40/58)
Fisher's exact test for biased strand distribution p-value = 7.45e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGAAGTCGCTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCG  >  NC_000913/1445746‑1446217

catttgaGTCGGTGCCCCTAAGGAACC‑GATCCCTATCAATATGAAAGTTAAAATTGGGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGaa                                                                                                                                                                                                                                                <  2:1352238/233‑1 (MQ=255)

acttGAAGGCGCTGCCCCAAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGGGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGaa                                                                                                                                                                                                                                                <  2:849858/237‑1 (MQ=255)

        ttgccgccccTAAGGAACC‑GATCCCTAACAATATGAACGTTAAAATTGGGCATGAAATCATCATCAAGGGAGACAAAGATGTCGAAGAGCTGCTGGAAACTCGTAGCGATCTTTTCTGAGGGGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTAAGGGGCGGGTA‑GGTGAAACATTTTa                                                                                                                                                                                                                                        <  2:1328360/231‑1 (MQ=255)

        tcccTGCCCCGAAGGAACCAGAT‑ACTATAAAATTGAACGTTAAAAGTGAGCTTGGGATCATCATAAAGGCAGACAAAGATGGCGAAGAGCTGCTTGAAACTCGTAGCGAGCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTa                                                                                                                                                                                                                                        <  2:1209923/236‑1 (MQ=255)

           cTGCCCTTAAGGACCC‑GATCACTTTAAATTGGAACGTTAAAATTGAGCATGAGACCAGCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACgg                                                                                                                                                                                                                                     <  2:378904/239‑1 (MQ=255)

           cTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGAAAAAATTTAAgg                                                                                                                                                                                                                                     >  2:670346/1‑239 (MQ=255)

           cTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGAAGAAA‑TGTAAGGGCATTGACTCTGGTTGAAGACATAACTCAGGGGCGGGGG‑GGGGAAAAAATTTAAgg                                                                                                                                                                                                                                     >  2:1365667/1‑239 (MQ=255)

                    aGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTAAGG‑GTATC‑TGGt                                                                                                                                                                                                                          >  1:221483/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGACAACATTACTTCTTTGTGCGTTCCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGAATGAAGCCAGAACTCAGGGGCGGGGA‑CGGGGAACAATTTCCTG‑TTTTT‑TGGTTTCGCTg                                                                                                                                                                                                                   >  1:933609/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑AGGGAAACATTTTACGG‑GTTTC‑TGGTTgcgccg                                                                                                                                                                                                                   >  1:596447/1‑237 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGGAACATTTTACGGTTTATC‑GGGTT‑CGCTg                                                                                                                                                                                                                   >  1:1263485/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCGGGTTGAAGCCTGACCTCAGGGGCGGGGA‑GGGGAAAAAATTTACGG‑TTTTA‑Tggggggggtg                                                                                                                                                                                                                   >  2:493863/1‑231 (MQ=255)

                             aTCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGGAAAACTTTTCGG‑GTATT‑TGGGTgggcggg                                                                                                                                                                                                                  >  1:1266558/1‑236 (MQ=255)

                                    tCAATATGAACGTTAAAATTGGGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGa                                                                                                                                                                                                           <  2:987609/239‑1 (MQ=255)

                                    tCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGT‑GGGGAAAAATTTTTTCG‑TTATCTTGTTTGCCTTGAAAACga                                                                                                                                                                                                            >  2:257992/1‑238 (MQ=255)

                                        tatGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGTGAAAAATTTTATGG‑GTATC‑TGGGTGGGCGGAAAACgggaggg                                                                                                                                                                                                       >  2:1218021/1‑234 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAAACA‑TTTACGG‑TTATC‑TGGCTGCGTCTAAAACggcggattgaa                                                                                                                                                                                                   >  1:717150/1‑230 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACAATTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGaaa                                                                                                                                                                                                   >  1:1008632/1‑239 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAAAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGAGGTTGaaa                                                                                                                                                                                                   >  1:720989/1‑239 (MQ=255)

                                              cGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACa                                                                                                                                                                                                 <  1:798329/239‑1 (MQ=255)

                                              cGTTAAAATTGAGCATGAGATCAACATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACa                                                                                                                                                                                                 <  1:493863/239‑1 (MQ=255)

                                                  aaaaTTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAAAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGAGGCTGAAAAAgaga                                                                                                                                                                                             >  1:195014/1‑239 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGGTGAAACAAAGAAATtagt                                                                                                                                                                                      >  1:675787/1‑239 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑AGTGGAACATTTTCCGG‑TTATT‑TGGTTGCGCTGGAAACGGCGGGTGAAAAAGGAAAATtatt                                                                                                                                                                                      >  1:468183/1‑237 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGACGGGGA‑GGTGAAACAATTTTGGG‑TTCTA‑TGGTTTCGCTGGAAACAGGTGCTGGAACAGAAAAATtaat                                                                                                                                                                                      >  2:641575/1‑237 (MQ=255)

                                                           cATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGGAACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACCGGTGCTGGAACAGAGAAATtaggag                                                                                                                                                                                    >  1:841768/1‑239 (MQ=255)

                                                           cATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGCA‑GGGGAAACATTTTACGG‑TTATC‑TGGTTGGCCTGAAAACGGATGCTGAAACAAGGAAAttatta                                                                                                                                                                                     >  1:1069276/1‑238 (MQ=255)

                                                             ttatatcttcatcAAGGCAAGATAAGGTGTGGAAGAGCTGCTTGAAACTCGTAGCGATTCTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                  <  2:118537/235‑1 (MQ=255)

                                                             tggattactaatcaAGGCAGACAAAGATTTGGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                  <  2:4335/229‑1 (MQ=255)

                                                             tgaggccattatcaAGGTCGAAAAAAGTGTCGAAGAGCTGCTTGAAACTCGGAGCTGTCTTTTGTGAGGTGATGGTGAAGAAAACATTGCTGCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                  <  2:1011458/233‑1 (MQ=255)

                                                                   catcatcaAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAAAATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACAtt                                                                                                                                                                            >  1:165324/1‑239 (MQ=255)

                                                                         caAGGCAGACAAAGATTGCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGcc                                                                                                                                                                      <  1:670346/239‑1 (MQ=255)

                                                                          aaGGCAGACAAAGTTTTCGAAGGGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGGGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCg                                                                                                                                                                     <  2:720989/239‑1 (MQ=255)

                                                                            ggCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACCTTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAACAGAGAAATTATTAGCCGACATTAATTCCGa                                                                                                                                                                    >  1:886816/1‑238 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGAAACTTTTTCGG‑GTATT‑TGGTTGCGGTGAAAACGGATGAGGGAAAAGAGAAATTGAAAGGCGACAATAATGCAGaaaggaa                                                                                                                                                              >  2:652871/1‑234 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGAAAATTTTACGG‑GTGTC‑TGGTTGGGTCTAAAAAGGTAGCTGAAAAAAAGAAATTATGAACCCACATTAAAGCCGaaagaaa                                                                                                                                                              >  2:405977/1‑234 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTTCGG‑TTATT‑TGGTTGGGCTGAAAACGGATGCTGAAACAAGGAAATTTGTAGCACACATTAATAGCGAACGTaa                                                                                                                                                              >  1:1025253/1‑239 (MQ=255)

                                                                                    aagttttCGAAAGGCTGCTTGAAACTCGTGGCGATTTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGc                                                                                                                                                           <  2:675787/233‑1 (MQ=255)

                                                                                        tGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAg                                                                                                                                                        <  1:1298300/238‑1 (MQ=255)

                                                                                        tGTCGAAGAGCTGCTGGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAATCATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGGCTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTTAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATTCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAg                                                                                                                                                        <  1:876687/238‑1 (MQ=255)

                                                                                           cGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGaa                                                                                                                                                                                                                                                >  2:835536/1‑149 (MQ=255)

                                                                                           cGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGaa                                                                                                                                                                                                                                                <  1:835536/149‑1 (MQ=255)

                                                                                                  ctgctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGGTGCTGAAACAGAGAAATTAGTAGCCGACCTTAATGCCGAACGGAAAGCGGG‑TTTCcaaaa                                                                                                                                              >  1:813536/1‑235 (MQ=255)

                                                                                                  ctgctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAAAATTTTTCGG‑TTATC‑GTGTTGGGGTGGAAACGGATGCTGAAAAAGAGAAATTAGTAGCCAGCATTAATTCCGAACGGAAAAGGAG‑GTAAcaaaaaa                                                                                                                                            >  2:220162/1‑238 (MQ=255)

                                                                                                    gctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGATGCGCTGAAAAAGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAAGGAG‑TTACCAACAACTg                                                                                                                                          >  1:508135/1‑239 (MQ=255)

                                                                                                     ctTTAAACTCGTTGGGATCTTTTCTGGGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTgg                                                                                                                                         <  2:1008632/239‑1 (MQ=255)

                                                                                                       tGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                       <  1:1271321/239‑1 (MQ=255)

                                                                                                         aaaCCCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgcaaa                                                                                                                                     <  1:120768/239‑1 (MQ=255)

                                                                                                            cTCGTAGTGATTTTTTCTGAGGTGATGATGAAGAAAACATTATTTCTGTGTGGGTTTCTTGTTGGTCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgcaaagca                                                                                                                                  <  2:645542/239‑1 (MQ=255)

                                                                                                               gTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAACAAAGAAATTAGTAGCCGGCATTAATGCCGAACGAAAAGCGAG‑TTTCCAACAAATGGCAAAGCaaaa                                                                                                                               >  1:1155766/1‑239 (MQ=255)

                                                                                                                   ggATCTTTTCTTGGGTGTTGGTGAAGAAAACATTCTTTTTTTGTGGGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAaataat                                                                                                                           <  2:877050/238‑1 (MQ=255)

                                                                                                                       tttttttGGGGTGATGATGAAGAAAACATTTCTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                       <  2:1069276/238‑1 (MQ=255)

                                                                                                                       cTTTTCTGAGGTTATGATTAAGAAAACATTACTTCCTTGGGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTTTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                       <  1:1255731/239‑1 (MQ=255)

                                                                                                                       cTTTTCTGAGGTTATGAGTAAAAAAACATTTCTTCTTTGTGCTTTTCTTGTTGGGGTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGAATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                       <  2:1132305/239‑1 (MQ=255)

                                                                                                                       cTTTTCTGAGGTGATGATGAAGAAAACATTCCTTCTTTGGTCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                       <  1:504649/239‑1 (MQ=255)

                                                                                                                         tttCTTAGGTGATTATGAAAAAAACTTTGCTGCTTTGGGTGTTTTTTGGTGGGCTGGTAAGCAGGAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACGCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAAGTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCg                                                                                                                     <  2:933609/239‑1 (MQ=255)

                                                                                                                              ggggTGAGGATGAAGAAAAACTTACTGCGTTGGGAGTTTCTTGGTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAGGTGTCGGTAga                                                                                                                <  1:962407/237‑1 (MQ=255)

                                                                                                                                   gaggatgaAGAAAACATTACTTTTTTGTGCGTTTCTTGGTGGGCTGGTAAGCAGCAA‑TGTAATGGCATGGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAgatgat                                                                                                            <  1:405977/235‑1 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAAAAAAGAAATTAATAGCCGACATTAAATCCCAAACTAAAGCGCG‑TTTCCCACAACTGGCAAAGCCAAATCATATGTCGGGAGATGata                                                                                                           >  1:895763/1‑238 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGAAAATTTTTCGG‑GTATA‑TGGTTGCGCGGGAAACGGGTGCGGAAAAAGGGAAATTAGAGACCTGAATTTATGCCCAAAGTAAA‑AGAGATTACAAACAACCGGGAAAAAAAAAAAATGAGTCTGAagaagaaaa                                                                                                          >  2:363972/1‑236 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTACGG‑TTTTA‑TGGGTGCGCTGGAAACGGATGCTGTAACAGAGAAATTAGTAGCCGACATTAATGCCTAACGTAAAGTGAG‑TTACCAACAACTGGCAAAAAAAAAAAAAGGGTGGGGagaagaaaa                                                                                                          >  2:286656/1‑236 (MQ=255)

                                                                                                                                       atgaAGAAAACATTTCTTCGTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑GGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCg                                                                                                                     <  1:767711/225‑1 (MQ=255)

                                                                                                                                                       ttttttGTGTTTTTTTTGTTGGGGTGTTAAGGAGCAA‑TTTAATGGAATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCaa                                                                                       <  2:165324/236‑1 (MQ=255)

                                                                                                                                                        tttttGTGGTTTTTTTTTTGGGCTGGTAAGCCGAAA‑TGTAATGGCATTGGCGTTGGATGAAGCCAGAACGCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCaaa                                                                                      <  2:657501/237‑1 (MQ=255)

                                                                                                                                                        tgtttGTGGTTTTCTTGTTGGGCTGGTAAGCGGAAA‑TTGAATGGAATTGACTATGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCaaa                                                                                      <  2:510758/237‑1 (MQ=255)

                                                                                                                                                           ttGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCCCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGGAAAGCGAG‑TTTCCAACAACTGGCAAAGCAAAAAAATGTGTCGGGAAGTGAAATCGCGAAAATCCCCGGGCcaaaa                                                                                    >  1:657501/1‑237 (MQ=255)

                                                                                                                                                            tgtgCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAAAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAAAGTAAAACGAG‑TTACCAACAAATGGCCAAGCAAAATAATGTGTCGGGAGATGATATCGCGAAAATCGCCCGGCAAAAGc                                                                                   >  1:1318953/1‑238 (MQ=255)

                                                                                                                                                                   ttCATGGTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCt                                                                          <  1:363972/239‑1 (MQ=255)

                                                                                                                                                                        gttgGGCGGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAACCAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGc                                                                    <  1:1149033/239‑1 (MQ=255)

                                                                                                                                                                            ggCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAg                                                                                                                                                        <  1:934260/154‑1 (MQ=255)

                                                                                                                                                                              cTGGTAAGC‑GCAATTGTAATGGCATTTGCTCTGGAAAAAGCCAGAAAGCAGGGGGGGGTA‑GGTGAAACATTTTACGG‑TTAAC‑TGGTTGAACTGAAAACGGATGCTGAAACAGAGAAAGTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAAc                                                               <  2:1169529/238‑1 (MQ=255)

                                                                                                                                                                                 gTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCggg                                                           <  1:220162/239‑1 (MQ=255)

                                                                                                                                                                                  tAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCggg                                                           <  1:257992/238‑1 (MQ=255)

                                                                                                                                                                                  tAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCggg                                                           <  1:652871/238‑1 (MQ=255)

                                                                                                                                                                                  gaaGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCggg                                                           <  1:42869/237‑1 (MQ=255)

                                                                                                                                                                                    aaggcgaaagttAATTGAATTGCTTCTGGGTGAATCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTACGG‑TTAAC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCggggg                                                         <  2:508135/228‑1 (MQ=255)

                                                                                                                                                                                      cagcaA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGAAACATTTTACGG‑GTATT‑TGGGTGCGCTGAAAAAGGAAGCTGAAACCGAAGAATTTGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAAAAAATGGCAAAGCAAAATAATGTGACGGGATAGGAGAAGGCGAAACACGCCAGGCAAAAACAC‑GCATGTC‑GGGCCAAAACggggggaa                                                      >  2:716215/1‑236 (MQ=255)

                                                                                                                                                                                       agcaA‑TGTAATGGCATTGACTCTGGACGAAGCCAGAACCCAGGGGCGGGTA‑GGTGAAACCTTTTACGG‑TTATC‑TGGGTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAat                                                      <  1:1365667/238‑1 (MQ=255)

                                                                                                                                                                                             tGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAat                                                      <  1:1275587/233‑1 (MQ=255)

                                                                                                                                                                                               tAAGGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACaa                                              <  1:1218021/239‑1 (MQ=255)

                                                                                                                                                                                                 aTGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTTCGG‑TTATC‑TGGTTGCGCTGAAAACGGGTGCTGGAACAAAGAAATTATTATCCGACATTTATGCCGAACCGAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAATTTTCGGTAGATGATATCGCGAAACTCGCCGGGGAAAAAACC‑GTAACGCGGGGCCAACCCGGGGGaaatagtaaaag                                             >  2:309841/1‑229 (MQ=255)

                                                                                                                                                                                                    gCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTTCGCTGGAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATTCCGAACGTAAAGCGGG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAAATGATAACGCGAAAACTGCCGGGCCAAAGCCC‑GTAGCTC‑GGGCCCAACCCGGGGGAAACGTACAAGAGAt                                         >  1:753167/1‑239 (MQ=255)

                                                                                                                                                                                                      gttGACTTTGGATTAAGCCCAAACTTCGGGGGGGGTG‑GGTGAAACCTTTTACGG‑TTAAC‑TGGTTGCGCTGAAAACGGATGCTGAAACAAAGAAATTAGTTGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAGCTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                                        <  2:596447/237‑1 (MQ=255)

                                                                                                                                                                                                      aTTGTCTCTGGATGAAGCCGGAACTCGGGGGGGGGTA‑GGTGAAACATTTTAAGG‑TTATC‑TGGTTGCGCTGAAAACGGATGTTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                                        <  2:195014/238‑1 (MQ=255)

                                                                                                                                                                                                      aTTGGCTGTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTGTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                                        <  1:641575/238‑1 (MQ=255)

                                                                                                                                                                                                      aTTGGCTCTGGATGAAGCCAGAACCAAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTTTA‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                                        <  2:895763/238‑1 (MQ=255)

                                                                                                                                                                                                      aTTGACTTTGGATGAGGCAAGAACTAAGGGGCGGGTA‑GGTGAAACTTTTTACGT‑TTATC‑TGTTGGCGCTGAAAACGGATGATTAAACAGATAAATTAGTAGCCGACATTAATGCCGAACGGAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGATC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                                        <  2:986427/238‑1 (MQ=255)

                                                                                                                                                                                                       tttcctttGGGTGAAGCCAGAACTCAGGGTCGGGTA‑GGGAAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATaaa                                      <  2:717150/232‑1 (MQ=255)

                                                                                                                                                                                                       tttaCTCTGGTTGAAGCCAGAATTCAGGGGGGGGTA‑GGTTAAGTATTTTACGG‑TTATC‑TGGTTGCGATGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAGTTGTCGGTAGATGATATCGCGAAACGCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATaaa                                      <  2:530154/236‑1 (MQ=255)

                                                                                                                                                                                                       ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGGAACAGAGAAAATAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTAACAACAACTGGCACAGCAAAATAATGTGTTGGGAGGTGGTATCGCGAAACTCACCGGGCAAAAACTCTGAAGCT‑‑GGGCCAAACCCGGGGAAAAATGACAAggggcgtg                                      >  2:1275587/1‑234 (MQ=255)

                                                                                                                                                                                                       ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                       <  2:1208942/144‑1 (MQ=255)

                                                                                                                                                                                                       ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAAATCGCCGGGCAAAAGCTC‑GTAGCCC‑GGGCCAAACCGGGGGAAAATGTACAAGGGATaaa                                      >  1:574236/1‑239 (MQ=255)

                                                                                                                                                                                                       ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                       >  1:1208942/1‑144 (MQ=255)

                                                                                                                                                                                                       ttGACTATGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATaaa                                      <  1:108859/239‑1 (MQ=255)

                                                                                                                                                                                                                        cAGAACTCAGGGGGGGGTA‑GGTTAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGGCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGaa                     <  1:1093486/239‑1 (MQ=255)

                                                                                                                                                                                                                                       ggTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCaa       <  1:1150833/238‑1 (MQ=255)

                                                                                                                                                                                                                                       ggTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCaa       <  1:340416/238‑1 (MQ=255)

                                                                                                                                                                                                                                            ggTTAAACCTTTTACGG‑TTTTC‑TGGTTGCGCTTAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCg  <  2:1155766/239‑1 (MQ=255)

CATTGAAGTCGCTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAA‑TGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAAACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCG  >  NC_000913/1445746‑1446217

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`