Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,726,600 G→C 7.2% A161A (GCC→GCG 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,6000GC 7.2% 219.9 / 6.1 69A161A (GCC→GCGyheSfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base C (3/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAGTTGTTCATTGCT  >  minE/2726531‑2726668
                                                                     |                                                                    
ggCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTggg                                                                     >  1:1460593/1‑71 (MQ=255)
 gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTggg                                                                     >  1:3071159/1‑70 (MQ=255)
 gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                                                    <  1:762173/71‑1 (MQ=255)
 gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                                                    <  1:3104206/71‑1 (MQ=255)
 gCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGc                                                                    <  1:1431423/71‑1 (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)
         atcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTaa                                                            <  1:1432369/71‑1 (MQ=255)
         atcgagGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTaa                                                            <  1:592740/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)
            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:2414877/71‑1 (MQ=255)
                                     agcagcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATc                                <  1:457950/71‑1 (MQ=255)
                                        agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt                             <  1:1044827/71‑1 (MQ=255)
                                        agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt                             <  1:3975827/71‑1 (MQ=255)
                                        agcaAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACt                             <  1:1470357/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:679759/39‑1 (MQ=255)
                                               cTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCa                                                      <  1:3165588/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:2622069/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:3044434/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:1168636/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:1850160/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:3510947/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:679216/1‑37 (MQ=255)
                                                   aCGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATa                                                    >  1:1746519/1‑37 (MQ=255)
                                                    cGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTc                 <  1:2746179/71‑1 (MQ=255)
                                                     ggCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTcc                <  1:840226/71‑1 (MQ=255)
                                                       cAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAg              <  1:3474123/71‑1 (MQ=255)
                                                          aTCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAgttg           <  1:3799255/71‑1 (MQ=255)
                                                          aTCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAg              >  1:3164248/1‑68 (MQ=255)
                                                            cAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCa               >  1:2301216/1‑65 (MQ=255)
                                                            cAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAgttgtt         >  1:2910513/1‑71 (MQ=255)
                                                             aGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCAc                              >  1:1307499/1‑49 (MQ=255)
                                                             aGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAgttgtt         <  1:369061/70‑1 (MQ=255)
                                                             aGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAgttgtt         <  1:3423204/70‑1 (MQ=255)
                                                              gcgcCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAGTTGTTCa       <  1:4100001/71‑1 (MQ=255)
                                                                   tGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTAc                          >  1:2411131/1‑47 (MQ=255)
                                                                   tGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTAc                          >  1:1202679/1‑47 (MQ=255)
                                                                   tGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTAc                          >  1:3218208/1‑47 (MQ=255)
                                                                     ggCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATc                                >  1:608534/1‑39 (MQ=255)
                                                                     ggCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAGTTGTTCATTGCt  <  1:2936268/69‑1 (MQ=255)
                                                                     |                                                                    
GGCATCGAGATCGAGGTGGTTAGTCGGTTCGTCGAGCAGCAGCAAGTCTGAACGGCAAATCAGCGCCTGGGCAAGGTTAAGACGCATACGCCAGCCCCCGGAAAAATCACTTACCGGGCGCTCCAGTTGTTCATTGCT  >  minE/2726531‑2726668

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: