Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,944 G→A 12.2% 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,9440GA12.2% 274.5 / 22.9 106intergenic (‑36/‑70)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/51);  new base A (7/6);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 122 total aligned reads displayed.

AAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAA  >  minE/217876‑218014
                                                                    |                                                                      
aaGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                                                                        >  1:4435033/1‑69 (MQ=255)
aaGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                                                                        >  1:3136085/1‑69 (MQ=255)
 aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                     >  1:1915426/1‑71 (MQ=255)
  gAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGt                                                                    >  1:1247148/1‑71 (MQ=255)
   aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                                   <  1:636679/71‑1 (MQ=255)
            cAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                          <  1:292028/71‑1 (MQ=255)
             aGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGc                                                          >  1:818974/1‑70 (MQ=255)
              gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCTCAGAGAAGTTCGCcac                                                             <  1:4690056/66‑1 (MQ=255)
              gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                        <  1:4370313/71‑1 (MQ=255)
              gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa                                                        <  1:264792/71‑1 (MQ=255)
                    gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                     <  1:4332326/68‑1 (MQ=255)
                    gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc                                                     <  1:3681812/68‑1 (MQ=255)
                    gCTCCGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:4354214/70‑1 (MQ=255)
                         gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAGGTTCGCCACAGCAATCCGAAACCa                                              >  1:3018440/1‑70 (MQ=255)
                         gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                             >  1:2903797/1‑71 (MQ=255)
                          aaCTAATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                              <  1:1460499/69‑1 (MQ=255)
                               ataataAATTCTATGGTTCTGCTGATTTAACCGCAgaga                                                                       >  1:5627028/1‑39 (MQ=255)
                               ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                      >  1:2025084/1‑40 (MQ=255)
                               ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                              >  1:1599342/1‑48 (MQ=255)
                               ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                              >  1:3649422/1‑48 (MQ=255)
                               ataataAATTCTATAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                              >  1:1536852/1‑48 (MQ=255)
                                taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGt                                                                    <  1:3306429/41‑1 (MQ=255)
                                  ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                    >  1:1605096/1‑71 (MQ=255)
                                   taaATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                      >  1:4129252/1‑36 (MQ=255)
                                   taaATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                                                                      >  1:3261225/1‑36 (MQ=255)
                                    aaaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                              >  1:1543275/1‑43 (MQ=255)
                                     aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTc                                                                  >  1:5571926/1‑38 (MQ=255)
                                     aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTc                                                                  >  1:4632689/1‑38 (MQ=255)
                                     aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                      <  1:4754874/66‑1 (MQ=255)
                                      aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:3646884/52‑1 (MQ=255)
                                      aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:1485011/52‑1 (MQ=255)
                                      aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:486304/52‑1 (MQ=255)
                                        tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                               >  1:1408437/1‑70 (MQ=255)
                                        tCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                               >  1:5215661/1‑70 (MQ=255)
                                           tcagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGATAATAAg                           <  1:914249/68‑1 (MQ=255)
                                           tcagTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAg                           <  1:701778/68‑1 (MQ=255)
                                           tcagTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                           >  1:1716862/4‑71 (MQ=255)
                                            tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                             <  1:2511962/68‑1 (MQ=255)
                                            cagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGa                                                   <  1:1747627/44‑1 (MQ=37)
                                            cagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                          >  1:1800812/3‑71 (MQ=255)
                                             agTTCTGCTGATTTAACCGCAGAGAAGATCGCCACAGCAATCCGAACCCACTGGCATGTCGAGAATAAGc                          >  1:2274820/2‑70 (MQ=255)
                                              gTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTg                        <  1:2634944/71‑1 (MQ=255)
                                              gTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCt                         >  1:3081042/1‑70 (MQ=255)
                                              gTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                        <  1:1500413/71‑1 (MQ=255)
                                              gTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                        >  1:5045809/1‑71 (MQ=255)
                                               ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                             <  1:2789026/49‑1 (MQ=255)
                                                tCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                          >  1:2127757/1‑67 (MQ=255)
                                                 ctgctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCAc                                             >  1:4101966/1‑47 (MQ=25)
                                                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                     >  1:116352/1‑71 (MQ=255)
                                                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                     >  1:346636/1‑71 (MQ=255)
                                                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                      <  1:1709404/70‑1 (MQ=255)
                                                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                      <  1:3278884/70‑1 (MQ=255)
                                                 ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                      <  1:2322516/70‑1 (MQ=255)
                                                 ctgctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcac                     >  1:775220/1‑71 (MQ=255)
                                                 ctgctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcac                     >  1:4592529/1‑71 (MQ=255)
                                                  tgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                        <  1:3364865/67‑1 (MQ=255)
                                                  tgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                        <  1:2292927/67‑1 (MQ=255)
                                                  tgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGCCATGTGGAGAATAAGCTg                        <  1:1386567/67‑1 (MQ=255)
                                                   gctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:1060539/39‑1 (MQ=255)
                                                   gctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGaa                               <  1:888866/59‑1 (MQ=255)
                                                    ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                          >  1:3438182/1‑47 (MQ=255)
                                                    ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTga                                          >  1:5483034/1‑46 (MQ=255)
                                                    ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                     <  1:2223114/68‑1 (MQ=255)
                                                    ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTAGCATGTGGAGAATAAGCTGCACTgg                  >  1:5356633/1‑71 (MQ=255)
                                                    ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCg                                                    <  1:5199602/37‑1 (MQ=255)
                                                    ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                   <  1:3693486/70‑3 (MQ=255)
                                                     tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                   <  1:2677556/37‑1 (MQ=255)
                                                      gATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcactg                   <  1:2472351/68‑3 (MQ=255)
                                                      gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCg                >  1:134022/1‑71 (MQ=255)
                                                       aTTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                           >  1:2922580/1‑59 (MQ=255)
                                                       aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                             >  1:1075564/1‑41 (MQ=255)
                                                       aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                          <  1:2494219/44‑1 (MQ=255)
                                                       aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGt               >  1:4721210/1‑71 (MQ=255)
                                                       aTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGc                 <  1:4158700/69‑1 (MQ=255)
                                                       aTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGt               >  1:902225/1‑71 (MQ=255)
                                                         ttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt             >  1:456631/1‑71 (MQ=255)
                                                         ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                      >  1:3839044/1‑62 (MQ=255)
                                                         ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt               >  1:3662270/1‑69 (MQ=255)
                                                          tAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt             <  1:1294206/70‑1 (MQ=255)
                                                          tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTg            >  1:2033247/1‑71 (MQ=255)
                                                          tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                          >  1:4052539/1‑57 (MQ=255)
                                                          tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                          >  1:3520978/1‑57 (MQ=255)
                                                          tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg            >  1:1009297/1‑71 (MQ=255)
                                                           aaCCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg            <  1:1689379/70‑1 (MQ=255)
                                                           aaCCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg            <  1:3691445/70‑1 (MQ=255)
                                                           aaCCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg            <  1:2028959/70‑1 (MQ=255)
                                                              cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg        <  1:4357549/71‑1 (MQ=255)
                                                              cGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgag                                 <  1:2625925/46‑1 (MQ=255)
                                                               gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt       <  1:4525289/71‑1 (MQ=255)
                                                               gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt       <  1:1187589/71‑1 (MQ=255)
                                                               gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt       <  1:320840/71‑1 (MQ=255)
                                                                ctgagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                           >  1:2269167/3‑50 (MQ=37)
                                                                ctgagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg      >  1:1464334/3‑71 (MQ=255)
                                                                cAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcactg                   <  1:5480294/58‑3 (MQ=255)
                                                                cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                      >  1:316088/1‑55 (MQ=255)
                                                                 agagAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                        <  1:4280936/52‑1 (MQ=255)
                                                                 agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg     <  1:5525763/71‑1 (MQ=255)
                                                                  gagaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt    <  1:2267324/71‑1 (MQ=255)
                                                                   agaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt    <  1:5014862/70‑1 (MQ=255)
                                                                    gaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa  >  1:2568315/1‑71 (MQ=255)
                                                                    |                                                                      
AAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAA  >  minE/217876‑218014

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: