Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,824,496 Δ1 bp coding (1081/1086 nt) mltB ← membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,824,4960G.97.7% 175.2 / ‑2.8 44coding (1081/1086 nt)mltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (20/23);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCGG  >  NC_000913/2824362‑2824615
                                                                                                                                      |                                                                                                                       
tcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTc                                                                                                                    >  2:236881/1‑139 (MQ=255)
           tGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                                         <  1:78001/139‑1 (MQ=255)
               ggACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACgg                                                                                                     <  1:167124/139‑1 (MQ=255)
                gACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGc                                                                                                    <  1:236881/139‑1 (MQ=255)
                gACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGc                                                                                                    <  2:28488/139‑1 (MQ=255)
                          cACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACgg                                                                                                     <  1:139009/128‑1 (MQ=255)
                          cACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACgg                                                                                                     >  2:139009/1‑128 (MQ=255)
                          cACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAAc                                                                                          >  1:460037/1‑139 (MQ=255)
                             tttCGTCAGGAAGAAGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGTAAGCTCACGCCCCCCTCTTGTAAATATTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGc                                                                                       >  2:740521/1‑139 (MQ=255)
                               tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCa                                                                                     >  1:254217/1‑139 (MQ=255)
                               tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCa                                                                                     >  2:151016/1‑139 (MQ=255)
                               tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCa                                                                                     >  2:1155047/1‑139 (MQ=255)
                                  tCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                  >  2:831354/1‑139 (MQ=255)
                                      gAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGcc                                                                              <  2:444595/139‑1 (MQ=255)
                                       aaGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCa                                                                             <  2:874095/139‑1 (MQ=255)
                                           tGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGc                                                                         <  2:624847/139‑1 (MQ=255)
                                            ggTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                  >  1:1015261/1‑129 (MQ=255)
                                            ggTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                  <  2:1015261/129‑1 (MQ=255)
                                                 tGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                                         >  1:1162527/1‑101 (MQ=255)
                                                 tGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                                         <  2:1162527/101‑1 (MQ=255)
                                                     gTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTg                                                               >  2:573629/1‑139 (MQ=255)
                                                              ttACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTg                                                      <  1:46615/139‑1 (MQ=255)
                                                                aCGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTa                                                    <  2:254217/139‑1 (MQ=255)
                                                                       gagacaggtaaaaaGCTGATTAGCCAGAGGCAAGATCAGGCCCCCCTCTTGTAAATAGTTACTGACTCGCGCCAGCGGCACGGCTTGTCCTAACTGCCAGACGGCCATTCCGTAATGGGAGCTGTGGTTGTAACGGGTg                                             <  2:813672/130‑1 (MQ=255)
                                                                        gaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGGAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                                            <  1:151016/139‑1 (MQ=255)
                                                                                  aaaGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCAt                                                                            >  2:372836/1‑97 (MQ=255)
                                                                                  aaaGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCAt                                                                            <  1:372836/97‑1 (MQ=255)
                                                                                         aTTAGCCAGAGGGCAACGCACGCCGCGCTTTTGCAAAGGGTTACTTACTCGCGCTAATACCCATGCTTGTCCTGACTGCAAGACCGCTATTTCGTAATGGGTGCTGTGGTTGTAACGGGTAATAGTGTAGaactgtctc                           >  1:813672/1‑132 (MQ=255)
                                                                                             gCCAGAGGTAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                       <  1:740521/139‑1 (MQ=255)
                                                                                             gCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                       <  1:573629/139‑1 (MQ=255)
                                                                                              ccAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAcc                      >  1:43102/1‑139 (MQ=255)
                                                                                                  agGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc                  <  2:43102/139‑1 (MQ=255)
                                                                                                  agGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc                  <  2:460037/139‑1 (MQ=255)
                                                                                                       aGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTa             >  1:597057/1‑139 (MQ=255)
                                                                                                       aGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTa             >  2:155/1‑139 (MQ=255)
                                                                                                        gCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                             >  2:508759/1‑90 (MQ=255)
                                                                                                        gCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                             <  1:508759/90‑1 (MQ=255)
                                                                                                            aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACgg                                                                                                     >  2:326987/1‑46 (MQ=255)
                                                                                                            aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACgg                                                                                                     <  1:326987/46‑1 (MQ=255)
                                                                                                            aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                                            >  2:599174/1‑103 (MQ=255)
                                                                                                            aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                                            <  1:599174/103‑1 (MQ=255)
                                                                                                                  cccTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCgg  <  1:155/139‑1 (MQ=255)
                                                                                                                            aaTAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAg               >  2:472418/1‑116 (MQ=255)
                                                                                                                            aaTAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAg               <  1:472418/116‑1 (MQ=255)
                                                                                                                                      |                                                                                                                       
TCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCGG  >  NC_000913/2824362‑2824615

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