Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,959 A→G 13.8% intergenic (‑51/‑55) 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,9590AG13.8% 263.8 / 30.6 101intergenic (‑51/‑55)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/46);  new base G (7/7);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGA  >  minE/217890‑218027
                                                                     |                                                                    
gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                                     <  1:3897863/71‑1 (MQ=255)
   ctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                                  >  1:773493/1‑71 (MQ=255)
      gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                                <  1:1239906/70‑1 (MQ=255)
       ctctGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                                              >  1:1249923/1‑71 (MQ=255)
          tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                             <  1:3178897/69‑1 (MQ=255)
           gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAcc                                                            <  1:3476126/69‑1 (MQ=255)
            aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                           <  1:3991144/69‑1 (MQ=255)
              cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                                       >  1:2931494/1‑71 (MQ=255)
               tCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                             <  1:394632/64‑1 (MQ=255)
                cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                                              >  1:2664722/1‑62 (MQ=255)
                cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                                     >  1:388499/1‑71 (MQ=255)
                cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                                     >  1:988695/1‑71 (MQ=255)
                 ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGc                                                      <  1:1423329/69‑1 (MQ=255)
                 ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                                     <  1:2323429/70‑1 (MQ=255)
                  taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                         <  1:2302957/65‑1 (MQ=255)
                  taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                   <  1:2386851/71‑1 (MQ=255)
                    ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                           >  1:3191486/1‑61 (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)
                          tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                                <  1:558680/50‑1 (MQ=255)
                          tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                            >  1:3966259/1‑70 (MQ=255)
                              tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATaa                                         <  1:975312/69‑1 (MQ=255)
                              tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                          <  1:245991/68‑1 (MQ=255)
                              cagTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa                                                              >  1:849788/3‑48 (MQ=255)
                               agTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                                       <  1:134594/69‑1 (MQ=255)
                                 ttCTGCTGATTTTACCGCAGAGAAGTTCGCCACAGCAAt                                                                    <  1:1140689/39‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAAt                                                                    <  1:2432146/39‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAAt                                                                    <  1:4005408/39‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:3905891/49‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:1022923/49‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:3856120/49‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                     <  1:2774097/70‑1 (MQ=255)
                                 ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                     <  1:1339770/70‑1 (MQ=255)
                                   ctgctgATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                                   >  1:118034/1‑38 (MQ=25)
                                   ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                        >  1:784773/1‑49 (MQ=255)
                                   ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAg                                        >  1:18459/1‑65 (MQ=255)
                                    tgctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaat                                                             >  1:3563026/1‑42 (MQ=25)
                                     gctgATTTAAGCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                 >  1:2471184/1‑54 (MQ=255)
                                      ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcact                                 <  1:4125316/69‑2 (MQ=255)
                                      ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                                <  1:2732487/38‑1 (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)
                                         aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAcc                                                            <  1:1639405/39‑1 (MQ=255)
                                         aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          >  1:618077/1‑41 (MQ=255)
                                         aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                         >  1:3464803/1‑42 (MQ=255)
                                         aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                         >  1:2921357/1‑42 (MQ=255)
                                          tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                           <  1:1901136/71‑1 (MQ=255)
                                          tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTc                           >  1:2025755/1‑71 (MQ=255)
                                          tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                 <  1:1959724/49‑1 (MQ=255)
                                           ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          >  1:2563916/1‑39 (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:184842/1‑60 (MQ=255)
                                                     agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                           <  1:3086913/60‑1 (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)
                                                                     |                                                                    
GATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGA  >  minE/217890‑218027

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: