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.5% 232.7 / ‑4.7 66N3S (AAC→AGC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (31/34);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

TCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGC  >  NC_000913/3078720‑3078982
                                                                                                                                    |                                                                                                                                  
tCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              <  1:1152314/139‑1 (MQ=255)
   ttAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCTCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataa                                                                                                                           <  1:1299386/139‑1 (MQ=255)
            cAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCt                                                                                                                  <  1:372850/139‑1 (MQ=255)
                 ggTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTcc                                                                                                                               <  2:460983/121‑1 (MQ=255)
                 ggTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTcc                                                                                                                               >  1:460983/1‑121 (MQ=255)
                    gggAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                          <  2:643351/139‑1 (MQ=255)
                             taaCGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                            >  2:973674/1‑128 (MQ=255)
                             taaCGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                            <  1:973674/128‑1 (MQ=255)
                                cGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                    <  2:491660/117‑1 (MQ=255)
                                cGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                    >  1:491660/1‑117 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataata                                                                                                                         <  1:1140728/107‑1 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataata                                                                                                                         >  2:1140728/1‑107 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              <  2:1464648/102‑1 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              >  1:1464648/1‑102 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATcc                                                                                            <  1:1529006/136‑1 (MQ=255)
                                     taataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATcc                                                                                            >  2:1529006/1‑136 (MQ=255)
                                               ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgagaga                                                                                                      >  2:318855/1‑116 (MQ=255)
                                               ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgagaga                                                                                                      <  1:318855/116‑1 (MQ=255)
                                               ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATg                                                                               >  1:920594/1‑139 (MQ=255)
                                               ccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATg                                                                               >  1:665199/1‑139 (MQ=255)
                                                  ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              >  2:800679/1‑89 (MQ=255)
                                                  ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              <  1:800679/89‑1 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              <  2:413520/78‑1 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              >  1:413520/1‑78 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              <  2:645687/78‑1 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                              >  1:645687/1‑78 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                          >  1:1235030/1‑98 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                          <  2:1235030/98‑1 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           <  1:1086470/113‑1 (MQ=255)
                                                             cACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           >  2:1086470/1‑113 (MQ=255)
                                                              aCCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGtttt                                                                >  1:1822546/1‑139 (MQ=255)
                                                                        gggATGACCATTGCAGTC‑AAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGc                                                     >  2:446542/1‑139 (MQ=255)
                                                                          gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTg                                                                                                                 <  1:363018/78‑1 (MQ=255)
                                                                          gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTg                                                                                                                 >  2:363018/1‑78 (MQ=255)
                                                                             gACCATTGCAGTCAAAAAGCCCCTAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCgg                                                                                        <  1:1317572/100‑1 (MQ=255)
                                                                             gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCgg                                                                                        >  2:1317572/1‑100 (MQ=255)
                                                                             gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGtt                                                 >  2:503970/1‑139 (MQ=255)
                                                                              aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTaga                                                                                                             >  1:1380553/1‑78 (MQ=255)
                                                                              aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTaga                                                                                                             <  2:1380553/78‑1 (MQ=255)
                                                                               ccATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                        >  1:632731/1‑82 (MQ=255)
                                                                               ccATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                        <  2:632731/82‑1 (MQ=255)
                                                                               ccATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTaa                                               >  1:876358/1‑139 (MQ=255)
                                                                               ccATTGCAGTC‑AAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAag                                              <  1:446542/139‑1 (MQ=255)
                                                                                cATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAag                                              <  1:1450962/139‑1 (MQ=255)
                                                                                  ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                    <  1:1437178/67‑1 (MQ=255)
                                                                                  ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                    >  2:1437178/1‑67 (MQ=255)
                                                                                        tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTc                                                                                                                                <  1:788326/49‑1 (MQ=255)
                                                                                        tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTc                                                                                                                                >  2:788326/1‑49 (MQ=255)
                                                                                        tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           <  2:700/86‑1 (MQ=255)
                                                                                        tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           >  1:700/1‑86 (MQ=255)
                                                                                        tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTc                                      >  2:1064381/1‑139 (MQ=255)
                                                                                          aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGa                                    >  2:190204/1‑139 (MQ=255)
                                                                                                cccccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTAttt                              <  1:503970/139‑1 (MQ=255)
                                                                                                cccccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTAttt                              <  2:920594/139‑1 (MQ=255)
                                                                                                            tGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGcaca                  >  1:55277/1‑139 (MQ=255)
                                                                                                             ggACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAg                 <  2:294020/139‑1 (MQ=255)
                                                                                                              gACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGt                >  2:806878/1‑139 (MQ=255)
                                                                                                                  tttGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAAtt            <  1:190204/139‑1 (MQ=255)
                                                                                                                    tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           >  2:1298628/1‑58 (MQ=255)
                                                                                                                    tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCg                                     <  2:1661537/112‑1 (MQ=255)
                                                                                                                    tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCg                                     >  1:1661537/1‑112 (MQ=255)
                                                                                                                    tGCACGAGCAGACATGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                           <  1:1298628/58‑1 (MQ=255)
                                                                                                                       aCGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt       <  2:657261/139‑1 (MQ=255)
                                                                                                                           gCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt       <  2:1468598/135‑1 (MQ=255)
                                                                                                                           gCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt       >  1:1468598/1‑135 (MQ=255)
                                                                                                                            cAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGc  <  1:1393623/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
TCGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGC  >  NC_000913/3078720‑3078982

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