Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,968 C→T 16.5% intergenic (‑60/‑46) 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,9680CT16.5% 261.5 / 37.4 103intergenic (‑60/‑46)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/40);  new base T (6/11);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 107 total aligned reads displayed.

TGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACA  >  minE/217900‑218032
                                                                    |                                                                
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)
      cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                                          >  1:388499/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)
          ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                      <  1:2953796/71‑1 (MQ=255)
            aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                      >  1:3918264/1‑69 (MQ=255)
            aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                      >  1:10688/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)
                       ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                               <  1:1022923/49‑1 (MQ=255)
                       ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                               <  1:3905891/49‑1 (MQ=255)
                       ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                               <  1:3856120/49‑1 (MQ=255)
                       ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                          <  1:1339770/70‑1 (MQ=255)
                       ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                          <  1:2774097/70‑1 (MQ=255)
                         ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                             >  1:784773/1‑49 (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)
                               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)
                                 ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGca                                        >  1:1873956/1‑62 (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)
                                                         gCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga       <  1:770845/71‑1 (MQ=255)
                                                           gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                          <  1:1590271/49‑1 (MQ=255)
                                                           gaTCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                          <  1:599508/49‑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:131304/1‑48 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTg                        >  1:1762635/1‑48 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                    >  1:9274/1‑52 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa  >  1:2686778/1‑70 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa  >  1:3864236/1‑70 (MQ=255)
                                                                 aaaTCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg              <  1:2334568/56‑1 (MQ=255)
                                                                   aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaat                   <  1:3152343/49‑1 (MQ=255)
                                                                    |                                                                
TGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACA  >  minE/217900‑218032

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: