Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,932 T→C 100% P33P (CCA→CCG)  ydgC ← GlpM family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,9320TC94.6% 233.7 / 1.8 74P33P (CCA→CCG) ydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (32/38);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681836‑1682027
                                                                                                |                                                                                               
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                              <  1:1893741/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             >  1:2078532/1‑100 (MQ=255)
    cacaTACTATAAATGATGGTTGCGTGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc                                                                                           <  2:2195366/99‑1 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                         >  1:990622/1‑100 (MQ=255)
      caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCTATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                         <  2:1315482/99‑1 (MQ=255)
      caTACTAAAAATGAGGGTTGCGCGTAAGGGTTCAATGCCGCGTGCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                        <  2:720151/100‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                      <  2:541316/100‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                      <  2:1067222/100‑1 (MQ=255)
         aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCg                                                                                     >  2:699529/1‑100 (MQ=255)
         aCTAAAAATGATAGTTGCGCATAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATTCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                      >  2:228756/1‑99 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                    <  1:1341317/100‑1 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                    >  1:830682/1‑100 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                    <  1:2010706/100‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                  <  1:1224899/100‑1 (MQ=255)
              aaaTGAAGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                                <  2:520824/100‑1 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                              <  2:1078557/100‑1 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                              >  1:2050363/1‑100 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCAatat                                                                              <  1:1387300/100‑1 (MQ=255)
                  gatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataa                                                                            >  1:709192/1‑100 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                         >  2:867530/1‑99 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                        >  1:2230721/1‑100 (MQ=255)
                         gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAtttt                                                                     <  1:1506787/100‑1 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                    >  1:251736/1‑100 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                    <  1:270263/100‑1 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                   >  1:2253853/1‑100 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAAtttt                                                                     >  2:2094713/1‑97 (MQ=255)
                             gTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                  <  1:1172540/99‑1 (MQ=255)
                                 ggCTTCTATGCCGCGTTCGCTGGCAACAATATAATGCGCGAGAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttt                                                               <  1:137607/98‑1 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTTGttttt                                                             >  2:197968/1‑100 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                             <  1:2098714/100‑1 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  2:1867513/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           >  2:1442332/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  1:1164185/100‑1 (MQ=255)
                                    ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                          <  2:1121016/100‑1 (MQ=255)
                                       aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       <  2:2180962/100‑1 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa                                                      >  1:624657/1‑100 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  2:2328327/100‑1 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCCAAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  2:940299/100‑1 (MQ=255)
                                         tGCCACGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     >  2:1862740/1‑100 (MQ=255)
                                          gCCGCGTTCGCTGGCAATAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGGCCGGCGATATAATAATTTTTCGTTTATGCTaaaaa                                                     >  2:1879215/1‑98 (MQ=255)
                                           ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                   >  1:1377860/1‑100 (MQ=255)
                                           ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                   >  1:1079108/1‑100 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                  >  2:2202186/1‑100 (MQ=255)
                                               gTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAAt                                               >  1:1045484/1‑100 (MQ=255)
                                               gTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAAt                                               >  2:893043/1‑100 (MQ=255)
                                                 tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                             <  2:2373489/100‑1 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                          <  2:1031316/100‑1 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                          <  2:1702507/100‑1 (MQ=255)
                                                      ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                        <  1:699537/100‑1 (MQ=255)
                                                           caaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTAATACACCAATCAAAAAAACTAc                                   >  2:1058759/1‑100 (MQ=255)
                                                           caaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAc                                   >  2:390830/1‑100 (MQ=255)
                                                            aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                  <  1:2288002/100‑1 (MQ=255)
                                                            aaTATAATGCGCCATAAGCGCAAAGGTCTGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                  <  2:970400/100‑1 (MQ=255)
                                                             atatAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCa                                 >  1:975553/1‑100 (MQ=255)
                                                                taATGCGCGATAAGCGCAAAGGTAGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                              >  2:60032/1‑100 (MQ=255)
                                                                 aaTGCGCGATAAGCGAAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                              <  1:1073025/99‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                            >  1:986827/1‑99 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           <  1:1522280/100‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           <  1:1530445/100‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           >  1:1538292/1‑100 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           >  1:872816/1‑100 (MQ=255)
                                                                       cgATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAg                       <  1:1889455/100‑1 (MQ=255)
                                                                         aTAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgcg                     >  2:957680/1‑100 (MQ=255)
                                                                                cAAAGGTCGGGAAGAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt              <  1:2276775/100‑1 (MQ=255)
                                                                                  aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa            <  2:664438/100‑1 (MQ=255)
                                                                                   aGGTCGAGAAACGCGGAATCAGTCCGGCGATATAATAATTTTTCGATTTTGCTAAAACACCTATCAACAGCACTTTCAGCGCACCGAGCGCGGCTTTg             <  2:240548/98‑1 (MQ=255)
                                                                                     gTCGGGAAAAGCGTAATCAGCACGGCGACGTAATAATTTTTCGTTTTTGCTACAACACAAATCAACAGCACTACCAGCGCACCGAGCGCGACTTTGa            >  2:1991509/1‑97 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc        <  2:1377846/100‑1 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAACACTACCAGCGCACCGAGCGCGGCTTTGATTAc        >  2:2016011/1‑100 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc       <  1:741532/100‑1 (MQ=255)
                                                                                        gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa      >  1:560153/1‑100 (MQ=255)
                                                                                         ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa      >  2:2290513/1‑99 (MQ=255)
                                                                                           aaaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCTGCGCACCGAGCGCGGCTTTGATTACCAGcc   <  2:1417678/100‑1 (MQ=255)
                                                                                            aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGATAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc  <  2:1679757/100‑1 (MQ=255)
                                                                                                |                                                                                               
CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681836‑1682027

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