Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,599,234 C→T G70S (GGT→AGT)  opgB ← OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,599,2340CT96.2% 154.9 / ‑1.6 52G70S (GGT→AGT) opgBOPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (23/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

AGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACC  >  NC_000913/4599104‑4599362
                                                                                                                                  |                                                                                                                                
agcagGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGc                                                                                                                          >  2:1082628/1‑139 (MQ=255)
       tGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGt                                                                                                                   <  2:1083352/139‑1 (MQ=255)
        gTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGtt                                                                                                                  <  2:915568/139‑1 (MQ=255)
         tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACACCAGCACTGGTCAGGCTGTTGGtt                                                                                                                  >  1:705812/1‑138 (MQ=255)
         tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACACCAGCACTGGTCAGGCTGTTGGtt                                                                                                                  <  2:705812/138‑1 (MQ=255)
           aCCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAg                                                                                                               <  2:897152/139‑1 (MQ=255)
                aaTGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGt                                                                                                                               >  1:109036/1‑118 (MQ=255)
                aaTGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGt                                                                                                                               <  2:109036/118‑1 (MQ=255)
                aaTGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGtata                                                                                                          <  1:1537813/139‑1 (MQ=255)
                     tGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGaa                                                                                                     >  1:1092930/1‑139 (MQ=255)
                      gCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAAc                                                                                                    >  1:727094/1‑139 (MQ=255)
                       cGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAAcc                                                                                                   <  2:727094/139‑1 (MQ=255)
                         gatgatGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgc                                                                                                 <  1:1082628/139‑1 (MQ=255)
                          atgatgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTa                                                                                                                 >  2:1333979/1‑122 (MQ=255)
                          atgatgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTa                                                                                                                 <  1:1333979/122‑1 (MQ=255)
                          atgatgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                                                                                                <  2:1293245/139‑1 (MQ=255)
                              tgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAg                                                                                                               <  2:452835/120‑1 (MQ=255)
                              tgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAg                                                                                                               >  1:452835/1‑120 (MQ=255)
                                   acgacgGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATa                                                                                       <  2:1579782/139‑1 (MQ=255)
                                       cgGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCAt                                                                                   >  2:919064/1‑139 (MQ=255)
                                                      cccagcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTa                                                                                                                 >  1:1558280/1‑94 (MQ=255)
                                                      cccagcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTa                                                                                                                 <  2:1558280/94‑1 (MQ=255)
                                                           cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTg                                                               >  1:296421/1‑139 (MQ=255)
                                                                cGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCaaa                                                          <  2:1079140/139‑1 (MQ=255)
                                                                     acCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAggg                                                     >  1:258500/1‑139 (MQ=255)
                                                                        gCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGa                                                  >  1:1006633/1‑139 (MQ=255)
                                                                              aGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTa                                            <  2:344533/139‑1 (MQ=255)
                                                                                    aGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTg                                                   <  2:1610232/126‑1 (MQ=255)
                                                                                    aGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTg                                                   >  1:1610232/1‑126 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGt                                                                                                                   >  1:76406/1‑55 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGt                                                                                                                   <  2:76406/55‑1 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAATGTATAGAGAACcgcg                                                                                                <  2:1157430/74‑1 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCaaa                                                          >  1:622735/1‑112 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCaaa                                                          <  2:622735/112‑1 (MQ=255)
                                                                                           tGCCGATACCCGGCAAAATATATTTGCTGGCGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                                                                                                >  1:1157430/1‑74 (MQ=255)
                                                                                              cGATACCCGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGccc                            <  2:1092930/139‑1 (MQ=255)
                                                                                                aTACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACa                                                        >  2:1521576/1‑109 (MQ=255)
                                                                                                aTACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACa                                                        <  1:1521576/109‑1 (MQ=255)
                                                                                                     cGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAg                                                                                                               >  2:303425/1‑49 (MQ=255)
                                                                                                     cGGCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAg                                                                                                               <  1:303425/49‑1 (MQ=255)
                                                                                                       gCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                                                                                                >  2:1537612/1‑62 (MQ=255)
                                                                                                       gCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                                                                                                <  1:1537612/62‑1 (MQ=255)
                                                                                                       gCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgc                                                                                                 >  1:258205/1‑61 (MQ=255)
                                                                                                       gCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACcgc                                                                                                 <  2:258205/61‑1 (MQ=255)
                                                                                                       gCAAAATGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGt                   >  2:1621862/1‑139 (MQ=255)
                                                                                                             tGTATTTGCTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGt             <  2:258500/139‑1 (MQ=255)
                                                                                                                     cTGACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCaa     >  2:66406/1‑139 (MQ=255)
                                                                                                                       gACGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAGAATAGTCGCTGGCaaa                                                          <  1:1168278/84‑1 (MQ=255)
                                                                                                                       gACGCCAGCACTCATCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCaaa                                                          >  2:1168278/1‑84 (MQ=255)
                                                                                                                        aCGCCAGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAacc  <  2:296421/139‑1 (MQ=255)
                                                                                                                             aGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTa                                            <  1:1282441/92‑1 (MQ=255)
                                                                                                                             aGCACTGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTa                                            >  2:1282441/1‑92 (MQ=255)
                                                                                                                                  |                                                                                                                                
AGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACC  >  NC_000913/4599104‑4599362

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: