Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,726,579 T→G 6.5% S168S (TCA→TCC yheS ← 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
*minE2,726,5790TG 6.5% 211.5 / 10.0 77S168S (TCA→TCCyheSfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (32/40);  new base G (3/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CCATTTTTCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACC  >  minE/2726510‑2726645
                                                                     |                                                                  
ccATTTTTCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:3323642/71‑1 (MQ=255)
     tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                              >  1:1953250/1‑71 (MQ=255)
     tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                              >  1:1165499/1‑71 (MQ=255)
     tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                              >  1:2155117/1‑71 (MQ=255)
     tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                              >  1:4025726/1‑71 (MQ=255)
      ttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCt                                                                    >  1:1342854/1‑64 (MQ=255)
      ttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCt                                                                    >  1:3216140/1‑64 (MQ=255)
            ccagATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCa                                                       <  1:1410610/71‑1 (MQ=255)
             ccgATACCGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCa                                                       <  1:1457253/68‑1 (MQ=255)
             cataTAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAg                                                      >  1:1769894/4‑71 (MQ=255)
             cagATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAg                                                      >  1:707835/1‑71 (MQ=255)
                aTAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:1192800/55‑1 (MQ=255)
                aTAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:3436228/55‑1 (MQ=255)
                aTAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:3739407/55‑1 (MQ=255)
                 tAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:20068/54‑1 (MQ=255)
                  aaCGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCt                                                 >  1:4290818/1‑71 (MQ=255)
                   aCGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTg                                                >  1:3580068/1‑71 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:21191/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:170956/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:4002062/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:3169035/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:2430409/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:3410204/1‑52 (MQ=255)
                     ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTggg                                              >  1:1460593/1‑71 (MQ=255)
                      gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTggg                                              >  1:3071159/1‑70 (MQ=255)
                      gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                             <  1:1431423/71‑1 (MQ=255)
                      gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                             <  1:3104206/71‑1 (MQ=255)
                      gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                             <  1:762173/71‑1 (MQ=255)
                         tcgagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTg                                                                   <  1:908211/46‑1 (MQ=255)
                           gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 >  1:123450/1‑46 (MQ=255)
                           gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCaa                                           <  1:577013/68‑1 (MQ=255)
                           gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCaa                                           <  1:1391996/68‑1 (MQ=255)
                           gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCaa                                           <  1:2919074/68‑1 (MQ=255)
                           gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGt                                        <  1:3586244/71‑1 (MQ=255)
                            agatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGtt                                       <  1:915776/71‑1 (MQ=255)
                            agatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGt                                        >  1:777171/1‑70 (MQ=255)
                            agatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGt                                        >  1:2544637/1‑70 (MQ=255)
                              atcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACg                                                               <  1:2719427/45‑1 (MQ=255)
                              atcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTaa                                     <  1:592740/71‑1 (MQ=255)
                              atcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTaa                                     <  1:1432369/71‑1 (MQ=255)
                               tcgagGTGGTTGGTCGGTTCGTCGAGCAACAGCAGGTCGGACGGGCACATCAGCGCCTGCGCCAGGtt                                       >  1:891876/1‑68 (MQ=12)
                               tcgagGTGGTTGGTCGGTTCGTCGAGCAACAGCAGGTCGGACGGGCACATCAGCGCCTGCGCCAGGtt                                       >  1:698128/1‑68 (MQ=12)
                               tcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAg                                    <  1:2313079/71‑1 (MQ=255)
                                cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                             >  1:1042042/1‑45 (MQ=255)
                                cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                             >  1:3591240/1‑45 (MQ=255)
                                cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                             >  1:3517874/1‑45 (MQ=255)
                                 gagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                   >  1:2821488/1‑70 (MQ=255)
                                 gagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                   <  1:3818389/70‑1 (MQ=255)
                                   ggtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGaa                                                                 <  1:1269891/38‑1 (MQ=255)
                                    gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:2489140/71‑1 (MQ=255)
                                    gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               >  1:3172757/1‑71 (MQ=255)
                                    gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:3978063/71‑1 (MQ=255)
                                      ggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                             >  1:1879184/1‑71 (MQ=255)
                                       gtTGGTCGGTTCGTCGAGCAGCAGCAGGTCGGAAGGGCACATCAGCGCCTGCGCCAGGTTCAGGCg                                 <  1:3835698/66‑1 (MQ=12)
                                       gtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATAc                            <  1:3482446/71‑1 (MQ=255)
                                          aGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                   >  1:658986/1‑61 (MQ=255)
                                          aGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                   >  1:2321576/1‑61 (MQ=255)
                                          aGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGc                          >  1:4308151/1‑70 (MQ=255)
                                                  cgtcgAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTaa                                     >  1:2524353/1‑51 (MQ=255)
                                                    tcgAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCgg                >  1:1448014/1‑70 (MQ=255)
                                                      gAGCAGCAGCAGGTCGGAAGGGCACATCAGCGCCTGCGCCAGGTTCAGGCGCATGCGCCAGCCGCCCGag              >  1:1745326/1‑69 (MQ=255)
                                                      gAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:53758/53‑1 (MQ=255)
                                                       agcagcagcAGGTCGGAAGGGCACATCAGCGCCTGCGCCAGGTTCAGGCGCATGCGCCAGCCGCCCGagaa            <  1:2994622/71‑1 (MQ=255)
                                                         cagcagcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAAt          <  1:3866905/71‑1 (MQ=255)
                                                          agcagcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATc         <  1:457950/71‑1 (MQ=255)
                                                          agcagcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATc         <  1:2414877/71‑1 (MQ=255)
                                                             agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt      <  1:1470357/71‑1 (MQ=255)
                                                             agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt      <  1:3975827/71‑1 (MQ=255)
                                                             agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt      <  1:1044827/71‑1 (MQ=255)
                                                              gcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACtt     <  1:1021921/71‑1 (MQ=255)
                                                                 aGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTAcc  <  1:2741484/71‑1 (MQ=255)
                                                                  gTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCgg                <  1:3578605/56‑1 (MQ=255)
                                                                   tCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACg                                 <  1:2686455/38‑1 (MQ=255)
                                                                   tCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACg                                 <  1:1957139/38‑1 (MQ=255)
                                                                    cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:3165588/39‑1 (MQ=255)
                                                                    cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:679759/39‑1 (MQ=255)
                                                                    cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               <  1:3807066/39‑1 (MQ=255)
                                                                    cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                             >  1:4431460/1‑41 (MQ=255)
                                                                    cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                             >  1:1769781/1‑41 (MQ=255)
                                                                     tGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                               >  1:2938612/1‑38 (MQ=255)
                                                                     |                                                                  
CCATTTTTCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACC  >  minE/2726510‑2726645

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: