Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 788,389 A→C 5.2% L150R (CTG→CGG)  galM ← aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913788,3890AC 5.2% 533.9 / 11.3 152L150R (CTG→CGG) galMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (72/71);  new base C (4/5);  total (76/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 154 total aligned reads displayed.

CTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTC  >  NC_000913/788117‑788660
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                               
cTTTGGGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATGTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                                                                                                                                               <  4:219833/275‑1 (MQ=255)
   tGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTGGTGATTCCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGCCGTGCTTAGTCATATTCACAGGGCAAGGTTTAGCAAATGTGGCGCGATAAGTAATGGAGATCCGGTTAACGCCGGTCAGACg                                                                                                                                                                                                                                                                             >  4:85162/1‑274 (MQ=255)
               cgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGt                                                                                                                                                                                                                                                                 >  2:7237/1‑274 (MQ=255)
                             tCACTGGGAATGATGTTGGTGCTGCGGAAATCGAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTAATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAc                                                                                                                                                                                                                                                   <  2:415912/274‑1 (MQ=255)
                             cgcgtgggcatgatTTTGGTGCGTGGGAAATAAAAAGCGTTGTCGGTCGCAGATTTAGGGCCTACGAGGGGAATGCCGCCTGCACCACCCGGGAGATATTCGTCCGCCAGAATCTGCAAGTTGGAATTGCCAACGTCAGACTGCTCGCCGTCAAGATGGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGAAACGCGATAAGTAGTGGAGGTACGGTGATGGGCGGTCAGACGATATTGCACCGACGCGCCGGGATTAc                                                                                                                                                                                                                                                   <  4:7237/265‑1 (MQ=255)
                                             tGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGCTTACCCGGGAAGCCCTgctg                                                                                                                                                                                                                                   >  1:9905/1‑271 (MQ=255)
                                                        aaaTAAACAGCCGTTCCGGGGNNCACATNTCCGGCCGTCGTCCGTAATGCCGCCTTCAGCACACGTCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGGCAGACTGCTCGCCGTCAAGAGTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCt                                                                                                                                                                                                                         <  3:9905/274‑1 (MQ=255)
                                                                      tCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACa                                                                                                                                                                                                         >  2:536135/1‑275 (MQ=255)
                                                                                         ggCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGAACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTc                                                                                                                                                                                        <  2:436407/275‑1 (MQ=255)
                                                                                                        tGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACCATCTGc                                                                                                                                                                         >  1:49416/1‑275 (MQ=255)
                                                                                                                    caACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                                                                                                                                               <  1:138768/161‑1 (MQ=255)
                                                                                                                                     tCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                                                                                                                                                 <  8:268085/142‑1 (MQ=255)
                                                                                                                                     tCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                                                                                                                                                 >  7:268085/1‑142 (MQ=255)
                                                                                                                                                        tgtgATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCCAACAGCACCTGACGATCGTTCTTGTTCACAATCTTCCAGCGACGTTTGTCGAACCCTTCAGGCCCGCCGTGCAGCTGGTTAACg                                                                                                                         >  2:486701/1‑275 (MQ=255)
                                                                                                                                                             ttGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                               >  8:803923/1‑136 (MQ=255)
                                                                                                                                                             ttGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                               <  7:803923/136‑1 (MQ=255)
                                                                                                                                                                 gcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGatat                                                                                                                                                                                                                                                                         >  5:139219/1‑122 (MQ=255)
                                                                                                                                                                 gcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGatat                                                                                                                                                                                                                                                                         <  6:139219/122‑1 (MQ=255)
                                                                                                                                                                  cACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                                >  7:729073/1‑114 (MQ=255)
                                                                                                                                                                  cACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                                <  5:492230/114‑1 (MQ=255)
                                                                                                                                                                  cACGTAAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                                >  6:492230/1‑114 (MQ=255)
                                                                                                                                                                      tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGc                                                                                                                                                                                                                                          >  6:959719/1‑148 (MQ=255)
                                                                                                                                                                      tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGc                                                                                                                                                                                                                                          >  8:506607/1‑148 (MQ=255)
                                                                                                                                                                       cAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGaa                                                                                                          <  1:9076/275‑1 (MQ=255)
                                                                                                                                                                          aCTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                               <  8:513558/123‑1 (MQ=255)
                                                                                                                                                                                cGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAc                                                                                                                                                                                                                                                   <  5:1102431/129‑1 (MQ=255)
                                                                                                                                                                                 gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAGGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTCGCGAAAGCGTCCCgg                                                                                                >  2:91525/1‑275 (MQ=255)
                                                                                                                                                                                 gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGc                                                                                                                                                                         >  2:554644/1‑202 (MQ=255)
                                                                                                                                                                                 gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACgg                                                                                                >  4:412424/1‑275 (MQ=255)
                                                                                                                                                                                      cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatc                                                                                                                                                                                                                             <  8:386543/145‑1 (MQ=255)
                                                                                                                                                                                        aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataa                                                                                                                                                                                                                                                                         >  8:625661/1‑98 (MQ=255)
                                                                                                                                                                                        aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      <  5:746327/134‑1 (MQ=255)
                                                                                                                                                                                        aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      >  6:746327/1‑134 (MQ=255)
                                                                                                                                                                                        aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAAAGGAGATACGGTTATCGTCTGACCGACGataa                                                                                                                                                                                                                                                                         <  7:625661/99‑2 (MQ=255)
                                                                                                                                                                                           ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGCCg                                                                                                                                                                                                                                                                             <  7:803666/92‑1 (MQ=255)
                                                                                                                                                                                           ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          <  5:161151/111‑1 (MQ=255)
                                                                                                                                                                                           ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          <  6:787790/111‑1 (MQ=255)
                                                                                                                                                                                           ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                                                                                                                                                     >  5:1019730/1‑148 (MQ=255)
                                                                                                                                                                                           tggAAAAAGCCGTGATTAGTGAAGTTCAACGGCACAAGTTTATGAACTGGGGCGGGAGAAGTAATAGAGAGTCAGTTATCGGCGGTCAGACGACAATGCACAGTGGGGCCGGGATTAACCGGTAAGCCCTGATCACCATCATCTGAACGCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTGCACCGTc                                                                                      <  4:361930/273‑1 (MQ=255)
                                                                                                                                                                                             gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                    <  2:342898/275‑1 (MQ=255)
                                                                                                                                                                                             gAAATAGACGTGATTAGTCATATTAACCGGGCAAGGTTTAGCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGGGATTACCCGGGAAGCCCTGATCACCATCATCTGAACGCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTGCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                    <  1:115203/275‑1 (MQ=255)
                                                                                                                                                                                              aaatggcagTGATTGGTCATCTTTACCGGGCAAGGTTTTTCAACGGTGGGGCGTTAAGTAATGGGGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaaa                                                                                   <  4:536134/267‑1 (MQ=255)
                                                                                                                                                                                               aaTAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                                                                                                                                                 <  6:1019730/148‑1 (MQ=255)
                                                                                                                                                                                                 tAGACTTGATTAGTAATATTCACCGGGCAAGGTTTATCAACTGTGGGGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCAGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAgg                                                                                 <  2:149437/274‑1 (MQ=255)
                                                                                                                                                                                                   gACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAAGGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCa                                                                                                                                                                                                             >  6:916314/1‑148 (MQ=255)
                                                                                                                                                                                                         aTTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGt                                                                                                                                       >  4:329274/1‑212 (MQ=255)
                                                                                                                                                                                                             gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAg                                                                                                                                   <  2:329275/212‑1 (MQ=255)
                                                                                                                                                                                                              tCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTg                                                                                                                                                                                                  <  7:506607/148‑1 (MQ=255)
                                                                                                                                                                                                                atatTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGata                                                                                                                                                                                                                                                                          <  5:333737/74‑1 (MQ=255)
                                                                                                                                                                                                                atatTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGata                                                                                                                                                                                                                                                                          >  6:333737/1‑74 (MQ=255)
                                                                                                                                                                                                                 tatTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGAt                                                                                                                                                                                                                                    >  7:247103/1‑111 (MQ=255)
                                                                                                                                                                                                                        cGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGTTTAACGCCCTGACTTTGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCCATTGGCGATActg                                                          >  2:202325/1‑272 (MQ=255)
                                                                                                                                                                                                                          ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTg                                                                                                                                                                                      <  7:553802/148‑1 (MQ=255)
                                                                                                                                                                                                                          ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACCGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCCGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTg                                                                                                                                                                                      <  5:916314/148‑1 (MQ=255)
                                                                                                                                                                                                                            cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   >  5:339578/1‑117 (MQ=255)
                                                                                                                                                                                                                             aaGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          >  7:1032092/1‑77 (MQ=255)
                                                                                                                                                                                                                             aaGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          <  8:1032092/77‑1 (MQ=255)
                                                                                                                                                                                                                              aGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                  >  7:212379/1‑100 (MQ=255)
                                                                                                                                                                                                                               ggTTTATCAACTGTGGCGCGATAAGTAATGGGGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGGGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATa                                                  <  2:263029/275‑1 (MQ=255)
                                                                                                                                                                                                                                gTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATGTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                  >  4:538817/1‑98 (MQ=255)
                                                                                                                                                                                                                                 tttATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                                                                                                                                                        >  5:270951/1‑91 (MQ=255)
                                                                                                                                                                                                                                    aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                  >  4:560550/1‑94 (MQ=255)
                                                                                                                                                                                                                                    aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCAc                                                                                          >  4:220411/1‑230 (MQ=255)
                                                                                                                                                                                                                                     tCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCt                                                                                                                                                                           >  5:464063/1‑148 (MQ=255)
                                                                                                                                                                                                                                      cAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTg                                                                                                                                                                          >  6:411154/1‑148 (MQ=255)
                                                                                                                                                                                                                                       aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACt                                                                                                                                                                                                                    <  5:432566/105‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCATACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                                                                                                                                              <  2:560551/94‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   >  6:221957/1‑105 (MQ=255)
                                                                                                                                                                                                                                        aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                          <  1:29619/114‑1 (MQ=255)
                                                                                                                                                                                                                                         cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                                                                                                                                                        >  6:1088476/1‑99 (MQ=255)
                                                                                                                                                                                                                                         cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                                                                                                                                                       >  8:682794/1‑100 (MQ=255)
                                                                                                                                                                                                                                             ggCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGa                                                                                                                                                                   <  6:464063/148‑1 (MQ=255)
                                                                                                                                                                                                                                              gcgcGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcaca                                                                                                                                                                               <  6:232686/135‑1 (MQ=255)
                                                                                                                                                                                                                                              gcgcGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcaca                                                                                                                                                                               >  5:232686/1‑135 (MQ=255)
                                                                                                                                                                                                                                                 cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                                                                                                                                                       >  8:1120775/1‑124 (MQ=255)
                                                                                                                                                                                                                                                 cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGCCTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGacc                                 >  2:230868/1‑274 (MQ=255)
                                                                                                                                                                                                                                                 cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGtt                                                                                                                                                               <  8:215417/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                                                                                                                                                 <  8:1015982/79‑1 (MQ=255)
                                                                                                                                                                                                                                                    tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAc                                                                                                                                                                                                                                 >  7:1015982/1‑79 (MQ=255)
                                                                                                                                                                                                                                                    tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACg                                                                                                                                                                                               <  6:302662/113‑1 (MQ=255)
                                                                                                                                                                                                                                                    tAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTTTGGTACACAATCTGCCAGCGACGTTTGt                                                                                                                                                            <  8:77863/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                            >  8:307424/1‑127 (MQ=255)
                                                                                                                                                                                                                                                         aaTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAAc                                                                                                                                                       >  7:328513/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          aTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                          >  7:195510/1‑96 (MQ=255)
                                                                                                                                                                                                                                                          aTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                          <  8:195510/96‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGcc                                                                                                                                                                        <  5:796162/129‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      >  7:447382/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                aTACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGccag                 >  2:210121/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                aTACGGTTATCGGCGGTCCAACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGTTCACCATCATCTGAACTCAGGTCAAACAGCACCTGACGATCGTGCTGGTTCACAAGCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcct                  <  4:486700/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc                                                                                                                                              <  6:138326/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCCTAACGACCAATCGAGGCCCCCAGAAACTCGGcctga                >  2:138027/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                  aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGAACACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTTGAACCTTTCCGGCCCGCCGTGCCGCTGGTTAACGCCCTGACTTGGCGAAAGCGGCACGGTTTCACCGTCAAATGTCTAACGGCTATTGGCGATACGGTTGGCag                                                   >  2:243104/1‑238 (MQ=255)
                                                                                                                                                                                                                                                                      ttATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCaa                                                                                                                                                                                                            <  6:429089/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                                                                                                                                                          <  7:758060/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                                                                                                                                                          >  8:758060/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                          cgtcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTg                                                                                                                                      >  8:375346/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                            tcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                            >  5:53567/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                            tcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                            <  6:53567/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ggTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCAt    >  2:378922/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                               gTCAGACGATATTGCACCGTCGCGCTGAGATTACCCGGGAAGCCCTGTTCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTc  <  2:400755/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                               
CTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTC  >  NC_000913/788117‑788660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: