Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,517,667 C→T 7.8% T286I (ACC→ATC)  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,6670CT 7.8% 143.8 / 2.2 51T286I (ACC→ATC) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (2/2);  total (26/25)
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

TCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTCAAGGTCACGCCGCAACTGCAGG  >  NC_002947/1517533‑1517813
                                                                                                                                       |                                                                                                                                                  
tcTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCAcc                                                                                                                                          <  8:189169/145‑1 (MQ=255)
tcTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGGCAcc                                                                                                                                          >  7:189169/1‑145 (MQ=255)
     cggcCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGt                                                                                                                                 <  3:84283/149‑1 (MQ=255)
     cggcCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGt                                                                                                                                 <  5:126021/149‑1 (MQ=255)
             ttCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAg                                                                                                                                   <  8:64424/139‑1 (MQ=255)
             ttCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAg                                                                                                                                   >  7:64424/1‑139 (MQ=255)
             ttCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTAt                                                                                                                         <  5:182615/149‑1 (MQ=255)
                gcgcgTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCtt                                                                                                                      >  4:314647/1‑149 (MQ=255)
                   cgTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTgcg                                                                                                                   <  6:161886/149‑1 (MQ=255)
                    gTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGCGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCACAAGGTCAGCGAATCTTgcg                                                                                                                   >  5:41332/1‑148 (MQ=255)
                       gCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGc                                                                                                               >  6:105775/1‑149 (MQ=255)
                          cAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCaa                                                                                                            <  3:314647/149‑1 (MQ=255)
                                  gACCATATCGAAAGCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCATGGTGAGCGTATCTTGCGCAGCCAACTGCTgg                                                                                                     <  6:161686/148‑1 (MQ=255)
                                      atatCGAAACCCTGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTAGCGCAGCCAACTGCTGGGCCGc                                                                                                <  7:206004/149‑1 (MQ=255)
                                         tCGAAACCCGGCT‑CTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTc                                                                                             >  7:211463/1‑148 (MQ=255)
                                              aCCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCt                                                                                        >  5:117301/1‑149 (MQ=255)
                                                ccGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCaa                                                                                           <  1:180113/144‑1 (MQ=255)
                                                ccGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCaa                                                                                           >  2:180113/1‑144 (MQ=255)
                                                          aGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGc                                                                             <  5:105775/148‑1 (MQ=255)
                                                              aCAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGcc                                                                        >  4:302944/1‑149 (MQ=255)
                                                               cAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTggg                                                                                                    >  3:132716/1‑120 (MQ=255)
                                                               cAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTggg                                                                                                    <  4:132716/120‑1 (MQ=255)
                                                                  ccTGGAAAGCCCAGTAGCGTCGTTGTGCTGCCGCGAAGCAGTATTTAGTGACGATGGGGTGCGCGC‑CACCCTCGTCACCGCCCAAGGGGAGCGTATCTTGCGCAGCCTACTTCTGGGCCGCTTCAACCTGAGCTACATGCTTGCCGCtg                                                                    <  6:41332/149‑1 (MQ=255)
                                                                  ccTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGc                                                                             <  1:325827/140‑1 (MQ=255)
                                                                  ccTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGc                                                                             >  2:325827/1‑140 (MQ=255)
                                                                        aaGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGa                                                              >  4:309466/1‑149 (MQ=255)
                                                                            ccAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCct                                                          >  8:178024/1‑149 (MQ=255)
                                                                             cAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTggg                                                                                                    <  1:61952/106‑1 (MQ=255)
                                                                             cAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTggg                                                                                                    >  2:61952/1‑106 (MQ=255)
                                                                              aGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCctgc                                                        <  3:302944/149‑1 (MQ=255)
                                                                                aaGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGc                                                                                    <  2:301340/119‑1 (MQ=255)
                                                                                aaGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGc                                                                                    >  1:301340/1‑119 (MQ=255)
                                                                                             cTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATg                                         >  7:270068/1‑149 (MQ=255)
                                                                                                 cgcgAAGCCGTTTTCAGTGACGATGGCGTGCGCGA‑AATACTCGTCACCGCCCTAGGTGAGCGTATCTTGCGCAGCAAACTGCTGGGACGCTTCAACCTGAGCTACAACCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATgc                                        >  3:64809/1‑146 (MQ=255)
                                                                                                 cgcgAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CATCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATgc                                        >  3:64812/1‑146 (MQ=255)
                                                                                                 cgcgAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CATCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATgc                                        <  4:64812/146‑1 (MQ=255)
                                                                                                 cgcgAAGCCGTATGCAGTCCCGATGGCGAGCGCGC‑CATCCTCGTCACCGCCCAACGTGAGCGTATCGTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGTGCAACATGATTGCCGCTGTGGCGACCCTGCTGGCGCTGGTCTAtca                                        <  4:64809/146‑3 (MQ=255)
                                                                                                  gcgAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTg                                    <  8:211463/149‑1 (MQ=255)
                                                                                                   cgAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACGGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTg                                    >  1:232441/1‑148 (MQ=255)
                                                                                                     aaGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGAc                                 <  3:309466/149‑1 (MQ=255)
                                                                                                        ccGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACgag                              <  2:232441/149‑1 (MQ=255)
                                                                                                           tATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGAt                            >  2:208487/1‑148 (MQ=255)
                                                                                                             ttCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATcc                          <  6:117301/148‑1 (MQ=255)
                                                                                                             ttCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCt                         >  2:86598/1‑149 (MQ=255)
                                                                                                                aGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGATCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTc                        >  6:278845/1‑147 (MQ=255)
                                                                                                                aGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGATCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTCaa                      <  5:278845/149‑1 (MQ=255)
                                                                                                                aGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTc                        >  4:216589/1‑147 (MQ=255)
                                                                                                                 gTGACGAGGGTG‑GCGCGCACGCCATCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTCAAg                     >  2:204976/1‑149 (MQ=255)
                                                                                                                  tGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATc                           >  2:147089/1‑142 (MQ=255)
                                                                                                                  tGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATc                           <  1:147089/142‑1 (MQ=255)
                                                                                                                                     cACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTCAAGGTCACGCCGCAACTGCAgg  >  8:150125/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
TCTGGCGGCCGACTTCGCGCGTCGCCCAAGTGTCGACCATATCGAAACCCGGCTGCTGAGCTACAGCCTGGAAAGCCCAGAAGCCTCGTTGTTCTGCCGCGAAGCCGTATTCAGTGACGATGGCGTGCGCGC‑CACCCTCGTCACCGCCCAAGGTGAGCGTATCTTGCGCAGCCAACTGCTGGGCCGCTTCAACCTGAGCAACATGCTTGCCGCTGTGGCGACCCTGCTGGCGCTGGACTATGCGCTGGACGAGATCCTCAAGGTCACGCCGCAACTGCAGG  >  NC_002947/1517533‑1517813

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