Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,324 G→T 5.8% T373T (ACC→ACA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,3240GT 5.8% 349.5 / 3.9 103T373T (ACC→ACAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (54/43);  new base T (2/4);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 107 total aligned reads displayed.

CTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCG  >  minE/1933257‑1933391
                                                                   |                                                                   
cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt                                                                   >  1:2090698/1‑70 (MQ=255)
cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCAAAGATCTCGGt                                                                   >  1:746323/1‑70 (MQ=255)
   ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAAc                                                                <  1:858309/70‑1 (MQ=255)
    cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:1175725/71‑1 (MQ=255)
    cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:2631428/71‑1 (MQ=255)
    cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:3512688/71‑1 (MQ=255)
     ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              >  1:2254113/1‑70 (MQ=255)
     ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              >  1:2501391/1‑70 (MQ=255)
     ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGAt                                                             <  1:358778/71‑1 (MQ=255)
       atatTCCGTAATAAAGGTAGGCTGAATCAGATGTGCTTCTGCCTCTTCATCAAAGATCTCTGTGACAATAc                                                           <  1:795829/71‑1 (MQ=255)
          ttCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGc                                                        >  1:2923334/1‑71 (MQ=255)
             cgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCa                                                     <  1:1707106/70‑1 (MQ=255)
             aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:2544326/62‑1 (MQ=255)
             aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa                                                     <  1:3744619/71‑1 (MQ=255)
               tAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGa                                                   >  1:2526580/1‑71 (MQ=255)
                 aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                                                                     >  1:1578023/1‑51 (MQ=255)
                 aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt                                                                   <  1:2412059/53‑1 (MQ=255)
                 aTGAAGGTAGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGc                                                        >  1:1472087/1‑64 (MQ=255)
                  taaaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAc                                                                <  1:4032444/53‑1 (MQ=255)
                  tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc                                                           <  1:2278648/60‑1 (MQ=255)
                   gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                                                                     >  1:764337/1‑49 (MQ=255)
                   gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                            >  1:849269/1‑58 (MQ=255)
                   gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                            >  1:2230538/1‑58 (MQ=255)
                   gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                >  1:3735313/1‑70 (MQ=255)
                   gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                >  1:3981425/1‑70 (MQ=255)
                     aGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt                                                                   <  1:223396/49‑1 (MQ=37)
                       gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc                                                           >  1:2142712/1‑55 (MQ=255)
                       gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                       >  1:3866125/1‑59 (MQ=255)
                       gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGc                                            >  1:3615909/1‑70 (MQ=255)
                           gCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtcttc                                       <  1:2266095/71‑1 (MQ=255)
                           gCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtctt                                        >  1:408421/1‑70 (MQ=255)
                             tGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtcttct                                      >  1:146083/1‑70 (MQ=255)
                              gAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAg                                             <  1:1381408/62‑1 (MQ=255)
                               aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   <  1:2300994/71‑1 (MQ=255)
                                aTCAGATGTGCTTCTGCCACTTCTTCGAAGCTCTCGGTAACGATACGGCCCAGACCCCAGct                                           >  1:809513/1‑62 (MQ=255)
                                aTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc                                                           >  1:2175364/1‑46 (MQ=255)
                                aTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAAc                                  >  1:845314/1‑71 (MQ=255)
                                 tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:2878206/42‑1 (MQ=255)
                                 ccAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:2488100/41‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAAc                                                                >  1:1353315/1‑38 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                 <  1:1832525/53‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGt                                <  1:2469012/70‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:3343223/71‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:1541852/71‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:3942020/71‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:3982203/71‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:4064049/71‑1 (MQ=255)
                                   aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               >  1:1269669/1‑71 (MQ=255)
                                    gATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              >  1:3235575/1‑39 (MQ=255)
                                    gATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTgg                              <  1:4045724/71‑1 (MQ=255)
                                     aTGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              <  1:1139404/38‑1 (MQ=255)
                                     aTGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                               >  1:3665814/1‑53 (MQ=255)
                                     aTGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGctct                                         <  1:3610464/59‑1 (MQ=255)
                                     aTGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                     >  1:1155314/1‑63 (MQ=255)
                                      tgtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                              >  1:1852934/1‑37 (MQ=255)
                                       gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGa                             >  1:337641/1‑69 (MQ=255)
                                        tgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATg                           <  1:4144918/70‑1 (MQ=255)
                                         gCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCa                                              >  1:3144897/1‑50 (MQ=255)
                                         gCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                         >  1:1222754/1‑71 (MQ=255)
                                          cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                       <  1:2185833/40‑1 (MQ=255)
                                          cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   >  1:2050678/1‑60 (MQ=255)
                                          cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                         <  1:1991402/70‑1 (MQ=255)
                                          cTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCa                                                     >  1:747605/1‑42 (MQ=16)
                                           ttCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                       >  1:1831073/1‑71 (MQ=255)
                                             cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   >  1:1304271/1‑57 (MQ=255)
                                             cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   >  1:1214733/1‑57 (MQ=255)
                                             cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGc                          >  1:2456128/1‑66 (MQ=255)
                                             cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                       <  1:72499/69‑1 (MQ=255)
                                             cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                       <  1:11766/69‑1 (MQ=255)
                                              tGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                       >  1:1364799/1‑68 (MQ=255)
                                               gCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                >  1:2856586/1‑42 (MQ=255)
                                                ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAg                    <  1:2902741/69‑1 (MQ=255)
                                                 cACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATa                     <  1:826830/67‑1 (MQ=255)
                                                  aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGa                             >  1:84642/1‑58 (MQ=255)
                                                  aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                 >  1:374953/1‑70 (MQ=255)
                                                  aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                 >  1:2993763/1‑70 (MQ=255)
                                                     tcttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAg              >  1:2160279/1‑70 (MQ=255)
                                                      cttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc             <  1:3599303/70‑1 (MQ=255)
                                                      cttcGAAGATCTCGGTAACGATACGGCCCAGACACCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc             >  1:2667430/1‑70 (MQ=255)
                                                       ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:4034701/51‑1 (MQ=255)
                                                       ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:7968/51‑1 (MQ=255)
                                                       ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGc                          >  1:1548436/1‑56 (MQ=255)
                                                       ttcGAAGATCTCGGTAACGATACGCCCCAGACCCCAGCTCTTCTCAACGt                                <  1:1103487/50‑1 (MQ=255)
                                                       ttcGAACATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCaa           >  1:4338/1‑71 (MQ=255)
                                                         cGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   >  1:44563/1‑45 (MQ=255)
                                                         cGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCa            <  1:2186850/68‑1 (MQ=255)
                                                         cGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCa            <  1:4113266/68‑1 (MQ=255)
                                                          gAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTg        <  1:1437023/71‑1 (MQ=255)
                                                            aGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                   >  1:1160441/1‑42 (MQ=255)
                                                             gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:12779/45‑1 (MQ=255)
                                                             gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt     >  1:3687825/1‑71 (MQ=255)
                                                             gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt     >  1:3092007/1‑71 (MQ=255)
                                                              aTCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                               <  1:200227/44‑1 (MQ=255)
                                                               tctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTg        >  1:704231/1‑66 (MQ=255)
                                                               tctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTg        >  1:71917/1‑66 (MQ=255)
                                                                ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTgg                              >  1:1969509/1‑43 (MQ=255)
                                                                ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGa                             >  1:3132611/1‑44 (MQ=255)
                                                                 tcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc             <  1:3741170/59‑1 (MQ=255)
                                                                   ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTc   <  1:3822186/67‑1 (MQ=255)
                                                                   ggTAACGATACGGCCCAGACCCCAGCTCTTATCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCg  >  1:3090198/1‑68 (MQ=255)
                                                                   |                                                                   
CTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCG  >  minE/1933257‑1933391

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: