Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,619,848 T→C 23.6% A297A (GCA→GCG 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,8480TC23.6% 94.1 / 23.9 76A297A (GCA→GCGzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (5/13);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

AGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCA  >  minE/2619786‑2619915
                                                              |                                                                   
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgggc                                                               >  1:1255469/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:1701686/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:1944671/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:1936985/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:40934/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:350950/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:1250536/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:2911044/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:1443920/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:2565282/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcgc                                                               >  1:3356719/1‑69 (MQ=255)
aGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTgcgcg                                                                >  1:743321/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCGCGTTGCGCGCTCCAc                                                          >  1:1790555/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:3348802/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:898804/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:1750718/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:148148/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:3417435/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:854604/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:291204/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCAc                                                          >  1:2578504/1‑68 (MQ=255)
      tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTCGCGCTCCAcc                                                         >  1:3214203/1‑68 (MQ=255)
         aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:1528751/68‑1 (MQ=255)
         aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:676952/68‑1 (MQ=255)
         aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:548657/68‑1 (MQ=255)
              cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1992614/69‑1 (MQ=255)
              cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1858112/69‑1 (MQ=255)
              cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1076221/69‑1 (MQ=255)
               aCCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCAGTTGCGCGCTCCACCgg                                                       <  1:1940958/62‑1 (MQ=255)
                  aGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:824394/59‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:2435162/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:1068257/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:689045/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:2913402/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:319047/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:3364528/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:501203/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:1574109/69‑1 (MQ=255)
                        tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCGCCACCGGAATGGATTCGCCGGTc                                       <  1:3148425/69‑1 (MQ=255)
                         cTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:1167199/68‑1 (MQ=255)
                                         agcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      <  1:1534129/69‑1 (MQ=255)
                                          gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:2040703/1‑62 (MQ=255)
                                          gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:1987590/1‑62 (MQ=255)
                                          gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:2844619/1‑62 (MQ=255)
                                          gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:972884/1‑62 (MQ=255)
                                          gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      <  1:2707471/68‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2370303/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:553/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:44217/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:1127068/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:1423827/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:1497406/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:3240369/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:1747799/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2722409/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2156101/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:797352/53‑1 (MQ=255)
                                              gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2428265/53‑1 (MQ=255)
                                               aaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:623287/53‑1 (MQ=255)
                                                   ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:171027/1‑59 (MQ=255)
                                                   ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:1971098/1‑59 (MQ=255)
                                                   ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:1117767/1‑59 (MQ=255)
                                                        gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:193431/1‑67 (MQ=255)
                                                        gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:2901632/1‑67 (MQ=255)
                                                        gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:1485231/1‑67 (MQ=255)
                                                        gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:2567298/1‑67 (MQ=255)
                                                        gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:378190/1‑67 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:88372/69‑1 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:2562924/69‑1 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:430527/69‑1 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:2902593/69‑1 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:2616733/69‑1 (MQ=255)
                                                           cGTTGCGCGCTCCACCGGAATGGATTCGCAGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:1515002/69‑1 (MQ=255)
                                                            gTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCGAACGGTGagag    <  1:995350/68‑1 (MQ=255)
                                                            gTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:500317/68‑1 (MQ=255)
                                                            gTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:1551781/68‑1 (MQ=255)
                                                             ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:1061026/69‑1 (MQ=255)
                                                             ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:2530106/69‑1 (MQ=255)
                                                             ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:1654969/69‑1 (MQ=255)
                                                             ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:1653409/69‑1 (MQ=255)
                                                             ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:1512228/69‑1 (MQ=255)
                                                              |                                                                   
AGCACTTCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCA  >  minE/2619786‑2619915

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: