Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,154,362 C→T A84V (GCC→GTC)  mltG → endolytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,154,3620CT96.2% 201.9 / ‑0.9 53A84V (GCC→GTC) mltGendolytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (26/25);  total (27/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 = 5.21e-01

TTTCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCT  >  NC_000913/1154316‑1154407
                                              |                                             
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            >  1:3754075/1‑50 (MQ=255)
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            >  1:3863859/1‑50 (MQ=255)
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGCCgg                                            >  1:4202782/1‑50 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            >  1:2391807/1‑49 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            >  1:1699418/1‑49 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            <  1:1908709/49‑1 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            >  1:3008908/1‑49 (MQ=255)
  tCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCggg                                           <  1:3857055/49‑1 (MQ=255)
    aaTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGAc                                         <  1:373497/49‑1 (MQ=255)
    aaTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACt                                        <  1:1139189/50‑1 (MQ=255)
     aTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACtt                                       <  1:3088157/50‑1 (MQ=255)
     aTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACt                                        >  1:2402165/1‑49 (MQ=255)
     aTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACt                                        >  1:3913953/1‑49 (MQ=255)
         ctgctgCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCg                                   <  1:4102671/50‑1 (MQ=255)
         ctgctgCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCg                                   <  1:4234725/50‑1 (MQ=255)
          tgctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGc                                  >  1:32834/1‑50 (MQ=255)
            ctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCtt                                <  1:3795652/50‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCttt                               <  1:1245685/48‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCttt                               <  1:2162588/48‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTac                             <  1:3567128/50‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTa                              <  1:4343833/49‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTac                             <  1:3734054/50‑1 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            >  1:2184960/1‑50 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            >  1:1553600/1‑50 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            >  1:639410/1‑50 (MQ=255)
                 tATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTacac                           <  1:1387480/50‑1 (MQ=255)
                   tCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACAcc                          <  1:2840324/49‑1 (MQ=255)
                   tCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACAcc                          <  1:3257603/49‑1 (MQ=255)
                       aCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGa                     <  1:1602072/50‑1 (MQ=255)
                        ccGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    <  1:4404059/50‑1 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    >  1:3944980/1‑49 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    >  1:2042810/1‑49 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    >  1:2654129/1‑49 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATg                   <  1:3681809/50‑1 (MQ=255)
                          ggATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGa                  >  1:2539568/1‑50 (MQ=255)
                           gATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAc                 <  1:3197434/50‑1 (MQ=255)
                            aTCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                >  1:1498428/1‑50 (MQ=255)
                            aTCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                >  1:3120634/1‑50 (MQ=255)
                            aTCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                >  1:1565420/1‑50 (MQ=255)
                              cTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGt              >  1:170301/1‑50 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                <  1:3148135/45‑1 (MQ=255)
                                    cACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:3479717/1‑48 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          <  1:3072545/47‑1 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCga        <  1:3276996/49‑1 (MQ=255)
                                      cTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:64960/1‑46 (MQ=255)
                                       tttAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGAt     >  1:1396812/1‑50 (MQ=255)
                                        ttAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATg    >  1:429904/1‑50 (MQ=255)
                                         tAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATg    <  1:629093/49‑1 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:2373235/1‑50 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:3233725/1‑50 (MQ=255)
                                           aaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:4289679/1‑49 (MQ=255)
                                           aaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  <  1:2734804/49‑1 (MQ=255)
                                           aaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  <  1:552404/49‑1 (MQ=255)
                                              |                                             
TTTCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCT  >  NC_000913/1154316‑1154407

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