Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,638,464 T→G 60.7% D152E (GAT→GAG ugpC → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,638,4640TG60.7% 12.3 / 20.1 28D152E (GAT→GAGugpCglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base G (7/10);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

CGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCA  >  minE/2638403‑2638523
                                                             |                                                           
cgcgcgAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCcggcgg                                                      <  1:2291219/69‑1 (MQ=255)
 gcgcgAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGAGCGCGATCcggcgg                                                      <  1:3152332/68‑1 (MQ=255)
                          cAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGctctct                            <  1:2921008/69‑1 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGctctct                            >  1:2012827/1‑68 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:72712/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:693940/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:3314261/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:3095589/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:2159279/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGGCGTGCGATTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:2815308/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTtg                                          >  1:1424517/1‑54 (MQ=255)
                           aGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:2756186/1‑69 (MQ=255)
                           aGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTa                           >  1:192388/1‑69 (MQ=255)
                               tgtgGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAAcc                        >  1:2319671/1‑68 (MQ=255)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:1299155/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:2770717/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:2381910/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:2861343/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:2980028/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:1993222/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:1656422/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:1636444/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:952664/68‑1 (MQ=25)
                                                   ttGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt     <  1:1945034/67‑1 (MQ=25)
                                                    tgtgCGCGATCCGGCTGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCa  <  1:2621889/69‑1 (MQ=255)
                                                    tgtgCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCa  <  1:2701989/69‑1 (MQ=255)
                                                    tgtgCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCa  <  1:3274769/69‑1 (MQ=255)
                                                    tgtgCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCa  <  1:383349/69‑1 (MQ=255)
                                                             |                                                           
CGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCA  >  minE/2638403‑2638523

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: