Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,503,335 C→T 18.9% S86S (TCC→TCT uhpT → hexose phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,503,3350CT18.9% 96.8 / 12.1 63S86S (TCC→TCTuhpThexose phosphate transporter
Reads supporting (aligned to +/- strand):  ref base C (35/16);  new base T (5/7);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ATCGGCCTGGGTTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTG  >  minE/2503285‑2503400
                                                  |                                                                 
atCGGCCTGGGTTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaa                                                 >  1:36433/1‑69 (MQ=255)
atCGGCCTGGGTTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaa                                                 >  1:1475302/1‑69 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:2965516/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:2566074/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:2635314/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:2213191/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:1791729/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:3190942/69‑1 (MQ=255)
atCGGCCTGGGCTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCTTACTACGCCGACGGCaaa                                                 <  1:1352500/69‑1 (MQ=255)
          gTTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:147030/1‑69 (MQ=255)
           tttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:1474364/1‑68 (MQ=255)
           tttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACa                                      >  1:170422/1‑69 (MQ=255)
           tttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACa                                      >  1:2994153/1‑69 (MQ=255)
           tttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACa                                      >  1:2051361/1‑69 (MQ=255)
           tttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACa                                      >  1:2169781/1‑69 (MQ=255)
            ttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACaa                                     <  1:573969/69‑1 (MQ=255)
            ttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACaa                                     <  1:2492326/69‑1 (MQ=255)
            ttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACaa                                     <  1:320647/69‑1 (MQ=255)
            ttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACaa                                     <  1:2192969/69‑1 (MQ=255)
            ttCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACaa                                     <  1:2002439/69‑1 (MQ=255)
               tcCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTc                                  >  1:3148975/1‑69 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaaaa                                               >  1:2739477/1‑50 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaaaa                                               >  1:593047/1‑50 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaaaa                                               >  1:2655966/1‑50 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCaaaaa                                               >  1:3005377/1‑50 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGcc                             >  1:631007/1‑68 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGcc                             >  1:3052161/1‑68 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGcc                             >  1:2578224/1‑68 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCg                            >  1:61609/1‑69 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCg                            >  1:864829/1‑69 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCg                            >  1:3054337/1‑69 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCg                            >  1:2631822/1‑69 (MQ=255)
                     aCTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCg                            >  1:1559957/1‑69 (MQ=255)
                         aTGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCa                        >  1:1402881/1‑69 (MQ=255)
                         aTGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCa                        >  1:206759/1‑69 (MQ=255)
                         aTGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCa                        >  1:16426/1‑69 (MQ=255)
                                 gttAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATc                >  1:2133929/3‑69 (MQ=255)
                                 ggTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATc                >  1:3420621/1‑69 (MQ=255)
                                 ggTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATc                >  1:184432/1‑69 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:2431637/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:1478733/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:1482545/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:1594428/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:2387517/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAAt                                    <  1:1092998/46‑1 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:1825736/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:742200/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:779059/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:2899145/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:2843689/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:516988/1‑62 (MQ=255)
                                    aaaaCGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCt                    >  1:2799907/1‑62 (MQ=255)
                                          cTGGTTTCTTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:472988/1‑37 (MQ=255)
                                          cTGGTTTCTTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:2615368/1‑37 (MQ=255)
                                          cTGGTTTCTTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:3158629/1‑37 (MQ=255)
                                          cTGGTTTCTTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:2921266/1‑37 (MQ=255)
                                          cTGGTTTCTTACTACGCCGACGGCAAAAACACCAAAc                                       >  1:191625/1‑37 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATg                      >  1:2797926/1‑49 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTg  <  1:3199656/69‑1 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTg  <  1:3079874/69‑1 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTg  <  1:2311137/69‑1 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTg  <  1:1590474/69‑1 (MQ=255)
                                               ttCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTg  <  1:1192727/69‑1 (MQ=255)
                                                  |                                                                 
ATCGGCCTGGGTTTCTCCATCACTTATGGCGTGGGTAAAACGCTGGTTTCCTACTACGCCGACGGCAAAAACACCAAACAATTCCTGCCGTTCATGCTGATCCTCTCTGCTATTTG  >  minE/2503285‑2503400

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: