Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,619,653 A→G 50.7% A362A (GCT→GCC zntA ← zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,619,6530AG50.7% ‑1.8 / 65.2 73A362A (GCT→GCCzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (31/6);  total (49/24)
Fisher's exact test for biased strand distribution p-value = 2.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AATCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATAAACCGCTCAAT  >  minE/2619596‑2619712
                                                         |                                                           
aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:1210718/69‑1 (MQ=255)
aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:533125/69‑1 (MQ=255)
aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:2764050/69‑1 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:948058/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2915643/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2088101/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2815065/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2228865/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:458685/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:15577/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:154882/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:150843/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:807664/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:837761/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:266928/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2641199/1‑69 (MQ=255)
        ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2509760/1‑69 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:846274/1‑69 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:19884/1‑69 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:1174669/1‑69 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:940490/1‑69 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2211539/69‑1 (MQ=255)
                 tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2207471/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:2815054/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:3213286/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:3202302/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:1607366/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:592156/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:135403/69‑1 (MQ=255)
                                      gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:1277104/69‑1 (MQ=255)
                                       cACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:209127/68‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:3069154/69‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:1094997/69‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:913888/69‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:2481791/69‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:2477753/69‑1 (MQ=255)
                                          cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:1403107/69‑1 (MQ=255)
                                           aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:153234/43‑1 (MQ=255)
                                           aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2458816/43‑1 (MQ=255)
                                           aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:1995341/43‑1 (MQ=255)
                                           aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:1396610/43‑1 (MQ=255)
                                           aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:3066179/43‑1 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2478615/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:327510/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2181769/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:1274166/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2321372/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2345198/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2763744/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2358850/3‑67 (MQ=255)
                                               ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:1068341/3‑67 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCaa   >  1:2386833/2‑68 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1324010/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1379873/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:814914/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1142794/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:71371/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:958432/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:578750/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1747421/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:534268/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:468730/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1559124/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:378010/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:3102736/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2985705/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:22133/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2904186/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2225678/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2562213/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:239925/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:240914/2‑69 (MQ=255)
                                                taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGATCAAt  >  1:981530/2‑69 (MQ=255)
                                                         |                                                           
AATCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATAAACCGCTCAAT  >  minE/2619596‑2619712

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: