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,3620CT95.2% 301.6 / 3.7 83A84V (GCC→GTC) mltGendolytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (33/46);  total (34/49)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TTTCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCT  >  NC_000913/1154316‑1154407
                                              |                                             
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCg                                             <  1:3256261/49‑1 (MQ=255)
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCg                                             <  1:1788496/49‑1 (MQ=255)
tttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCg                                             <  1:3922740/49‑1 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            <  1:772408/49‑1 (MQ=255)
 ttCAATGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCgg                                            <  1:5279297/49‑1 (MQ=255)
    aaTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGAc                                         <  1:1215017/49‑1 (MQ=255)
    aaTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACt                                        <  1:2846129/50‑1 (MQ=255)
     aTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACtt                                       <  1:2071222/50‑1 (MQ=255)
     aTGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACtt                                       <  1:5151334/50‑1 (MQ=255)
      tGGCTGCTGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTa                                      >  1:5930623/1‑50 (MQ=255)
        gctgctGCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTAc                                     <  1:252875/49‑1 (MQ=255)
         ctgctgCGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCg                                   <  1:1181441/50‑1 (MQ=255)
          tgctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGc                                  >  1:490707/1‑50 (MQ=255)
          tgctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGc                                  >  1:4058837/1‑50 (MQ=255)
           gctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCt                                 >  1:5885196/1‑50 (MQ=255)
           gctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCt                                 >  1:1934673/1‑50 (MQ=255)
            ctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGc                                  <  1:1364510/48‑1 (MQ=255)
            ctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCtt                                >  1:4143349/1‑50 (MQ=255)
            ctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCt                                 <  1:4692519/49‑1 (MQ=255)
            ctgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCt                                 <  1:3115874/49‑1 (MQ=255)
             tgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCtt                                >  1:2531284/1‑49 (MQ=255)
             tgcGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCtt                                >  1:1106985/1‑49 (MQ=255)
               cGTATCGAACCGGATCTTTCTCGCTTTAAAGTCGGGACTTACCGCTTTa                              <  1:1160988/49‑1 (MQ=38)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCttt                               <  1:1217970/48‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCttt                               <  1:3580933/48‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTa                              <  1:1905310/49‑1 (MQ=255)
               cGTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTa                              <  1:4327672/49‑1 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            <  1:50096/50‑1 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            >  1:5774303/1‑50 (MQ=255)
                gTATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTaca                            <  1:420450/50‑1 (MQ=255)
                 tATCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTacac                           <  1:4343202/50‑1 (MQ=255)
                  aTCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACAcc                          >  1:1667534/1‑50 (MQ=255)
                   tCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACAcc                          <  1:2039875/49‑1 (MQ=255)
                   tCGAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         <  1:4980382/50‑1 (MQ=255)
                     gAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:5126323/1‑48 (MQ=255)
                     gAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:3489604/1‑48 (MQ=255)
                     gAACCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCa                       <  1:5539395/50‑1 (MQ=255)
                      aaCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:1175259/1‑47 (MQ=255)
                      aaCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:2381956/1‑47 (MQ=255)
                      aaCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAg                      <  1:120376/50‑1 (MQ=255)
                      aaCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAg                      <  1:5162383/50‑1 (MQ=255)
                       aCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:268237/1‑46 (MQ=255)
                       aCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCg                         >  1:4444087/1‑46 (MQ=255)
                       aCCGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGa                     <  1:3958206/50‑1 (MQ=255)
                        ccGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    <  1:2701289/50‑1 (MQ=255)
                        ccGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    <  1:2182045/50‑1 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGAt                    >  1:155734/1‑49 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATg                   <  1:5248278/50‑1 (MQ=255)
                         cGGATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATg                   <  1:459677/50‑1 (MQ=255)
                          ggATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGa                  >  1:148040/1‑50 (MQ=255)
                          ggATCTTTCTCACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGa                  <  1:551316/50‑1 (MQ=255)
                           gATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAc                 >  1:1803562/1‑50 (MQ=255)
                           gATCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAc                 >  1:2866931/1‑50 (MQ=255)
                            aTCTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                >  1:3255398/1‑50 (MQ=255)
                              cTTTCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGt              >  1:2729999/1‑50 (MQ=255)
                               tttCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTg             >  1:4257733/1‑50 (MQ=255)
                                ttCTCACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgc            >  1:3473583/1‑50 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                <  1:4875863/45‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGAcc                <  1:2111144/45‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCg               <  1:2878286/46‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCg               <  1:618345/46‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGt              <  1:5529512/47‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcg           >  1:383056/1‑50 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgc            <  1:2071024/49‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgc            <  1:1691220/49‑1 (MQ=255)
                                 tctcACTTTAAAGTCGGGACTTACCGCTTTACACAGCAGATGACCg               <  1:4491209/46‑1 (MQ=255)
                                 tctcACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGACCGt              <  1:2771266/47‑1 (MQ=255)
                                  ctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:3428608/1‑50 (MQ=255)
                                  ctcACTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcg           <  1:4616308/49‑1 (MQ=255)
                                    cACTTTAAAGCCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:199859/1‑48 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          <  1:3981128/47‑1 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCgag       <  1:2710741/50‑1 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCgag       <  1:2488039/50‑1 (MQ=255)
                                     aCTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCga        <  1:2220212/49‑1 (MQ=255)
                                      cTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:5732742/1‑46 (MQ=255)
                                      cTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTgcgc          >  1:109241/1‑46 (MQ=255)
                                      cTTTAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCgaga      <  1:2804730/50‑1 (MQ=255)
                                         tAAAGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGc   >  1:2933618/1‑50 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:1656597/1‑50 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:347503/1‑50 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:3460984/1‑50 (MQ=255)
                                          aaaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  >  1:179630/1‑50 (MQ=255)
                                           aaGTCGGGACTTACCGCTTTACACCGCAGATGACCGTGCGCGAGATGCt  <  1:4418672/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: