Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,220,238 C→T 9.4% 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,2380CT 9.4% 147.6 / 4.8 53intergenic (+498/‑480)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base C (18/30);  new base T (3/2);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCC  >  minE/2220168‑2220305
                                                                      |                                                                   
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:3357673/71‑1 (MQ=255)
         aTGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAg                                                            <  1:3578032/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)
                        tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTg                                                               <  1:3474209/53‑1 (MQ=255)
                        tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCa                                             <  1:2276963/71‑1 (MQ=255)
                           cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAgg                                           <  1:2946723/70‑1 (MQ=255)
                             aTCCTGGCATCCCTACATTATTTGTGGTTTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                        >  1:2966677/1‑71 (MQ=255)
                             aTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                        >  1:3882101/1‑71 (MQ=255)
                              tCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctg                                       <  1:2441963/71‑1 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttt                                  <  1:2145023/71‑1 (MQ=255)
                                   gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttt                                  <  1:2678425/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)
                                                            tgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaag         <  1:2770125/71‑1 (MQ=255)
                                                              aaGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttt                                  <  1:2309994/44‑1 (MQ=255)
                                                              aaGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaagaa       <  1:3737436/71‑1 (MQ=255)
                                                              aaGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACCTTATaaga        <  1:3274971/70‑1 (MQ=255)
                                                               aGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttca                       <  1:3621684/54‑1 (MQ=255)
                                                                 aTTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                        >  1:1956931/1‑35 (MQ=255)
                                                                  ttATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGtt    <  1:306202/70‑1 (MQ=255)
                                                                   tATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTcc  >  1:1724198/1‑71 (MQ=255)
                                                                    aTCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTcc  <  1:555951/70‑1 (MQ=255)
                                                                      |                                                                   
CGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCC  >  minE/2220168‑2220305

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: