Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,726,583 C→G 7.0% R167P (CGT→CCT)  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,5830CG 7.0% 228.6 / 5.7 71R167P (CGT→CCT) yheSfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base G (3/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

TTTCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTC  >  minE/2726515‑2726653
                                                                    |                                                                      
tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                                      >  1:1953250/1‑71 (MQ=255)
tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                                      >  1:4025726/1‑71 (MQ=255)
tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                                      >  1:1165499/1‑71 (MQ=255)
tttCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACgg                                                                      >  1:2155117/1‑71 (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)
             aaCGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCt                                                         >  1:4290818/1‑71 (MQ=255)
              aCGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTg                                                        >  1:3580068/1‑71 (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:762173/71‑1 (MQ=255)
                 gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                                     <  1:3104206/71‑1 (MQ=255)
                      gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCaa                                                   <  1:1391996/68‑1 (MQ=255)
                      gagatcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCaa                                                   <  1:577013/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:698128/1‑68 (MQ=12)
                          tcgagGTGGTTGGTCGGTTCGTCGAGCAACAGCAGGTCGGACGGGCACATCAGCGCCTGCGCCAGGtt                                               >  1:891876/1‑68 (MQ=12)
                          tcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAg                                            <  1:2313079/71‑1 (MQ=255)
                           cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                                     >  1:3591240/1‑45 (MQ=255)
                           cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                                     >  1:3517874/1‑45 (MQ=255)
                           cgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGc                                                                     >  1:1042042/1‑45 (MQ=255)
                            gagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                           <  1:3818389/70‑1 (MQ=255)
                            gagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGa                                           >  1:2821488/1‑70 (MQ=255)
                               gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                                       >  1:3172757/1‑71 (MQ=255)
                               gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                                       <  1:3978063/71‑1 (MQ=255)
                               gtggtTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                                       <  1:2489140/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:3975827/71‑1 (MQ=255)
                                                        agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt              <  1:1470357/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:1957139/38‑1 (MQ=255)
                                                              tCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACg                                         <  1:2686455/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)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:3510947/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:1746519/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:2622069/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:679216/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:3044434/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:1850160/1‑37 (MQ=255)
                                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                     >  1:1168636/1‑37 (MQ=255)
                                                                    cGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTc  <  1:2746179/71‑1 (MQ=255)
                                                                    |                                                                      
TTTCCAGCCAGATAACGGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTC  >  minE/2726515‑2726653

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: