Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,078,852 T→C N3S (AAC→AGC)  cmtA ← putative mannitol‑specific PTS IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,8520TC98.4% 220.9 / NA 64N3S (AAC→AGC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/35);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTT  >  NC_000913/3078718‑3078977
                                                                                                                                      |                                                                                                                             
gTTCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTc                                                                                                                           <  2:999127/139‑1 (MQ=255)
  tCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                         <  2:582626/139‑1 (MQ=255)
  tCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                         <  1:250204/139‑1 (MQ=255)
          cAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAAtatta                                                                                                                 <  2:700449/139‑1 (MQ=255)
           aGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtattat                                                                                                                >  1:432012/1‑139 (MQ=255)
            accAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               <  1:161343/138‑1 (MQ=255)
             ccAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATcc                                                                                                              >  1:309736/1‑139 (MQ=255)
                    gTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                       <  2:283809/139‑1 (MQ=255)
                      gggAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                     <  1:1059222/139‑1 (MQ=255)
                              ataaCGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTaa                                                                                             <  2:549573/139‑1 (MQ=255)
                               taaCGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAAt                                                                                            >  2:285802/1‑139 (MQ=255)
                                  cGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataata                                                                                                                    >  1:270638/1‑112 (MQ=255)
                                  cGCAGTAATAAAACCCCAGGCAATAAAAGCACCAAGATTGGGGATGACCTTTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataata                                                                                                                    <  2:270638/112‑1 (MQ=255)
                                   gCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACGATTGCAGTGAAAAAGCCCACAAAAGCGTGGAACTTTGAACGAGCAGAATTGCTTTACATAATATTATAGTGTAGAGGAGAGAGTAATTATc                                                                                        >  2:303489/1‑139 (MQ=255)
                                   gCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATc                                                                                        >  1:34583/1‑139 (MQ=255)
                                     aGTAATAAAACCCCTGGCAATAAAAGCACCAATATTGGGGAAGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTCTCCATAATATTATc                                                                                                               <  1:654267/114‑1 (MQ=255)
                                     aGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               >  2:654267/1‑114 (MQ=255)
                                     aGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                      <  1:231011/139‑1 (MQ=255)
                                       taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               <  1:346796/112‑1 (MQ=255)
                                       taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               >  2:346796/1‑112 (MQ=255)
                                       taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                     <  1:399162/122‑1 (MQ=255)
                                       taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                     >  2:399162/1‑122 (MQ=255)
                                       taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcgcg                                                                                    >  1:14492/1‑139 (MQ=255)
                                               ccccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATaa                                                                            <  1:177986/139‑1 (MQ=255)
                                                 ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               <  2:527080/102‑1 (MQ=255)
                                                 ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               >  1:527080/1‑102 (MQ=255)
                                                   aGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               <  1:873253/100‑1 (MQ=255)
                                                   aGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                               >  2:873253/1‑100 (MQ=255)
                                                   aGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTc                                                                        >  2:523676/1‑139 (MQ=255)
                                                     gCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcgc                                                                                     <  1:576389/124‑1 (MQ=255)
                                                     gCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcgc                                                                                     >  2:576389/1‑124 (MQ=255)
                                                     gCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCg                                                                      <  2:432012/139‑1 (MQ=255)
                                                          aaaaGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAAt                                                                 >  2:420912/1‑139 (MQ=255)
                                                                          gggATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTc                                                                                                                           >  1:996692/1‑65 (MQ=255)
                                                                          gggATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTc                                                                                                                           <  2:996692/65‑1 (MQ=255)
                                                                          gggATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                 <  1:523676/139‑1 (MQ=255)
                                                                            gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                       >  2:167228/1‑83 (MQ=255)
                                                                            gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                       <  1:167228/83‑1 (MQ=255)
                                                                            gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCt                                               >  2:743899/1‑139 (MQ=255)
                                                                              tGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGt                                             <  1:743899/139‑1 (MQ=255)
                                                                               gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                 >  1:527870/1‑70 (MQ=255)
                                                                               gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                 <  2:527870/70‑1 (MQ=255)
                                                                                aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                     <  1:469717/81‑1 (MQ=255)
                                                                                aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                     >  2:469717/1‑81 (MQ=255)
                                                                                    ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATc                                                                                        <  1:887224/90‑1 (MQ=255)
                                                                                    ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATc                                                                                        >  2:887224/1‑90 (MQ=255)
                                                                                        aGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAgttgt                                   <  1:420912/139‑1 (MQ=255)
                                                                                             aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAg                              <  1:930680/139‑1 (MQ=255)
                                                                                                  cccccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                   >  1:228030/1‑65 (MQ=255)
                                                                                                  cccccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                   <  2:228030/65‑1 (MQ=255)
                                                                                                    cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                 <  1:1013433/49‑1 (MQ=255)
                                                                                                    cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                 >  2:1013433/1‑49 (MQ=255)
                                                                                                                   cTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcgcg                                                                                    >  1:277035/1‑63 (MQ=255)
                                                                                                                   cTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcgcg                                                                                    <  2:277035/63‑1 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTa                                                                                                          >  1:172485/1‑38 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTa                                                                                                          <  2:172485/38‑1 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCgg                                                                                   >  2:551667/1‑61 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCgg                                                                                   <  1:551667/61‑1 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                 <  2:554742/95‑1 (MQ=255)
                                                                                                                      tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                 >  1:554742/1‑95 (MQ=255)
                                                                                                                         gagAGCAGACTTGCTTTCGATTGTTTTTTGCTGTAGTGGAGAGAGTAATTATCACCGGCTGATAATGTGCGATAGTTGTTCTTCTGATGATGCTGTAAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTGACtt  <  2:213277/137‑1 (MQ=255)
                                                                                                                         aCGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt  <  2:34583/139‑1 (MQ=255)
                                                                                                                         aCGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt  <  2:309736/139‑1 (MQ=255)
                                                                                                                                gACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTgatgat                                                  <  2:868951/84‑1 (MQ=255)
                                                                                                                                gACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTgatgat                                                  >  1:868951/1‑84 (MQ=255)
                                                                                                                                      |                                                                                                                             
GTTCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTT  >  NC_000913/3078718‑3078977

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