BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,445,981	T→G	24.1%	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	24.1%	189.1 / 11.8	107	G34G (GGT→GGG)	ydbL	DUF1318 family protein
Reads supporting (aligned to +/- strand): ref base T (18/58); new base G (27/3); total (46/61)
Fisher's exact test for biased strand distribution p-value = 2.78e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CATTGAAGTCGCTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCG  >  NC_000913/1445746‑1446217

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

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

        ttgccgccccTAAGGAACC‑GATCCCTAACAATATGAACGTTAAAATTGGGCATGAAATCATCATCAAGGGAGACAAAGATGTCGAAGAGCTGCTGGAAACTCGTAGCGATCTTTTCTGAGGGGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTAAGGGGCGGGTA‑GGTGAA‑ACATTTTa                                                                                                                                                                                                                                           <  2:1328360/231‑1 (MQ=255)

        tcccTGCCCCGAAGGAACCAGAT‑ACTATAAAATTGAACGTTAAAAGTGAGCTTGGGATCATCATAAAGGCAGACAAAGATGGCGAAGAGCTGCTTGAAACTCGTAGCGAGCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTa                                                                                                                                                                                                                                           <  2:1209923/236‑1 (MQ=255)

           cTGCCCTTAAGGACCC‑GATCACTTTAAATTGGAACGTTAAAATTGAGCATGAGACCAGCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACgg                                                                                                                                                                                                                                        <  2:378904/239‑1 (MQ=255)

           cTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGAA‑AAAATTTAAgg                                                                                                                                                                                                                                        >  2:670346/1‑239 (MQ=255)

           cTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGAAGAAATGTAAGGGCATTGACTCTGGTTGAAGACATAACTCAGGGGCGGGGG‑GGGGAA‑AAAATTTAAgg                                                                                                                                                                                                                                        >  2:1365667/1‑239 (MQ=255)

                    aGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTAAGG‑GTATC‑TGGt                                                                                                                                                                                                                             >  1:221483/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGCTGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGACTTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAACACATTTTAGGG‑taacttgttggccct                                                                                                                                                                                                                        >  1:68187/1‑225 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGACAACATTACTTCTTTGTGCGTTCCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGAATGAAGCCAGAACTCAGGGGCGGGGA‑CGGGGA‑ACAATTTCCTG‑TTTTT‑TGGTTTCGCTg                                                                                                                                                                                                                      >  1:933609/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑AGGGAA‑ACATTTTACGG‑GTTTC‑TGGTTgcgccg                                                                                                                                                                                                                      >  1:596447/1‑237 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGGA‑ACATTTTACGGTTTATC‑GGGTT‑CGCTg                                                                                                                                                                                                                      >  1:1263485/1‑239 (MQ=255)

                            gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCGGGTTGAAGCCTGACCTCAGGGGCGGGGA‑GGGGAA‑AAAATTTACGG‑TTTTA‑Tggggggggtg                                                                                                                                                                                                                      >  2:493863/1‑231 (MQ=255)

                             aTCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGGA‑AAACTTTTCGG‑GTATT‑TGGGTgggcggg                                                                                                                                                                                                                     >  1:1266558/1‑236 (MQ=255)

                                    tCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGT‑GGGGAA‑AAATTTTTTCG‑TTATCTTGTTTGCCTTGAAAACga                                                                                                                                                                                                               >  2:257992/1‑238 (MQ=255)

                                        tatGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGTGAA‑AAATTTTATGG‑GTATC‑TGGGTGGGCGGAAAACgggaggg                                                                                                                                                                                                          >  2:1218021/1‑234 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAA‑ACA‑TTTACGG‑TTATC‑TGGCTGCGTCTAAAACggcggattgaa                                                                                                                                                                                                      >  1:717150/1‑230 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACAATTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGaaa                                                                                                                                                                                                      >  1:1008632/1‑239 (MQ=255)

                                            aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑AAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGAGGTTGaaa                                                                                                                                                                                                      >  1:720989/1‑239 (MQ=255)

                                              cGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACa                                                                                                                                                                                                    <  1:798329/239‑1 (MQ=255)

                                              cGTTAAAATTGAGCATGAGATCAACATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACa                                                                                                                                                                                                    <  1:493863/239‑1 (MQ=255)

                                                  aaaaTTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑AAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGAGGCTGAAAAAgaga                                                                                                                                                                                                >  1:195014/1‑239 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGGTGAAACAAAGAAATtagt                                                                                                                                                                                         >  1:675787/1‑239 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑AGTGGA‑ACATTTTCCGG‑TTATT‑TGGTTGCGCTGGAAACGGCGGGTGAAAAAGGAAAATtatt                                                                                                                                                                                         >  1:468183/1‑237 (MQ=255)

                                                         aGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGACGGGGA‑GGTGAA‑ACAATTTTGGG‑TTCTA‑TGGTTTCGCTGGAAACAGGTGCTGGAACAGAAAAATtaat                                                                                                                                                                                         >  2:641575/1‑237 (MQ=255)

                                                           cATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGGA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACCGGTGCTGGAACAGAGAAATtaggag                                                                                                                                                                                       >  1:841768/1‑239 (MQ=255)

                                                           cATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGCA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTGGCCTGAAAACGGATGCTGAAACAAGGAAAttatta                                                                                                                                                                                        >  1:1069276/1‑238 (MQ=255)

                                                             ttatatcttcatcAAGGCAAGATAAGGTGTGGAAGAGCTGCTTGAAACTCGTAGCGATTCTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                     <  2:118537/235‑1 (MQ=255)

                                                             tggattactaatcaAGGCAGACAAAGATTTGGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                     <  2:4335/229‑1 (MQ=255)

                                                             tgaggccattatcaAGGTCGAAAAAAGTGTCGAAGAGCTGCTTGAAACTCGGAGCTGTCTTTTGTGAGGTGATGGTGAAGAAAACATTGCTGCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                                                                     <  2:1011458/233‑1 (MQ=255)

                                                                   catcatcaAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑AAATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑tt                                                                                                                                                                              >  1:165324/1‑239 (MQ=255)

                                                                         caAGGCAGACAAAGATTGCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGcc                                                                                                                                                                        <  1:670346/239‑1 (MQ=255)

                                                                          aaGGCAGACAAAGTTTTCGAAGGGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGGGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCg                                                                                                                                                                       <  2:720989/239‑1 (MQ=255)

                                                                            ggCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACCTTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAACAGAGAAATTATTAGCCGACA‑TTAATTCCGa                                                                                                                                                                      >  1:886816/1‑238 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGA‑AACTTTTTCGG‑GTATT‑TGGTTGCGGTGAAAACGGATGAGGGAAAAGAGAAATTGAAAGGCGACA‑ATAATGCAGaaaggaa                                                                                                                                                                >  2:652871/1‑234 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGA‑AAATTTTACGG‑GTGTC‑TGGTTGGGTCTAAAAAGGTAGCTGAAAAAAAGAAATTATGAACCCACA‑TTAAAGCCGaaagaaa                                                                                                                                                                >  2:405977/1‑234 (MQ=255)

                                                                                 aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTTCGG‑TTATT‑TGGTTGGGCTGAAAACGGATGCTGAAACAAGGAAATTTGTAGCACACA‑TTAATAGCGAACGTaa                                                                                                                                                                >  1:1025253/1‑239 (MQ=255)

                                                                                    aagttttCGAAAGGCTGCTTGAAACTCGTGGCGATTTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGc                                                                                                                                                             <  2:675787/233‑1 (MQ=255)

                                                                                        tGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAg                                                                                                                                                          <  1:1298300/238‑1 (MQ=255)

                                                                                        tGTCGAAGAGCTGCTGGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAATCATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGGCTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTTAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATTCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAg                                                                                                                                                          <  1:876687/238‑1 (MQ=255)

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

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

                                                                                                  ctgctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑AAATTTTTCGG‑TTATC‑GTGTTGGGGTGGAAACGGATGCTGAAAAAGAGAAATTAGTAGCCAGCA‑TTAATTCCGAACGGAAAAGGAG‑GTAAcaaaaaa                                                                                                                                              >  2:220162/1‑238 (MQ=255)

                                                                                                    gctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGATGCGCTGAAAAAGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAAGGAG‑TTACCAACAACTg                                                                                                                                            >  1:508135/1‑239 (MQ=255)

                                                                                                     ctTTAAACTCGTTGGGATCTTTTCTGGGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTgg                                                                                                                                           <  2:1008632/239‑1 (MQ=255)

                                                                                                       ttaaactcgtgggttttttttttGGGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                         <  2:221483/222‑1 (MQ=255)

                                                                                                       tgaacttcgtgtggATCTTTTTTGTGTGGATGATGAAGAAAACATTACTTCTTTGTGCTTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                         <  2:1266558/226‑1 (MQ=255)

                                                                                                       tGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                         <  1:1271321/239‑1 (MQ=255)

                                                                                                         aaaCCCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgcaaa                                                                                                                                       <  1:120768/239‑1 (MQ=255)

                                                                                                            cTCGTAGTGATTTTTTCTGAGGTGATGATGAAGAAAACATTATTTCTGTGTGGGTTTCTTGTTGGTCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgcaaagca                                                                                                                                    <  2:645542/239‑1 (MQ=255)

                                                                                                               gTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAACAAAGAAATTAGTAGCCGGCA‑TTAATGCCGAACGAAAAGCGAG‑TTTCCAACAAATGGCAAAGCaaaa                                                                                                                                 >  1:1155766/1‑239 (MQ=255)

                                                                                                                   ggATCTTTTCTTGGGTGTTGGTGAAGAAAACATTCTTTTTTTGTGGGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAaataat                                                                                                                             <  2:877050/238‑1 (MQ=255)

                                                                                                                       tttttttGGGGTGATGATGAAGAAAACATTTCTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                         <  2:1069276/238‑1 (MQ=255)

                                                                                                                       cTTTTCTGAGGTTATGATTAAGAAAACATTACTTCCTTGGGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTTTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                         <  1:1255731/239‑1 (MQ=255)

                                                                                                                       cTTTTCTGAGGTTATGAGTAAAAAAACATTTCTTCTTTGTGCTTTTCTTGTTGGGGTGGTAAGCAGCAATGTAATGGCATTGACTCTGAATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATgtgt                                                                                                                         <  2:1132305/239‑1 (MQ=255)

                                                                                                                         tttCTTAGGTGATTATGAAAAAAACTTTGCTGCTTTGGGTGTTTTTTGGTGGGCTGGTAAGCAGGAATGTAATGGCATTGACTCTGGATGAAGCCAGAACGCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAAGTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCg                                                                                                                       <  2:933609/239‑1 (MQ=255)

                                                                                                                              ggggTGAGGATGAAGAAAAACTTACTGCGTTGGGAGTTTCTTGGTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAGGTGTCGGTAga                                                                                                                  <  1:962407/237‑1 (MQ=255)

                                                                                                                               aGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGC‑GGAGAA‑ACATTTTCTGG‑TGATC‑TTGGTGCACTGAACACGGAACCTGAAACAGAGAAATTAGTAACGGCCC‑TTAATGCCGAAAGTAAATCGAG‑GTACCAACAAATGGCAAAACAaaaaaaaaagtggggagaa                                                                                                                 >  2:668884/1‑226 (MQ=255)

                                                                                                                                   gaggatgaAGAAAACATTACTTTTTTGTGCGTTTCTTGGTGGGCTGGTAAGCAGCAATGTAATGGCATGGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAgat‑gat                                                                                                             <  1:405977/235‑1 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGGAAACGGATGCTGAAAAAAAGAAATTAATAGCCGACA‑TTAAATCCCAAACTAAAGCGCG‑TTTCCCACAACTGGCAAAGCCAAATCATATGTCGGGAGAT‑Gata                                                                                                            >  1:895763/1‑238 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGGA‑AAATTTTTCGG‑GTATA‑TGGTTGCGCGGGAAACGGGTGCGGAAAAAGGGAAATTAGAGACCTGAA‑TTTATGCCCAAAGTAAA‑AGAGATTACAAACAACCGGGAAAAAAAAAAAATGAGTCTGAagaa‑gaaaa                                                                                                           >  2:363972/1‑236 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTTGGCTGGAAACGGATGCTGAAACAGAGAAATTAGTATCCGCCA‑TTAATGCCGAACGTAAAGCGAG‑TTACCCAAAAATGGCCAAACAAAAAAAAgtgtgtgaagaagaaaa                                                                                                            >  2:1150833/1‑226 (MQ=255)

                                                                                                                                    atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAAAAGAGAGATAATTAACCCGCATTTAATG‑CGAAGCGAAAGCCGG‑TTAACAAAAAATGGCACAGaaaaaaaaaagagaggggaagaaaa                                                                                                             >  2:340416/1‑221 (MQ=255)

                                                                                                                                       atgaAGAAAACATTTCTTCGTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAAGGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCg                                                                                                                       <  1:767711/225‑1 (MQ=255)

                                                                                                                                                 cATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑AACTTTTACGG‑TTATC‑TGGTTTCGCTGGAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAACGGG‑TTACCAACAAACGGCCAAGCAAAAAAAATGGTCGGAAgag‑gaaaacagtccctttaaa                                                                                              >  2:767711/1‑223 (MQ=255)

                                                                                                                                                       ttttttGTGTTTTTTTTGTTGGGGTGTTAAGGAGCAATTTAATGGAATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCaa                                                                                        <  2:165324/236‑1 (MQ=255)

                                                                                                                                                        tgtttGTGGTTTTCTTGTTGGGCTGGTAAGCGGAAATTGAATGGAATTGACTATGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCaaa                                                                                       <  2:510758/237‑1 (MQ=255)

                                                                                                                                                           ttGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCCCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGGAAAGCGAG‑TTTCCAACAACTGGCAAAGCAAAAAAATGTGTCGGGAAGT‑GAAATCGCGAAAATCCCCGGGCcaaaa                                                                                     >  1:657501/1‑237 (MQ=255)

                                                                                                                                                            tgtgCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑AAATTTTACGG‑GTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAAAGTAAAACGAG‑TTACCAACAAATGGCCAAGCAAAATAATGTGTCGGGAGAT‑GATATCGCGAAAATCGCCCGGCAAAAGc                                                                                    >  1:1318953/1‑238 (MQ=255)

                                                                                                                                                                   ttCATGGTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCt                                                                           <  1:363972/239‑1 (MQ=255)

                                                                                                                                                                        tgttgggtgttagggaggattGTATTGGGATTGACCCGGGATGAAGCCAGAACTTAGGGGGGGGTA‑GGTGAA‑ACTTGTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGc                                                                     <  2:841768/219‑1 (MQ=255)

                                                                                                                                                                        gttgGGCGGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAACCAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGc                                                                     <  1:1149033/239‑1 (MQ=255)

                                                                                                                                                                            ggCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAg                                                                                                                                                          <  1:934260/154‑1 (MQ=255)

                                                                                                                                                                                 gTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑ggg                                                           <  1:220162/239‑1 (MQ=255)

                                                                                                                                                                                  tAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑ggg                                                           <  1:257992/238‑1 (MQ=255)

                                                                                                                                                                                  tAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑ggg                                                           <  1:652871/238‑1 (MQ=255)

                                                                                                                                                                                  gaaGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑ggg                                                           <  1:42869/237‑1 (MQ=255)

                                                                                                                                                                                   aaggcgaaagttAATTGAATTGCTTCTGGGTGAATCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTAAC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑ggggg                                                         <  2:508135/228‑1 (MQ=255)

                                                                                                                                                                                      cagcaATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGA‑GGGGAA‑ACATTTTACGG‑GTATT‑TGGGTGCGCTGAAAAAGGAAGCTGAAACCGAAGAATTTGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAAAAAATGGCAAAGCAAAATAATGTGACGGGATAG‑GAGAAGGCGAAACACGCCAGGCAAAAACAC‑GCATGTC‑GGGCCAAAAC‑ggggggaa                                                      >  2:716215/1‑236 (MQ=255)

                                                                                                                                                                                       agcaATGTAATGGCATTGACTCTGGACGAAGCCAGAACCCAGGGGCGGGTA‑GGTGAA‑ACCTTTTACGG‑TTATC‑TGGGTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAat                                                      <  1:1365667/238‑1 (MQ=255)

                                                                                                                                                                                              tAAGGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACaa                                              <  1:1218021/239‑1 (MQ=255)

                                                                                                                                                                                                aTGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTTCGG‑TTATC‑TGGTTGCGCTGAAAACGGGTGCTGGAACAAAGAAATTATTATCCGACA‑TTTATGCCGAACCGAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAATTTTCGGTAGAT‑GATATCGCGAAACTCGCCGGGGAAAAAACC‑GTAACGCGGGGCCAACCC‑GGGGGaaatagtaaaag                                             >  2:309841/1‑229 (MQ=255)

                                                                                                                                                                                                   gCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTTCGCTGGAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATTCCGAACGTAAAGCGGG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAAAT‑GATAACGCGAAAACTGCCGGGCCAAAGCCC‑GTAGCTC‑GGGCCCAACC‑CGGGGGAAACGTACAAGAGAt                                         >  1:753167/1‑239 (MQ=255)

                                                                                                                                                                                                     gttGACTTTGGATTAAGCCCAAACTTCGGGGGGGGTG‑GGTGAA‑ACCTTTTACGG‑TTAAC‑TGGTTGCGCTGAAAACGGATGCTGAAACAAAGAAATTAGTTGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAGCTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATa                                        <  2:596447/237‑1 (MQ=255)

                                                                                                                                                                                                     aTTGTCTCTGGATGAAGCCGGAACTCGGGGGGGGGTA‑GGTGAA‑ACATTTTAAGG‑TTATC‑TGGTTGCGCTGAAAACGGATGTTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATa                                        <  2:195014/238‑1 (MQ=255)

                                                                                                                                                                                                     aTTGGCTGTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTGTACGG‑TTATC‑GGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATa                                        <  1:641575/238‑1 (MQ=255)

                                                                                                                                                                                                     aTTGGCTCTGGATGAAGCCAGAACCAAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTTTA‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATa                                        <  2:895763/238‑1 (MQ=255)

                                                                                                                                                                                                     aTTGACTTTGGATGAGGCAAGAACTAAGGGGCGGGTA‑GGTGAA‑ACTTTTTACGT‑TTATC‑TGTTGGCGCTGAAAACGGATGATTAAACAGATAAATTAGTAGCCGACA‑TTAATGCCGAACGGAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGATC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATa                                        <  2:986427/238‑1 (MQ=255)

                                                                                                                                                                                                      tttcctttGGGTGAAGCCAGAACTCAGGGTCGGGTA‑GGGAAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATaaa                                      <  2:717150/232‑1 (MQ=255)

                                                                                                                                                                                                      tttaCTCTGGTTGAAGCCAGAATTCAGGGGGGGGTA‑GGTTAA‑GTATTTTACGG‑TTATC‑TGGTTGCGATGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAAGTTGTCGGTAGAT‑GATATCGCGAAACGCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATaaa                                      <  2:530154/236‑1 (MQ=255)

                                                                                                                                                                                                      ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGGAACAGAGAAAATAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTAACAACAACTGGCACAGCAAAATAATGTGTTGGGAGGT‑GGTATCGCGAAACTCACCGGGCAAAAACTCTGAAGCT‑‑GGGCCAAACC‑CGGGGAAAAATGACAAggggcgtg                                      >  2:1275587/1‑234 (MQ=255)

                                                                                                                                                                                                      ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                         <  2:1208942/144‑1 (MQ=255)

                                                                                                                                                                                                      ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAAATCGCCGGGCAAAAGCTC‑GTAGCCC‑GGGCCAAACC‑GGGGGAAAATGTACAAGGGATaaa                                      >  1:574236/1‑239 (MQ=255)

                                                                                                                                                                                                      ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGgca                                                                                                                                         >  1:1208942/1‑144 (MQ=255)

                                                                                                                                                                                                      ttGACTATGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATaaa                                      <  1:108859/239‑1 (MQ=255)

                                                                                                                                                                                                             tGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGGGAA‑ACATTTTAAGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTAACAACAAACGGGAAAGCAAAATAATTTGTCGGGATATGGATA‑CGCCAAACTCCCCGGGGAAAACATA‑CGAGCCC‑GGGCCAAACCAGGGGGaaatagtaaaggggaaaaagggaaa                                >  2:1076010/1‑216 (MQ=255)

                                                                                                                                                                                                                       cAGAACTCAGGGGGGGGTA‑GGTTAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGGCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGaa                     <  1:1093486/239‑1 (MQ=255)

                                                                                                                                                                                                                                      ggTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCaa       <  1:340416/238‑1 (MQ=255)

                                                                                                                                                                                                                                      ggTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCaa       <  1:1150833/238‑1 (MQ=255)

                                                                                                                                                                                                                                           ggTTAA‑ACCTTTTACGG‑TTTTC‑TGGTTGCGCTTAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCg  <  2:1155766/239‑1 (MQ=255)

CATTGAAGTCGCTGCACCTAAGGAACC‑GATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTA‑GGTGAA‑ACATTTTACGG‑TTATC‑TGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACA‑TTAATGCCGAACGTAAAGCGAG‑TTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGAT‑GATATCGCGAAACTCGCCGGGCAAAAGCTC‑GTAGCTC‑GGGCCAAACC‑GGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATAGCAAAACCG  >  NC_000913/1445746‑1446217

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