Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,324 G→T 14.0% 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,3240GT14.0% 233.7 / 15.8 78T373T (ACC→ACAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base T (5/6);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCA  >  minE/1933254‑1933393
                                                                      |                                                                     
cTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                                                                       <  1:677601/71‑1 (MQ=255)
 ttCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCgg                                                                      >  1:546764/1‑71 (MQ=255)
   cTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt                                                                     <  1:3167784/70‑1 (MQ=255)
   cTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt                                                                     >  1:1528363/1‑70 (MQ=255)
   cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt                                                                     >  1:324878/1‑70 (MQ=255)
    tGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTaa                                                                   >  1:1404275/1‑71 (MQ=255)
      ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTa                                                                    <  1:2387211/68‑1 (MQ=255)
      ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTa                                                                    <  1:1810719/68‑1 (MQ=255)
       cGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAca                                                                 <  1:2562765/70‑2 (MQ=255)
       cGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAc                                                                  <  1:3765587/69‑1 (MQ=255)
       cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGACACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:831199/1‑71 (MQ=255)
        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:3495522/1‑70 (MQ=255)
        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:1592644/1‑70 (MQ=255)
         gATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                                                              <  1:3900767/71‑1 (MQ=255)
           tatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACg                                                            >  1:436776/1‑71 (MQ=255)
              tCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                         >  1:3371050/1‑71 (MQ=255)
               ccgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATa                                                              >  1:489336/3‑65 (MQ=255)
               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:1451068/1‑63 (MQ=255)
                   aaTAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCaacc                                                    >  1:522311/1‑68 (MQ=255)
                     tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                  <  1:3098484/71‑1 (MQ=255)
                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                  >  1:2167486/1‑70 (MQ=255)
                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                  >  1:2327411/1‑70 (MQ=255)
                       aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGtga                                                                   <  1:1770045/52‑3 (MQ=255)
                       aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAAcccc                                                 >  1:815524/1‑69 (MQ=255)
                       aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                                 <  1:3802313/70‑1 (MQ=255)
                        aGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                                                                       >  1:783357/1‑47 (MQ=255)
                        aGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                                                                       >  1:1087968/1‑47 (MQ=255)
                           ttggCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGctt                                            <  1:3338806/69‑2 (MQ=255)
                                  aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:854070/69‑1 (MQ=255)
                                  aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:2603798/69‑1 (MQ=255)
                                  aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCa                                      <  1:999637/70‑1 (MQ=255)
                                  aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCa                                      <  1:2218973/70‑1 (MQ=255)
                                   aTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAAc                                    <  1:2021833/71‑1 (MQ=255)
                                    tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGct                                             <  1:1377821/61‑1 (MQ=255)
                                    tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACg                                   <  1:934048/71‑1 (MQ=255)
                                    tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACg                                   <  1:952308/71‑1 (MQ=255)
                                     cAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCaa                                     >  1:1336831/1‑68 (MQ=255)
                                      aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAg                                               >  1:2657762/1‑57 (MQ=255)
                                      aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGctct                                           >  1:958536/1‑61 (MQ=255)
                                      aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGt                                  <  1:1039940/70‑1 (MQ=255)
                                      aGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                 >  1:634405/1‑71 (MQ=255)
                                      aGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGCTTTTCTCTAccg                                  <  1:3288293/70‑3 (MQ=255)
                                        aTGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc                                                         <  1:91851/45‑1 (MQ=255)
                                          gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:1051086/1‑36 (MQ=255)
                                          gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGAt                                                               <  1:3281668/37‑1 (MQ=255)
                                          gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATg                             >  1:3084274/1‑71 (MQ=255)
                                          gtgCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATg                             >  1:3077592/1‑71 (MQ=255)
                                             cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGTCCCAGACCCCAGCTCTTCTCAACGTGGATGCCg                          <  1:2445001/71‑1 (MQ=255)
                                             cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGcc                           >  1:830991/1‑70 (MQ=255)
                                             cTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCg                          <  1:2270213/71‑1 (MQ=255)
                                              ttCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                  >  1:2435912/1‑46 (MQ=255)
                                               tCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                         >  1:3042568/1‑70 (MQ=255)
                                               tCTGCCACTTATTCGAAGATCTCGGTAACGATACGGCCCAGAcc                                                   >  1:136297/1‑44 (MQ=255)
                                                cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCtctt                                          >  1:12180/1‑52 (MQ=255)
                                                cTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                         >  1:375866/1‑69 (MQ=255)
                                                  gCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                 >  1:644045/1‑59 (MQ=255)
                                                   ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:755818/52‑1 (MQ=255)
                                                   ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:1017272/52‑1 (MQ=255)
                                                   ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCg                          <  1:398359/65‑1 (MQ=255)
                                                   ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAg                      <  1:2631763/69‑1 (MQ=255)
                                                   ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAt                    <  1:3622973/71‑1 (MQ=255)
                                                    cACTTCTTCGAAGATCTCGGTAACTATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                   <  1:2911318/71‑1 (MQ=255)
                                                    cACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAAc                                    >  1:899521/1‑54 (MQ=255)
                                                    cACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGa                         >  1:1345345/1‑65 (MQ=255)
                                                     aCTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAtt                   >  1:2278324/1‑70 (MQ=255)
                                                      cttcttCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAAc                                    >  1:3102110/1‑52 (MQ=255)
                                                         cttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGAt                    <  1:1293342/65‑1 (MQ=255)
                                                         cttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGc               >  1:1524940/1‑70 (MQ=255)
                                                          ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGa                     <  1:2562539/63‑1 (MQ=255)
                                                          ttcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCa              >  1:1968865/1‑70 (MQ=255)
                                                           tcGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTc                  >  1:872003/1‑65 (MQ=255)
                                                            cGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTg                                 >  1:3761292/1‑49 (MQ=255)
                                                             gAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAAtt           >  1:3411652/1‑70 (MQ=255)
                                                              aaGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGc         >  1:1393313/1‑71 (MQ=255)
                                                                gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt       <  1:3576120/71‑1 (MQ=255)
                                                                gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt       <  1:2615521/71‑1 (MQ=255)
                                                                 aTCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCa                 >  1:1632895/1‑60 (MQ=255)
                                                                   ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACg                                   <  1:1922717/40‑1 (MQ=255)
                                                                   ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACg                                   <  1:2088080/40‑1 (MQ=255)
                                                                   ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGAt                        >  1:3912951/1‑51 (MQ=255)
                                                                   ctcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCaa             >  1:2139955/1‑62 (MQ=255)
                                                                    tcGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCg    >  1:1590480/1‑70 (MQ=255)
                                                                      ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTc     <  1:2143175/67‑1 (MQ=255)
                                                                      ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCg    >  1:1344165/1‑68 (MQ=255)
                                                                      ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGc   <  1:505269/69‑1 (MQ=255)
                                                                      ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCa  <  1:920000/70‑1 (MQ=255)
                                                                      ggTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCa  <  1:924726/70‑1 (MQ=255)
                                                                      |                                                                     
CTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCA  >  minE/1933254‑1933393

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: