Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,919 T→C 7.6% Y346H (TAC→CAC)  ttuD → hydroxypyruvate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,9190TC 7.6% 211.5 / 6.6 82Y346H (TAC→CAC) ttuDhydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base C (3/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCC  >  NC_002947/4891773‑4892068
                                                                                                                                                   |                                                                                                                                                     
ttCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACg                                                                                                                                                    >  1:1893619/1‑150 (MQ=255)
   aaCACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTTCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGccc                                                                                                                                                 >  2:6842/1‑150 (MQ=255)
   aaCACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGccc                                                                                                                                                 >  1:1530672/1‑150 (MQ=255)
   aaCACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGGCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTCACGCCGACTGCCTCCGCAGCCTCACCGAAAACCTCAAAGCCCTGCCGGGCGCGGCCGccc                                                                                                                                                 >  2:337450/1‑150 (MQ=255)
    acacGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCAGGGCGTACACGccc                                                                                                                                                 >  1:655612/1‑150 (MQ=255)
      acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggc                                                                                                                                              >  2:1333405/1‑147 (MQ=255)
      acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggc                                                                                                                                              >  2:1333523/1‑147 (MQ=255)
        ggCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCNGCCTCANCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggnag                                                                                                                                            <  1:1333405/150‑6 (MQ=255)
        ggCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggcag                                                                                                                                            <  1:1333523/150‑6 (MQ=255)
         gCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCAGGGCGTACACGcccggcag                                                                                                                                            <  2:655612/150‑6 (MQ=255)
         gCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCg                                                                                                                                           <  2:60908/150‑1 (MQ=255)
         gCCAGCCGCGCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCg                                                                                                                                           <  1:443817/150‑1 (MQ=255)
           cAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGt                                                                                                                                         <  1:840939/150‑1 (MQ=255)
           cAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCGGCTCAGCCTCACCGACAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCGGGt                                                                                                                                         <  2:274579/150‑1 (MQ=255)
               cGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGaca                                                                                                                                     <  1:1653197/150‑1 (MQ=255)
                gCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGacac                                                                                                                                    <  2:34722/150‑1 (MQ=255)
                 cTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACAcc                                                                                                                                   >  2:1244387/1‑150 (MQ=255)
                   cAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGa                                                                                                                                 <  1:1244387/150‑1 (MQ=255)
                    aaaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGAc                                                                                                                                <  2:970594/150‑1 (MQ=255)
                     aaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACg                                                                                                                               <  2:58073/150‑1 (MQ=255)
                     aaGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTACCGGGCGTGTACGCCCTGGCCGGTGACACCGACg                                                                                                                               <  1:340959/150‑1 (MQ=255)
                           ccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTAc                                                                                                                                                     >  2:476257/1‑122 (MQ=255)
                           ccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTAc                                                                                                                                                     <  1:476257/122‑1 (MQ=255)
                           ccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCg                                                                                                                         >  2:1132813/1‑150 (MQ=255)
                           ccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGCGTACGCCCTGGCCGGTGACACCGACGGCATCg                                                                                                                         >  2:1318313/1‑150 (MQ=255)
                              gCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATg                                                                                                                      <  2:1854430/150‑1 (MQ=255)
                                   aTCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGacg                                                                                                                      >  1:1137507/1‑143 (MQ=255)
                                   aTCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGacg                                                                                                                      <  2:1137507/145‑3 (MQ=255)
                                   aTCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCg                                                                                                                 >  1:1247080/1‑150 (MQ=255)
                                    tCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCgg                                                                                                                <  2:403485/150‑1 (MQ=255)
                                     ccTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGga                                                                                                               <  1:1919443/150‑1 (MQ=255)
                                     ccTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGga                                                                                                               <  1:579310/150‑1 (MQ=255)
                                      cTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGTAAGTCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGgaa                                                                                                              >  1:1755579/1‑150 (MQ=255)
                                      cTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGgaa                                                                                                              >  1:1740268/1‑150 (MQ=255)
                                      cTGTCGGGTGGCGAAACCACCGTGACCGTGCGAGGGAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAGATCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGgaa                                                                                                              >  1:151907/1‑150 (MQ=255)
                                                gCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCg                                                                                                    >  2:399962/1‑150 (MQ=255)
                                                gCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCg                                                                                                    >  1:474792/1‑150 (MQ=255)
                                                gCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCg                                                                                                    >  1:1764880/1‑150 (MQ=255)
                                                  gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCAGAACGTGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                                  <  1:1060085/150‑1 (MQ=255)
                                                  gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGAc                                                                                                                                <  2:1451314/120‑1 (MQ=255)
                                                  gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGAc                                                                                                                                >  1:1451314/1‑120 (MQ=255)
                                                  gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                                  >  1:1085527/1‑150 (MQ=255)
                                                       caccGTGACCGTGCGCGGCTATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGCAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCtt                                                                                             >  2:337649/1‑150 (MQ=255)
                                                          cgtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCAt                                                                                          >  2:1774592/1‑150 (MQ=255)
                                                               ccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGAccc                                                                                     <  2:1530672/150‑1 (MQ=255)
                                                                       ggCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGc                                                                             >  2:727875/1‑150 (MQ=255)
                                                                          aaTGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTAc                                                                          >  2:281885/1‑150 (MQ=255)
                                                                            tGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTg                                                                                                                                               <  1:1655305/79‑1 (MQ=255)
                                                                            tGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTg                                                                                                                                               >  2:1655305/1‑79 (MQ=255)
                                                                              gCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCa                                                                      <  1:399962/150‑1 (MQ=255)
                                                                               ccGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAg                                                                     >  1:1167453/1‑150 (MQ=255)
                                                                                      gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGaa                                                              >  1:791422/1‑150 (MQ=255)
                                                                                            aaCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTg                                                        >  2:1517773/1‑150 (MQ=255)
                                                                                             aCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTgg                                                       <  2:1105069/150‑1 (MQ=255)
                                                                                             aCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTgg                                                       <  2:474792/150‑1 (MQ=255)
                                                                                               gcagcgTTCCTGCTAAGCCTGACCGATAGCCTGAAATTCCGTCCGGGGGGGTACGCCCTGGCCGGAGACACCGACGGAATCGATCGCCCGGACGAACGCGCCGGTGGCTTTCTGACCCCCGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                     >  2:1859086/6‑150 (MQ=255)
                                                                                               gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                                  >  1:783198/1‑105 (MQ=255)
                                                                                               gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGt                                                                                                  <  2:783198/105‑1 (MQ=255)
                                                                                               aacGAGTTCCTCCTCAGCTTCACCGAAATCCGGAAGCGCCGGCCGGGCCTGTACGTCCTGGCCGGTGACACCGACAACATCGATGGCGCGGAAGAAAACGNCGGTNCCTTCATGACCCCGNCAAGCTACGCCAGCGCCGAAGCCCTGGGc                                                     <  1:57874/148‑1 (MQ=255)
                                                                                                  gAGTTCCTGCTCAGGCCAACCGAAAGACTGAAAGGCCTGCCGGGCCTGACCGTCCTGAACGGTGACACCGATGGCATCGATGGCCCGGGAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTg                                                  >  2:828450/1‑150 (MQ=255)
                                                                                                    gTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTTGACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                                <  1:108039/150‑1 (MQ=255)
                                                                                                     ttCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCg                                               <  1:1114794/150‑1 (MQ=255)
                                                                                                     ttCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCg                                               >  2:520355/1‑150 (MQ=255)
                                                                                                           cTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTCCGCCCTGGCAGGCAAAACCGACGGCATCGATGTGTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGc                                         >  1:1673420/1‑150 (MQ=255)
                                                                                                              aGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACATCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgac                                      <  2:56449/150‑1 (MQ=255)
                                                                                                                gttCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGAGGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgacga                                    >  2:675534/3‑150 (MQ=255)
                                                                                                                 cTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                                   >  1:524131/1‑150 (MQ=255)
                                                                                                                  tgaCCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGCAGAAAACGCCGGTGCCTTCATGACCCCTGCCAGCTACTCCTGCGCCTAAGCCTTGGGCCTGTCGGCCAGcgacgacc                                  >  2:1534056/3‑148 (MQ=255)
                                                                                                                  tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTg                                                        <  2:72604/128‑1 (MQ=255)
                                                                                                                  tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTg                                                        >  1:72604/1‑128 (MQ=255)
                                                                                                                  tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGc                                  >  2:1445679/1‑150 (MQ=255)
                                                                                                                  tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGc                                  <  1:727875/150‑1 (MQ=255)
                                                                                                                  acACCGAAAGACTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAACCGCCGGTGCCTTCATGCCCATGGCCAGCTACGTCAGCGCCGAACCCTTGTGCTTGTCGGCCATCGACGAGc                                  <  2:1755579/149‑1 (MQ=255)
                                                                                                                    cccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                                <  1:1528242/149‑1 (MQ=255)
                                                                                                                    aCCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGTGAGCTACGGGAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                                <  1:1774592/150‑1 (MQ=255)
                                                                                                                              tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTc                                                                                            >  2:1429846/1‑80 (MQ=255)
                                                                                                                              tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTc                                                                                            <  1:1429846/80‑1 (MQ=255)
                                                                                                                                   ggCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGcc                                                                 >  2:1771547/1‑102 (MQ=255)
                                                                                                                                   ggCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGcc                                                                 <  1:1771547/102‑1 (MQ=255)
                                                                                                                                    gCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggct                >  2:1445074/1‑150 (MQ=255)
                                                                                                                                         ccGGGCGTGTACGCCCTGGCCGGTGCCACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCCATCGATGTAACCTTGCGGCCAGCGACGAGCTGGACAACAACAacggctacggc           <  1:799672/150‑1 (MQ=255)
                                                                                                                                            gggCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAAcggctacggcta         >  2:1557144/1‑150 (MQ=255)
                                                                                                                                                 tgtACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTcg    <  1:339162/150‑1 (MQ=255)
                                                                                                                                                   tACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTcgcc  >  2:139578/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGCC  >  NC_002947/4891773‑4892068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: