Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,330,467 A→G 48.1% F373F (TTT→TTC)  cpsG ← phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,330,4670AG48.1% ‑2.5 / 46.4 54F373F (TTT→TTC) cpsGphosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (17/11);  new base G (14/12);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGTTCCACGCGGTTAATCGCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTC  >  minE/1330399‑1330529
                                                                    |                                                              
tGCTGTTCCACGCGGTTAATCGCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGa                                                                <  1:3234209/69‑1 (MQ=255)
                tAATCGCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCg                                                  >  1:1478773/1‑67 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2047239/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:3164265/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2771218/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:3395409/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2651155/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:3413344/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2393435/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2337217/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:2085712/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:431936/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:1435839/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:1497946/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:80729/1‑68 (MQ=255)
                     gCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGGAATGCCGCCATCCGGTcgcgc                                           >  1:1174611/1‑68 (MQ=255)
                                ttGCGCGAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCa                                       >  1:438162/1‑62 (MQ=255)
                                ttGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCa                                       >  1:375788/1‑62 (MQ=255)
                                ttGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCa                                       >  1:479278/1‑62 (MQ=255)
                                ttGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAg                                      >  1:2617209/1‑63 (MQ=255)
                                              tGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:1405103/68‑1 (MQ=255)
                                              tGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:2334857/68‑1 (MQ=255)
                                              tGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:1966921/68‑1 (MQ=255)
                                              tGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:1945338/68‑1 (MQ=255)
                                              tGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:1708360/68‑1 (MQ=255)
                                              tGTTGATCTCACCGCTGGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCt                   <  1:1781619/68‑1 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTcg                                              >  1:423435/1‑37 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCg                                  >  1:2158272/1‑49 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCg                                  >  1:873513/1‑49 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGcc                                >  1:1630606/1‑51 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGt                          >  1:2839974/1‑57 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCa               >  1:2979062/1‑68 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCa               >  1:25926/1‑68 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCa               >  1:1234924/1‑68 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGGACCAGTTCg                                  >  1:2024209/1‑49 (MQ=255)
                                                  gATCTCACCGCTTGCCGGAAACGCC‑CCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAgg             >  1:1716538/1‑69 (MQ=255)
                                                  gATCTCACCGCTTGCCGCAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCa               >  1:169730/1‑68 (MQ=255)
                                                        aCCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGcac          >  1:2814904/1‑67 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACTCCAGtt   <  1:2487322/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:3260176/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:957831/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:934541/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:3123186/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:2717983/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:2450410/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:2287193/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:2240056/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:1559185/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:1211806/69‑1 (MQ=255)
                                                             ttGCCGGGAATGCCGCCATCCGGTCGCGCACCAGTTCGCCCAGCGTTTTTCCTTTCAGACACACCAGtt   <  1:1113034/69‑1 (MQ=255)
                                                              tGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTc  <  1:3352073/69‑1 (MQ=255)
                                                              tGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTc  <  1:191940/69‑1 (MQ=255)
                                                              tGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTc  <  1:1604668/69‑1 (MQ=255)
                                                              tGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTc  <  1:1351308/69‑1 (MQ=255)
                                                              tGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTc  <  1:1031247/69‑1 (MQ=255)
                                                                    |                                                              
TGCTGTTCCACGCGGTTAATCGCCTCAACGGGTTGCGCCAGTTTGCTGTTGATCTCACCGCTTGCCGGAAACGCCGCCATCCGGTCGCGTACCAGTTCGCCCAGCGTTTTATCTTTCAGGCACACCAGTTC  >  minE/1330399‑1330529

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