Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 788,387 T→A 5.2% T151S (ACC→TCC) ‡ 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,3870TA 5.2% 547.5 / 13.3 152T151S (ACC→TCC) ‡galMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base T (73/71);  new base A (4/4);  total (77/75)
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.15e-01

Only 100 of 157 total aligned reads displayed.

TCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTC  >  NC_000913/788114‑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)
                                                  ggccTGCGGAAACCAAACGCCGTTCCGGCGACAGTTTTAAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCAAGAAGCTGCAACTGGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAcc                                                                                                                                                                                                                                <  4:518388/272‑1 (MQ=255)
                                                           aaaTAAACAGCCGTTCCGGGGNNCACATNTCCGGCCGTCGTCCGTAATGCCGCCTTCAGCACACGTCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGGCAGACTGCTCGCCGTCAAGAGTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCt                                                                                                                                                                                                                         <  3:9905/274‑1 (MQ=255)
                                                                         tCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACa                                                                                                                                                                                                         >  2:536135/1‑275 (MQ=255)
                                                                                            ggCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGCCGATCGtt                                                                                                                                                                                         >  2:261234/1‑274 (MQ=255)
                                                                                                           tGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACCATCTGc                                                                                                                                                                         >  1:49416/1‑275 (MQ=255)
                                                                                                                     atcaACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGAATACCCGGGAAGCCCTGATCAACATCATCTGAACTCAGGGCAAACAGCCACTGACGACCGTTCGTGTTCCCAATCGGcacgccacct                                                                                                                                                                >  4:400754/1‑265 (MQ=255)
                                                                                                                       caACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                                                                                                                                               <  1:138768/161‑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                                                                                                                                                                                                                                                                         <  6:139219/122‑1 (MQ=255)
                                                                                                                                                                     cACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                                <  5:492230/114‑1 (MQ=255)
                                                                                                                                                                     cACGTAAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAg                                                                                                                                                                                                                                                                                >  6:492230/1‑114 (MQ=255)
                                                                                                                                                                         tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGc                                                                                                                                                                                                                                          >  8:506607/1‑148 (MQ=255)
                                                                                                                                                                          cAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGaa                                                                                                          <  1:9076/275‑1 (MQ=255)
                                                                                                                                                                             aCTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                               <  8:513558/123‑1 (MQ=255)
                                                                                                                                                                             aCTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                                                                                                                                               >  7:513558/1‑123 (MQ=255)
                                                                                                                                                                                   cGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTAc                                                                                                                                                                                                                                                   <  5:1102431/129‑1 (MQ=255)
                                                                                                                                                                                    gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAGGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTCGCGAAAGCGTCCCgg                                                                                                >  2:91525/1‑275 (MQ=255)
                                                                                                                                                                                    gCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACgg                                                                                                >  4:412424/1‑275 (MQ=255)
                                                                                                                                                                                        naaaGATNANAATAGACGTGATTAGTCATATTCACCGGGCAAGGNTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGc                                                                                                                                                                     <  4:554643/200‑1 (MQ=255)
                                                                                                                                                                                         cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatc                                                                                                                                                                                                                             <  8:386543/145‑1 (MQ=255)
                                                                                                                                                                                         cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatc                                                                                                                                                                                                                             >  7:386543/1‑145 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataa                                                                                                                                                                                                                                                                         >  8:625661/1‑98 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      >  6:746327/1‑134 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      <  5:746327/134‑1 (MQ=255)
                                                                                                                                                                                           aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAAAGGAGATACGGTTATCGTCTGACCGACGataa                                                                                                                                                                                                                                                                         <  7:625661/99‑2 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          <  5:161151/111‑1 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                                                                                                                                                     >  5:1019730/1‑148 (MQ=255)
                                                                                                                                                                                              ttGAAATAGACGTAATTAGTCATATTCACCGGGCAAGTTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                                                                                                                                                     >  7:77863/1‑148 (MQ=255)
                                                                                                                                                                                                gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                    <  2:342898/275‑1 (MQ=255)
                                                                                                                                                                                                gAAATAGACGTGATTAGTCATATTAACCGGGCAAGGTTTAGCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGGGATTACCCGGGAAGCCCTGATCACCATCATCTGAACGCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTGCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCaa                                                                                    <  1:115203/275‑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)
                                                                                                                                                                                                                   atatTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGata                                                                                                                                                                                                                                                                          <  5:333737/74‑1 (MQ=255)
                                                                                                                                                                                                                    tatTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGAt                                                                                                                                                                                                                                    <  8:247103/111‑1 (MQ=255)
                                                                                                                                                                                                                           cGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGTTTAACGCCCTGACTTTGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCCATTGGCGATActg                                                          >  2:202325/1‑272 (MQ=255)
                                                                                                                                                                                                                             ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACCGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCCGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTg                                                                                                                                                                                      <  5:916314/148‑1 (MQ=255)
                                                                                                                                                                                                                               cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   >  5:339578/1‑117 (MQ=255)
                                                                                                                                                                                                                               cAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   <  6:339578/117‑1 (MQ=255)
                                                                                                                                                                                                                                aaGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCg                                                                                                                                                                                                                                                          <  8:1032092/77‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                  >  7:212379/1‑100 (MQ=255)
                                                                                                                                                                                                                                  ggTTTATCAACTGTGGCGCGATAAGTAATGGGGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGGGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATa                                                  <  2:263029/275‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                                                                                                                                                 <  6:709385/116‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATa                                                  <  2:412425/275‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATGTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                                                                                                                                                  >  4:538817/1‑98 (MQ=255)
                                                                                                                                                                                                                                   gTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTcac                                                                                                                                                                                >  6:963320/1‑148 (MQ=255)
                                                                                                                                                                                                                                    tttATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                                                                                                                                                        <  6:270951/91‑1 (MQ=255)
                                                                                                                                                                                                                                    tttATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                                                                                                                                                        >  5:270951/1‑91 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGt                                                                                                                                                                                                                                                                                   <  4:526267/45‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATGTTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                                                                                                                                              <  2:538818/98‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGc                                                                                                                                                                                                              >  3:29619/1‑114 (MQ=255)
                                                                                                                                                                                                                                       aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCAc                                                                                          >  4:220411/1‑230 (MQ=255)
                                                                                                                                                                                                                                        tCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCt                                                                                                                                                                           >  5:464063/1‑148 (MQ=255)
                                                                                                                                                                                                                                         cAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTg                                                                                                                                                                          >  6:411154/1‑148 (MQ=255)
                                                                                                                                                                                                                                          aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACt                                                                                                                                                                                                                    >  6:432566/1‑105 (MQ=255)
                                                                                                                                                                                                                                          aaCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACt                                                                                                                                                                                                                    <  5:432566/105‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCATACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                                                                                                                                              <  2:560551/94‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   <  5:221957/105‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTc                                                                                                                                                                                                                   >  6:221957/1‑105 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                          <  1:29619/114‑1 (MQ=255)
                                                                                                                                                                                                                                           aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGAGCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTc                                                                                      <  2:220412/230‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                                                                                                                                                        >  6:1088476/1‑99 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                                                                                                                                                       >  8:682794/1‑100 (MQ=255)
                                                                                                                                                                                                                                            cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGa                                                                                                                                                                                                                       <  7:682794/100‑1 (MQ=255)
                                                                                                                                                                                                                                               tgGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATagag                                     <  2:465890/275‑1 (MQ=255)
                                                                                                                                                                                                                                                ggCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGa                                                                                                                                                                   <  6:464063/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgcgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACg                                                                                                                                                                 >  7:215417/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                                                                                                                                                       >  8:1120775/1‑124 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                                                                                                                                                       <  7:1120775/124‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGCCTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGacc                                 >  2:230868/1‑274 (MQ=255)
                                                                                                                                                                                                                                                    cgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGtt                                                                                                                                                               <  8:215417/148‑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                                                                                                                                                                            <  7:307424/127‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                            >  8:307424/1‑127 (MQ=255)
                                                                                                                                                                                                                                                             aTGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAAc                                                                                                                                                                                                          <  8:195510/96‑1 (MQ=255)
                                                                                                                                                                                                                                                              tGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGcc                                                                                                                                                                        >  6:796162/1‑129 (MQ=255)
                                                                                                                                                                                                                                                              tGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGcc                                                                                                                                                                        <  5:796162/129‑1 (MQ=255)
                                                                                                                                                                                                                                                               ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTg                                                                                                                                                                                                                                      >  7:447382/1‑66 (MQ=255)
                                                                                                                                                                                                                                                               ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCt                                                                                                                                                    >  5:138326/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc                                                                                                                                              <  6:138326/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCCTAACGACCAATCGAGGCCCCCAGAAACTCGGcctga                >  2:138027/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          tATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                                                                                                                                                          <  7:758060/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGt                                                                                                                                                                                          >  8:758060/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                               tcgGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATc                                                                                                                                                                            >  5:53567/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                                 ggTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCAt    >  2:378922/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                   tCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTg                                                                                                                                <  7:826096/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   tCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTg                                                                                                                                <  8:328513/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                 
TCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAATGCCGCCTTCATCAACCGGCAGATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTC  >  NC_000913/788114‑788660

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