Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 863,792 T→A 6.3% Y31* (TAT→TAA) ‡ yciQ → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE863,7920TA 6.3% 236.8 / 6.7 79Y31* (TAT→TAA) ‡yciQpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/43);  new base A (3/2);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAA  >  minE/863722‑863861
                                                                      |                                                                     
ttttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTAt                                                                       >  1:1785746/1‑71 (MQ=255)
ttttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTAt                                                                       >  1:435405/1‑71 (MQ=255)
  ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      <  1:4531384/70‑1 (MQ=255)
  ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:1769590/1‑70 (MQ=255)
   tGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:215396/1‑69 (MQ=255)
    gCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:1225798/1‑68 (MQ=255)
     cTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaag                                                                   <  1:285349/70‑1 (MQ=255)
        ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAgg                                                               >  1:1696408/1‑71 (MQ=255)
        ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAgg                                                               >  1:2065495/1‑71 (MQ=255)
            tAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtc                                                           >  1:1590703/1‑71 (MQ=255)
              gCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                         >  1:2046183/1‑71 (MQ=255)
               cAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                         >  1:220302/1‑70 (MQ=255)
               cAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        >  1:4117445/1‑71 (MQ=255)
                aGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        <  1:2007772/70‑1 (MQ=255)
                aGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        <  1:1889471/70‑1 (MQ=255)
                       tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGaagaag                                                                <  1:3138625/55‑1 (MQ=255)
                       tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGaagaag                                                                <  1:4417079/55‑1 (MQ=255)
                       tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTg                                                 <  1:1280187/70‑1 (MQ=255)
                         tGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATg                                              >  1:3393488/1‑71 (MQ=255)
                           tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGtt                                            <  1:306970/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAgg                                                               <  1:3494380/50‑1 (MQ=25)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAgg                                                               <  1:1781082/50‑1 (MQ=25)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                           <  1:4247866/70‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                           <  1:574056/70‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:3635176/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:3781643/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1603042/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:4410732/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:699198/71‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1019086/71‑1 (MQ=255)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:270836/1‑50 (MQ=21)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:3087692/1‑50 (MQ=21)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:829399/1‑50 (MQ=21)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAg                                        >  1:1666873/1‑71 (MQ=255)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAg                                        >  1:4423694/1‑71 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:4261111/69‑1 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:2587104/69‑1 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:2632304/69‑1 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:1884001/69‑1 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   >  1:2234809/1‑71 (MQ=255)
                                      aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGa                                                >  1:2750890/1‑56 (MQ=255)
                                      aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGa                                                >  1:2995220/1‑56 (MQ=255)
                                      aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCActct                                 <  1:4695120/71‑1 (MQ=255)
                                      aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   <  1:640971/69‑1 (MQ=255)
                                      aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCAActct                                 >  1:2307132/1‑71 (MQ=255)
                                       gTTATGCAGAAAACACGGAGATCCCTTCTTATgaaga                                                                  <  1:4407412/37‑1 (MQ=255)
                                        ttATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTg                                >  1:2616481/1‑70 (MQ=255)
                                         tATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGc                               >  1:3133113/1‑70 (MQ=255)
                                         tATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGc                               >  1:1385211/1‑70 (MQ=255)
                                          aTGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCaa                             >  1:4431598/1‑71 (MQ=255)
                                          aTGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCGTTGATGTTGAAGCCACTCTGCaa                             >  1:1155152/1‑71 (MQ=255)
                                             cAGAAAACACGGAGATCCCTTCTTATGTAGAAGGGATCTCGCTCTTTGATGTTGa                                          >  1:2020278/1‑55 (MQ=255)
                                             cAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1108434/71‑1 (MQ=255)
                                              aGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:3976209/71‑1 (MQ=255)
                                               gAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:4383889/70‑1 (MQ=255)
                                                  aaCACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                     >  1:2178086/1‑71 (MQ=255)
                                                    cacGGAGATCCCTTCTTATTAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1272243/48‑1 (MQ=255)
                                                    cacGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1640396/48‑1 (MQ=255)
                                                        gagaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   <  1:1833316/51‑1 (MQ=255)
                                                           aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1289934/57‑1 (MQ=255)
                                                           aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:3627075/57‑1 (MQ=255)
                                                           aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1506010/57‑1 (MQ=255)
                                                            tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc           <  1:369722/71‑1 (MQ=255)
                                                            tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc           <  1:2320173/71‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:2176745/55‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:28689/55‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:2210202/55‑1 (MQ=255)
                                                              ccTTCTTATGAAGAAGTTATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:262598/54‑1 (MQ=255)
                                                              ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:3163966/55‑1 (MQ=255)
                                                                ttcttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                    <  1:1653596/58‑1 (MQ=255)
                                                                ttcttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATc        >  1:1456205/1‑70 (MQ=255)
                                                                  cttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCa       >  1:3547894/1‑69 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:1652286/63‑1 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:3327521/63‑1 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            >  1:4398449/1‑63 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:3824416/63‑1 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            >  1:4738369/1‑63 (MQ=255)
                                                                     aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                      >  1:1463483/1‑51 (MQ=255)
                                                                     aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                      >  1:4130459/1‑51 (MQ=255)
                                                                     aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAg    >  1:1100416/1‑69 (MQ=255)
                                                                     aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa   >  1:3336499/1‑70 (MQ=255)
                                                                      tGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGaa  <  1:518666/70‑1 (MQ=255)
                                                                      tGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGaa  <  1:3871669/70‑1 (MQ=255)
                                                                      |                                                                     
TTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAA  >  minE/863722‑863861

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: