Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,700 T→C L323P (CTT→CCT)  slt → lytic murein transglycosylase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,7000TC95.6% 224.7 / 0.4 68L323P (CTT→CCT) sltlytic murein transglycosylase, soluble
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (37/28);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAT  >  NC_000913/4631603‑4631794
                                                                                                 |                                                                                              
ggCAACGATGTCACCGACGAGCATGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTg                                                                                              <  2:427008/100‑1 (MQ=255)
 gCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTgg                                                                                             >  1:1115032/1‑100 (MQ=255)
     cGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCAcc                                                                                         <  1:901114/100‑1 (MQ=255)
       aTGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCgg                                                                                       <  2:708860/100‑1 (MQ=255)
        tGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGc                                                                                      >  1:867595/1‑100 (MQ=255)
         gTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCg                                                                                     >  1:1084504/1‑100 (MQ=255)
         gTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCg                                                                                     >  2:512796/1‑100 (MQ=255)
         gTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCg                                                                                     >  2:400129/1‑100 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                    >  2:868052/1‑100 (MQ=255)
          tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGa                                                                                    >  1:1524144/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGAt                                                                                   >  1:526493/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                   >  1:1525971/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                   >  1:735920/1‑100 (MQ=255)
           cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAt                                                                                   >  1:660478/1‑100 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                  >  2:481254/1‑100 (MQ=255)
            aCCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtc                                                                                  <  1:726264/100‑1 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                 <  2:1185263/100‑1 (MQ=255)
             ccGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGAtcg                                                                                 >  1:1390491/1‑100 (MQ=255)
                      aGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCt                                                                        >  1:861894/1‑100 (MQ=255)
                       ggCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTg                                                                       <  2:277618/100‑1 (MQ=255)
                        gCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGa                                                                      <  2:107740/100‑1 (MQ=255)
                         cGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGaa                                                                     <  2:388881/100‑1 (MQ=255)
                          gAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAAt                                                                    <  1:512797/100‑1 (MQ=255)
                           aaaTGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATa                                                                   <  2:1524135/100‑1 (MQ=255)
                                   cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACctggctg                                                           >  1:871854/1‑100 (MQ=255)
                                   cgATGACGCCATTATGCGCTCGCAATCTACTTCGCTCATTGTACGCCGTGTACGAATCGCGCCTGGCACCTGCGATCGTCGCGGCCTGAATACctggctg                                                           >  2:161428/1‑100 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                      <  1:1580608/100‑1 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                      >  1:1232933/1‑100 (MQ=255)
                                        aCGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGcgcg                                                      <  1:139725/100‑1 (MQ=255)
                                              ttATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGcc                                                <  2:689555/100‑1 (MQ=255)
                                                aTGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTa                                              <  2:136348/100‑1 (MQ=255)
                                                  gcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                            >  1:217475/1‑100 (MQ=255)
                                                  gcgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATg                                            >  1:1093777/1‑100 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                           <  2:1277881/100‑1 (MQ=255)
                                                   cgcTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATgg                                           >  2:1093007/1‑100 (MQ=255)
                                                     cTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGaa                                         >  1:1077192/1‑100 (MQ=255)
                                                     cTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGaa                                         >  2:1468547/1‑100 (MQ=255)
                                                        gCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCg                                      <  1:1379499/100‑1 (MQ=255)
                                                           aTCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaa                                   <  2:63826/100‑1 (MQ=255)
                                                             cTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaaga                                 >  2:1331575/1‑100 (MQ=255)
                                                               aCTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCTTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagag                                <  2:865112/99‑1 (MQ=255)
                                                               aCTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagag                                >  2:1185777/1‑99 (MQ=255)
                                                                cTTCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaa                              >  1:1409287/1‑100 (MQ=255)
                                                                 ttCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaa                             <  2:340051/100‑1 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                            <  2:1466313/100‑1 (MQ=255)
                                                                  tCGCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                            >  2:463102/1‑100 (MQ=255)
                                                                   ngCTTATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCgaaagagaaag                            >  2:1276060/2‑99 (MQ=255)
                                                                       tATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGaa                       >  1:1263294/1‑100 (MQ=255)
                                                                       tATTGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGaa                       >  1:133553/1‑100 (MQ=255)
                                                                        aTTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAAt                      >  2:740331/1‑100 (MQ=255)
                                                                         ttGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                     <  2:595538/100‑1 (MQ=255)
                                                                         ttGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATg                     >  1:415019/1‑100 (MQ=255)
                                                                          tGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                    <  2:1420020/100‑1 (MQ=255)
                                                                          tGAACGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATgg                    >  1:1337229/1‑100 (MQ=255)
                                                                           gAACGCCGTGTACGAATGGGGCCTGGCACTGGCGATCGTCGCGGCCTGAATAGCTGGCTGGGGCGTCTGCTTATGGAAGCGAAAGAGAAAGATGAATGGc                   >  2:1314296/1‑100 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt                 >  1:280217/1‑100 (MQ=255)
                                                                             aCGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGt                 >  1:326633/1‑100 (MQ=255)
                                                                              cGCCGTGTACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGtt                <  2:1525962/100‑1 (MQ=255)
                                                                                  gtgtACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg            <  2:526492/100‑1 (MQ=255)
                                                                                   tgtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTg            <  2:380887/99‑1 (MQ=255)
                                                                                    gtACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGc          <  2:55834/100‑1 (MQ=255)
                                                                                     tACGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAGTGGCGTTACTGGCa         <  1:1468555/100‑1 (MQ=255)
                                                                                      aCGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAg        >  2:1121640/1‑100 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg       >  1:1498416/1‑100 (MQ=255)
                                                                                       cGAATGGCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAgg       >  1:284378/1‑100 (MQ=255)
                                                                                            ggCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt  <  1:1500301/100‑1 (MQ=255)
                                                                                            ggCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt  >  1:133963/1‑100 (MQ=255)
                                                                                            ggCGCCTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAt  <  2:153149/100‑1 (MQ=255)
                                                                                                 |                                                                                              
GGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAT  >  NC_000913/4631603‑4631794

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