Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,721,332 T→C 55.0% L83L (CTA→CTG) ‡ yhfK ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,721,3320TC55.0% 8.5 / 38.6 51L83L (CTA→CTG) ‡yhfKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/3);  new base C (14/14);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 7.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGTT  >  minE/2721269‑2721399
                                                               |                                                                   
gACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCtg                                                                >  1:113147/1‑69 (MQ=255)
gACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCtg                                                                >  1:240255/1‑69 (MQ=255)
gACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCtg                                                                >  1:2224734/1‑69 (MQ=255)
gACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCtg                                                                >  1:2137574/1‑69 (MQ=255)
 aCGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCtgt                                                               <  1:1257077/69‑1 (MQ=255)
         tACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCa                                                        <  1:218316/68‑1 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:1634803/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:263691/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:818271/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:2726655/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:3166770/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:3353646/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:1364474/1‑69 (MQ=255)
           ccAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGc                                                     >  1:931216/1‑69 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:3138316/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:2949960/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:2894284/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:2307544/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:3233153/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:2718088/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:612600/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:1041790/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:1894116/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:152775/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:164596/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:1846930/69‑1 (MQ=255)
                 ttAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:1774813/69‑1 (MQ=255)
                   aaTCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAgg                                               <  1:2015292/67‑1 (MQ=255)
                       cGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCa                                          >  1:3374293/1‑68 (MQ=255)
                       cGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCa                                          >  1:910196/1‑68 (MQ=255)
                       cGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCa                                          >  1:1548599/1‑68 (MQ=255)
                       cGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCa                                          >  1:409884/1‑68 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCt                                                    >  1:350621/1‑51 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:565821/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:1387750/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:3174887/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:220833/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:2609987/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:1685294/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:3827/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:1662805/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:3389072/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:67188/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:155524/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:2894738/1‑69 (MQ=255)
                              cAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGCACCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGa                                  >  1:1218085/1‑69 (MQ=255)
                                     ggCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCa                           <  1:1901879/69‑1 (MQ=255)
                                                 tCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGt               >  1:692391/1‑69 (MQ=255)
                                                 tCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGt               >  1:912595/1‑69 (MQ=255)
                                                 tCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGt               >  1:521745/1‑69 (MQ=255)
                                                 tCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGAGt               >  1:255072/1‑69 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGAc                                 >  1:260933/1‑39 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:369306/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:982476/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:2927146/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:1314864/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:1506016/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:2275184/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:2518440/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACg    >  1:3136741/1‑68 (MQ=255)
                                                             agtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACAAATGATTAAGCGTTTGAAAAAACg    >  1:1289435/1‑68 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:257974/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:2510601/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:245895/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:2659971/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:1900487/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  <  1:1650177/69‑1 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:875199/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:272042/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:1495820/1‑69 (MQ=255)
                                                              gtagCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGtt  >  1:947535/1‑69 (MQ=255)
                                                               |                                                                   
GACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAGTAGTAGCTGTGTCAGCAAGCTACAGGTGGCAAACAGCGACGCACCAATGATTAAGCGTTTGAAAAAACGTT  >  minE/2721269‑2721399

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: