Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,247 C→G 11.9% F168L (TTC→TTG) ‡ murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,521,2470CG11.9% 118.8 / 8.9 42F168L (TTC→TTG) ‡murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (17/20);  new base G (3/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

AACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGG  >  NC_002947/1521101‑1521389
                                                                                                                                                  |                                                                                                                                              
aaCGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccc                                                                                                                                              >  8:58997/1‑148 (MQ=255)
    ccAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTg                                                                                                                                          <  5:112978/149‑1 (MQ=255)
     cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTg                                                                                                                                          >  7:68080/1‑148 (MQ=255)
           cACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                   <  8:68080/149‑1 (MQ=255)
           cACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                   <  8:255708/149‑1 (MQ=255)
                  acgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagct                                                                                                                                   >  3:112919/1‑139 (MQ=255)
                  acgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagct                                                                                                                                   >  3:117562/1‑139 (MQ=255)
                  acgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagct                                                                                                                                   <  4:112919/142‑4 (MQ=255)
                  acgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagct                                                                                                                                   <  4:117562/142‑4 (MQ=255)
                    gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGc                                                                                                                          >  3:19916/1‑149 (MQ=255)
                       gCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                        >  6:8323/1‑148 (MQ=255)
                       gCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                        >  4:186158/1‑148 (MQ=255)
                                  aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTg                                                                                                            <  8:6795/149‑1 (MQ=255)
                                    aTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGc                                                                                                          <  5:23099/149‑1 (MQ=255)
                                     tGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCg                                                                                                         <  3:300429/149‑1 (MQ=255)
                                         cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                      >  6:238554/1‑116 (MQ=255)
                                                gCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGaa                                                                                              <  5:8323/149‑1 (MQ=255)
                                                gCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGaa                                                                                              <  4:203733/149‑1 (MQ=255)
                                                        aCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGc                                                                                       >  6:153796/1‑148 (MQ=255)
                                                                 cGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACcaca                                                                             >  4:76056/1‑149 (MQ=255)
                                                                   tGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATg                                                                           <  3:76056/149‑1 (MQ=255)
                                                                              cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGa                                                                         <  4:298138/140‑1 (MQ=255)
                                                                              cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGa                                                                         >  3:298138/1‑140 (MQ=255)
                                                                               tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGc                                                               >  7:275066/1‑149 (MQ=255)
                                                                                 gggaCCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCGGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             <  4:140957/146‑1 (MQ=255)
                                                                                 ggTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             <  4:19916/149‑1 (MQ=255)
                                                                                   tACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGcc                                                           >  4:149787/1‑149 (MQ=255)
                                                                                       ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTcc                                                       <  3:149787/149‑1 (MQ=255)
                                                                                         ggCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCgg                                                     <  5:153796/149‑1 (MQ=255)
                                                                                          gCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGc                                                    >  5:36501/1‑149 (MQ=255)
                                                                                                           cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCaaa                                    >  7:192757/1‑148 (MQ=255)
                                                                                                            gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCg                                                              <  5:194993/121‑1 (MQ=255)
                                                                                                            gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCg                                                              >  6:194993/1‑121 (MQ=255)
                                                                                                               gACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             >  1:207323/1‑119 (MQ=255)
                                                                                                               gACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             <  2:207323/119‑1 (MQ=255)
                                                                                                                     gAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtttt                         <  6:36501/149‑1 (MQ=255)
                                                                                                                                 cTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGccc              <  8:192757/148‑1 (MQ=255)
                                                                                                                                  tGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtata                                                                 <  2:114764/96‑1 (MQ=255)
                                                                                                                                  tGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtata                                                                 >  1:114764/1‑96 (MQ=255)
                                                                                                                                    gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc            >  3:253113/1‑147 (MQ=255)
                                                                                                                                    gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCGCAACGCCGAAGTGGCGCCCGTTCTGAACATCAGCGCAGCCCACATGGACCGCTATAGCGGCCGTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGccc          >  5:138791/1‑148 (MQ=255)
                                                                                                                                     aGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAg         <  8:275066/149‑1 (MQ=255)
                                                                                                                                        tttCCAGCTTCCAGCTGGAAAGCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg      >  2:302563/1‑149 (MQ=255)
                                                                                                                                            cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtgg  <  3:186158/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGG  >  NC_002947/1521101‑1521389

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: ‑