Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 863,791 A→C 6.5% Y31S (TAT→TCT) ‡ yciQ → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE863,7910AC 6.5% 217.2 / 9.4 77Y31S (TAT→TCT) ‡yciQpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/43);  new base C (3/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGA  >  minE/863722‑863860
                                                                     |                                                                     
ttttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTAt                                                                      >  1:1785746/1‑71 (MQ=255)
ttttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTAt                                                                      >  1:435405/1‑71 (MQ=255)
  ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                     >  1:1769590/1‑70 (MQ=255)
  ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                     <  1:4531384/70‑1 (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:1781082/50‑1 (MQ=25)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAgg                                                              <  1:3494380/50‑1 (MQ=25)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                          <  1:4247866/70‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                          <  1:574056/70‑1 (MQ=255)
                             tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                         <  1:1019086/71‑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)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                            >  1:3087692/1‑50 (MQ=21)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                            >  1:829399/1‑50 (MQ=21)
                               atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                            >  1:270836/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:2632304/69‑1 (MQ=255)
                                    tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                    <  1:2587104/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:1506010/57‑1 (MQ=255)
                                                           aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                         <  1:3627075/57‑1 (MQ=255)
                                                            tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc          <  1:369722/71‑1 (MQ=255)
                                                            tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc          <  1:2320173/71‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                         <  1:28689/55‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                         <  1:2210202/55‑1 (MQ=255)
                                                             cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                         <  1:2176745/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:4398449/1‑63 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa           <  1:1652286/63‑1 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa           <  1:3824416/63‑1 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa           >  1:4738369/1‑63 (MQ=255)
                                                                   ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa           <  1:3327521/63‑1 (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)
                                                                     |                                                                     
TTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGA  >  minE/863722‑863860

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: