Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,266,180 A→G 9.6% Y88Y (TAT→TAC 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,1800AG 9.6% 237.2 / 14.7 83Y88Y (TAT→TACtufBprotein 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.38e-01

GCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACG  >  minE/2266111‑2266221
                                                                     |                                         
gCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACAt                                          <  1:957302/71‑1 (MQ=34)
    tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      >  1:4802934/1‑71 (MQ=34)
    tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      <  1:4269337/71‑1 (MQ=34)
    tCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTc                                      <  1:3212939/71‑1 (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:1301986/71‑1 (MQ=34)
              caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            <  1:525607/71‑1 (MQ=34)
              caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            <  1:2159842/71‑1 (MQ=34)
              caACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATCTTTTTAACATAGTCGGCg                                  >  1:362186/1‑65 (MQ=255)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:2835901/1‑71 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:3835567/1‑71 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:4277662/71‑1 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:3137970/71‑1 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           >  1:4448339/1‑71 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:6012114/71‑1 (MQ=34)
               aaCTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCg                           <  1:1346157/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:6441363/70‑1 (MQ=35)
                            gcgcCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACgt               <  1:4563692/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:2619866/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:6564143/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:3479845/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:2165152/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:2701333/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGt      >  1:4500411/1‑71 (MQ=34)
                                    cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:3641480/1‑48 (MQ=2)
                                    cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:6080001/1‑48 (MQ=2)
                                    cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:2428734/1‑48 (MQ=2)
                                    cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:2632474/1‑48 (MQ=2)
                                    cATCTGCGCAGCACCGGTAATCATGTTCTTGATGTAGTCGGCGTGgcc                             >  1:4271202/1‑48 (MQ=2)
                                        tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGc                                   >  1:1620691/1‑38 (MQ=40)
                                        tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGc                                   >  1:4184696/1‑38 (MQ=40)
                                        tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGcccc                            >  1:4639724/1‑45 (MQ=38)
                                        tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:16282/56‑1 (MQ=35)
                                        tGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTAc                 <  1:2273957/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)
                                                                     |                                         
GCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGCGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACG  >  minE/2266111‑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: