Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,199 C→T 15.8% intergenic (+459/‑519) yhdH → / → accB predicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,220,1990CT15.8% 150.4 / 17.9 63intergenic (+459/‑519)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (5/5);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

CAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTG  >  minE/2220131‑2220268
                                                                    |                                                                     
cAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCt                                                                     <  1:1963329/71‑1 (MQ=255)
 aGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                                                    >  1:1632614/1‑71 (MQ=255)
   aCGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGc                                                                  <  1:1735119/71‑1 (MQ=255)
    cGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCa                                                                 >  1:3223493/1‑71 (MQ=255)
     gATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGc                                                                  >  1:3273377/1‑69 (MQ=255)
       tGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCATGTCGCGATCCTGGCATcc                                                              >  1:725309/1‑71 (MQ=255)
       tGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGc                                                                  >  1:2068044/1‑67 (MQ=255)
         cAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                                                             >  1:2237226/1‑70 (MQ=255)
              aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACat                                                        >  1:1959895/1‑70 (MQ=255)
               aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                                                      <  1:537918/71‑1 (MQ=255)
               aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACa                                                         >  1:3355265/1‑68 (MQ=255)
                 aTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                                                    >  1:362339/1‑71 (MQ=255)
                 aTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                                                    >  1:1730527/1‑71 (MQ=255)
                    gggTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtg                                                  <  1:3593622/70‑1 (MQ=255)
                     ggTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                                                   <  1:3064823/68‑1 (MQ=255)
                       tAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTgg                                               >  1:899234/1‑70 (MQ=255)
                        aaTCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTc                                             >  1:3277027/1‑71 (MQ=255)
                        aaTCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTc                                             >  1:1333357/1‑71 (MQ=255)
                          tCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCt                                            >  1:1493443/1‑70 (MQ=255)
                            cccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTa                                                           >  1:3750518/1‑53 (MQ=255)
                            cccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg                                         <  1:2076917/71‑1 (MQ=255)
                                gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                                                      <  1:3631409/38‑1 (MQ=255)
                                gCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCcgt                                          >  1:746696/1‑66 (MQ=255)
                                gCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAg                                      <  1:2025432/70‑1 (MQ=255)
                                  aaacGCTGCGTAATGGCGTCTATCTTGTCGCGATCCTGGCTTGCCTACATTATTTGTGGTCTTTGAAGAtt                                   <  1:932014/68‑1 (MQ=255)
                                  aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                                                       >  1:1554008/1‑35 (MQ=255)
                                  aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt                                   <  1:3621841/71‑1 (MQ=255)
                                  aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt                                   <  1:1592181/71‑1 (MQ=255)
                                   aCCGCTGCGTAATGTCGTCTATCTTGTCGCGATcg                                                                      <  1:114431/35‑2 (MQ=255)
                                     cGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtc                                <  1:560546/71‑1 (MQ=255)
                                      gCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACa                                                         >  1:181925/1‑45 (MQ=255)
                                      gCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACa                                                         >  1:3457964/1‑45 (MQ=255)
                                         gcatAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGGCCGTGAAGATTATTTCa                             <  1:2331853/67‑1 (MQ=255)
                                         gCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCAt                            >  1:3208725/1‑71 (MQ=255)
                                          cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                             <  1:1021865/69‑1 (MQ=255)
                                          cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                             <  1:3613836/69‑1 (MQ=255)
                                          cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCAtt                           <  1:1665413/71‑1 (MQ=255)
                                            tAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                                                             <  1:2851323/35‑1 (MQ=255)
                                              aTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAg                       <  1:3578032/71‑1 (MQ=255)
                                              aTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAg                       <  1:3357673/71‑1 (MQ=255)
                                               tGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                                                     <  1:2685490/40‑1 (MQ=255)
                                                  cgtcTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGcccc                   >  1:43452/1‑71 (MQ=255)
                                                   gtcTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGcccat                  >  1:2706800/1‑69 (MQ=255)
                                                   gtcTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCt                  >  1:1692820/1‑71 (MQ=255)
                                                     cTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCa                >  1:263531/1‑71 (MQ=255)
                                                      tATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAt               >  1:1218499/1‑71 (MQ=255)
                                                      tATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGAAGCCGCTCa                >  1:1786018/1‑70 (MQ=255)
                                                      tATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCAt               >  1:2893/1‑71 (MQ=255)
                                                       aTCTTGTCGCGATTCTGGCCCCAATACATTATCTGTgg                                               <  1:113801/38‑1 (MQ=255)
                                                       aTCTTGTCGCGATTCTGGCCCCAATACATTATCTGTgg                                               <  1:3500207/38‑1 (MQ=255)
                                                       aTCTTGTCGCGATTCTGGCCCCAATACATTATCTGTgg                                               <  1:1514676/38‑1 (MQ=255)
                                                       aTCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATc              >  1:3737366/1‑71 (MQ=255)
                                                        tCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCt             >  1:1963595/1‑71 (MQ=255)
                                                         cTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGt                                              >  1:719520/1‑37 (MQ=255)
                                                            gTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTc           <  1:2189772/69‑1 (MQ=255)
                                                             tCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgt                                          <  1:1790647/37‑1 (MQ=255)
                                                             tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTg                          <  1:3474209/53‑1 (MQ=255)
                                                             tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCa        <  1:2276963/71‑1 (MQ=255)
                                                                cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGa                                        >  1:1734788/1‑36 (MQ=255)
                                                                cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGa                                        >  1:2477830/1‑36 (MQ=255)
                                                                cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGa                                        >  1:2482351/1‑36 (MQ=255)
                                                                cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAgg      <  1:2946723/70‑1 (MQ=255)
                                                                  aTCCTGGCATCCCTACATTATTTGTGGTTTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct   >  1:2966677/1‑71 (MQ=255)
                                                                  aTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct   >  1:3882101/1‑71 (MQ=255)
                                                                   tCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTa                                  <  1:3620748/39‑1 (MQ=255)
                                                                   tCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctg  <  1:2441963/71‑1 (MQ=255)
                                                                    |                                                                     
CAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTG  >  minE/2220131‑2220268

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: