Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,959 A→G 14.1% 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,9590AG14.1% 263.5 / 29.8 106intergenic (‑51/‑55)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/54);  new base G (8/7);  total (45/61)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGA  >  minE/217889‑218027
                                                                      |                                                                    
aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCa                                                                      <  1:1364410/71‑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)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                           <  1:2888069/69‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:2663366/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:3761061/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                          <  1:424125/70‑1 (MQ=255)
             aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                         >  1:724093/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)
                     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)
                                ggTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                            <  1:550551/65‑1 (MQ=255)
                                agTTCTGCTGATTTAACCGCTGAGAAATTCGCCAGAGCGATCCGAAATCACTGGCACg                                                   <  1:3873982/57‑1 (MQ=255)
                                agTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAtgtg                                                 >  1:1451689/2‑60 (MQ=255)
                                  ttCTGCTGATTTAACCGCAGAGAATTTCGCCACAGCaa                                                                     <  1:778613/38‑1 (MQ=255)
                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                                <  1:1956657/43‑1 (MQ=255)
                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                             >  1:2779023/1‑46 (MQ=255)
                                  ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGc                                    <  1:3013101/71‑1 (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)
                                        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:3661896/1‑70 (MQ=255)
                                              aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                         >  1:3740812/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:2536068/1‑71 (MQ=255)
                                                ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                      >  1:1126576/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)
                                                   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)
                                                     gaaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 <  1:1641344/68‑1 (MQ=255)
                                                      agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                  >  1:2536499/1‑69 (MQ=255)
                                                      aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                     <  1:294466/48‑1 (MQ=255)
                                                      aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                     <  1:2232146/48‑1 (MQ=255)
                                                      aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 >  1:1028133/3‑70 (MQ=255)
                                                      aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 >  1:570190/3‑70 (MQ=255)
                                                         aattCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg             <  1:1675544/69‑1 (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)
                                                          attCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg             <  1:3447269/69‑1 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg                   >  1:3804863/1‑61 (MQ=38)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat          >  1:752732/1‑70 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                >  1:1765062/1‑48 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                     <  1:3808671/59‑1 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                     <  1:1816668/59‑1 (MQ=255)
                                                             cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAACCTGCACAGGCGTCTGGACGTGGTaatgaat          >  1:3490900/1‑70 (MQ=255)
                                                              gCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg         >  1:1425715/1‑70 (MQ=255)
                                                              gCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAatgaatga        <  1:522549/71‑1 (MQ=255)
                                                               ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 <  1:607213/61‑1 (MQ=38)
                                                               ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 <  1:2608476/61‑1 (MQ=38)
                                                               ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg         <  1:3806731/69‑1 (MQ=39)
                                                               ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat          <  1:618958/68‑1 (MQ=39)
                                                               ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa       <  1:3257089/71‑1 (MQ=255)
                                                                cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg      <  1:3904199/71‑1 (MQ=255)
                                                                cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg      <  1:1179097/71‑1 (MQ=255)
                                                                cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa                       >  1:876085/1‑54 (MQ=255)
                                                                 acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGCAGTGGTAATGAATGAAga     >  1:3267468/1‑71 (MQ=255)
                                                                 acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                 <  1:2976693/59‑1 (MQ=255)
                                                                  caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg                               >  1:281643/1‑44 (MQ=38)
                                                                  caGCAATCGGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac    >  1:392748/1‑71 (MQ=255)
                                                                  caGCAATCCGAAACCCCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac    >  1:1126792/1‑71 (MQ=255)
                                                                  caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                      <  1:252004/53‑1 (MQ=38)
                                                                  caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac    >  1:794802/1‑71 (MQ=39)
                                                                  caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat          >  1:337863/1‑65 (MQ=255)
                                                                  caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa               <  1:1757602/60‑1 (MQ=255)
                                                                  caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa               <  1:329316/60‑1 (MQ=255)
                                                                   aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac    >  1:1256877/1‑70 (MQ=255)
                                                                     cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga  <  1:569411/70‑1 (MQ=255)
                                                                     cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga  <  1:341886/70‑1 (MQ=255)
                                                                     cAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTa                >  1:945283/1‑56 (MQ=255)
                                                                      |                                                                    
AGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGA  >  minE/217889‑218027

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: