Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,168 C→A 11.8% intergenic (+428/‑550) 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,1680CA11.8% 280.7 / 14.8 93intergenic (+428/‑550)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base A (7/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

TTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATC  >  minE/2220104‑2220238
                                                                |                                                                      
ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGc                                                                   >  1:3907912/1‑70 (MQ=255)
ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCg                                                                  >  1:5294878/1‑71 (MQ=255)
ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCg                                                                  >  1:48508/1‑71 (MQ=255)
  gCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                 <  1:3470197/70‑1 (MQ=255)
  gCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                 <  1:1463045/70‑1 (MQ=255)
    gTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                 <  1:1023665/68‑1 (MQ=255)
     tGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTa                                                                <  1:1816896/68‑1 (MQ=255)
     tGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATg                                                             <  1:918950/71‑1 (MQ=255)
      gCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgt                                                            >  1:1633661/1‑71 (MQ=255)
      gCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgt                                                            >  1:481108/1‑71 (MQ=255)
       cTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtc                                                           >  1:4470631/1‑71 (MQ=255)
          tAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtcgt                                                         <  1:3038778/70‑1 (MQ=255)
              aTTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATc                                                    >  1:2328295/1‑71 (MQ=255)
                  acaTCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTa                                                                >  1:1958122/1‑55 (MQ=255)
                  acaTCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTa                                                                >  1:502112/1‑55 (MQ=255)
                     tCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                 <  1:5001991/51‑1 (MQ=255)
                     tCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCt                                                   >  1:725991/1‑65 (MQ=255)
                      cTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgc                                             >  1:3169424/1‑70 (MQ=255)
                       aacccAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                                           <  1:4306231/68‑1 (MQ=255)
                       aacccAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                                           <  1:4653670/68‑1 (MQ=255)
                        aTCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGa                                           >  1:468223/1‑70 (MQ=255)
                         tCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgc                                             <  1:1828754/67‑1 (MQ=255)
                         tCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGa                                           <  1:5020788/69‑1 (MQ=255)
                         tCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:4505061/71‑1 (MQ=255)
                           cAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                        <  1:4283239/70‑1 (MQ=255)
                            aggcGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTg                                      >  1:2603526/4‑71 (MQ=255)
                            aggcGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTg                                      <  1:3139613/68‑1 (MQ=255)
                                gATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt                                  >  1:2346852/1‑71 (MQ=255)
                                gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                                           >  1:1835493/1‑62 (MQ=255)
                                gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                                           >  1:1681589/1‑62 (MQ=255)
                                gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                                           >  1:3365152/1‑62 (MQ=255)
                                   gCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                               >  1:5396563/1‑71 (MQ=255)
                                     aGTGAAAATTGGGTAATCCCCAGCTACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTa                             <  1:475777/71‑1 (MQ=255)
                                     aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                                               >  1:439939/1‑37 (MQ=255)
                                     aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                                               >  1:5144609/1‑37 (MQ=255)
                                     aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                                               >  1:3432889/1‑37 (MQ=255)
                                     aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTc                                               >  1:4740542/1‑53 (MQ=255)
                                      gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgcg                                            <  1:2699941/55‑1 (MQ=255)
                                      gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt                                  <  1:3595115/65‑1 (MQ=255)
                                      gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt                                  <  1:1836574/65‑1 (MQ=255)
                                       tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                        <  1:4539228/58‑1 (MQ=255)
                                       tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                      <  1:3516440/60‑1 (MQ=255)
                                       tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                      <  1:4295675/60‑1 (MQ=255)
                                       tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                      <  1:1083546/60‑1 (MQ=255)
                                        gAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                  >  1:1414804/1‑47 (MQ=255)
                                        gAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         >  1:2118056/1‑56 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATgtcg                                                          >  1:2629802/1‑38 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTAtaa                                                   >  1:1496903/1‑43 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTc                                               >  1:3933827/1‑49 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                      >  1:925135/1‑58 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                      >  1:5627965/1‑58 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACat                          >  1:267293/1‑70 (MQ=255)
                                         aaaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGAACCTGGCATCCCTACat                          >  1:4004785/1‑70 (MQ=255)
                                          aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:3311709/54‑1 (MQ=255)
                                          aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:577148/54‑1 (MQ=255)
                                          aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:3975953/54‑1 (MQ=255)
                                           aaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                       >  1:3353499/1‑71 (MQ=255)
                                            aTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                      <  1:3386335/71‑1 (MQ=255)
                                            aTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                      <  1:5607574/71‑1 (MQ=255)
                                             ttGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                      >  1:2948697/1‑70 (MQ=255)
                                             ttGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                      >  1:2122867/1‑70 (MQ=255)
                                             ttGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                     <  1:2925540/71‑1 (MQ=255)
                                             ttGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                     <  1:162078/71‑1 (MQ=255)
                                              tGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                  <  1:4014550/41‑1 (MQ=255)
                                              tGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                  <  1:5149069/41‑1 (MQ=255)
                                                  tAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:2871434/46‑1 (MQ=255)
                                                  tAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                                         <  1:1427390/46‑1 (MQ=255)
                                                   aaTCCCCAGCAACCGCTGCGTAATGTCGTGGATCTTGTCGCGTTCCTGGCATCCCTACATTTTTTGTgg                 >  1:4540535/1‑69 (MQ=255)
                                                   aaTCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTAc                            <  1:3924974/58‑1 (MQ=255)
                                                    aTCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCt                                       >  1:5601656/1‑46 (MQ=255)
                                                    aTCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTc               <  1:1818535/70‑1 (MQ=255)
                                                      ccccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgt            >  1:627322/1‑71 (MQ=255)
                                                       cccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg           <  1:2005162/71‑1 (MQ=255)
                                                       cccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg           <  1:1736173/71‑1 (MQ=255)
                                                        ctggCAACAGTTGCATAATTTCGTCTATCTTGTGGCGATTCTGGCCCCAATACATTATCTGTGGGCCGTGa          <  1:4910396/68‑1 (MQ=255)
                                                        ctggCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCcgtg           >  1:4042929/4‑70 (MQ=255)
                                                        ccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCt                                       <  1:1148480/42‑1 (MQ=255)
                                                         cAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg           >  1:3741447/1‑69 (MQ=255)
                                                          ggCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCcgtg           >  1:3262699/2‑68 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                        >  1:1471143/1‑38 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                        >  1:2532016/1‑38 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                               <  1:2398680/47‑1 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                               <  1:3870149/47‑1 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                               <  1:4257752/47‑1 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                               <  1:4311493/47‑1 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                >  1:673363/1‑62 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                >  1:5648965/1‑62 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                >  1:3018610/1‑62 (MQ=255)
                                                           gCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGa       <  1:1257394/71‑1 (MQ=255)
                                                           gCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCcgtg           >  1:1397027/1‑67 (MQ=255)
                                                             aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                        <  1:204342/52‑1 (MQ=255)
                                                             aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                        <  1:2460804/52‑1 (MQ=255)
                                                             aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt     <  1:5545810/71‑1 (MQ=255)
                                                                cGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGc                                    >  1:4422049/1‑37 (MQ=255)
                                                                cGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtc  >  1:2923897/1‑71 (MQ=255)
                                                                cGCTGCGTAATGTCGTCTAACTTGTCGCGATCCTGGc                                    >  1:5650794/1‑37 (MQ=255)
                                                                |                                                                      
TTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATC  >  minE/2220104‑2220238

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: