Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,330 G→A 7.7% I371I (ATC→ATT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,3300GA 7.7% 231.6 / 7.1 78I371I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base A (4/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCG  >  minE/1933261‑1933399
                                                                     |                                                                     
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)
                aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAAcccc                                                       >  1:815524/1‑69 (MQ=255)
                aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc                                                       <  1:3802313/70‑1 (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: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: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)
                                                                gTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGc                                  <  1:1315475/43‑1 (MQ=255)
                                                                gTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTCCTTTCGCa        >  1:3709521/1‑69 (MQ=255)
                                                                   aCGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTc   >  1:3278154/1‑71 (MQ=255)
                                                                    cGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCg  <  1:179982/71‑1 (MQ=255)
                                                                     gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCaga      >  1:989844/1‑66 (MQ=255)
                                                                     gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCg  >  1:1725207/1‑70 (MQ=255)
                                                                     gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCg  >  1:1737968/1‑70 (MQ=255)
                                                                     gATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCg  >  1:2230439/1‑70 (MQ=255)
                                                                     |                                                                     
CGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGCTCTTCTCAACGTGGATGCCGATAGATTCAGCAATTGCTTTCGCAGAGTCG  >  minE/1933261‑1933399

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: