Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,151 G→C 7.4% intergenic (+411/‑567) 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,1510GC 7.4% 277.6 / 9.4 94intergenic (+411/‑567)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base C (4/3);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTT  >  minE/2220081‑2220219
                                                                      |                                                                    
tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTg                                                                      >  1:4081791/1‑71 (MQ=255)
tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTg                                                                      >  1:240438/1‑71 (MQ=255)
tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTg                                                                      >  1:5572976/1‑71 (MQ=255)
   gCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTggg                                                                    >  1:1358528/1‑70 (MQ=255)
   gCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTggg                                                                    >  1:4081054/1‑70 (MQ=255)
     aTTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGt                                                                   <  1:1154700/69‑1 (MQ=255)
     aTTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTaa                                                                 <  1:4585115/71‑1 (MQ=255)
        gTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                              >  1:4129086/1‑71 (MQ=255)
        gTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                              >  1:3215650/1‑71 (MQ=255)
        gTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                              >  1:3547341/1‑71 (MQ=255)
         tCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATccc                                                             <  1:1612330/71‑1 (MQ=255)
         tCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATccc                                                             <  1:1493587/71‑1 (MQ=255)
         tCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                              <  1:3755849/70‑1 (MQ=255)
          cAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATccc                                                             >  1:4383629/1‑70 (MQ=255)
          cAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATccc                                                             >  1:3634342/1‑70 (MQ=255)
           aaTCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcccc                                                            <  1:1749445/70‑1 (MQ=255)
            aTCGGTTATTTTTGGGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                              <  1:2301052/67‑1 (MQ=255)
             tCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATc                                                               <  1:2258334/65‑1 (MQ=255)
             tCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGc                                                         <  1:3751093/71‑1 (MQ=255)
              cGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGc                                                         >  1:4576284/1‑70 (MQ=255)
              cGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCAAGACGATGCAGTGAAAATTGGGTAATCCCCAGc                                                         >  1:1518973/1‑70 (MQ=255)
               ggTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCaa                                                       >  1:1695143/1‑71 (MQ=255)
                 ttAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTa                                                                  <  1:5219739/58‑1 (MQ=255)
                   aaTTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTaa                                                                 >  1:5341459/1‑57 (MQ=255)
                   aaTTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAAt                                                                >  1:66742/1‑58 (MQ=255)
                       ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATccc                                                             >  1:4346937/1‑57 (MQ=255)
                       ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGc                                                >  1:3907912/1‑70 (MQ=255)
                       ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCg                                               >  1:5294878/1‑71 (MQ=255)
                       ttGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCg                                               >  1:48508/1‑71 (MQ=255)
                        tGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTaa                                                                 >  1:5617110/1‑52 (MQ=255)
                        tGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTaa                                                                 >  1:3501263/1‑52 (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:481108/1‑71 (MQ=255)
                             gCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgt                                         >  1:1633661/1‑71 (MQ=255)
                              cTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAc                                                      >  1:4526042/1‑57 (MQ=255)
                              cTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtc                                        >  1:4470631/1‑71 (MQ=255)
                                 tAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtcgt                                      <  1:3038778/70‑1 (MQ=255)
                                    cATTCACATCTATCCAGACGATGCAGTGAAAATTgg                                                                     >  1:2314996/1‑36 (MQ=255)
                                     aTTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATc                                 >  1:2328295/1‑71 (MQ=255)
                                      ttCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCaa                                                       >  1:850230/1‑48 (MQ=255)
                                      ttCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCa                                                        <  1:5426005/47‑1 (MQ=255)
                                        cacaTCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAc                                                      >  1:3265866/1‑47 (MQ=255)
                                        cacaTCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAc                                                      >  1:1742366/1‑47 (MQ=255)
                                        cacaTCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAc                                                      >  1:2252702/1‑47 (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:4653670/68‑1 (MQ=255)
                                              aacccAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                        <  1:4306231/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:3139613/68‑1 (MQ=255)
                                                   aggcGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTg                   >  1:2603526/4‑71 (MQ=255)
                                                       gATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt               >  1:2346852/1‑71 (MQ=255)
                                                       gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                        >  1:3365152/1‑62 (MQ=255)
                                                       gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                        >  1:1681589/1‑62 (MQ=255)
                                                       gATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGa                        >  1:1835493/1‑62 (MQ=255)
                                                          gCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc            >  1:5396563/1‑71 (MQ=255)
                                                            aGTGAAAATTGGGTAATCCCCAGCTACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTa          <  1:475777/71‑1 (MQ=255)
                                                            aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                            >  1:3432889/1‑37 (MQ=255)
                                                            aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                            >  1:439939/1‑37 (MQ=255)
                                                            aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTaa                                            >  1:5144609/1‑37 (MQ=255)
                                                            aGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTc                            >  1:4740542/1‑53 (MQ=255)
                                                             gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgcg                         <  1:2699941/55‑1 (MQ=255)
                                                             gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt               <  1:1836574/65‑1 (MQ=255)
                                                             gTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCAt               <  1:3595115/65‑1 (MQ=255)
                                                              tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                     <  1:4539228/58‑1 (MQ=255)
                                                              tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                   <  1:1083546/60‑1 (MQ=255)
                                                              tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                   <  1:3516440/60‑1 (MQ=255)
                                                              tGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                   <  1:4295675/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:3975953/54‑1 (MQ=255)
                                                                 aaaTTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATc                      <  1:577148/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:5149069/41‑1 (MQ=255)
                                                                     tGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                               <  1:4014550/41‑1 (MQ=255)
                                                                      |                                                                    
TGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTT  >  minE/2220081‑2220219

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: