Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 668,428 C→A 24.9% T589K (ACG→AAG)  uup → fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE668,4280CA24.9% 82.9 / 10.5 32T589K (ACG→AAG) uupfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (18/6);  new base A (4/4);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATATGGCT  >  minE/668361‑668498
                                                                   |                                                                      
cAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGc                                                                      <  1:1388872/70‑1 (MQ=255)
            gAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGc                                                         >  1:118231/1‑71 (MQ=255)
                agcTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGcttcc                                                     >  1:409157/1‑71 (MQ=255)
                      cGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCttcttc                                               <  1:1685028/71‑1 (MQ=255)
                        cAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTtcag                                             <  1:1427468/71‑1 (MQ=255)
                         aaTTGCTCGAAGATCTGGAGGCTAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagt                                            <  1:1864924/71‑1 (MQ=255)
                         aaTTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagt                                            <  1:865083/71‑1 (MQ=255)
                         aaTTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagt                                            <  1:1033926/71‑1 (MQ=255)
                         aaTTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTtcag                                             >  1:944954/1‑70 (MQ=255)
                          aTTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagt                                            >  1:378317/1‑70 (MQ=255)
                          aTTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagt                                            >  1:1754535/1‑70 (MQ=255)
                             gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagtcag                                         >  1:1925383/1‑70 (MQ=255)
                             gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTtcagtc                                           >  1:1396776/1‑68 (MQ=255)
                             gCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCATCTtcagtc                                           >  1:135874/1‑68 (MQ=255)
                                   aGATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCAt                                   >  1:1013690/1‑70 (MQ=255)
                                    gATCTGGAGGCGAAGCTGGAAGCCCTACAGAAGCAAg                                                                   >  1:1019039/1‑37 (MQ=255)
                                     aTCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTg                                                                 <  1:1214783/38‑1 (MQ=255)
                                       cTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGc                                                               >  1:1004637/1‑38 (MQ=255)
                                       cTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGc                                                               >  1:663244/1‑38 (MQ=255)
                                                aaGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                      >  1:1841556/1‑70 (MQ=255)
                                                aaGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                      >  1:1136541/1‑70 (MQ=255)
                                                  gCTGGAAGCCCTACAGACGCAAGTGGCGTATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAg                   >  1:790415/1‑71 (MQ=255)
                                                        aGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgcttg             >  1:1117503/1‑71 (MQ=255)
                                                        aGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgcttg             >  1:913697/1‑71 (MQ=255)
                                                        aGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgctag             >  1:513007/1‑69 (MQ=255)
                                                            cTACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                      <  1:1199347/58‑1 (MQ=255)
                                                             tACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCaaa                      <  1:1241580/57‑1 (MQ=255)
                                                             tACAGAAGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgcttgct           <  1:2120293/68‑1 (MQ=255)
                                                              aCAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGtgcttgc            >  1:1701681/1‑66 (MQ=255)
                                                              aCAGAAGCAAGTGGCGGATGCTTCCTTCTtcagtc                                           >  1:1986522/1‑35 (MQ=255)
                                                               cAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGatat      >  1:923069/1‑71 (MQ=255)
                                                               cAGACGCAAGTGGCCGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGatat      >  1:1296871/1‑71 (MQ=255)
                                                                   cGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATATGgct  >  1:104825/1‑71 (MQ=255)
                                                                   cGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATATGgct  >  1:55944/1‑71 (MQ=255)
                                                                   |                                                                      
CAGCGCGAACTGGAGCAGCTACCGCAATTGCTCGAAGATCTGGAGGCGAAGCTGGAAGCCCTACAGACGCAAGTGGCGGATGCTTCCTTCTTCAGTCAGCCGCATGAGCAGACGCAAAAAGTGCTTGCTGATATGGCT  >  minE/668361‑668498

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: