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,8520TC100.0% 256.3 / NA 71N3S (AAC→AGC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/37);  total (34/37)

TTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAA  >  NC_000913/3078723‑3078979
                                                                                                                                 |                                                                                                                               
ttAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataa                                                                                                                        <  1:268690/139‑1 (MQ=255)
            ccGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTa                                                                                                            <  1:753730/139‑1 (MQ=255)
              ggTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCt                                                                                                               >  1:1473263/1‑134 (MQ=255)
              ggTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCt                                                                                                               <  2:1473263/134‑1 (MQ=255)
               gTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                         <  1:54776/139‑1 (MQ=255)
                             cGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTAt                                                                                           >  2:678149/1‑139 (MQ=255)
                                    ataaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGCGAGAGTAATTATCCGCGGc                                                                                    <  1:678149/139‑1 (MQ=255)
                                     taaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           <  2:226002/99‑1 (MQ=255)
                                     taaAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           >  1:226002/1‑99 (MQ=255)
                                           cccAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataata                                                                                                                      <  1:333318/98‑1 (MQ=255)
                                           cacAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACATTTGAACGAGCAGACTTGCTTTCCataata                                                                                                                      >  2:333318/3‑98 (MQ=255)
                                               ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           >  2:240347/1‑89 (MQ=255)
                                               ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           <  1:240347/89‑1 (MQ=255)
                                               ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                   <  2:154762/97‑1 (MQ=255)
                                               ggCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                   >  1:154762/1‑97 (MQ=255)
                                                    tAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAAtatta                                                                                                                   >  1:836599/1‑92 (MQ=255)
                                                    tAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGAATTGCTTTCCATAAtatta                                                                                                                   <  2:836599/92‑1 (MQ=255)
                                                                     gggATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   <  1:314758/139‑1 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                         >  1:1260107/1‑83 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTagag                                                                                                         <  2:1260107/83‑1 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   >  1:19257/1‑137 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   <  2:19257/137‑1 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   >  1:14344/1‑137 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   <  2:14344/137‑1 (MQ=255)
                                                                       gATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCt                                                 >  1:158407/1‑139 (MQ=255)
                                                                          gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTa                                                                                                            <  1:865629/77‑1 (MQ=255)
                                                                          gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTa                                                                                                            >  2:865629/1‑77 (MQ=255)
                                                                          gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                       >  2:230009/1‑82 (MQ=255)
                                                                          gACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAgagga                                                                                                       <  1:230009/82‑1 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           <  1:572744/61‑1 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           >  2:572744/1‑61 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 <  1:544235/71‑1 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 >  2:544235/1‑71 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                        <  2:191109/96‑1 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                        <  1:865434/96‑1 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                        >  1:191109/1‑96 (MQ=255)
                                                                           aCCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCATGTAGAGGAGAGAGTAATTATCcg                                                                                        >  2:865434/1‑96 (MQ=255)
                                                                            ccATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 <  1:597988/70‑1 (MQ=255)
                                                                            ccATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 >  2:597988/1‑70 (MQ=255)
                                                                               ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 <  1:1140788/67‑1 (MQ=255)
                                                                               ttGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATc                                                                                                                 >  2:1140788/1‑67 (MQ=255)
                                                                                tGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           <  1:610258/56‑1 (MQ=255)
                                                                                tGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCa                                                                                                                           >  2:610258/1‑56 (MQ=255)
                                                                                tGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCt                                                 >  2:1368948/1‑130 (MQ=255)
                                                                                tGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCt                                                 <  1:1368948/130‑1 (MQ=255)
                                                                                tGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAgt                                        >  1:962668/1‑139 (MQ=255)
                                                                                 gCAGTCAAAAAGCACCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATCATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAgtt                                       >  2:1017563/1‑139 (MQ=255)
                                                                                     tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCtgatg                                                      <  2:1493344/120‑1 (MQ=255)
                                                                                     tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCtgatg                                                      >  1:1493344/1‑120 (MQ=255)
                                                                                      cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataa                                                                                                                        >  1:343997/1‑53 (MQ=255)
                                                                                      cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCataa                                                                                                                        <  2:343997/53‑1 (MQ=255)
                                                                                      cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   >  2:749546/1‑122 (MQ=255)
                                                                                      cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGatgatg                                                   <  1:749546/122‑1 (MQ=255)
                                                                                       aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGa                                 >  1:856585/1‑139 (MQ=255)
                                                                                       aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGa                                 >  2:767537/1‑139 (MQ=255)
                                                                                               cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGAAAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttct                         >  2:269565/1‑139 (MQ=255)
                                                                                                  aaaaGCCTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttct                      <  2:856585/139‑1 (MQ=255)
                                                                                                       ccTGGACCTTTGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGca                 <  1:269565/139‑1 (MQ=255)
                                                                                                                ttGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATc                                                                                          <  2:873775/57‑1 (MQ=255)
                                                                                                                ttGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATc                                                                                          >  1:873775/1‑57 (MQ=255)
                                                                                                                 tGCATGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                        >  1:799165/1‑58 (MQ=255)
                                                                                                                 tGCATGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCcg                                                                                        <  2:799165/58‑1 (MQ=255)
                                                                                                                 tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                     <  1:1269101/45‑1 (MQ=255)
                                                                                                                 tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGgaga                                                                                                     >  2:1269101/1‑45 (MQ=255)
                                                                                                                 tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAAtt         <  1:1283154/137‑1 (MQ=255)
                                                                                                                 tGCACGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAAtt         >  2:1283154/1‑137 (MQ=255)
                                                                                                                    aCGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt    <  1:767537/139‑1 (MQ=255)
                                                                                                                    aCGAGCAGACTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt    <  2:158407/139‑1 (MQ=255)
                                                                                                                      gAGCAGCCTTGCTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTaa  <  1:1017563/139‑1 (MQ=255)
                                                                                                                        gCAGACTTGCTTTCCATAATATTTTCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCg                                  <  1:676270/105‑1 (MQ=255)
                                                                                                                        gCAGACTTGCTTTCCATAATATTTTCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCg                                  >  2:676270/1‑105 (MQ=255)
                                                                                                                                 |                                                                                                                               
TTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCATTGCAGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAA  >  NC_000913/3078723‑3078979

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