Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,932 T→C 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,9320TC98.4% 219.1 / ‑4.3 62P33P (CCA→CCG) ydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (26/35);  total (26/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 = 7.74e-01

CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681836‑1682028
                                                                                                |                                                                                                
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                               <  2:1158526/100‑1 (MQ=255)
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                               >  1:142829/1‑100 (MQ=255)
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                               <  1:1017634/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              >  2:578635/1‑100 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              <  2:753907/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              <  2:420066/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              >  1:1034335/1‑100 (MQ=255)
   ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc                                                                                            <  1:764450/100‑1 (MQ=255)
      caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTTGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                         <  2:265201/100‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       <  2:296117/100‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       >  2:990044/1‑100 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       >  1:934834/1‑100 (MQ=255)
         aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                      >  1:389486/1‑100 (MQ=255)
          ctcaaaaTGATCGTTGCGCGCAACCCTTCAGTGCCGCGTTCGCGTGCAACAATACAATGCGCGATAAGCGCAAAGGTCTGGAAAAGCGGAATCAGCCCg                                                                                      <  2:823030/96‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    >  1:943438/1‑99 (MQ=255)
              aaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                                 <  2:1324487/100‑1 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                     >  1:1386418/1‑100 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                     <  1:1082952/100‑1 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                    <  2:1225779/100‑1 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAAGATAATTCGCGATAAGCGCAAAGGTCGTGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                    <  2:75999/100‑1 (MQ=255)
                              tAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtt                                                                 >  1:1030012/1‑100 (MQ=255)
                                aGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt                                                               <  1:449095/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATTAGCGATACGCGCAACGGTCGGGCAAAGCGGAGTCAGAACGGCGATATAATAATTTTTCGTTTTTGc                                                            >  2:543133/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             >  2:5733/1‑99 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            >  1:1464286/1‑100 (MQ=255)
                                     tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGAGAAAAGCGGGATCAGCCCGGCGATATAATAATTATTCGTTTTTGCTa                                                          >  2:1066244/1‑100 (MQ=255)
                                     gatctGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                          <  2:28968/96‑1 (MQ=255)
                                      cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                         <  1:578637/100‑1 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATTCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  1:518939/100‑1 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  1:1428002/100‑1 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  2:882122/100‑1 (MQ=255)
                                          gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaagatc                                                     >  1:28968/1‑96 (MQ=255)
                                           ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                    >  2:808612/1‑100 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   >  1:407904/1‑100 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   >  1:1223543/1‑100 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   <  2:1219736/100‑1 (MQ=255)
                                             gcgTTCGCTGGCAACAATAGAATGCGGGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCGGGCGATATAATAATTTTTCGTGTTTGCTAAAACACCa                                                  <  2:1084025/100‑1 (MQ=255)
                                              cgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCaa                                                 >  2:539087/1‑100 (MQ=255)
                                                ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                               <  2:1122239/100‑1 (MQ=255)
                                                   gcgggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            <  2:1364646/97‑1 (MQ=255)
                                                   gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            <  2:1148053/100‑1 (MQ=255)
                                                   gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            >  1:537025/1‑100 (MQ=255)
                                                    cTGGCAACAATATAATGTGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           >  2:995558/1‑100 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           <  2:450255/100‑1 (MQ=255)
                                                     tGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAACTTTTCGTTTTTGCTAAACCACCAATCAACAg                                          <  1:1062489/100‑1 (MQ=255)
                                                      ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                         <  1:808617/100‑1 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAACCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                       <  1:394955/100‑1 (MQ=255)
                                                          acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTa                                     >  1:777364/1‑100 (MQ=255)
                                                            aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                   <  2:619844/100‑1 (MQ=255)
                                                            aaTATAAGGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                   >  1:831684/1‑100 (MQ=255)
                                                               ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                                <  2:610471/100‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                            <  1:1347819/100‑1 (MQ=255)
                                                                          tAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgcgc                     >  1:121415/1‑100 (MQ=255)
                                                                             gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                  <  1:478196/100‑1 (MQ=255)
                                                                               gcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCtt                >  1:122396/1‑100 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATTAACAGCACTACCAGCGCACCGAGCGCGGCTTTg              <  2:1030006/100‑1 (MQ=255)
                                                                                     gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa          <  2:327588/100‑1 (MQ=255)
                                                                                     gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa          >  1:586333/1‑100 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt           <  2:138039/98‑1 (MQ=255)
                                                                                        tggAAAGGCGGAGGCAGCCCGGCGATATAATAATTTTTGGTATTTGCTAAAACACCAATCAGCAGCACTACCAGCGCACCTAGCGCGGCTTTGATTACCa       >  2:861423/2‑100 (MQ=255)
                                                                                        gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa       <  2:1117435/100‑1 (MQ=255)
                                                                                             aaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcccc  <  2:862038/100‑1 (MQ=255)
                                                                                                |                                                                                                
CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681836‑1682028

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: ‑