Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,266,177 A→G 9.6% V89V (GTT→GTC) ‡ tufB ← protein chain elongation factor EF‑Tu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,266,1770AG 9.6% 237.6 / 14.7 83V89V (GTT→GTC) ‡tufBprotein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base G (5/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACG  >  minE/2266109‑2266221
                                                                    |                                            
gggCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTaa                                             >  1:4176876/1‑70 (MQ=35)
gggCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTaa                                             >  1:5788848/1‑70 (MQ=35)
  gCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACAt                                          <  1:957302/71‑1 (MQ=34)
      tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      <  1:3212939/71‑1 (MQ=34)
      tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      <  1:4269337/71‑1 (MQ=34)
      tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      >  1:4802934/1‑71 (MQ=34)
      tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      <  1:1591960/71‑1 (MQ=34)
        aGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCgg                                    <  1:3002446/71‑1 (MQ=34)
        aGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCgg                                    >  1:6749481/1‑71 (MQ=34)
        aGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCgg                                    >  1:1463028/1‑71 (MQ=34)
        aGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     >  1:2431770/1‑70 (MQ=35)
        aGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:470721/70‑1 (MQ=35)
         gTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGc                                   <  1:849833/71‑1 (MQ=34)
          tCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTGGGCg                                  <  1:4715803/71‑1 (MQ=21)
          tCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCg                                  <  1:2916673/71‑1 (MQ=34)
             cagcaACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGc                               <  1:1619526/71‑1 (MQ=34)
                caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            <  1:2159842/71‑1 (MQ=34)
                caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            <  1:1301986/71‑1 (MQ=34)
                caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            <  1:525607/71‑1 (MQ=34)
                caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATCTTTTTAACATAGTCGGCg                                  >  1:362186/1‑65 (MQ=255)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:2835901/1‑71 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:1346157/71‑1 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:3137970/71‑1 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:3835567/1‑71 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:6012114/71‑1 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:4448339/1‑71 (MQ=34)
                 aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:4277662/71‑1 (MQ=34)
                    tACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                         >  1:358940/1‑70 (MQ=35)
                    tACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCgg                          >  1:5942885/1‑69 (MQ=35)
                     aCCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCggg                         >  1:6132928/1‑69 (MQ=35)
                       cAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGt                     <  1:2206827/71‑1 (MQ=34)
                         ggATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTc                    >  1:6338468/1‑70 (MQ=35)
                         ggATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTc                    >  1:1037028/1‑70 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:4717679/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:6752560/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:919969/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:1192888/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:1951630/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:2083631/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:2229981/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:2256709/52‑1 (MQ=35)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTa                  >  1:4562119/1‑71 (MQ=34)
                          gATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTa                  <  1:4440998/71‑1 (MQ=34)
                           aTCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:3780647/51‑1 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 >  1:1884785/1‑70 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 >  1:3839918/1‑70 (MQ=35)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                >  1:3976009/1‑71 (MQ=34)
                            tCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACg                >  1:4389025/1‑71 (MQ=34)
                              gcgcCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACgt               <  1:4563692/70‑1 (MQ=35)
                              gcgcCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACgt               <  1:6441363/70‑1 (MQ=35)
                              gcgcCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACgt               <  1:1420581/70‑1 (MQ=11)
                              gcgcCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACgt               <  1:6204939/70‑1 (MQ=11)
                              gcgcCGTCCATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGGCCTGGGCAGTCAACgt               <  1:5514204/70‑1 (MQ=11)
                               cgcCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACgt               <  1:136568/69‑1 (MQ=35)
                                 ccGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGt                                 >  1:2047949/1‑49 (MQ=38)
                                 ccGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGc           >  1:6018642/1‑71 (MQ=34)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAg                      <  1:4517093/59‑1 (MQ=35)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAg                      <  1:3644896/59‑1 (MQ=35)
                                  cGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAg                      <  1:283892/59‑1 (MQ=35)
                                    tCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCgg                                    <  1:1776275/43‑1 (MQ=38)
                                    tCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCg                                     <  1:351835/42‑1 (MQ=38)
                                    tCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTc                    >  1:6106337/1‑59 (MQ=35)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcc                              >  1:2915887/1‑47 (MQ=38)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:3479845/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:6564143/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:2165152/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:4500411/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:2619866/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:2701333/1‑71 (MQ=34)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:2428734/1‑48 (MQ=2)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:3641480/1‑48 (MQ=2)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:6080001/1‑48 (MQ=2)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:4271202/1‑48 (MQ=2)
                                      cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:2632474/1‑48 (MQ=2)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGc                                   >  1:4184696/1‑38 (MQ=40)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGc                                   >  1:1620691/1‑38 (MQ=40)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            >  1:4639724/1‑45 (MQ=38)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:2273957/56‑1 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:16282/56‑1 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:4300367/56‑1 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:2815004/56‑1 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACgt               >  1:632122/1‑58 (MQ=35)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACg  >  1:6172737/1‑71 (MQ=34)
                                          tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACg  >  1:1313985/1‑71 (MQ=34)
                                            cgcAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGa    <  1:6288857/67‑1 (MQ=31)
                                                                    |                                            
GGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACG  >  minE/2266109‑2266221

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: