Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,266 T→C 10.4% E59E (GAA→GAG) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2660TC10.4% 209.5 / 13.0 77E59E (GAA→GAG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (3/5);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGG  >  minE/1934196‑1934334
                                                                      |                                                                    
cTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttct                                                                      <  1:3938307/71‑1 (MQ=255)
  ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                                                    >  1:3348784/1‑71 (MQ=255)
      ccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt                                                                     <  1:3335348/66‑1 (MQ=255)
       cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                                                    <  1:1460137/66‑1 (MQ=255)
         tgataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcc                                                             <  1:2818886/69‑2 (MQ=255)
         taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                             <  1:2329692/71‑1 (MQ=255)
         taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                             <  1:170121/71‑1 (MQ=255)
         taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                             <  1:425590/71‑1 (MQ=255)
          aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            >  1:3586116/1‑71 (MQ=255)
          aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            >  1:3785940/1‑71 (MQ=255)
          aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            >  1:37839/1‑71 (MQ=255)
           ataCGACGGGTCATCATTCGGCCAGCAACCGAGCCTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcct                                                            <  1:860711/70‑3 (MQ=255)
            taCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                                                           <  1:2698314/70‑1 (MQ=255)
            taCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            <  1:190178/69‑1 (MQ=255)
                acgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            <  1:760319/65‑1 (MQ=255)
                acgGGTCATCATGCGGCCAGCAACGGAGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                            <  1:2024128/65‑1 (MQ=255)
                 cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                                                           >  1:3508625/1‑65 (MQ=255)
                  gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                    <  1:3623730/71‑1 (MQ=255)
                    gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                  <  1:1904250/71‑1 (MQ=255)
                     tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                                              >  1:509796/1‑58 (MQ=255)
                     tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                 >  1:2079312/1‑71 (MQ=255)
                      catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGt                                                                  >  1:3870703/1‑53 (MQ=255)
                      catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                 >  1:2449958/1‑70 (MQ=255)
                       atcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                             <  1:2672363/57‑1 (MQ=255)
                         catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                                                           >  1:1086959/1‑57 (MQ=255)
                         catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                                              <  1:3895473/70‑1 (MQ=255)
                         catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                                             >  1:1004369/1‑71 (MQ=255)
                          atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc                                            >  1:3580720/1‑71 (MQ=255)
                               gCCAGCAACGGCGACTTCTATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                       <  1:3632756/71‑1 (MQ=255)
                               gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                  >  1:430752/1‑60 (MQ=255)
                                ccAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                 <  1:590094/60‑1 (MQ=255)
                                 cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt                                     >  1:797984/1‑71 (MQ=255)
                                  aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                >  1:989762/1‑59 (MQ=255)
                                  aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                    >  1:56651/1‑71 (MQ=255)
                                  aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                    <  1:2212748/71‑1 (MQ=255)
                                   gCAACGGCGACTTCGATGTTCAGTGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                   <  1:3438780/71‑1 (MQ=255)
                                   gCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                   <  1:994345/71‑1 (MQ=255)
                                    cAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                  >  1:1360918/1‑55 (MQ=255)
                                      aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                       <  1:3577816/64‑1 (MQ=255)
                                      aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                   <  1:2617371/68‑1 (MQ=255)
                                      aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTc                                <  1:2290633/71‑1 (MQ=255)
                                           gACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGcc                                                        <  1:238993/42‑1 (MQ=255)
                                            aCTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAgg                           <  1:1001321/70‑1 (MQ=255)
                                             cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGt                                                      >  1:413024/1‑42 (MQ=255)
                                             cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                          >  1:1893917/1‑54 (MQ=255)
                                             cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                          >  1:3077641/1‑54 (MQ=255)
                                             cTTCGATGTTCAGCGCTTACAGTTCTTCGTTCTCTTTGCCg                                                       >  1:2315693/1‑41 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                 <  1:1798553/44‑1 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTg                                         >  1:1663257/1‑52 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGt                                 >  1:3374137/1‑60 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTAt                        >  1:2489960/1‑69 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATg                       <  1:19321/70‑1 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                      >  1:2930972/1‑71 (MQ=255)
                                                cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                      >  1:874366/1‑71 (MQ=255)
                                                 gATGTTCAGCGCTTCCAGTTCTTCGGTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                      <  1:1010110/70‑1 (MQ=255)
                                                  aTGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCaga                             >  1:1201739/1‑62 (MQ=255)
                                                   tGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                  >  1:3409677/1‑40 (MQ=255)
                                                     ttCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTc                                               <  1:2001492/41‑1 (MQ=255)
                                                      tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                    <  1:810349/35‑1 (MQ=255)
                                                      tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                  >  1:177321/1‑69 (MQ=255)
                                                       cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                >  1:2138374/1‑38 (MQ=255)
                                                       cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc               >  1:3699312/1‑71 (MQ=255)
                                                       cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc               >  1:3419140/1‑71 (MQ=255)
                                                       cAGAGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                >  1:1711027/1‑38 (MQ=255)
                                                        aGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACg              >  1:2465454/1‑71 (MQ=255)
                                                         gcgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACgg             <  1:249032/71‑1 (MQ=255)
                                                          cgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGa            <  1:640171/71‑1 (MQ=255)
                                                           gcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGaa           <  1:95399/71‑1 (MQ=255)
                                                           gattCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaa           <  1:3352812/69‑1 (MQ=255)
                                                            cTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTTGTCAGAGGTATGATCGCGACGGaa           <  1:1306342/70‑1 (MQ=255)
                                                             ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaaa          >  1:3717359/1‑70 (MQ=255)
                                                             ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt         >  1:1094387/1‑71 (MQ=255)
                                                              tCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATc        >  1:1628646/1‑71 (MQ=255)
                                                              tCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt         >  1:1441068/1‑70 (MQ=255)
                                                               ccAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAAt         <  1:914869/69‑1 (MQ=255)
                                                                cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGa                >  1:135864/1‑61 (MQ=255)
                                                                cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCg       <  1:277073/70‑1 (MQ=255)
                                                                cAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATca       <  1:2108073/70‑2 (MQ=255)
                                                                cAGTTATTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                   >  1:1832902/1‑42 (MQ=255)
                                                                 aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt     >  1:2157153/1‑71 (MQ=255)
                                                                 aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt     >  1:1236295/1‑71 (MQ=255)
                                                                   ttcttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCg   <  1:3722214/71‑1 (MQ=255)
                                                                     cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                               <  1:514162/41‑1 (MQ=255)
                                                                     cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                               <  1:3893408/41‑1 (MQ=255)
                                                                     cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg  >  1:1925090/1‑70 (MQ=255)
                                                                     cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg  >  1:1519415/1‑70 (MQ=255)
                                                                      ttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCagag                            <  1:1117770/43‑1 (MQ=255)
                                                                      |                                                                    
CTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGG  >  minE/1934196‑1934334

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: