Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,517,528 C→T 7.8% R240C (CGT→TGT)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,517,5280CT 7.8% 141.6 / 2.2 51R240C (CGT→TGT) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/2);  total (27/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.24e-01

GGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGT  >  NC_002947/1517380‑1517673
                                                                                                                                                     |                                                                                                                                                 
ggAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGgtc                                                                                                                                                   >  4:100531/1‑149 (MQ=255)
   aTTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGgtcgtc                                                                                                                                                <  7:315198/149‑1 (MQ=255)
           tCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTAAGAGGCCGCCAAGGCCAAGCTGTCACGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTggc                                                                                                                                         <  1:272268/149‑1 (MQ=255)
             gcgGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCCTCTggcggc                                                                                                                                      >  4:93415/1‑149 (MQ=255)
              cgGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGcc                                                                                                                                     >  4:5657/1‑149 (MQ=255)
                      aCCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgc                                                                                                                             <  5:272404/149‑1 (MQ=255)
                      aCCAACCTGTCCCGCGACCACCTCGACTATCACCGCAGCATGGAAGCCTATGATGCCGCCAAGGCCAAGCTGTTAGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgc                                                                                                                             >  6:192429/1‑149 (MQ=255)
                        cAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTcgcgc                                                                                                                           >  1:34496/1‑149 (MQ=255)
                        cAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGACTGGCGGCCGACTTcgcgc                                                                                                                           >  3:45819/1‑149 (MQ=255)
                                 ccGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCa                                                                                                                   <  3:100531/148‑1 (MQ=255)
                                             cGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTGGAccttt                                                                                                      >  8:254443/1‑146 (MQ=255)
                                             cGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCata                                                                                                      >  8:259139/1‑149 (MQ=255)
                                               aCTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc                                                                                                    >  1:206864/1‑149 (MQ=255)
                                                 tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCg                                                                                                   >  4:166978/1‑148 (MQ=255)
                                                 tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGCCTGGCGGCCGACTTCGCGCGTCGCCCaa                                                                                                                  >  4:116925/1‑133 (MQ=255)
                                                 tATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGCCTGGCGGCCGACTTCGCGCGTCGCCCaa                                                                                                                  <  3:116925/133‑1 (MQ=255)
                                                   tCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGaaa                                                                                                >  3:208142/1‑149 (MQ=255)
                                                    cACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAAc                                                                                               >  8:101005/1‑149 (MQ=255)
                                                    cACGGCAGCATGGAAGCCTATGAGGCCGCCAACGCCAAGCAGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCCTATCGAAAc                                                                                               >  7:211654/1‑149 (MQ=255)
                                                             aTGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGCCTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGgctgct                                                                                      >  7:289465/1‑149 (MQ=255)
                                                             aTGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGgctgct                                                                                      >  4:63251/1‑149 (MQ=255)
                                                                gAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGa                                                                                    <  7:101005/148‑1 (MQ=255)
                                                                 aaGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGc                                                                                  >  7:117756/1‑149 (MQ=255)
                                                                 aaGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGc                                                                                  <  4:208142/149‑1 (MQ=255)
                                                                   gCCTCTGCGGCCGCCAAGGCCATGCTGTTCGCCTGGCCGTGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTa                                                                                <  8:211654/149‑1 (MQ=255)
                                                                     cTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc                                                                                                    >  8:206996/1‑127 (MQ=255)
                                                                     cTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATc                                                                                                    <  7:206996/127‑1 (MQ=255)
                                                                     cTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACa                                                                              <  3:63251/149‑1 (MQ=255)
                                                                     cTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACa                                                                              <  3:166978/149‑1 (MQ=255)
                                                                      tATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAg                                                                             >  5:277037/1‑149 (MQ=255)
                                                                          aGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCt                                                                          <  2:206864/148‑1 (MQ=255)
                                                                          aGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTg                                                                         <  3:5657/149‑1 (MQ=255)
                                                                                       caagCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTGCAGCCTGGAAAGCCCAGAAg                                                            >  4:83339/1‑149 (MQ=255)
                                                                                           cTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTc                                                        <  7:259139/149‑1 (MQ=255)
                                                                                           cTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTTGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTc                                                        <  7:254443/149‑1 (MQ=255)
                                                                                              ttCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgtt                                                     <  2:34496/149‑1 (MQ=255)
                                                                                              ttCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgt                                                      <  3:93415/148‑1 (MQ=255)
                                                                                                cagcttGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCgttg                                                    <  7:201665/144‑1 (MQ=255)
                                                                                                   ccTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTc                                                 <  8:117756/149‑1 (MQ=255)
                                                                                                       gCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGcc                                             <  5:192429/149‑1 (MQ=255)
                                                                                                             ccTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAg                                       >  7:231511/1‑149 (MQ=255)
                                                                                                             ccTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAg                                       >  8:141063/1‑149 (MQ=255)
                                                                                                               tGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCTTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCAAGAAGCCTAGTTGTTCTGCCGAGAAGcc                                     >  5:319907/1‑149 (MQ=255)
                                                                                                                 cGTTGCCAGTTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCg                                    <  8:231511/148‑1 (MQ=255)
                                                                                                                      ccAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTAtt                                >  2:341352/1‑147 (MQ=255)
                                                                                                                             gtCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGa                       >  6:37248/1‑149 (MQ=255)
                                                                                                                                aaCCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATg                     <  1:341352/148‑1 (MQ=255)
                                                                                                                                  ccTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCg                  >  2:165737/1‑149 (MQ=255)
                                                                                                                                     ggATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTgc               >  7:257526/1‑149 (MQ=255)
                                                                                                                                       atgatgACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTgcgc             <  6:277037/149‑1 (MQ=255)
                                                                                                                                        tgatgaCTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCg                                                       <  6:111898/106‑1 (MQ=255)
                                                                                                                                        cgatgaCTTGGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCg                                                       <  6:111899/105‑1 (MQ=255)
                                                                                                                                               tttGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCt     >  8:257170/1‑149 (MQ=255)
                                                                                                                                                   gtcgtcgtcTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGt  <  7:141063/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTCGACTATCACGGCAGCATGGAAGCCTATGAGGCCGCCAAGGCCAAGCTGTTCGCCTGGCCGAGCCTGCGTTGCCAGGTGGTCAACCTGGATGATGACTTTGGTCGTCGTCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGCCACCCTCGT  >  NC_002947/1517380‑1517673

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