Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,619,872 G→T 49.2% G289G (GGC→GGA 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,8720GT49.2% 10.1 / 74.3 69G289G (GGC→GGAzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (9/26);  new base T (5/29);  total (14/55)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCA  >  minE/2619810‑2619933
                                                              |                                                             
tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                                         <  1:1068257/69‑1 (MQ=255)
tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                                         <  1:2913402/69‑1 (MQ=255)
tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                                         <  1:2435162/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:689045/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:2844619/1‑62 (MQ=255)
                  gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                              >  1:2040703/1‑62 (MQ=255)
                  gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                              >  1:1987590/1‑62 (MQ=255)
                  gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                              >  1:972884/1‑62 (MQ=255)
                  gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                        <  1:2707471/68‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:553/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:2722409/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:1127068/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:2428265/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:3240369/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:1423827/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:2370303/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:2156101/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:1497406/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:797352/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:1747799/53‑1 (MQ=255)
                      gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:44217/53‑1 (MQ=255)
                       aaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                 <  1:623287/53‑1 (MQ=255)
                           ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                        >  1:1971098/1‑59 (MQ=255)
                           ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                        >  1:171027/1‑59 (MQ=255)
                           ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                        >  1:1117767/1‑59 (MQ=255)
                                gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                           >  1:378190/1‑67 (MQ=255)
                                gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                           >  1:2901632/1‑67 (MQ=255)
                                gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                           >  1:2567298/1‑67 (MQ=255)
                                gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                           >  1:193431/1‑67 (MQ=255)
                                gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                           >  1:1485231/1‑67 (MQ=255)
                                   cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                      <  1:88372/69‑1 (MQ=255)
                                   cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                      <  1:2902593/69‑1 (MQ=255)
                                   cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                      <  1:430527/69‑1 (MQ=255)
                                   cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                      <  1:2562924/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:1551781/68‑1 (MQ=255)
                                    gTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                      <  1:500317/68‑1 (MQ=255)
                                     ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                    <  1:2530106/69‑1 (MQ=255)
                                     ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                    <  1:1512228/69‑1 (MQ=255)
                                     ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                    <  1:1061026/69‑1 (MQ=255)
                                     ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                    <  1:1654969/69‑1 (MQ=255)
                                     ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                    <  1:1653409/69‑1 (MQ=255)
                                                 ccGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACg                             >  1:422820/1‑48 (MQ=255)
                                                   ggAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggca      >  1:203154/1‑69 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:2793601/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:1182814/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:709151/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:614518/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:430511/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:213215/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:2385166/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:3362658/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:2387657/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:320614/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:2400269/69‑1 (MQ=255)
                                                       tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTAGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:760532/69‑1 (MQ=255)
                                                        ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:1634369/68‑1 (MQ=255)
                                                        ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:565594/68‑1 (MQ=255)
                                                        ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:75477/68‑1 (MQ=255)
                                                        ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca  <  1:1079775/68‑1 (MQ=255)
                                                              |                                                             
TCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCA  >  minE/2619810‑2619933

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: