Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 79,533 T→C 54.9% F48F (TTT→TTC murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE79,5330TC54.9% 2.2 / 21.6 31F48F (TTT→TTCmurFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base C (2/15);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

GCTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGCAGGCGCACTA  >  minE/79465‑79584
                                                                    |                                                   
gCTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:679330/69‑1 (MQ=255)
 cTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:88475/68‑1 (MQ=255)
 cTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:3136844/68‑1 (MQ=255)
 cTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:346035/68‑1 (MQ=255)
 cTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:202082/68‑1 (MQ=255)
 cTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:204241/68‑1 (MQ=255)
  tGTAACCACTGACACGCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGa                                                   >  1:1121976/1‑69 (MQ=255)
  tGTAACCACTGACACGCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGa                                                   >  1:763409/1‑69 (MQ=255)
   gTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGAt                                                  >  1:562131/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:1093593/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:3130828/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:2256065/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:1813022/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:708952/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGAcc                               >  1:971296/1‑69 (MQ=255)
                      aaCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGGTTTGATGCCCACGATTTTGCCGAcc                               >  1:3278991/1‑69 (MQ=255)
                       aCTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCa                              <  1:692611/69‑1 (MQ=255)
                           aCGCCGGGCTGGCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:3010220/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATTATTTTg                                     <  1:3131007/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:561631/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:1103409/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:829489/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:1570661/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:1969102/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:2003847/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:3367727/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:3302978/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:2154420/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:237902/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:2475636/58‑1 (MQ=255)
                           aCGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:3124411/58‑1 (MQ=255)
                            cGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTCGATGCTCATGATTTTg                                     <  1:1887062/57‑1 (MQ=255)
                              ccGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:786266/39‑1 (MQ=255)
                               cGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGtttt                                                     <  1:2907841/38‑1 (MQ=255)
                                    tGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCa                              <  1:2122265/56‑1 (MQ=255)
                                                   cTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGCAGGCGCacta  <  1:677724/69‑1 (MQ=255)
                                                   cTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGCAGGCGCacta  <  1:1443795/69‑1 (MQ=255)
                                                   cTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGCAGGCGCacta  <  1:2596198/69‑1 (MQ=255)
                                                   cTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGAAGGCGCacta  <  1:1208579/69‑1 (MQ=255)
                                                                    |                                                   
GCTGTAACCACTGATACCCGAAAACTGACGCCGGGCTGCCTGTTTGTTGCCCTGAAAGGCGAACGTTTTGATGCCCACGATTTTGCCGACCAGGCGAAAGCTGGCGGCGCAGGCGCACTA  >  minE/79465‑79584

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