Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,619,764 C→T 51.4% E4K (GAA→AAA)  eutI ← predicted phosphotransacetylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,619,7640CT51.4% ‑2.5 / 57.2 65E4K (GAA→AAA) eutIpredicted phosphotransacetylase subunit
Reads supporting (aligned to +/- strand):  ref base C (21/10);  new base T (23/10);  total (45/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

CGTTGGTCTAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTC  >  minE/1619697‑1619824
                                                                   |                                                            
cGTTGGTCTAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTc                                                              >  1:55290/1‑68 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCa                                                      >  1:170155/1‑68 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:83505/1‑69 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:2436956/1‑69 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:814201/1‑69 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:2522258/1‑69 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:808646/1‑69 (MQ=255)
        tAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCAt                                                     >  1:166041/1‑69 (MQ=255)
                 cGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTct                                            >  1:2658405/1‑68 (MQ=255)
                 cGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTcc                                            >  1:2325388/1‑69 (MQ=255)
                 cGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTcc                                            >  1:2321470/1‑69 (MQ=255)
                 cGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTcc                                            >  1:230625/1‑69 (MQ=255)
                        aCCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGt                                      >  1:2566387/1‑68 (MQ=255)
                        aCCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGt                                      >  1:816210/1‑68 (MQ=255)
                        aCCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTc                                     >  1:3092969/1‑69 (MQ=255)
                        aCCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTc                                     >  1:508298/1‑69 (MQ=255)
                          cACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCgt                                   >  1:1095444/1‑69 (MQ=255)
                          cACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCgt                                   >  1:2360071/1‑69 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:2130220/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:3335605/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:3182965/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:3156240/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:2801519/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:1890847/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:2634814/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:1920895/66‑1 (MQ=255)
                                   gggcgCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:828370/66‑1 (MQ=255)
                                              aCGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTgg                             <  1:2869757/55‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:749619/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:638972/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:3292201/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:3187067/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:3122055/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:1683687/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:2497068/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:2557092/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:2546905/69‑1 (MQ=255)
                                                gCCAGTTCACGACAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTa             <  1:2521781/69‑1 (MQ=255)
                                                      tCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCgttgt                                >  1:1551174/1‑44 (MQ=255)
                                                      tCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCtt                     >  1:2049372/1‑55 (MQ=255)
                                                      tCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCgg         >  1:2073878/1‑67 (MQ=255)
                                                      tCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGgct       >  1:2466591/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCgttgt                                >  1:199845/1‑39 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCgttgt                                >  1:2805339/1‑39 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCgttgt                                >  1:258583/1‑39 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCgtt                                  >  1:2472485/1‑37 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCgt                                   >  1:2036748/1‑36 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGttgttg                              >  1:1951865/1‑41 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGa                            >  1:1583762/1‑43 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGAt                           >  1:2922875/1‑44 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTg                        >  1:3192478/1‑47 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTg                        >  1:2303241/1‑47 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTg                        >  1:1975812/1‑47 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCAcc                 >  1:2579861/1‑54 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCAcc                 >  1:2907775/1‑54 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAg            >  1:456883/1‑59 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAg            >  1:2049847/1‑59 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGgc        >  1:777136/1‑63 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:1458817/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:3415807/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:2704669/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:2507453/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:2406855/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTc  >  1:1948022/1‑69 (MQ=255)
                                                           aCAACGTTTAATAATCATTGCGTTTCCCCCAGTCGTGGTTggt                            >  1:2350056/1‑42 (MQ=38)
                                                           aCAACGTTAAATAATCATTGCGTTTCCCCCAGTCGttgttg                              >  1:989338/1‑41 (MQ=255)
                                                                   |                                                            
CGTTGGTCTAACGCATCCGGAAAAACCACTCTGGCGGGCGCTCGCAACGCCAGTTCACGACAACGTTCAATAATCATTGCGTTTCCCCCAGTCGTTGTTGGATCTGCTTCACCGTTAGCGGCTGCTTC  >  minE/1619697‑1619824

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