Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,977 C→T 8.3% intergenic (‑69/‑37) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE217,9770CT 8.3% 302.8 / 12.3 96intergenic (‑69/‑37)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base T (6/2);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGA  >  minE/217908‑218046
                                                                     |                                                                     
taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                                      <  1:2386851/71‑1 (MQ=255)
  ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                                    <  1:2953796/71‑1 (MQ=255)
    aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                                    >  1:10688/1‑69 (MQ=255)
    aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                                    >  1:3918264/1‑69 (MQ=255)
      aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgaga                                                                <  1:2564520/71‑1 (MQ=255)
       ttCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAc                                                                       <  1:1307870/63‑1 (MQ=255)
        tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                                               >  1:3966259/1‑70 (MQ=255)
            tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATaa                                                            <  1:975312/69‑1 (MQ=255)
            tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                                             <  1:245991/68‑1 (MQ=255)
             agTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                                                          <  1:134594/69‑1 (MQ=255)
               ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                                        <  1:2774097/70‑1 (MQ=255)
               ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                                        <  1:1339770/70‑1 (MQ=255)
                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAg                                                           >  1:18459/1‑65 (MQ=255)
                   gctgATTTAAGCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                                    >  1:2471184/1‑54 (MQ=255)
                    ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcact                                                    <  1:4125316/69‑2 (MQ=255)
                    ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                                   <  1:3776470/70‑1 (MQ=255)
                      gATTTAACCGCAGAGAAGTTCGCCAGAGCAATCCGAAACCACTGGCACGt                                                                     <  1:1817019/50‑1 (MQ=255)
                      gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg                                                >  1:326242/1‑71 (MQ=255)
                       aTTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGt                                               >  1:512236/1‑71 (MQ=255)
                       aTTTAACCGCTGAGAAATCCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGc                                                       <  1:3992780/63‑1 (MQ=255)
                        tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                                              >  1:2025755/1‑71 (MQ=255)
                        tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                                              <  1:1901136/71‑1 (MQ=255)
                        tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                                    <  1:1959724/49‑1 (MQ=255)
                         ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGca                                                      >  1:1873956/1‑62 (MQ=255)
                         ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCt                                             >  1:496185/1‑71 (MQ=255)
                         ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt                                             >  1:215941/1‑71 (MQ=255)
                         ttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                             >  1:2440319/1‑71 (MQ=255)
                            aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGga                                                                  >  1:3381260/1‑47 (MQ=255)
                            aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                            <  1:1041706/69‑1 (MQ=255)
                            aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa                                          <  1:1979510/71‑1 (MQ=255)
                            aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTCCACAGGCGTCTgg                                           <  1:1929880/70‑1 (MQ=255)
                                cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg                                      >  1:3055515/1‑71 (MQ=255)
                                  gagaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGt                                    <  1:2527608/71‑1 (MQ=255)
                                  gagaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGt                                    <  1:2411306/71‑1 (MQ=255)
                                   agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                                              <  1:3086913/60‑1 (MQ=255)
                                   agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                                              >  1:184842/1‑60 (MQ=255)
                                   agaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTa                                   <  1:1513302/71‑1 (MQ=255)
                                    gaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                                  >  1:1978075/1‑71 (MQ=255)
                                    gaAATTCGCCACAGCGATACGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                                  >  1:4088144/1‑71 (MQ=255)
                                    aaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                        >  1:1383737/2‑65 (MQ=255)
                                     aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCt                                                         >  1:1863280/1‑47 (MQ=255)
                                     aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                  <  1:625707/70‑1 (MQ=255)
                                     aaGTTCACCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                                                          >  1:534716/1‑46 (MQ=255)
                                        ttCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                                <  1:198313/69‑1 (MQ=39)
                                        ttCGCCACAGCAATCCGAAACCACTGGCACGTGgag                                                                 <  1:3291827/36‑1 (MQ=255)
                                        ttCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGt                                    <  1:2250256/65‑1 (MQ=255)
                                          cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                              >  1:3391866/1‑69 (MQ=255)
                                          cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                              >  1:1405483/1‑69 (MQ=255)
                                           gCCACAGCGATCTGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                            <  1:3161061/70‑1 (MQ=255)
                                           gCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                            <  1:3823459/70‑1 (MQ=39)
                                           gCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                            <  1:1720862/70‑1 (MQ=255)
                                           gCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaat                                 <  1:3222011/65‑1 (MQ=255)
                                           gCCACAGAAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                            <  1:425924/70‑1 (MQ=255)
                                            ccACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAAtgaatgaa                          <  1:1327535/71‑1 (MQ=255)
                                             cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                                              >  1:3019099/1‑50 (MQ=255)
                                             cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                            >  1:2275848/1‑68 (MQ=255)
                                             cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                  >  1:583595/1‑62 (MQ=255)
                                              acaGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAga                        >  1:1267931/1‑71 (MQ=39)
                                              acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAga                        >  1:745010/1‑71 (MQ=255)
                                              acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                        >  1:618963/1‑71 (MQ=255)
                                               caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                            >  1:553838/1‑50 (MQ=255)
                                               caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatga                               >  1:3796999/1‑63 (MQ=255)
                                               caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                       <  1:742332/71‑1 (MQ=255)
                                               caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                       <  1:46986/71‑1 (MQ=255)
                                                aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                       >  1:2079355/1‑70 (MQ=255)
                                                aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                       <  1:3563126/70‑1 (MQ=255)
                                                aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                       >  1:3847868/1‑70 (MQ=255)
                                                 gCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                     <  1:770845/71‑1 (MQ=255)
                                                   gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                        <  1:1590271/49‑1 (MQ=255)
                                                   gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                        <  1:599508/49‑1 (MQ=255)
                                                    aTCCGAAACCACTGGCACGTGGAGAATAAGCTGca                                                      <  1:2053383/35‑1 (MQ=255)
                                                     tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg                                     >  1:4095779/1‑51 (MQ=255)
                                                     tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg                                     >  1:1520299/1‑51 (MQ=255)
                                                     tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                  >  1:2276241/1‑54 (MQ=255)
                                                      ccGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                    <  1:3402317/51‑1 (MQ=255)
                                                      ccGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                                  >  1:2670122/1‑53 (MQ=255)
                                                      ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                                 <  1:3433863/38‑1 (MQ=255)
                                                      ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                                 <  1:3298503/38‑1 (MQ=255)
                                                       cGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAatgaatga                           >  1:3173880/1‑59 (MQ=38)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                            >  1:261192/1‑42 (MQ=255)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTg                                      >  1:1762635/1‑48 (MQ=255)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTg                                      >  1:131304/1‑48 (MQ=255)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                  >  1:9274/1‑52 (MQ=255)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa                >  1:3864236/1‑70 (MQ=255)
                                                       cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa                >  1:2686778/1‑70 (MQ=255)
                                                         aaaTCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                            <  1:2334568/56‑1 (MQ=255)
                                                           aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaat                                 <  1:3152343/49‑1 (MQ=255)
                                                             cACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgacga                     >  1:1360817/1‑59 (MQ=255)
                                                             cACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt            >  1:3404292/1‑68 (MQ=255)
                                                              aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgacga                     >  1:3222757/1‑58 (MQ=255)
                                                               cTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg                         >  1:1464668/1‑53 (MQ=255)
                                                                tGGCACGTGGAGAATAAGCTGCACAGTCGTCTGGACGTGGTaatg                                <  1:211378/45‑1 (MQ=255)
                                                                  gCACGTGGAGAATAAGCTGCACTGGCGTCTGGAGGTGGTAATTAATGAAgacga                     <  1:3326906/54‑1 (MQ=255)
                                                                  gCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaag      <  1:427443/69‑1 (MQ=255)
                                                                  gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                                        <  1:1911939/35‑1 (MQ=255)
                                                                   cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                              >  1:1340826/1‑44 (MQ=255)
                                                                    aCGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa  <  1:3254373/71‑1 (MQ=255)
                                                                     cGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa  >  1:868611/1‑70 (MQ=255)
                                                                     |                                                                     
TAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGA  >  minE/217908‑218046

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: