Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,327 A→G 10.2% V372V (GTT→GTC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,3270AG10.2% 210.2 / 13.0 78V372V (GTT→GTClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base G (4/4);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCA  >  minE/1933258‑1933393
                                                                     |                                                                  
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:1592644/1‑70 (MQ=255)
    ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                                                                >  1:3495522/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:2327411/1‑70 (MQ=255)
                  gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc                                                  >  1:2167486/1‑70 (MQ=255)
                   aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAAcccc                                                 >  1:815524/1‑69 (MQ=255)
                   aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                                 <  1:3802313/70‑1 (MQ=255)
                       ttggCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAACCCCCAGctt                                            <  1:3338806/69‑2 (MQ=255)
                              aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:2603798/69‑1 (MQ=255)
                              aaTCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:854070/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:952308/71‑1 (MQ=255)
                                tCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACg                                   <  1:934048/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:1017272/52‑1 (MQ=255)
                                               ccACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTtctc                                       <  1:755818/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:2615521/71‑1 (MQ=255)
                                                            gATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCtt       <  1:3576120/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)
                                                                   gTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGc                            <  1:1315475/43‑1 (MQ=255)
                                                                   gTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTCCTTTCGCa  >  1:3709521/1‑69 (MQ=255)
                                                                     |                                                                  
TGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCA  >  minE/1933258‑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: