Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,619,860 C→A 30.5% P293P (CCG→CCT 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,8600CA30.5% 73.2 / 33.6 59P293P (CCG→CCTzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base C (7/34);  new base A (5/13);  total (12/47)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

TCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCA  >  minE/2619792‑2619915
                                                                    |                                                       
tCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTCGCGCTCCAcc                                                         >  1:3214203/1‑68 (MQ=255)
   aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:548657/68‑1 (MQ=255)
   aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:676952/68‑1 (MQ=255)
   aaCGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:1528751/68‑1 (MQ=255)
        cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1076221/69‑1 (MQ=255)
        cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1992614/69‑1 (MQ=255)
        cACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGa                                                 <  1:1858112/69‑1 (MQ=255)
         aCCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCAGTTGCGCGCTCCACCgg                                                       <  1:1940958/62‑1 (MQ=255)
            aGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCgg                                                       <  1:824394/59‑1 (MQ=255)
                  tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:1574109/69‑1 (MQ=255)
                  tCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTc                                       <  1:2435162/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:1068257/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:972884/1‑62 (MQ=255)
                                    gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:2844619/1‑62 (MQ=255)
                                    gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:1987590/1‑62 (MQ=255)
                                    gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTtc                            >  1:2040703/1‑62 (MQ=255)
                                    gcagGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      <  1:2707471/68‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:553/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:3240369/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:44217/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2722409/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2370303/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:797352/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:1747799/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2428265/53‑1 (MQ=255)
                                        gaaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:2156101/53‑1 (MQ=255)
                                         aaCTTTATCGCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACt                               <  1:623287/53‑1 (MQ=255)
                                             ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:1971098/1‑59 (MQ=255)
                                             ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:1117767/1‑59 (MQ=255)
                                             ttATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCaaa                      >  1:171027/1‑59 (MQ=255)
                                                  gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:193431/1‑67 (MQ=255)
                                                  gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:1485231/1‑67 (MQ=255)
                                                  gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:378190/1‑67 (MQ=255)
                                                  gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:2567298/1‑67 (MQ=255)
                                                  gCCCGTCGCGCGCTCCACAGGAATGGATTCTCCGGTCAGGGCACTTTCATCAAAACTGGCAAACGGt         >  1:2901632/1‑67 (MQ=255)
                                                     cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:88372/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)
                                                     cGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGagag    <  1:2902593/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: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)
                                                       ttGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCa  <  1:1061026/69‑1 (MQ=255)
                                                                   ccGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACg           >  1:422820/1‑48 (MQ=255)
                                                                    |                                                       
TCCAACGTCACCAGACGGTCTACGCTGGTGGCACCAGCAGGGACTTTATCGCCCGTTGCGCGCTCCACCGGAATGGATTCGCCGGTCAGGGCGCTTTCATCAAAACTGGCAAACGGTGAGAGCA  >  minE/2619792‑2619915

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: