Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,619,884 G→A 61.8% S285S (AGC→AGT 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,8840GA61.8% 40.1 / 52.6 81S285S (AGC→AGTzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base A (20/30);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCT  >  minE/2619827‑2619944
                                                         |                                                            
agcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                                   <  1:1534129/69‑1 (MQ=255)
 gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                                         >  1:972884/1‑62 (MQ=255)
 gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                                         >  1:2040703/1‑62 (MQ=255)
 gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                                         >  1:1987590/1‑62 (MQ=255)
 gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                                                         >  1:2844619/1‑62 (MQ=255)
 gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                                   <  1:2707471/68‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:2722409/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:3240369/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:2156101/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:1747799/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:2370303/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:44217/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:2428265/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:1497406/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:797352/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:1423827/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:553/53‑1 (MQ=255)
     gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:1127068/53‑1 (MQ=255)
      aaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                                                            <  1:623287/53‑1 (MQ=255)
          ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                                   >  1:1117767/1‑59 (MQ=255)
          ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                                   >  1:1971098/1‑59 (MQ=255)
          ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                                                   >  1:171027/1‑59 (MQ=255)
               gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                                      >  1:2567298/1‑67 (MQ=255)
               gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                                      >  1:193431/1‑67 (MQ=255)
               gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                                      >  1:378190/1‑67 (MQ=255)
               gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                                      >  1:2901632/1‑67 (MQ=255)
               gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt                                      >  1:1485231/1‑67 (MQ=255)
                  cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:2902593/69‑1 (MQ=255)
                  cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:430527/69‑1 (MQ=255)
                  cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:2616733/69‑1 (MQ=255)
                  cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:2562924/69‑1 (MQ=255)
                  cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:88372/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:1654969/69‑1 (MQ=255)
                    ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                               <  1:1061026/69‑1 (MQ=255)
                    ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                               <  1:1653409/69‑1 (MQ=255)
                    ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                               <  1:2530106/69‑1 (MQ=255)
                    ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa                               <  1:1512228/69‑1 (MQ=255)
                                ccGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACg                                        >  1:422820/1‑48 (MQ=255)
                                  ggAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggca                 >  1:203154/1‑69 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:3362658/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:320614/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:430511/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:2793601/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:614518/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:709151/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:2385166/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:2387657/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:1182814/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:213215/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:2400269/69‑1 (MQ=255)
                                      tGGATTCTCCGGTCAGGGCACTTTCATCAAAACTAGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:760532/69‑1 (MQ=255)
                                       ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:75477/68‑1 (MQ=255)
                                       ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:565594/68‑1 (MQ=255)
                                       ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:1634369/68‑1 (MQ=255)
                                       ggATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:1079775/68‑1 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGccc      >  1:229192/1‑67 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:1062995/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:135956/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:1635375/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:1080289/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:38383/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:3267233/1‑69 (MQ=255)
                                               cGGTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCAc    >  1:1069988/1‑69 (MQ=255)
                                               cGGTCAGGGCACTATCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCggcaggca             <  1:794802/60‑1 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAAGTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:2209190/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:1757069/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:2045856/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:49466/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:3245816/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:211045/1‑69 (MQ=255)
                                                 gTCAGGGCACTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCt  >  1:1095312/1‑69 (MQ=255)
                                                   cAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag                                 <  1:1168963/36‑1 (MQ=255)
                                                   cAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGttt                           <  1:1661538/42‑1 (MQ=255)
                                                   cAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGttt                           <  1:1760319/42‑1 (MQ=255)
                                                   cAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGttt                           <  1:1862803/42‑1 (MQ=255)
                                                    aGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAAc          >  1:3191833/1‑58 (MQ=255)
                                                    aGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAAc          >  1:2829089/1‑58 (MQ=255)
                                                         |                                                            
AGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCAGTTTACCGTCGGCAGGCAAACGCCCACCT  >  minE/2619827‑2619944

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: