Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 586,692 G→A 59.2% L324L (CTG→TTG)  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,6920GA59.2% 25.7 / 49.9 71L324L (CTG→TTG) cydDfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (18/24);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGC  >  minE/586630‑586756
                                                              |                                                                
cGATGCTAATTCCGCCTCACCACGTTGCGGATGTGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                            <  1:470628/69‑1 (MQ=255)
cGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                            <  1:2168189/69‑1 (MQ=255)
cGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                            <  1:1553081/69‑1 (MQ=255)
    gCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGc                                                         >  1:361501/1‑68 (MQ=255)
    gCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCt                                                        >  1:1066927/1‑69 (MQ=255)
           ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                   >  1:981469/1‑67 (MQ=255)
           ccGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGcc                                                   >  1:2949387/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)
                tCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTc                                                            <  1:2604136/53‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGAt             >  1:1492334/1‑64 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:2854522/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:1071013/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:891135/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:577131/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:521668/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:1123559/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:346822/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:3319821/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:3231661/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:3118006/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:30083/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:298850/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:1127607/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:2608468/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        >  1:2050419/1‑69 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:2014008/69‑1 (MQ=255)
                                                    aaCGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGt        <  1:1647431/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCTAGATCGCGTAATg    <  1:3191617/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:331982/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2958824/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:1786315/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:1294858/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:3383940/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:1674325/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:1081647/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:1528427/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:641896/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:681730/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:842066/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:121523/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2553191/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2488531/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2367632/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2685717/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2222194/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2754871/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:1416464/69‑1 (MQ=255)
                                                        ttttCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    <  1:2940688/69‑1 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATgg   >  1:2214079/1‑67 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATg    >  1:2395320/1‑66 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1514330/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1429158/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:970383/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1384445/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:2483380/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1342572/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:256210/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1744073/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1192614/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:3308650/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:3202086/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:3150955/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:1886626/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:2751521/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:2814456/1‑68 (MQ=255)
                                                           tCAAACTGTCAGCTGCACCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGc  >  1:2908742/1‑68 (MQ=255)
                                                              |                                                                
CGATGCTAATTCCGCCTCACCACGTTGCGGATGGGCGAGCGGGGTTTCCATAAACGTTTTCAGACTGTCAGCTGCGCCAACAGCCTGGGCTTTAGCATGATAAAACGTACCGAGATCGCGTAATGGC  >  minE/586630‑586756

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: