Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 679,872 C→T 55.6% T237T (ACG→ACA yccS ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE679,8720CT55.6% 9.2 / 49.4 63T237T (ACG→ACAyccSpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (20/8);  new base T (17/18);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 7.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGT  >  minE/679805‑679934
                                                                   |                                                              
tttGATACTGAATATGAGAAGAGCTGGCACGCTCGTTAATATCCTGTGCGACAAAGTAATAATGCAGCg                                                               >  1:1412909/1‑69 (MQ=255)
tttGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCg                                                               >  1:366786/1‑69 (MQ=255)
    aTACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGc                                                            <  1:791502/68‑1 (MQ=255)
    aTACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGc                                                            <  1:2037241/68‑1 (MQ=255)
    aTACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGc                                                            <  1:168687/68‑1 (MQ=255)
    aTACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGc                                                            <  1:1773401/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2953000/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:3433005/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:338587/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2218590/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:3112861/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2291669/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:8887/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2322992/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2560421/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2606599/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2626916/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:2798362/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:850877/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:1651254/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:633611/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:1761831/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:1809576/68‑1 (MQ=255)
                                    gAATATCCTGTGCGACAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGt                            <  1:675762/68‑1 (MQ=255)
                                         tCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCgg                      >  1:385435/1‑69 (MQ=255)
                                         tCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCgg                      >  1:265008/1‑69 (MQ=255)
                                         tCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCgg                      >  1:1738057/1‑69 (MQ=255)
                                         tCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCgg                      >  1:2014042/1‑69 (MQ=255)
                                         tCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCgg                      >  1:575254/1‑69 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:297995/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:3261784/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:750789/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:492838/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:420257/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:941154/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:1473273/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2124592/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2120896/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2129038/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2138954/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:1869079/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:1773355/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2141504/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2172457/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2822891/1‑66 (MQ=255)
                                                   cAAAGTAATAATGCAGTGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTcagca               >  1:2293233/1‑66 (MQ=255)
                                                     aaGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCgag          <  1:142651/69‑1 (MQ=255)
                                                     aaGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCgag          <  1:2909367/69‑1 (MQ=255)
                                                     aaGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCgag          <  1:2127996/69‑1 (MQ=255)
                                                     aaGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCgag          <  1:2892984/69‑1 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATc                                   >  1:2260343/1‑42 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCg                       >  1:2857847/1‑54 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:1101822/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:1426164/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:1243999/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:36418/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:362421/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:2250185/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:3253464/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:2292678/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGagag        >  1:2954949/1‑69 (MQ=255)
                                                       gTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACCTAAGCGGGTCAGCAGCGagag        >  1:46371/1‑69 (MQ=255)
                                                             atGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGt  >  1:3372933/1‑69 (MQ=255)
                                                                   |                                                              
TTTGATACTGAATATGAGAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGT  >  minE/679805‑679934

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: