Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,620,337 T→A 35.7% L134L (CTA→CTT zntA ← zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,620,3370TA35.7% 83.9 / 77.1 95L134L (CTA→CTTzntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base T (47/14);  new base A (19/15);  total (66/29)
Fisher's exact test for biased strand distribution p-value = 3.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AACGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTGTTCA  >  minE/2620272‑2620400
                                                                 |                                                               
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGCCCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:3355206/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTa                                                                >  1:3287456/1‑67 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:2093967/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:3168259/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:3089738/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:3000377/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:3301573/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:341299/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:224054/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:2058161/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:528100/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:2014849/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:598314/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:1587299/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:1269234/1‑68 (MQ=255)
aaCGCCAGTTGCCCGAACGGATGATTGAACTGCTCAAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:733713/1‑68 (MQ=255)
aaCGCCAGTTGCCCCAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAg                                                               >  1:2252322/1‑68 (MQ=255)
  cGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGt                                                            <  1:2081111/69‑1 (MQ=255)
   gCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTa                                                           <  1:2414297/69‑1 (MQ=255)
   gCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTa                                                           <  1:2368494/69‑1 (MQ=255)
   gCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTa                                                           <  1:1004543/69‑1 (MQ=255)
   gCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTa                                                           <  1:3005268/69‑1 (MQ=255)
   gCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTa                                                           <  1:830944/69‑1 (MQ=255)
          gCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGt                                                            >  1:2226627/1‑61 (MQ=255)
          gCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGt                                                            >  1:2824499/1‑61 (MQ=255)
          gCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGt                                                            >  1:1842910/1‑61 (MQ=255)
          gCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCAACACGAtt                                                                 >  1:1993048/1‑56 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:3311995/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:2669697/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:2957655/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:1668761/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:174241/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:2707514/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:272143/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:2234014/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:1576992/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:1175781/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:2791877/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:1576279/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:1392888/57‑1 (MQ=255)
              aaaCGGATGATTGAACTGCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:908439/57‑1 (MQ=255)
              aaaCGGATGATTGAACTCCTCCAGCCCCCAACTGATTGCCATCATCACGATAAGCGTa                                                           <  1:3266032/57‑1 (MQ=37)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCa                                                       >  1:442350/1‑59 (MQ=255)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGAtt                                             >  1:843288/1‑69 (MQ=255)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGAtt                                             >  1:2081798/1‑69 (MQ=255)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGAtt                                             >  1:319026/1‑69 (MQ=255)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGAt                                              >  1:925739/1‑68 (MQ=255)
                 cGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGAt                                              >  1:3381202/1‑68 (MQ=255)
                     tgatTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGAGTAATCAGCGGCAGATtctc                                          >  1:3355191/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCtgcggc         <  1:702755/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1244044/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1042423/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:2800136/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1271170/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1325681/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1873377/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCggcggc         <  1:1624022/65‑1 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCt       >  1:3400343/1‑67 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:3432237/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:429859/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:510822/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:998817/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:761222/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:3065266/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:139569/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:1749609/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:2320557/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:2323086/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:2449545/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:2853280/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:3114968/1‑68 (MQ=255)
                                                         atcaCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTg      >  1:301191/1‑68 (MQ=255)
                                                            aCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATg                            >  1:1394485/1‑43 (MQ=255)
                                                            aCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTGTTCa  >  1:2735764/1‑69 (MQ=255)
                                                            aCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTGTTCa  >  1:1869448/1‑69 (MQ=255)
                                                            aCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTGTTCa  >  1:588459/1‑69 (MQ=255)
                                                            aCGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCATGTTCa  >  1:2916101/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTgtg                      >  1:1246175/1‑47 (MQ=38)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGttcg                 >  1:2461339/1‑53 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:711056/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:1110103/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:332859/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:890658/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:870076/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:1495998/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:3238426/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:1781794/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:371857/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:2481014/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:1970568/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:2033795/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:515081/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:2207966/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:2948386/1‑69 (MQ=255)
                                                            aCGATAAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTAggaggc         >  1:2705350/1‑62 (MQ=255)
                                                            aCGATAAGCGTAATCAGAGGCAGATTCTCTTTCAGGCGTGATTCTTGTGGTTCGTCGGAGGCCTGTTCa  >  1:3010782/1‑69 (MQ=255)
                                                                 |                                                               
AACGCCAGTTGCCCGAACGGATGATTGAACTGCTCCAGACCCCAGCTGATTGCCATCATCACGATTAGCGTAATCAGCGGCAGATTCTCTTTCAGGCGTGATGCTTGCGGTTCTTCGGCGGCCTGTTCA  >  minE/2620272‑2620400

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: