Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,947 G→A 11.8% 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,9470GA11.8% 228.4 / 17.9 85intergenic (‑39/‑67)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/46);  new base A (5/5);  total (33/52)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATG  >  minE/217877‑218016
                                                                      |                                                                     
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                       >  1:2097651/1‑71 (MQ=255)
 gAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGt                                                                      <  1:3020782/71‑1 (MQ=255)
  aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                                     >  1:1594437/1‑71 (MQ=255)
   aaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTc                                                                    >  1:3441981/1‑71 (MQ=255)
      ttAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                                                                  >  1:3787525/1‑70 (MQ=255)
        aaCCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                                <  1:2568774/70‑1 (MQ=255)
           cAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                            <  1:1350411/71‑1 (MQ=255)
           cAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                            <  1:530907/71‑1 (MQ=255)
            aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                            >  1:1300896/1‑70 (MQ=255)
            aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCa                                                           <  1:1364410/71‑1 (MQ=255)
                 tcGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                            <  1:1332634/65‑1 (MQ=255)
                   gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATc                                                        <  1:1384266/67‑1 (MQ=255)
                   gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:498420/70‑1 (MQ=255)
                    ctctGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:1755118/69‑1 (MQ=255)
                       tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                          <  1:717035/61‑1 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCg                                                      >  1:3522061/1‑63 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                <  1:2888069/69‑1 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               <  1:2663366/70‑1 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               <  1:424125/70‑1 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               <  1:3761061/70‑1 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                              >  1:724093/1‑71 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                              >  1:2458313/1‑71 (MQ=255)
                           cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                            <  1:1390659/71‑1 (MQ=255)
                           cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                            <  1:2077113/71‑1 (MQ=255)
                             cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               <  1:245371/66‑1 (MQ=255)
                                 ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGt                                       <  1:2888312/70‑1 (MQ=255)
                                  taaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                     <  1:2942122/71‑1 (MQ=255)
                                     aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:2464969/52‑1 (MQ=255)
                                          tcagTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                             >  1:1686334/4‑71 (MQ=255)
                                           tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                               <  1:2987154/68‑1 (MQ=255)
                                           cagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgag                                   <  1:1418205/62‑1 (MQ=255)
                                            ggTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                 <  1:550551/65‑1 (MQ=255)
                                            agTTCTGCTGATTTAACCGCTGAGAAATTCGCCAGAGCGATCCGAAATCACTGGCACg                                        <  1:3873982/57‑1 (MQ=255)
                                            agTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAtgtg                                      >  1:1451689/2‑60 (MQ=255)
                                            agTTCTGCTGATTTAACCGCAGAAAAGTTCGCcac                                                               <  1:3461769/34‑1 (MQ=38)
                                              ttCTGCTGATTTAACCGCAGAGAATTTCGCCACAGCaa                                                          <  1:778613/38‑1 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                            <  1:2991690/36‑1 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:1956657/43‑1 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                  >  1:2779023/1‑46 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGc                         <  1:3013101/71‑1 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAAAAGTTCGCCACAg                                                             >  1:1578345/1‑35 (MQ=255)
                                              ttCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTgg                                            <  1:3736631/52‑1 (MQ=255)
                                               tCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGca                        >  1:2583378/1‑71 (MQ=255)
                                                ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:3199143/41‑1 (MQ=255)
                                                ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                       >  1:1189466/1‑71 (MQ=255)
                                                  gctgTTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGc                                           <  1:2825032/49‑1 (MQ=255)
                                                   ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcactg                     <  1:918029/70‑3 (MQ=255)
                                                   ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg                    >  1:2072619/1‑71 (MQ=255)
                                                   ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                     <  1:2612697/70‑3 (MQ=255)
                                                   ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                     <  1:696136/70‑3 (MQ=255)
                                                    tgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                   >  1:2602757/1‑71 (MQ=255)
                                                    tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                >  1:920853/1‑42 (MQ=255)
                                                     gATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                                  <  1:1378176/55‑1 (MQ=255)
                                                       tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                   <  1:410811/68‑1 (MQ=255)
                                                       tttAACCGCAGAGAAGTTCGCCACAGCAATCTGAAACCACTGGCATGTGGAGAATAAGCTg                          >  1:464516/1‑61 (MQ=255)
                                                       tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTc                <  1:3655820/71‑1 (MQ=255)
                                                        ttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt               >  1:964310/1‑71 (MQ=255)
                                                         tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGga                                    <  1:2104179/49‑1 (MQ=255)
                                                          aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg              >  1:3740812/1‑70 (MQ=255)
                                                          aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg              >  1:3661896/1‑70 (MQ=255)
                                                           aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAc                                         >  1:2507046/1‑42 (MQ=255)
                                                           aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa            <  1:1274479/71‑1 (MQ=255)
                                                           aCCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                     >  1:2422969/1‑46 (MQ=255)
                                                            ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc           >  1:1126576/1‑71 (MQ=255)
                                                            ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc           <  1:2699247/71‑1 (MQ=255)
                                                            ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc           >  1:2536068/1‑71 (MQ=255)
                                                            ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGga                                    >  1:1471115/1‑46 (MQ=255)
                                                            ccGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGAc           <  1:906744/71‑1 (MQ=255)
                                                             cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACg          <  1:3368884/71‑1 (MQ=255)
                                                              gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                          <  1:2992773/54‑1 (MQ=255)
                                                               ctgagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAAt                                >  1:2266602/3‑47 (MQ=37)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTgg                                     <  1:1674402/41‑1 (MQ=255)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg       >  1:473284/1‑71 (MQ=255)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                        <  1:1954758/38‑1 (MQ=255)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        >  1:989865/1‑54 (MQ=255)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        >  1:2096358/1‑54 (MQ=255)
                                                                 gagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTCCACTGGCGTCTGGACGTgg       >  1:3931825/1‑70 (MQ=255)
                                                                 gaaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      <  1:1641344/68‑1 (MQ=255)
                                                                  agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg       >  1:2536499/1‑69 (MQ=255)
                                                                  aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                          <  1:2232146/48‑1 (MQ=255)
                                                                  aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                          <  1:294466/48‑1 (MQ=255)
                                                                  aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      >  1:570190/3‑70 (MQ=255)
                                                                  aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      >  1:1028133/3‑70 (MQ=255)
                                                                     aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg  <  1:1511678/71‑1 (MQ=255)
                                                                     aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg  <  1:3521627/71‑1 (MQ=255)
                                                                     aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg  <  1:2857703/71‑1 (MQ=255)
                                                                      gTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg  >  1:124471/1‑70 (MQ=39)
                                                                      gTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg  >  1:2895760/1‑70 (MQ=255)
                                                                      |                                                                     
AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATG  >  minE/217877‑218016

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: