Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,238 C→T 29.2% intergenic (+498/‑480) 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,2380CT29.2% 97.4 / 43.8 65intergenic (+498/‑480)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (5/14);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTT  >  minE/2220168‑2220303
                                                                      |                                                                 
cGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtc                                                                   >  1:2923897/1‑71 (MQ=255)
 gCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtct                                                                  >  1:4713829/1‑71 (MQ=255)
  ttGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTAtttc                                                                 <  1:2915009/70‑1 (MQ=255)
  ttGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTAtttc                                                                 <  1:431565/70‑1 (MQ=255)
  cTGCGTAATGTCGTCTATCTTGTCGCGCTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtct                                                                  >  1:4510203/1‑70 (MQ=255)
   tGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtctc                                                                 <  1:1517157/70‑1 (MQ=255)
   tGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCa                                                                >  1:5541425/1‑71 (MQ=255)
     cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                <  1:1744959/69‑1 (MQ=255)
     cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                <  1:1943736/69‑1 (MQ=255)
     cGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCAt                                                               >  1:5020110/1‑70 (MQ=255)
      gTAATGTCGTCTCTCTTGTCGCGATCCTGGCATCCCTACATTATTTGTTGTCTGTGAAGATTATCTCATTg                                                             <  1:1870458/71‑1 (MQ=255)
      gTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTg                                                             <  1:3085671/71‑1 (MQ=255)
        aaTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTg                                                             >  1:917296/1‑69 (MQ=255)
            tcgtcTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGccg                                                       >  1:2110651/1‑70 (MQ=255)
            tcgtcTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGccc                                                       >  1:3997921/1‑71 (MQ=255)
               tcTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTc                                                    <  1:4592050/71‑1 (MQ=255)
                 tATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtct                                                                  <  1:1077216/55‑1 (MQ=255)
                 tATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCCCAt                                                  >  1:5414755/1‑71 (MQ=255)
                   tCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCa                                                                <  1:1933036/55‑1 (MQ=255)
                   tCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCa                                                                <  1:1778168/55‑1 (MQ=255)
                   tCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATc                                                 <  1:3981274/70‑1 (MQ=255)
                   tCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                >  1:675890/1‑55 (MQ=255)
                      tGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGccc                                                       >  1:4032344/1‑61 (MQ=255)
                      tGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGccc                                                       >  1:5464369/1‑61 (MQ=255)
                       gTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCt                                                <  1:465123/67‑1 (MQ=255)
                       gTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGc                                            >  1:1666235/1‑71 (MQ=255)
                       gTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCtt                                               <  1:483433/68‑1 (MQ=255)
                       gTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCg                                             <  1:876015/70‑1 (MQ=255)
                        tCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGct                                           >  1:610291/1‑70 (MQ=255)
                        tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGc                                            <  1:5574589/70‑1 (MQ=255)
                          gcgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAg                                          <  1:2060615/70‑1 (MQ=255)
                           cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTc                                                    <  1:3710448/59‑1 (MQ=255)
                           cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTc                                                    <  1:3938216/59‑1 (MQ=255)
                           cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTc                                                    <  1:2089726/59‑1 (MQ=255)
                             aTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                      <  1:5189383/71‑1 (MQ=255)
                                 tGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtct                                                                  <  1:141270/39‑1 (MQ=255)
                                 tGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAtct                                                                  <  1:968893/39‑1 (MQ=255)
                                  ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAg                                          <  1:2581014/62‑1 (MQ=255)
                                  ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctggcgt                                 >  1:61763/1‑69 (MQ=255)
                                  ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctt                                 >  1:4919699/1‑71 (MQ=255)
                                  ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctt                                 >  1:68418/1‑71 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                >  1:985965/1‑39 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                >  1:4842262/1‑39 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCa                                                                >  1:598487/1‑39 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctggct                                  <  1:3455878/69‑1 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctggct                                  <  1:84931/69‑1 (MQ=255)
                                     aTCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctttt                               >  1:4336204/1‑70 (MQ=255)
                                        ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagc                           >  1:873650/1‑71 (MQ=255)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGccg                                                       <  1:905693/41‑2 (MQ=16)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGccg                                                       <  1:5139796/41‑2 (MQ=16)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGccg                                                       <  1:5179624/41‑2 (MQ=16)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTc                                                    <  1:2956533/44‑1 (MQ=16)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTc                                                    <  1:685174/44‑1 (MQ=16)
                                          tACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTc                                                    <  1:5512540/44‑1 (MQ=16)
                                            cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                       >  1:4523401/1‑71 (MQ=255)
                                            cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                       >  1:4992352/1‑71 (MQ=255)
                                            cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                       >  1:4123580/1‑71 (MQ=255)
                                            cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctt                         <  1:1391394/69‑1 (MQ=255)
                                            cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctt                         <  1:5504834/69‑1 (MQ=255)
                                            cATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGgctgtactgctttt                       <  1:3746792/71‑1 (MQ=255)
                                             attattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttc                      <  1:2441561/71‑1 (MQ=255)
                                             attattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttc                      <  1:4611400/71‑1 (MQ=255)
                                             attattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctt                         <  1:4298870/68‑1 (MQ=255)
                                              ttatcTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCtac                                              <  1:3175017/46‑3 (MQ=16)
                                                attTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCa                                           >  1:2353023/1‑47 (MQ=255)
                                                attTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCAcc                   >  1:606583/1‑71 (MQ=255)
                                                    gtgGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTAc               >  1:1561470/1‑71 (MQ=255)
                                                    gtgGTCTGTGAAGATTATCCCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTAc               >  1:2110177/1‑71 (MQ=255)
                                                    gtgGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGgctgtactgctt                         >  1:1615761/1‑61 (MQ=255)
                                                        tccgtgAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGta           <  1:714237/68‑1 (MQ=255)
                                                          tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat          >  1:3065471/1‑70 (MQ=255)
                                                            tgAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTa                <  1:364042/62‑1 (MQ=255)
                                                            tgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaag       <  1:3122521/71‑1 (MQ=255)
                                                                gATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGt   <  1:1220387/71‑1 (MQ=255)
                                                                 attattTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGtt  <  1:4528504/71‑1 (MQ=255)
                                                                      |                                                                 
CGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTT  >  minE/2220168‑2220303

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: