Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,680,429 A→G 36.7% D1515D (GAT→GAC yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,680,4290AG36.7% 45.2 / 48.3 79D1515D (GAT→GACyfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/19);  new base G (8/21);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 4.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACCTGCAACCACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCG  >  minE/1680363‑1680493
                                                                  |                                                                
ttttACCTGCAACCATACCAGCACCAGATCTCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATc                                                                <  1:353064/69‑1 (MQ=255)
ttttACCTGCAACCACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATc                                                                <  1:157423/69‑1 (MQ=255)
ttttACCTGCAACCACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATc                                                                <  1:225942/69‑1 (MQ=255)
ttttACCTGCAACCACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATc                                                                <  1:2626522/69‑1 (MQ=255)
              acacCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAatat                                                  >  1:2230593/1‑69 (MQ=255)
              acacCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAatat                                                  >  1:1360146/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:869739/1‑57 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAatat                                                  >  1:2361200/1‑60 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:2647119/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:2080480/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:368050/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:1916803/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:2361780/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                         >  1:615176/1‑69 (MQ=255)
                       ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGAAGCTCGa                                         >  1:1291055/1‑69 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:2992885/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:961182/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:228796/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:232692/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:1067007/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:786924/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:1332502/1‑44 (MQ=255)
                                    aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                     >  1:85235/1‑44 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:3257791/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2977438/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:3168181/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2354138/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:3006385/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2634297/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2299635/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:209842/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1168551/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1402823/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:761137/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1475999/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:701273/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:630421/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1544247/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1606064/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:431392/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2041503/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:2191917/69‑1 (MQ=255)
                                        aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCACCACATTGCTc                        <  1:2069956/69‑1 (MQ=255)
                                              tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTc                        >  1:1783113/1‑63 (MQ=255)
                                              tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCgg                   >  1:3340585/1‑68 (MQ=255)
                                              tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCgg                   >  1:193906/1‑68 (MQ=255)
                                              tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCgg                   >  1:2632810/1‑68 (MQ=255)
                                                cAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTaa                <  1:2477789/69‑1 (MQ=255)
                                                   cGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                        <  1:1829940/58‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:1119687/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:1493398/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:3416577/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:2419497/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:1417996/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:2445363/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:3254433/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:1305825/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:266390/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:2978992/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:2801221/69‑1 (MQ=255)
                                                       ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:126743/69‑1 (MQ=255)
                                                        tGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg         <  1:1487190/68‑1 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAgg             >  1:786730/1‑58 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:861807/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:541742/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:2401994/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:518988/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:2021227/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:2580015/1‑69 (MQ=255)
                                                              tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg  >  1:3134919/1‑69 (MQ=255)
                                                               aCCATCAGTACCAATAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAgg             <  1:3285114/57‑1 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGCCGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:1790379/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGc                      >  1:2876206/1‑48 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:478502/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:2842619/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:3429185/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:2628371/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg  >  1:2956425/1‑68 (MQ=255)
                                                               aCCATCAGTACCAAGAAAATGACGCTCGATTTGCAGCACATAGTTCGCCGGTAAAGGtgcggattgcg  >  1:1343807/1‑68 (MQ=255)
                                                                  |                                                                
TTTTACCTGCAACCACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCG  >  minE/1680363‑1680493

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: