Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 855,061 A→G G228G (GGT→GGC)  ldtB ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913855,0610AG93.5% 91.5 / ‑3.0 31G228G (GGT→GGC) ldtBL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (15/14);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTACCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGTGCATACAGCCCC  >  NC_000913/854916‑855205
                                                                                                                                                 |                                                                                                                                                
cGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTgccgac                                                                                                                                             <  2:721449/151‑1 (MQ=255)
 gCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTgccgacc                                                                                                                                            <  2:759015/151‑1 (MQ=255)
   tctTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCgg                                                                                                                                          >  2:566104/1‑151 (MQ=255)
        gTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACtt                                                                                                                                     <  2:806025/151‑1 (MQ=255)
          cAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTtc                                                                                                                                   <  1:611254/151‑1 (MQ=255)
                    ggCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAgg                                                                                                                          <  1:708220/150‑1 (MQ=255)
                         aaTTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTg                                                                                                                    <  1:835248/151‑1 (MQ=255)
                           tttCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGAt                                                                                                                  >  2:604314/1‑151 (MQ=255)
                            ttCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATg                                                                                                                 >  1:592900/1‑151 (MQ=255)
                                 gACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTc                                                                                                            <  1:145074/151‑1 (MQ=255)
                                             gggCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACg                                                                                                <  2:592938/151‑1 (MQ=255)
                                                     gtggtAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAAc                                                                                        >  1:521799/1‑151 (MQ=255)
                                                      tggtAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAAcc                                                                                       >  1:693834/1‑151 (MQ=255)
                                                       ggtAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCa                                                                                      <  2:754253/151‑1 (MQ=255)
                                                        gtAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCAt                                                                                     >  2:699058/1‑151 (MQ=255)
                                                            aCAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACt                                                                                 <  2:306340/151‑1 (MQ=255)
                                                                        ggACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGcc                                                                     <  2:235954/151‑1 (MQ=255)
                                                                          aCTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCg                                                                    >  2:702206/1‑150 (MQ=255)
                                                                                 tataACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGa                                                             >  2:166633/1‑150 (MQ=255)
                                                                                  ataACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGaa                                                            <  2:330899/150‑1 (MQ=255)
                                                                                    aaCGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGtt                                                         <  2:521837/151‑1 (MQ=255)
                                                                                    aaCGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGtt                                                         <  2:169376/151‑1 (MQ=255)
                                                                                        gCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCg                                                     <  1:566066/151‑1 (MQ=255)
                                                                                                      cGGTGGTCGCTTTTACCTGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGAt                                       >  1:121260/1‑151 (MQ=255)
                                                                                                                  ttACCGGCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGc                            >  1:295443/1‑150 (MQ=255)
                                                                                                                     ccGGCTCATCAATAAACTGTACGCGGGGGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGAcc                        >  2:622199/1‑151 (MQ=255)
                                                                                                                        gCTCATCAATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGat                     >  2:363333/1‑151 (MQ=255)
                                                                                                                              caATAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATaga                <  2:738438/150‑1 (MQ=255)
                                                                                                                               aaTAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGt              >  2:346846/1‑151 (MQ=255)
                                                                                                                                aTAAACTGTACGCGGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGTg             >  1:786251/1‑151 (MQ=255)
                                                                                                                                  aaaCTGTACGCGGGTACCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGTGc            >  2:328520/1‑150 (MQ=255)
                                                                                                                                   caCTATACGCGTGAACGGCCAGGTAATTTCTCGAACAGGAATTTGATGTCTTCGTCACGCCGACGCGCACAACCCCGCCTTACACGCAGGCCGATAGCGCACCTGGCGTTGGTGCCATGGATAGAATACATAAGACAGatatctcgtcca           >  1:671439/2‑142 (MQ=255)
                                                                                                                                            cgGGTGCCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGTGCATACAGcccc  >  2:240386/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTTACCGGCTCATCAATAAACTGTACGCGGGTACCGACCGGTACTTTCTCGAACAGGAATTTGATGTCTTCGTTACGCAGACGCACACAACCATGACTTACACGCAGGCCGATACCGAAGTTGGCGTTGGTGCCATGGATAGCATACAGGCGACCGATATAGAGTGCATACAGCCCC  >  NC_000913/854916‑855205

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