Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 584,839 A→G 38.2% S352S (TCT→TCC cydC ← 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
*minE584,8390AG38.2% 42.4 / 59.0 89S352S (TCT→TCCcydCfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (36/19);  new base G (18/16);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAAC  >  minE/584779‑584898
                                                            |                                                           
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:2464840/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:482215/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:2036439/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:2716144/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:1539882/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:220931/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:1299601/69‑1 (MQ=255)
gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:324056/69‑1 (MQ=255)
       gcaatatTTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc                                                     <  1:3115626/59‑1 (MQ=255)
             tGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGtgaa                                         <  1:263834/68‑1 (MQ=255)
             tGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGtgaa                                         <  1:1208402/68‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:3163524/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:2097561/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:404235/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:724596/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:227732/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:518216/69‑1 (MQ=255)
                     ggCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTCCGGATAAGTGAACTGAACATc                                <  1:2039717/69‑1 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:2550959/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:2550794/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          <  1:61750/69‑1 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:895444/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:3309389/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          <  1:3265829/69‑1 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:941322/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          >  1:3124747/1‑69 (MQ=255)
                           tACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTaa                          <  1:1737159/69‑1 (MQ=255)
                                    agaAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGa                 <  1:3238394/69‑1 (MQ=255)
                                        aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAAc              >  1:493412/1‑68 (MQ=255)
                                        aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAAc              >  1:3199331/1‑68 (MQ=255)
                                        aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAAc              >  1:1857169/1‑68 (MQ=255)
                                                   tGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1796532/1‑69 (MQ=255)
                                                   tGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:222156/1‑69 (MQ=255)
                                                   tGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:3274249/1‑69 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:3409127/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:450213/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:2924009/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:2528717/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:3028661/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1108280/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:2392354/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1431432/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1989579/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1631836/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAac  >  1:1768354/1‑67 (MQ=255)
                                                     ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGAAACGTCAGCGAAACGCGATCGGCAac  >  1:3084928/1‑67 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:1071878/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:1220084/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:138500/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:1411340/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:1551348/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:408028/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2456968/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:1698418/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:586930/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:487865/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2043214/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:526145/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2696670/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2102815/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2212053/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCaa   >  1:2396272/1‑63 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATCGGCa    >  1:1954972/1‑62 (MQ=255)
                                                        gTTGGGATTGCTCCGGATAAGTGAACTGAACATCCCGTAATGTCAGCGAAACGCGATAGGCaa   >  1:3117574/1‑63 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:1013300/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:997635/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:351394/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:355390/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:714001/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:2069367/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:101955/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1505981/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1521637/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1610876/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1643479/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1741924/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:1841536/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:336184/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:2745525/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:2761244/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:2776107/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:3009131/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:3111685/1‑60 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:3127842/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:3140723/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:3145250/60‑1 (MQ=255)
                                                         ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:3307128/60‑1 (MQ=255)
                                                         ttGAGATTGCTACGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     <  1:808642/60‑1 (MQ=255)
                                                         ttGAGATTGCTACGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGc     >  1:826112/1‑60 (MQ=255)
                                                            |                                                           
GAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAAC  >  minE/584779‑584898

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: