Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,947 G→A 13.7% 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,9470GA13.7% 329.9 / 34.0 131intergenic (‑39/‑67)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (50/63);  new base A (7/11);  total (57/74)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 132 total aligned reads displayed.

AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATG  >  minE/217877‑218016
                                                                      |                                                                     
aGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAg                                                                       >  1:1915426/1‑71 (MQ=255)
  aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGtt                                                                     <  1:636679/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)
                        gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               >  1:2903797/1‑71 (MQ=255)
                         aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                               <  1:2669981/70‑1 (MQ=255)
                         aaCTAATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                <  1:1460499/69‑1 (MQ=255)
                              ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                                >  1:3649422/1‑48 (MQ=255)
                              ataataAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                                >  1:1599342/1‑48 (MQ=255)
                              ataataAATTCTATAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCca                                                                >  1:1536852/1‑48 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGt                                                                      <  1:2218769/41‑1 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGt                                                                      <  1:3306429/41‑1 (MQ=255)
                               taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                                                                  <  1:2927701/45‑1 (MQ=255)
                                 ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                      <  1:1512108/71‑1 (MQ=255)
                                 ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                      >  1:1605096/1‑71 (MQ=255)
                                    aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTc                                                                    >  1:4632689/1‑38 (MQ=255)
                                    aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTc                                                                    >  1:5571926/1‑38 (MQ=255)
                                    aaTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                        <  1:4754874/66‑1 (MQ=255)
                                     aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:486304/52‑1 (MQ=255)
                                     aTTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGa                                                     <  1:3646884/52‑1 (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)
                                           cagTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                             >  1:2283110/3‑54 (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)
                                             gTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATaa                              <  1:4777835/67‑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)
                                                ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        <  1:1709404/70‑1 (MQ=255)
                                                ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        <  1:2322516/70‑1 (MQ=255)
                                                ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        <  1:3278884/70‑1 (MQ=255)
                                                ctgctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcac                       >  1:4592529/1‑71 (MQ=255)
                                                ctgctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcac                       >  1:775220/1‑71 (MQ=255)
                                                 tgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                          <  1:2292927/67‑1 (MQ=255)
                                                 tgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTg                          <  1:3371385/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)
                                                   ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                >  1:5249360/1‑43 (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)
                                                     gATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcactg                     <  1:2472351/68‑3 (MQ=255)
                                                     gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCa                                          <  1:4359350/47‑1 (MQ=255)
                                                     gATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCg                  >  1:134022/1‑71 (MQ=255)
                                                      aTTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                             >  1:2922580/1‑59 (MQ=255)
                                                      aTTTAACCGCAGAAAAGTTCGCCACAGGAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGc                   <  1:2536454/69‑1 (MQ=255)
                                                      aTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGc                   <  1:4158700/69‑1 (MQ=255)
                                                      aTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGt                 >  1:902225/1‑71 (MQ=255)
                                                        ttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACg                                        >  1:1863310/1‑46 (MQ=25)
                                                        ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        >  1:3839044/1‑62 (MQ=255)
                                                        ttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt                 >  1:3662270/1‑69 (MQ=255)
                                                         tAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt               <  1:1294206/70‑1 (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:2028959/70‑1 (MQ=255)
                                                          aaCCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg              <  1:3691445/70‑1 (MQ=255)
                                                           aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTg              <  1:1220340/69‑1 (MQ=255)
                                                            ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                       <  1:789490/59‑1 (MQ=255)
                                                             cGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgag                                   <  1:3727569/46‑1 (MQ=255)
                                                              gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt         <  1:1187589/71‑1 (MQ=255)
                                                              gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt         <  1:320840/71‑1 (MQ=255)
                                                              gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGt         <  1:4525289/71‑1 (MQ=255)
                                                               ctgagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg        >  1:1464334/3‑71 (MQ=255)
                                                               cAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcactg                     <  1:5480294/58‑3 (MQ=255)
                                                               cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                        >  1:316088/1‑55 (MQ=255)
                                                                agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg       <  1:5525763/71‑1 (MQ=255)
                                                                 gagaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      <  1:2267324/71‑1 (MQ=255)
                                                                 gaaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      <  1:1425037/68‑1 (MQ=255)
                                                                  agaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      <  1:4600385/70‑1 (MQ=255)
                                                                  agaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt      <  1:5014862/70‑1 (MQ=255)
                                                                  aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATaa                              <  1:3403752/44‑1 (MQ=255)
                                                                   gaAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa    >  1:2568315/1‑71 (MQ=255)
                                                                   gaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa    <  1:4244904/71‑1 (MQ=255)
                                                                   aaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                             >  1:4014494/2‑46 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTg              <  1:1060063/60‑1 (MQ=38)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACg          <  1:2818678/64‑1 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaat   >  1:2452531/1‑71 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaat   >  1:3358243/1‑71 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa    >  1:296888/1‑70 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                            >  1:3090713/1‑46 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGc                            >  1:1007590/1‑46 (MQ=255)
                                                                    aaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa    >  1:4839950/1‑70 (MQ=255)
                                                                     aGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaat   <  1:4325223/70‑1 (MQ=255)
                                                                      |                                                                     
AGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATG  >  minE/217877‑218016

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: