Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,947 G→A 12.9% intergenic (‑39/‑67) 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,9470GA12.9% 239.8 / 21.1 85intergenic (‑39/‑67)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (6/5);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAA  >  minE/217877‑218014
                                                                      |                                                                   
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     <  1:4112173/71‑1 (MQ=255)
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     >  1:273513/1‑71 (MQ=255)
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     <  1:3752923/71‑1 (MQ=255)
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     <  1:2941092/71‑1 (MQ=255)
  aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                                   >  1:421166/1‑71 (MQ=255)
    aaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCg                                                                 >  1:4090969/1‑71 (MQ=255)
     aTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                                                                <  1:812661/71‑1 (MQ=255)
          ccAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAg                                                           >  1:1662292/1‑71 (MQ=255)
             gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                        <  1:3897863/71‑1 (MQ=255)
               tctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     >  1:3820880/1‑56 (MQ=255)
                ctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGcc                                                               >  1:2854829/1‑61 (MQ=255)
                ctcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                     >  1:773493/1‑71 (MQ=255)
                   gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   >  1:988750/1‑70 (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)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                              >  1:573681/1‑53 (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:988695/1‑71 (MQ=255)
                             cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                        >  1:388499/1‑71 (MQ=255)
                              ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                                   >  1:933284/1‑43 (MQ=255)
                              ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGc                                         <  1:1423329/69‑1 (MQ=255)
                              ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                        <  1:2323429/70‑1 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCg                                                                 <  1:1708743/44‑1 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                            <  1:2302957/65‑1 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                      <  1:2386851/71‑1 (MQ=255)
                                 ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     <  1:166028/38‑1 (MQ=255)
                                 ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                              >  1:3191486/1‑61 (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)
                                       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:4005408/39‑1 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAAt                                                       <  1:2432146/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:2921357/1‑42 (MQ=255)
                                                      aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                            >  1:3464803/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:2411306/71‑1 (MQ=255)
                                                                 gagaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGt    <  1:2527608/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)
                                                                      |                                                                   
AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAA  >  minE/217877‑218014

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: