Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,824,491 C→T *362* (TAG→TAA 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,4910CT98.1% 161.5 / ‑5.0 52*362* (TAG→TAAmltBmembrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/26);  total (26/26)
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

GTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCGGTGCCAACATC  >  NC_000913/2824357‑2824625
                                                                                                                                      |                                                                                                                                      
gtcagtcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACt                                                                                                                                    >  1:612529/1‑139 (MQ=255)
  cagtcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACtta                                                                                                                                  <  1:374398/139‑3 (MQ=255)
   agtcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACttac                                                                                                                                 >  2:1352196/1‑136 (MQ=255)
                tGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTTCGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                                                   <  1:1484880/139‑1 (MQ=255)
                     gACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGc                                                                                                              <  2:1332637/139‑1 (MQ=255)
                               cACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAAc                                                                                                    >  2:1332761/1‑139 (MQ=255)
                                    tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGc                                                                                                              >  1:223569/1‑124 (MQ=255)
                                    tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGc                                                                                                              <  2:223569/124‑1 (MQ=255)
                                    tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCa                                                                                               >  1:1478582/1‑139 (MQ=255)
                                    tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCACCAGCGCCACGGCTTGTCCTAACTGCCa                                                                                               >  1:416551/1‑139 (MQ=255)
                                     cGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAg                                                                                              >  2:296657/1‑139 (MQ=255)
                                       tCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                            <  2:494/139‑1 (MQ=255)
                                              gATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                            <  2:340092/132‑1 (MQ=255)
                                              gATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                            >  1:340092/1‑132 (MQ=255)
                                                   ggtTGGCGTTGCAATTTTACGATGAGGATTCAGAAAAAAGCTGATTATCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCGCGGCTTGTCCTAACTGCCAGACCGCCATTGCGTa                                                                                <  1:442747/138‑1 (MQ=255)
                                                      tGGTGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATg                                                                             >  1:1080464/1‑139 (MQ=255)
                                                         cGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTg                                                                                                           >  2:107576/1‑106 (MQ=255)
                                                         cGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTg                                                                                                           <  1:107576/106‑1 (MQ=255)
                                                           ttGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGc                                                                        <  1:296657/139‑1 (MQ=255)
                                                              cAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCt                                                                                                       <  1:978575/105‑1 (MQ=255)
                                                              cAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTTACTCGCGCCAGCGCCACGGCTAGTCCt                                                                                                       >  2:978575/1‑105 (MQ=255)
                                                               aaTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCACCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg                                                                    <  2:416551/139‑1 (MQ=255)
                                                                    tACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCa                                                                                       <  1:1016699/115‑1 (MQ=255)
                                                                    tACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCa                                                                                       >  2:1016699/1‑115 (MQ=255)
                                                                           aggaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGt                                                        >  1:1266792/1‑139 (MQ=255)
                                                                             gaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTc                                                                                                         >  2:908499/1‑88 (MQ=255)
                                                                             gaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTc                                                                                                         <  1:908499/88‑1 (MQ=255)
                                                                                  aGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGgtg                                                 >  2:479036/1‑139 (MQ=255)
                                                                                         aGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCt                                                                                                             <  2:955729/72‑1 (MQ=255)
                                                                                         aGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCt                                                                                                             >  1:955729/1‑72 (MQ=255)
                                                                                             gATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCgg                                      <  1:1542912/139‑1 (MQ=255)
                                                                                                tAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAg                                   >  2:133457/1‑139 (MQ=255)
                                                                                                  gCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                                 <  1:274431/139‑1 (MQ=255)
                                                                                                     agagGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAAt                                                                              >  1:1540661/1‑91 (MQ=255)
                                                                                                     agagGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAAt                                                                              <  2:1540661/91‑1 (MQ=255)
                                                                                                     agagGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGt                              <  1:133457/139‑1 (MQ=255)
                                                                                                      gagGGAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa                             <  2:612529/139‑1 (MQ=255)
                                                                                                         ggAAGCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCa                          <  2:1367415/139‑1 (MQ=255)
                                                                                                             gCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                                 <  2:1238658/128‑1 (MQ=255)
                                                                                                             gCTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                                 >  1:1238658/1‑128 (MQ=255)
                                                                                                              cTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                                 >  2:1529446/1‑127 (MQ=255)
                                                                                                              cTCACGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc                                 <  1:1529446/127‑1 (MQ=255)
                                                                                                                 aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                            <  2:1268857/65‑1 (MQ=255)
                                                                                                                 aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGAc                                                                                            >  1:1268857/1‑65 (MQ=255)
                                                                                                                 aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                                                      >  2:1106050/1‑103 (MQ=255)
                                                                                                                 aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                                                      <  1:1106050/103‑1 (MQ=255)
                                                                                                                 aCGCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGt                  >  2:1352624/1‑139 (MQ=255)
                                                                                                                   gCCCCCCTCTTGTAAATAGTTACTTACTCgcgccag                                                                                                                       <  2:49710/36‑1 (MQ=255)
                                                                                                                   gCCCCCCTCTTGTAAATAGTTACTTACTCgcgccag                                                                                                                       >  1:49710/1‑36 (MQ=255)
                                                                                                                                 aaTAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                                       >  2:173399/1‑70 (MQ=255)
                                                                                                                                 aaTAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                                       <  1:173399/70‑1 (MQ=255)
                                                                                                                                 aaTAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCGGTGCCAACATc  >  1:427740/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAGTACTGGTAGCCGGTGCCAACATC  >  NC_000913/2824357‑2824625

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: