Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,282 C→G 11.8% W31C (TGG→TGC)  ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,2820CG11.8% 120.0 / 7.7 39W31C (TGG→TGC) ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (15/19);  new base G (2/3);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-02

CCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGT  >  NC_000913/1170067‑1170542
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                    
cagtgggggcgtgcATCTGGAGGTGGGGGTGGGAAAGGGGTGGTTAAAAGGCGTCCCGGTAGCGGTACGTAAGGGTCAACGCCGTGGGTAGTCTGTAACAGGGCGGGAAGTCCGACGTTGTATTTATTAGCGATGGCGTCCAGAGAACCACCGTCATTTTCCACCACATGATAGTTGTTCTCGCCAACCAGACGACTGCCCGCGGGCGGCAGCGGCCAGGTGTTGGATTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGt                                                                                                                                                                                                           <  4:2254/263‑1 (MQ=255)
 cTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGNCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCNCGCTGGCGGCACCGGCCAGGTGGTGGCTTTTGCCGGTAGCGCCAGNCNCCCGGCAGGGGCGACCGTaaaaaaacgt                                                                                                                                                                                                          >  3:48229/1‑271 (MQ=255)
   tCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAg                                                                                                                                                                                                        >  2:159461/1‑275 (MQ=255)
                      gTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGtt                                                                                                                                                                                     <  2:224655/275‑1 (MQ=255)
                          ggTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGACATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTNGGCTTTTGCCGGTAGCGCCATCGCCACGGCAGCGGCGAACGTAAAACAGGTAATCAAGAGAGAAAAACCAGttttta                                                                                                                                                                                 >  3:129731/1‑273 (MQ=255)
                              tGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCAGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGGTAGCCAGCGGAAAAAAcc                                                                                                                                                                                         >  4:224654/1‑262 (MQ=255)
                                     ggATCTGTAACCCGCGGCCCGCGGGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTGTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaa                                                                                                                                                                       <  4:70606/274‑1 (MQ=255)
                                               cacGGTGCCTGCGGGTGGTTCGTAAGGATCAACGCCGGGGTTAGCCTGTAAAAGAGCGAGAAAGCCTACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGACAGCCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                            <  4:39985/275‑1 (MQ=255)
                                               cacGCTGCCCTCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAGATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                            <  1:53515/275‑1 (MQ=255)
                                               cacGCTGCCCGCGGGGGGTACGTAATGATCAACGCCGGGGTTAGCCTGTAACAGGGCGTGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCAGTTTCCACCACATGAAATTTGTTGTCGCCAACCAGACGACTGCCCGCGGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCAACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataa                                                                                                                                                             <  4:371272/274‑1 (MQ=255)
                                                            cgcgGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGATAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTa                                                                                                                                               <  2:30243/275‑1 (MQ=255)
                                                                     tAAGGATCAACGCGGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                                      <  1:129731/275‑1 (MQ=255)
                                                                     tAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                                      <  4:159460/275‑1 (MQ=255)
                                                                           tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt                                                                                                                                >  4:450810/1‑275 (MQ=255)
                                                                                gCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGc                                                                                                                            <  2:534362/274‑1 (MQ=255)
                                                                                                     agCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAATTTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCGGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATACAAAATAACTATAAGGTTATAGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAgt                                                                                                      >  2:316019/1‑275 (MQ=255)
                                                                                                        gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTACAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCa                                                                                                                                                                                                                                                                    >  8:1078807/1‑114 (MQ=255)
                                                                                                        gagaAAGCCGACGTTGTATTTTTTGGCGATAGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCa                                                                                                                                                                                                                                                                    <  7:1078807/114‑1 (MQ=255)
                                                                                                           aaaGCGGACATGGTATTTTTTTGCGATGGCGTCCAGAGAACCATCGTCATTTTCCACCACATGAAATTTGGTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCACGGCCACGGCAGCGGCTAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCAACAATAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                                <  4:316018/275‑1 (MQ=255)
                                                                                                           aaaGCCGACGTTGTTTTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                                <  1:48229/275‑1 (MQ=255)
                                                                                                                       gTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGa                                                                                     <  2:450811/274‑1 (MQ=255)
                                                                                                                          ttttttGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                                                                                                                                                                                     >  7:425644/1‑95 (MQ=255)
                                                                                                                          ttttttGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                                                                                                                                                                                     <  8:425644/95‑1 (MQ=255)
                                                                                                                                  cGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGcc                                                                                                                                                                                                                                                    >  7:462695/1‑104 (MQ=255)
                                                                                                                                  cGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGcc                                                                                                                                                                                                                                                    <  8:462695/104‑1 (MQ=255)
                                                                                                                                    aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCa                                                                                                                                                                                                                               <  7:32366/123‑1 (MQ=255)
                                                                                                                                    aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCa                                                                                                                                                                                                                               >  8:32366/1‑123 (MQ=255)
                                                                                                                                                     accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                                      <  7:565889/67‑1 (MQ=255)
                                                                                                                                                     accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                                      >  8:565889/1‑67 (MQ=255)
                                                                                                                                                                   ccacATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGc                                        >  2:228770/1‑275 (MQ=255)
                                                                                                                                                                                                  ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGaa                                                                                                                                                                                                                       <  7:333472/68‑1 (MQ=25)
                                                                                                                                                                                                  ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGaa                                                                                                                                                                                                                       >  8:333472/2‑69 (MQ=25)
                                                                                                                                                                                                  acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGt                                                                                                                                                                                      <  6:717697/102‑1 (MQ=255)
                                                                                                                                                                                                  acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGt                                                                                                                                                                                      >  5:717697/1‑102 (MQ=255)
                                                                                                                                                                                                      gccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                                                                                                                                                                                                >  6:957125/2‑88 (MQ=255)
                                                                                                                                                                                                      gccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                                                                                                                                                                                                <  5:957125/87‑1 (MQ=255)
                                                                                                                                                                                                      cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                                  <  5:124565/102‑1 (MQ=255)
                                                                                                                                                                                                      cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                                  >  6:124565/1‑102 (MQ=255)
                                                                                                                                                                                                          cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGt  >  2:313950/1‑274 (MQ=255)
                                                                                                                                                                                                            ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaaa                                                                                                                                                                    <  7:632266/110‑1 (MQ=255)
                                                                                                                                                                                                            ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaaa                                                                                                                                                                    >  8:632266/1‑110 (MQ=255)
                                                                                                                                                                                                               ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGtt                                                                                                                                                                                     <  8:329281/90‑1 (MQ=255)
                                                                                                                                                                                                               ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGtt                                                                                                                                                                                     >  7:329281/1‑90 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                    
CCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTTCTCTTCCAGT  >  NC_000913/1170067‑1170542

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