Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,638,470 G→A 67.8% A154A (GCG→GCA ugpC → glycerol‑3‑phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,638,4700GA67.8% 14.8 / 9.0 26A154A (GCG→GCAugpCglycerol‑3‑phosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (7/10);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTGA  >  minE/2638403‑2638535
                                                                   |                                                                 
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:2815308/1‑69 (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)
                           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:2980028/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:2861343/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:1299155/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:952664/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:2770717/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:1636444/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:1656422/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:2381910/68‑1 (MQ=25)
                                                  aTTGTGCGTGAGCCGGCAGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:1993222/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:383349/69‑1 (MQ=255)
                                                    tgtgCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCa              <  1:3274769/69‑1 (MQ=255)
                                                               cggcggTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTg   <  1:2437847/69‑1 (MQ=255)
                                                               cggcggTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTg   <  1:2919821/69‑1 (MQ=255)
                                                               cggcggTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTg   <  1:1389545/69‑1 (MQ=255)
                                                               cggcagTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGTTGCgtgt                 <  1:3282093/55‑1 (MQ=25)
                                                                 gcggTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTGa  <  1:1827238/68‑1 (MQ=255)
                                                                 gcggTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTGa  <  1:2907542/68‑1 (MQ=255)
                                                                 gcggTGGTCCTGTTTGATGAGCCGCGCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTGa  <  1:2889836/68‑1 (MQ=255)
                                                                   |                                                                 
CGCGCGAGCTTTCCGGCGGTCAGCGCCAGCGTGTGGCGATGGGCCGCGCGATTGTGCGCGATCCGGCGGTGTTCCTGTTTGATGAGCCGCTCTCTAACCTCGATGCCAAGCTGCGCGTGCAGATGCGTCTTGA  >  minE/2638403‑2638535

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: