Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 586,705 G→A 65.9% G319G (GGC→GGT cydD ← fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE586,7050GA65.9% 53.4 / 44.0 79G319G (GGC→GGTcydDfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (19/8);  new base A (28/24);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

CCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCTGGAAA  >  minE/586641‑586762
                                                                |                                                         
ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                         >  1:2949387/1‑67 (MQ=255)
ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                         >  1:981469/1‑67 (MQ=255)
ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                         >  1:2230740/1‑67 (MQ=255)
ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                         >  1:1357107/1‑67 (MQ=255)
ccGCCTCACCACGTTGGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAAc                                                      >  1:2304959/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGAt                   >  1:1492334/1‑64 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:1647431/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:2608468/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:2854522/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:298850/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:2050419/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:2014008/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:30083/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:3118006/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:3231661/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:1071013/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:346822/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:1123559/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:521668/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:577131/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              <  1:1127607/69‑1 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:3319821/1‑69 (MQ=255)
                                         aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt              >  1:891135/1‑69 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCTAGATCGCGTAATg          <  1:3191617/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:641896/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2685717/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2754871/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2940688/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:842066/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2958824/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:681730/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2990887/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:3000196/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:30899/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1081647/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:3383940/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:331982/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2553191/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1674325/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1528427/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1786315/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1416464/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2222194/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2367632/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:1294858/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:121523/69‑1 (MQ=255)
                                             ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          <  1:2488531/69‑1 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATgg         >  1:2214079/1‑67 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg          >  1:2395320/1‑66 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1342572/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1384445/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1192614/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1429158/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1514330/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:970383/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:3308650/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:3150955/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:3202086/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:256210/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:2751521/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:2814456/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1744073/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:2483380/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:2908742/1‑68 (MQ=255)
                                                tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:1886626/1‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc        >  1:2616810/4‑67 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:3159692/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:923529/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:2487988/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:385224/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:3009865/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:3360979/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:3116177/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:1642233/4‑68 (MQ=255)
                                                 caaaCTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCt       >  1:1653295/4‑68 (MQ=255)
                                                     cTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCTGGaaa  >  1:1600022/1‑69 (MQ=255)
                                                     cTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCTGGaaa  >  1:3144158/1‑69 (MQ=255)
                                                     cTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCTGGaaa  >  1:1181419/1‑69 (MQ=255)
                                                     cTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACAGAGATCGCGTAATGGCTg      >  1:2043667/1‑65 (MQ=255)
                                                                |                                                         
CCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGCTGGAAA  >  minE/586641‑586762

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: