Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,944 G→A 14.3% intergenic (‑36/‑70) 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,9440GA14.3% 175.8 / 18.4 70intergenic (‑36/‑70)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/38);  new base A (4/6);  total (26/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

AAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGG  >  minE/217876‑218011
                                                                    |                                                                   
aaGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                                                                     >  1:3516243/1‑69 (MQ=255)
 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:424125/70‑1 (MQ=255)
                          aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                          <  1:2663366/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:2077113/71‑1 (MQ=255)
                            cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                       <  1:1390659/71‑1 (MQ=255)
                              cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                          <  1:245371/66‑1 (MQ=255)
                                taataaATTCTATGGTTCTGCTGATTTAACCGCagag                                                                     <  1:3390617/37‑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:2536068/1‑71 (MQ=255)
                                                             ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc      <  1:2699247/71‑1 (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:2096358/1‑54 (MQ=255)
                                                                 agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                   >  1:989865/1‑54 (MQ=255)
                                                                  gagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTCCACTGGCGTCTGGACGTgg  >  1:3931825/1‑70 (MQ=255)
                                                                   agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg  >  1:2536499/1‑69 (MQ=255)
                                                                    |                                                                   
AAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGG  >  minE/217876‑218011

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: