Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,902,161 T→C L89S (TTA→TCA)  ygcE → putative sugar kinase YgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,1610TC96.4% 182.7 / ‑2.1 56L89S (TTA→TCA) ygcEputative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (29/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCTGACGCATCGCTTAAC  >  NC_000913/2902031‑2902296
                                                                                                                                  |                                                                                                                                       
gCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGc                                                                                                                                >  1:444518/1‑140 (MQ=255)
 cAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGcc                                                                                                                               >  1:971142/1‑140 (MQ=255)
 cAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGcc                                                                                                                               >  1:361862/1‑140 (MQ=255)
  aGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCg                                                                                                                              <  2:652830/140‑1 (MQ=255)
   gAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGa                                                                                                                             >  1:83443/1‑140 (MQ=255)
      cATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTTACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATgg                                                                                                                          <  2:615722/140‑1 (MQ=255)
         ccTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCAc                                                                                                                       <  2:83443/140‑1 (MQ=255)
                acgaTTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCctgct                                                                                                                >  1:74069/1‑140 (MQ=255)
                            catcatTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGAtt                                                                                                    >  1:903828/1‑140 (MQ=255)
                              tcatTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAg                                                                                                  <  2:971142/140‑1 (MQ=255)
                               catTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGc                                                                                                 >  2:389579/1‑140 (MQ=255)
                               catTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGa                                                                                                      >  1:487222/1‑135 (MQ=255)
                               catTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGa                                                                                                      <  2:487222/135‑1 (MQ=255)
                                atTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCAGATGGCACGCCTGCTGCGCCGTTGATTAGCt                                                                                                <  1:362947/140‑1 (MQ=255)
                                     gTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAg                                                                                           >  2:659719/1‑140 (MQ=255)
                                      ttttGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAgg                                                                                          <  2:444518/140‑1 (MQ=255)
                                         tGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATg                                                                                       >  2:218205/1‑140 (MQ=255)
                                         tGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGAAGGATg                                                                                       >  2:268835/1‑140 (MQ=255)
                                           ccGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCa                                                                                     <  1:218205/140‑1 (MQ=255)
                                             ggTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcg                                                                                   <  1:659719/140‑1 (MQ=255)
                                             ggTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGAAGGATGCAcg                                                                                   <  1:268835/140‑1 (MQ=255)
                                                    aTTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTaca                                                                            <  2:903828/140‑1 (MQ=255)
                                                    aTTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTaca                                                                            <  2:74069/140‑1 (MQ=255)
                                                    aTTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATACACGCGTTaca                                                                            <  2:651319/140‑1 (MQ=255)
                                                        gatgaGTCAGTTTTCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGcc                                                                        <  1:538256/140‑1 (MQ=255)
                                                                aGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGaa                                                                <  1:389579/140‑1 (MQ=255)
                                                                    tGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacaca                                                            >  1:607693/1‑140 (MQ=255)
                                                                     gCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacacac                                                           <  2:798671/140‑1 (MQ=255)
                                                                        gggAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTaa                                                        >  1:937015/1‑140 (MQ=255)
                                                                            aTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACg                                                                  >  2:28429/1‑126 (MQ=255)
                                                                            aTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACg                                                                  <  1:28429/126‑1 (MQ=255)
                                                                            aTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACg                                                                  >  1:436955/1‑126 (MQ=255)
                                                                            aTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACg                                                                  <  2:436955/126‑1 (MQ=255)
                                                                                     atatTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCa                                           >  1:693747/1‑140 (MQ=255)
                                                                                       atTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAccc                                                     <  1:490150/128‑1 (MQ=255)
                                                                                       atTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAccc                                                     >  2:490150/1‑128 (MQ=255)
                                                                                       atTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACg                                                >  2:791729/1‑133 (MQ=255)
                                                                                       atTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACg                                                <  1:791729/133‑1 (MQ=255)
                                                                                         tGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATg                                       >  1:86136/1‑140 (MQ=255)
                                                                                               tATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTgg                                             <  2:974007/128‑1 (MQ=255)
                                                                                               tATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTgg                                             >  1:974007/1‑128 (MQ=255)
                                                                                               tATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCa                                           >  2:417660/1‑130 (MQ=255)
                                                                                               tATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCa                                           <  1:417660/130‑1 (MQ=255)
                                                                                                    gTCTGGGATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCtttt                            >  1:832769/1‑140 (MQ=255)
                                                                                                         ggATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCggg                       <  1:82754/140‑1 (MQ=255)
                                                                                                         ggATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCggg                       <  2:937015/140‑1 (MQ=255)
                                                                                                         ggATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCggg                       <  1:42659/140‑1 (MQ=255)
                                                                                                         ggATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCggg                       >  2:359325/1‑140 (MQ=255)
                                                                                                         ggATCCATCCGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCggg                       >  2:324946/1‑140 (MQ=255)
                                                                                                              catccGTTGCTGCCGTGCGTCATTGAAAGCCGATGGTACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATc                  >  2:464800/1‑140 (MQ=255)
                                                                                                              catccGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATc                  >  2:213212/1‑140 (MQ=255)
                                                                                                              catccGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATc                  >  2:703956/1‑140 (MQ=255)
                                                                                                               atccGTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCt                 >  1:848407/1‑140 (MQ=255)
                                                                                                                   gTTGCTGCCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCTGACg             <  2:693747/140‑1 (MQ=255)
                                                                                                                        tgcCGTGCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCTGACGCATCg        >  2:651447/1‑140 (MQ=255)
                                                                                                                              gCGTCATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCTGACGCATCGCTTAAc  >  1:251447/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                       
GCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCTCTTTTTCGGGTTATCTGACGCATCGCTTAAC  >  NC_000913/2902031‑2902296

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