Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,996 A→T 41.5% intergenic (‑88/‑18) 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,9960AT41.5% 56.7 / 118.7 94intergenic (‑88/‑18)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (21/18);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAAT  >  minE/217928‑218059
                                                                    |                                                               
ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                                                <  1:3776470/70‑1 (MQ=255)
  gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg                                                             >  1:326242/1‑71 (MQ=255)
   aTTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGt                                                            >  1:512236/1‑71 (MQ=255)
    tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                                                           >  1:2025755/1‑71 (MQ=255)
    tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                                                           <  1:1901136/71‑1 (MQ=255)
     ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCt                                                          >  1:496185/1‑71 (MQ=255)
     ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt                                                          >  1:215941/1‑71 (MQ=255)
     ttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                                          >  1:2440319/1‑71 (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)
                 aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                               <  1:625707/70‑1 (MQ=255)
                    ttCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                                             <  1:198313/69‑1 (MQ=39)
                    ttCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGt                                                 <  1:2250256/65‑1 (MQ=255)
                      cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                                           >  1:1405483/1‑69 (MQ=255)
                      cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                                           >  1:3391866/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:46986/71‑1 (MQ=255)
                           caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                    <  1:742332/71‑1 (MQ=255)
                            aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                    >  1:2079355/1‑70 (MQ=255)
                            aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                    >  1:3847868/1‑70 (MQ=255)
                            aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                    <  1:3563126/70‑1 (MQ=255)
                             gCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                                  <  1:770845/71‑1 (MQ=255)
                               gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                                     <  1:1590271/49‑1 (MQ=255)
                               gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                                     <  1:599508/49‑1 (MQ=255)
                                 tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg                                                  >  1:1520299/1‑51 (MQ=255)
                                 tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg                                                  >  1:4095779/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:131304/1‑48 (MQ=255)
                                   cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTg                                                   >  1:1762635/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)
                                                   tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAATATaaga                     <  1:77/62‑1 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                                             <  1:3415811/38‑1 (MQ=255)
                                                    ggAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACg           >  1:3494015/1‑71 (MQ=255)
                                                     gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaatg           >  1:1985133/1‑68 (MQ=255)
                                                       gaATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                                         <  1:3401481/38‑1 (MQ=255)
                                                       gaATAAGCTGCACAGGCGTCTGGACGTGGTaatga                                            <  1:3400197/35‑1 (MQ=255)
                                                       gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag        <  1:1101612/71‑1 (MQ=255)
                                                        aaTAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                                         <  1:827227/37‑1 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTaatgaa                                           <  1:4109538/35‑1 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                         <  1:2377378/53‑1 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag        >  1:4056465/1‑70 (MQ=255)
                                                         aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACt                                <  1:3445503/45‑1 (MQ=255)
                                                         aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAAc            <  1:3589637/65‑1 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaa                       >  1:233314/1‑53 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag        >  1:484529/1‑68 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag        >  1:1863877/1‑68 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag     <  1:327962/71‑1 (MQ=255)
                                                           aaGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag     >  1:3768469/1‑70 (MQ=255)
                                                              cTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAat  >  1:2275791/1‑70 (MQ=255)
                                                              cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag     <  1:247099/67‑1 (MQ=255)
                                                                 cactGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaaa             >  1:3968423/1‑56 (MQ=255)
                                                                 cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGa    >  1:970008/1‑65 (MQ=255)
                                                                    |                                                               
CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAAT  >  minE/217928‑218059

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: