Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,966,716 T→C E573E (GAA→GAG ptsP ← phosphoenolpyruvate‑protein phosphotransferase PtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,7160TC95.7% 158.3 / NA 47E573E (GAA→GAGptsPphosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (21/24);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGA  >  NC_000913/2966581‑2966845
                                                                                                                                       |                                                                                                                                 
gCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCaa                                                                                                                               >  1:735626/1‑140 (MQ=255)
   aaCAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCa                                                                                                                            >  2:36650/1‑140 (MQ=255)
        agaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAggcgg                                                                                                                       >  2:559451/1‑140 (MQ=255)
        agaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAggcgg                                                                                                                       >  1:513637/1‑140 (MQ=255)
         gaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATGTCCTCGAGTTCACGACCGGCATGCTCAATCAggcggc                                                                                                                      >  1:303277/1‑140 (MQ=255)
         gaTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAggcggc                                                                                                                      >  2:16769/1‑140 (MQ=255)
          aTGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGgcggcg                                                                                                                     <  2:371324/140‑1 (MQ=255)
           tGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGt                                                                                                                    <  2:727809/140‑1 (MQ=255)
              aaTCGACCCGCTTTGACAGATGAGGTAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCt                                                                                                                 <  1:338981/140‑1 (MQ=255)
               aTCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCtt                                                                                                                >  1:699378/1‑140 (MQ=255)
                 cGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTtca                                                                                                              >  1:550746/1‑140 (MQ=255)
                   aCCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTtcatc                                                                                                            >  1:599920/1‑140 (MQ=255)
                       gCTTTGCCAGATGCGGCCGCATAAATACCATTGACCGCAATTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACt                                                                                                        <  2:303277/140‑1 (MQ=255)
                         tttGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTc                                                                                                      >  2:366933/1‑140 (MQ=255)
                             ccAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCg                                                                                                  <  1:53504/140‑1 (MQ=255)
                             ccAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCg                                                                                                  <  2:436534/140‑1 (MQ=255)
                              cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGa                                                                                                 <  1:366933/140‑1 (MQ=255)
                              cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGa                                                                                                 <  1:309526/140‑1 (MQ=255)
                                           aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCAt                                                                                    <  2:735626/140‑1 (MQ=255)
                                            tAAATACCATTGACGGCACTTCGAGCATGATCCCGCTACCTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATc                                                                                   <  2:226997/140‑1 (MQ=255)
                                            tAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATc                                                                                   <  2:513637/140‑1 (MQ=255)
                                              aaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCgg                                                                                 <  2:599920/140‑1 (MQ=255)
                                               aTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGc                                                                                >  2:406353/1‑140 (MQ=255)
                                                tACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCa                                                                               <  2:670896/140‑1 (MQ=255)
                                                tACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCa                                                                               <  2:699378/140‑1 (MQ=255)
                                                 aCCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCaa                                                                              >  1:377655/1‑140 (MQ=255)
                                                   cATTGACGGCACTTCCAGCATGTTGCCGATCGGAGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACa                                                                            <  2:339854/140‑1 (MQ=255)
                                                   cATTGACGGCACTTCCAGCATGCTGCCGATACGTGTTTTGGGAATTTCCTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACa                                                                            <  2:314055/140‑1 (MQ=255)
                                                    aTTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTGCATCAACGTCATCGGGGCTTGAGACCATCGGCAACAg                                                                           >  1:210363/1‑140 (MQ=255)
                                                      tGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTAATCCAGGCTTGTGTCCATCGGCAACAGaa                                                                         >  1:262651/1‑140 (MQ=255)
                                                        aCGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAata                                                                       >  1:92590/1‑140 (MQ=255)
                                                          ggCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCa                                                                               >  2:92447/1‑130 (MQ=255)
                                                          ggCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCa                                                                               <  1:92447/130‑1 (MQ=255)
                                                               ttCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGt                                                                >  1:586593/1‑140 (MQ=255)
                                                                         atgCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGc                                                      >  1:696287/1‑140 (MQ=255)
                                                                         atgCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGc                                                      >  2:505294/1‑140 (MQ=255)
                                                                           gCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCgg                                                    <  2:35460/140‑1 (MQ=255)
                                                                              gATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCAt                                                 <  2:550746/140‑1 (MQ=255)
                                                                                               tttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGcc                                <  1:505294/140‑1 (MQ=255)
                                                                                               tttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGcc                                <  1:436610/140‑1 (MQ=255)
                                                                                               tttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGcc                                <  2:74623/140‑1 (MQ=255)
                                                                                                    tAACCGATCATCTCCTCGACTTTACGTCCGGTACGTTCAATCAGGCGGTGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCAc                           <  2:267949/140‑1 (MQ=255)
                                                                                                        cGATCATCTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTgg                       <  1:730176/140‑1 (MQ=255)
                                                                                                          atcatcTCCTCGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGAt                     >  2:615956/1‑140 (MQ=255)
                                                                                                              tctcCTCGACTTCACCTCCGGCGCGTTCATTCAGGCGACGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaag                 <  1:283165/140‑1 (MQ=255)
                                                                                                               ctcctcGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaaga                <  1:151564/140‑1 (MQ=255)
                                                                                                                    cGACTTCACGTCCGGCACGCTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtc           >  2:499669/1‑140 (MQ=255)
                                                                                                                          cACGTCCGGCACGCTCAATCA‑‑‑GGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGa  >  1:29069/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGA  >  NC_000913/2966581‑2966845

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