Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,226 T→C 11.8% intergenic (+486/‑492) 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,2260TC11.8% 131.5 / 9.0 51intergenic (+486/‑492)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base C (3/3);  total (22/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 = 9.79e-01

TCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTAT  >  minE/2220157‑2220295
                                                                     |                                                                     
tCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCt                                                                       >  1:1493443/1‑70 (MQ=255)
  cccAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg                                                                    <  1:2076917/71‑1 (MQ=255)
      gCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCcgt                                                                     >  1:746696/1‑66 (MQ=255)
      gCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAg                                                                 <  1:2025432/70‑1 (MQ=255)
        aaacGCTGCGTAATGGCGTCTATCTTGTCGCGATCCTGGCTTGCCTACATTATTTGTGGTCTTTGAAGAtt                                                              <  1:932014/68‑1 (MQ=255)
        aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt                                                              <  1:1592181/71‑1 (MQ=255)
        aaCCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt                                                              <  1:3621841/71‑1 (MQ=255)
           cGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtc                                                           <  1:560546/71‑1 (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)
                    aTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAg                                                  <  1:3578032/71‑1 (MQ=255)
                    aTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAg                                                  <  1:3357673/71‑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)
                             aTCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATc                                         >  1:3737366/1‑71 (MQ=255)
                              tCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCt                                        >  1:1963595/1‑71 (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:2482351/1‑36 (MQ=255)
                                      cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGa                                                                   >  1:2477830/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)
                                              gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttt                        <  1:2678425/71‑1 (MQ=255)
                                              gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttt                        <  1:2145023/71‑1 (MQ=255)
                                               cATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAg                                  <  1:1018372/60‑1 (MQ=255)
                                                aTCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGcctttt                       >  1:3007249/1‑70 (MQ=255)
                                                   ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGc                                                    >  1:2367900/1‑38 (MQ=255)
                                                   ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagc                   <  1:874545/71‑1 (MQ=255)
                                                   ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCATCATCTTCGCAGGGCTGgcttttcagc                   <  1:3012277/71‑1 (MQ=255)
                                                    cTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagct                  >  1:85787/1‑71 (MQ=255)
                                                    cTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagct                  >  1:3951021/1‑71 (MQ=255)
                                                      aCATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttt                >  1:230707/1‑71 (MQ=255)
                                                       cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGc                                    >  1:2659388/1‑50 (MQ=255)
                                                       cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttt                <  1:1748476/70‑1 (MQ=255)
                                                         ttattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt               <  1:1858827/69‑1 (MQ=255)
                                                          tattTGTGGTCCGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt               <  1:3894412/68‑1 (MQ=255)
                                                            tttGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCtt                                       <  1:2077782/42‑1 (MQ=255)
                                                             gattgGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt         <  1:1793023/68‑1 (MQ=255)
                                                              tgtgGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt         <  1:365151/70‑1 (MQ=255)
                                                              tgtgGGCTGTGAAGATTATCTCATTGCAGCCCCTCATCtt                                       <  1:3173215/40‑1 (MQ=255)
                                                                tgGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACg      >  1:1629281/1‑71 (MQ=255)
                                                                     tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttc              <  1:2495725/58‑1 (MQ=255)
                                                                     tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat  <  1:341319/70‑1 (MQ=255)
                                                                     |                                                                     
TCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTAT  >  minE/2220157‑2220295

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: