Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,329,107 T→C 19.6% intergenic (‑24/+125) yiiM ← / ← kdgT conserved hypothetical protein/2‑keto‑3‑deoxy‑D‑gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,329,1070TC19.6% 103.3 / 20.4 51intergenic (‑24/+125)yiiM/kdgTconserved hypothetical protein/2‑keto‑3‑deoxy‑D‑gluconate transporter
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/8);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

AAGCCTGAATCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCCAGGCAGGAACTC  >  minE/2329041‑2329175
                                                                  |                                                                    
aaGCCTGAATCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTccc                                                                   >  1:3637418/1‑70 (MQ=255)
   ccTGAATCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCAtt                                                               <  1:3687289/71‑1 (MQ=255)
    cTGAATCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACCAGAAAttt                                                                      <  1:3946676/63‑1 (MQ=255)
       aaTCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGAc                                                           >  1:1434257/1‑71 (MQ=255)
             gCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCat                                                        <  1:2512501/68‑1 (MQ=255)
             gCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTa                                                              >  1:260438/1‑62 (MQ=255)
              ccTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCAtattt                                                    <  1:3384385/71‑1 (MQ=255)
               cTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCAtattta                                                   <  1:2804273/71‑1 (MQ=255)
                 gtgtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTccc                                                                   <  1:631761/53‑1 (MQ=255)
                  tgtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATatttattt                                                <  1:1924296/71‑1 (MQ=255)
                  tgtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCATatttattt                                                <  1:3513773/71‑1 (MQ=255)
                  tgtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCATatttattt                                                <  1:2849653/71‑1 (MQ=255)
                  tgtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCATatttattt                                                <  1:936263/71‑1 (MQ=255)
                   gtATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATatttattt                                                >  1:1302484/1‑70 (MQ=255)
                     ataCATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTccc                                                                    >  1:1783400/1‑48 (MQ=255)
                      taCATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATAtttattttat                                             <  1:1627835/70‑1 (MQ=255)
                      taCATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATg                                            >  1:3679198/1‑71 (MQ=255)
                        cATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGt                                           <  1:3151320/70‑1 (MQ=255)
                        cATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGt                                           <  1:3015439/70‑1 (MQ=255)
                         aTCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGt                                           >  1:1482082/1‑69 (MQ=255)
                         aTCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTc                                          >  1:2302332/1‑70 (MQ=255)
                            aaCCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGc                                      >  1:3049056/1‑71 (MQ=255)
                              ccGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCAtatt                                                     <  1:2163599/54‑1 (MQ=255)
                              ccGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCAtatt                                                     <  1:3498402/54‑1 (MQ=255)
                              ccGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCAtatt                                                     <  1:997621/54‑1 (MQ=255)
                                ggATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGc                                  <  1:2533819/71‑1 (MQ=255)
                                  atatCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTg                                >  1:518141/1‑71 (MQ=255)
                                   tatCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGt                               >  1:3114685/1‑71 (MQ=255)
                                   tatCGCATCTGATCCCTCTCCACGAGAAATTTCCCCATTAGACCAtatt                                                     <  1:1205477/49‑1 (MQ=255)
                                     tCGCATCTGTTCCCTCTCCACGAGAAATTCCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGtt                              <  1:943311/70‑1 (MQ=255)
                                       gCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTc                           <  1:1729129/71‑1 (MQ=255)
                                         aTCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTcgc                         >  1:3896046/1‑71 (MQ=255)
                                          tCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTcgcg                        >  1:2531009/1‑71 (MQ=255)
                                           cTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATAtttattttat                                             <  1:2865118/49‑1 (MQ=255)
                                            tGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGGGTGCTGTTTTcgcg                        >  1:1861221/1‑69 (MQ=255)
                                            tGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCgg                       >  1:3933049/1‑70 (MQ=255)
                                               tCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCgg                       >  1:1707359/1‑67 (MQ=255)
                                               tCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAg                   >  1:2808028/1‑71 (MQ=255)
                                                cccTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTc                                          <  1:3410916/47‑1 (MQ=255)
                                                cccTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCa                                         >  1:1218090/1‑48 (MQ=255)
                                                 ccTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCCGGc                                      <  1:2936218/50‑1 (MQ=255)
                                                    ctcCACGAGAAATTTCCCCATTAGACCAtatttat                                                  >  1:2363141/1‑35 (MQ=255)
                                                    ctcCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATc               >  1:503038/1‑70 (MQ=255)
                                                     tcCACGAGAAATTCCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATc               <  1:342920/69‑1 (MQ=255)
                                                      ccACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTcg                          <  1:2791751/57‑1 (MQ=255)
                                                       cACGAGAAATTTCCCCATTAGACCATAtttatttta                                              >  1:2988843/1‑36 (MQ=255)
                                                       cACGAGAAATTCCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcag            <  1:2321988/70‑1 (MQ=255)
                                                          gagaAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcag        >  1:1449888/1‑71 (MQ=255)
                                                          gagaAATTCCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcag        >  1:82456/1‑71 (MQ=255)
                                                            gaAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTAt                <  1:3390562/61‑1 (MQ=255)
                                                            gaAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTAt                <  1:769391/61‑1 (MQ=255)
                                                            gaAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcagg       >  1:2427541/1‑70 (MQ=255)
                                                               aTTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCCAGGCAGGAACt   <  1:994103/71‑1 (MQ=255)
                                                                tttCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCCAGGCAGGAACTc  <  1:1628931/71‑1 (MQ=255)
                                                                 ttCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcagg       >  1:3600548/1‑65 (MQ=255)
                                                                 ttCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcagg       >  1:80125/1‑65 (MQ=255)
                                                                 ttCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCcaggcagg       >  1:2082868/1‑65 (MQ=255)
                                                                  |                                                                    
AAGCCTGAATCTTGCCTGTGTATACATCAACCGGATATCGCATCTGTTCCCTCTCCACGAGAAATTTCCCCATTAGACCATATTTATTTTATGTCAGGCGTGCTGTTTTCGCGGTAAGTATCCAGGCAGGAACTC  >  minE/2329041‑2329175

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: