Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 882,856 C→A R178L (CGT→CTT)  ybjG ← undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913882,8560CA98.4% 224.2 / ‑3.4 63R178L (CGT→CTT) ybjGundecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (27/35);  total (27/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 = 9.17e-01

AGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGG  >  NC_000913/882761‑882950
                                                                                               |                                                                                              
aGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTg                                                                                            <  1:146441/100‑1 (MQ=255)
 gCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGa                                                                                           <  2:784336/100‑1 (MQ=255)
 gCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGa                                                                                           <  1:1901547/100‑1 (MQ=255)
  ccGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGAt                                                                                          <  2:1073314/100‑1 (MQ=255)
  ccGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGAt                                                                                          <  1:1265553/100‑1 (MQ=255)
  ccGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGAt                                                                                          >  1:1849534/1‑100 (MQ=255)
  ccGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGAt                                                                                          <  1:1732912/100‑1 (MQ=255)
   cGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGAtn                                                                                         <  2:149330/100‑2 (MQ=255)
    gCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATag                                                                                        <  2:1177846/100‑1 (MQ=255)
        aaTACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGtt                                                                                    <  2:369984/100‑1 (MQ=255)
        aaTACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGtn                                                                                    <  1:880805/100‑2 (MQ=255)
            ctacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATg                                                                                >  2:1653573/1‑100 (MQ=255)
              acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGcc                                                                              >  2:116339/1‑100 (MQ=255)
              acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGcc                                                                              <  1:1231865/100‑1 (MQ=255)
               cacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGccc                                                                             >  2:481798/1‑100 (MQ=255)
                  tttttAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCAt                                                                          >  1:1254801/1‑100 (MQ=255)
                         cAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGc                                                                   >  2:773721/1‑100 (MQ=255)
                         cAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGc                                                                   <  2:1849518/100‑1 (MQ=255)
                          agagATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGcc                                                                  >  1:309638/1‑100 (MQ=255)
                          agagATCAGTCACGCACCCAGCCTCTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGcc                                                                  <  1:1541571/100‑1 (MQ=255)
                           gagaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCa                                                                 <  2:1306052/100‑1 (MQ=255)
                           gagaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCa                                                                 >  2:887569/1‑100 (MQ=255)
                           gagaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCa                                                                 <  1:1906749/100‑1 (MQ=255)
                            agaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAg                                                                <  1:428659/100‑1 (MQ=255)
                             gaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGa                                                               >  2:1383625/1‑100 (MQ=255)
                              aTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGAt                                                              >  1:1235781/1‑100 (MQ=255)
                              aTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGAt                                                              <  1:481799/100‑1 (MQ=255)
                               tCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATa                                                             <  1:1901703/100‑1 (MQ=255)
                                   tCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCt                                                         <  1:1292261/100‑1 (MQ=255)
                                    cACGCACCCAGCCTTTGTGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTg                                                        >  2:1719017/1‑100 (MQ=255)
                                    cACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTg                                                        <  2:1051558/100‑1 (MQ=255)
                                        cACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCa                                                    >  1:995209/1‑100 (MQ=255)
                                          cccAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACt                                                  >  1:1882143/1‑100 (MQ=255)
                                           ccAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTa                                                 <  1:1634255/100‑1 (MQ=255)
                                            cAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTaa                                                >  2:1503553/1‑100 (MQ=255)
                                               ccTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGc                                             <  2:971926/100‑1 (MQ=255)
                                                  ttGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGc                                          <  1:1451120/100‑1 (MQ=255)
                                                  ttGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGc                                          <  2:1892400/100‑1 (MQ=255)
                                                   tGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGcc                                         <  1:1503562/100‑1 (MQ=255)
                                                       gATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATc                                     >  2:616428/1‑100 (MQ=255)
                                                       gATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATc                                     <  2:1235775/100‑1 (MQ=255)
                                                       gATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATc                                     <  2:6687/100‑1 (MQ=255)
                                                          cGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCAt                                   <  2:1894598/99‑1 (MQ=255)
                                                          cGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATa                                  <  1:1841745/100‑1 (MQ=255)
                                                              aaTGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTg                              >  2:1620292/1‑100 (MQ=255)
                                                              aaTGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTg                              >  2:177124/1‑100 (MQ=255)
                                                                tGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGcc                            <  2:1127961/100‑1 (MQ=255)
                                                                 gCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCa                           >  2:162987/1‑100 (MQ=255)
                                                                  cAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAg                          >  2:773339/1‑100 (MQ=255)
                                                                  cAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAg                          <  2:816039/100‑1 (MQ=255)
                                                                   aaaaCAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGc                         >  1:399605/1‑100 (MQ=255)
                                                                    aaaCAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCa                        >  2:815786/1‑100 (MQ=255)
                                                                      aCAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAAt                      <  1:1620176/100‑1 (MQ=255)
                                                                       cAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGTCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATc                     <  1:1520279/100‑1 (MQ=255)
                                                                           cgcgATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCAcc                 >  2:1059304/1‑100 (MQ=255)
                                                                           cgcgATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCAcc                 <  1:2040859/100‑1 (MQ=255)
                                                                             cgATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGa               >  2:2028565/1‑100 (MQ=255)
                                                                                  ccACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCata          >  1:856544/1‑100 (MQ=255)
                                                                                  ccACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCat           >  2:1866796/1‑99 (MQ=255)
                                                                                    aCGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATc        <  2:995205/100‑1 (MQ=255)
                                                                                        ttGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGc    >  2:234420/1‑100 (MQ=255)
                                                                                        ttGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGc    <  1:670442/100‑1 (MQ=255)
                                                                                          gCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCgg  >  2:1761668/1‑100 (MQ=255)
                                                                                               |                                                                                              
AGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGG  >  NC_000913/882761‑882950

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