Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,895,677 G→A 29.9% A237A (GCG→GCA yjfF → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,895,6770GA29.9% 58.6 / 24.5 60A237A (GCG→GCAyjfFpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (5/37);  new base A (8/10);  total (13/47)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCT  >  minE/2895610‑2895736
                                                                   |                                                           
tGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATgcg                                                             >  1:2334921/1‑68 (MQ=255)
  ctctCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTg                                                          >  1:1201823/1‑69 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:298444/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:2175298/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:392964/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:1680916/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:285180/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:1531941/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcggg                                                     >  1:1142218/1‑67 (MQ=255)
         ccGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTggcgg                                                      >  1:2895436/1‑66 (MQ=255)
             aCTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAgg                                               <  1:1232680/69‑1 (MQ=255)
                 gCAATGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGgtgt                                            >  1:3124268/1‑68 (MQ=255)
                 gCAACGCTGGCTGGGATTGTCTTCTCGATTTATACACAGGCCGGATATgcgc                                                            >  1:781624/1‑52 (MQ=255)
                 gCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGgtgt                                            >  1:2288721/1‑68 (MQ=255)
                 gCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGgtgt                                            >  1:734933/1‑68 (MQ=255)
                 gCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGgtgt                                            >  1:2594282/1‑68 (MQ=255)
                   aaCGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGa                                         <  1:2861024/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:2460705/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:1232739/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:2044959/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:579794/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:1659966/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:635469/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:1541379/69‑1 (MQ=255)
                       cTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTg                                     <  1:1450392/69‑1 (MQ=255)
                             gggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:1038713/69‑1 (MQ=255)
                             gggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:3142446/69‑1 (MQ=255)
                             gggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:403894/69‑1 (MQ=255)
                             gggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:1569678/69‑1 (MQ=255)
                             gggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:847379/69‑1 (MQ=255)
                              ggATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:3081023/68‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:47707/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:1554451/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:286971/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:3287232/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:355991/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:2843585/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:1068374/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:1114820/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:2486681/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:2424244/69‑1 (MQ=255)
                               gATTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:1642826/69‑1 (MQ=255)
                                aTTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:2940305/68‑1 (MQ=255)
                                aTTGTCTTCTCGATTTATACCCAGGCCGGATATGCACTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:323173/68‑1 (MQ=255)
                                   gTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAgg                                               <  1:2553923/47‑1 (MQ=255)
                                   gTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAgg                                               <  1:1881509/47‑1 (MQ=255)
                                             tttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTa                <  1:1666624/68‑1 (MQ=255)
                                             tttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTa                <  1:1230612/68‑1 (MQ=255)
                                             tttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTa                <  1:2356214/68‑1 (MQ=255)
                                             tttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTa                <  1:214442/68‑1 (MQ=255)
                                             tttATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTa                <  1:632669/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:688614/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:932185/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:2913298/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:3222114/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:3212880/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:2895991/68‑1 (MQ=255)
                                                     cAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCgg        <  1:2923530/68‑1 (MQ=255)
                                                          cGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCt  <  1:426572/69‑1 (MQ=255)
                                                          cGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCt  <  1:926733/69‑1 (MQ=255)
                                                           ggATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTAt                             <  1:1027745/41‑1 (MQ=255)
                                                               atGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCt                               <  1:2044579/35‑1 (MQ=255)
                                                                   |                                                           
TGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTTATACCCAGGCCGGATATGCGCTGGCGGGCGTAGGTGTGGAACTGGATGCTATCGCCTCAGTGGTAATTGGCGGTACGCT  >  minE/2895610‑2895736

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: