Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,945 G→A 26.6% R166R (CGC→CGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9450GA26.6% 128.6 / 49.4 83R166R (CGC→CGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (23/37);  new base A (11/11);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933875‑1934014
                                                                      |                                                                     
tCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAg                                                                       >  1:651522/1‑71 (MQ=255)
tCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATAGAGATAACGCTGACGATAg                                                                       >  1:2970932/1‑71 (MQ=255)
   gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcg                                                                     >  1:607602/1‑70 (MQ=255)
              ttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATc                                                         >  1:2134591/1‑71 (MQ=255)
                 aaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTg                                                      <  1:2590196/71‑1 (MQ=255)
                  aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                                     <  1:1453479/71‑1 (MQ=255)
                   ggTGTTGCGGGATTCATCGTTGGAGATTAGATCGAGATAACGCTGACGATAGCGCGCTTcc                                                              <  1:4083618/61‑1 (MQ=255)
                   ggTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTcc                                                              <  1:785053/61‑1 (MQ=255)
                   ggTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCa                                                    <  1:2689667/71‑1 (MQ=255)
                    gtgtTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCa                                                    <  1:24950/70‑1 (MQ=255)
                     tgtTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTg                                                            >  1:749180/1‑61 (MQ=255)
                       ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGc                                                 <  1:2095220/70‑1 (MQ=255)
                       ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGc                                                 <  1:4792362/70‑1 (MQ=255)
                         gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                              >  1:3136849/1‑71 (MQ=255)
                         gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                              >  1:1526006/1‑71 (MQ=255)
                         gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                              >  1:1441388/1‑71 (MQ=255)
                         gCGGGATTCATCGTTGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTg                                             >  1:5055691/1‑71 (MQ=255)
                           gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                            <  1:4292694/71‑1 (MQ=255)
                           gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                            <  1:2935649/71‑1 (MQ=255)
                           gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                            <  1:397563/71‑1 (MQ=255)
                            ggATTTATCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGa                                           <  1:94396/71‑1 (MQ=255)
                            ggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:5435180/71‑1 (MQ=255)
                             gATTTATCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTg                                                      >  1:2460759/1‑59 (MQ=255)
                                 taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       >  1:5341320/2‑70 (MQ=255)
                                 taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       <  1:1652139/69‑1 (MQ=255)
                                 cATCGTAGGAGATGACATCGAGATAACGCTGACGATAGCGCGCTTcc                                                              >  1:3249342/1‑47 (MQ=38)
                                  aTCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTAt                                     <  1:2347465/71‑1 (MQ=255)
                                   tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                    >  1:5022455/1‑71 (MQ=255)
                                   tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                    <  1:5118446/71‑1 (MQ=255)
                                   tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                    >  1:1704494/1‑71 (MQ=255)
                                   tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                    <  1:962843/71‑1 (MQ=255)
                                    cGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                    <  1:4640566/70‑1 (MQ=255)
                                     gTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCg                                  >  1:2002469/1‑71 (MQ=255)
                                     gTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCg                                  >  1:4588203/1‑71 (MQ=255)
                                      ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                 <  1:3679940/71‑1 (MQ=255)
                                        ggAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                 >  1:2110155/1‑69 (MQ=255)
                                           gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:3141770/71‑1 (MQ=255)
                                           gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:468999/71‑1 (MQ=255)
                                           gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:4376784/71‑1 (MQ=255)
                                           gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                            <  1:1631048/71‑1 (MQ=255)
                                              gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                          >  1:4377217/1‑70 (MQ=255)
                                              gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                          >  1:3111569/1‑70 (MQ=255)
                                              gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         >  1:2724307/1‑71 (MQ=255)
                                              gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         >  1:2646964/1‑71 (MQ=255)
                                               agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                         <  1:570914/70‑1 (MQ=255)
                                                 atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTg                                             >  1:1639718/1‑48 (MQ=255)
                                                 atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCa                               >  1:1219607/1‑62 (MQ=255)
                                                 atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCa                               <  1:605297/62‑1 (MQ=255)
                                                 atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCa                               <  1:2942886/62‑1 (MQ=255)
                                                 atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGc                       >  1:3726896/1‑70 (MQ=255)
                                                   cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                              >  1:5324149/1‑45 (MQ=255)
                                                   cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                    >  1:462623/1‑71 (MQ=255)
                                                   cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                    >  1:5245815/1‑71 (MQ=255)
                                                   cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                    >  1:2749457/1‑71 (MQ=255)
                                                   cGAGATAACGCTGACGAT‑‑CGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                  >  1:421839/1‑71 (MQ=255)
                                                     agaTAACGCTGTCGATAGCGCGCTTCCTGATCCTGCaa                                                   <  1:136216/38‑1 (MQ=255)
                                                     agaTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                  <  1:2135753/71‑1 (MQ=255)
                                                       aTAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACgt                       <  1:5205120/64‑2 (MQ=255)
                                                       aTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgta                                           <  1:1407849/44‑3 (MQ=255)
                                                       aTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:892463/44‑1 (MQ=255)
                                                        tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:4025334/43‑1 (MQ=255)
                                                        tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTAt                                     <  1:1386521/49‑1 (MQ=255)
                                                         aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGaa                                          >  1:22073/1‑43 (MQ=255)
                                                         aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtt                >  1:4850361/1‑69 (MQ=255)
                                                         aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtt                >  1:3160493/1‑69 (MQ=255)
                                                          aCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtttag             <  1:3315533/71‑3 (MQ=255)
                                                          aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGcc                                                >  1:1559433/1‑36 (MQ=255)
                                                          aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATcc                                   >  1:2591524/1‑49 (MQ=255)
                                                            gttGACGATAACGGACTTCCTGATCCTGCAGACCATGGGATTTATCTGGTAAAGGACGTAGTGCTTTAGTc           <  1:880468/69‑1 (MQ=255)
                                                            gCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                <  1:3626070/66‑1 (MQ=255)
                                                             ttGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCttta              <  1:4037612/66‑2 (MQ=255)
                                                             cTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                <  1:4503008/65‑1 (MQ=255)
                                                              tGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:1782687/37‑1 (MQ=255)
                                                              tGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                                           <  1:3266450/37‑1 (MQ=255)
                                                               gACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       >  1:4853913/1‑40 (MQ=255)
                                                               gACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                       >  1:243451/1‑40 (MQ=255)
                                                               gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag         >  1:1287451/1‑70 (MQ=255)
                                                               gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag         <  1:408624/70‑1 (MQ=255)
                                                               gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag         >  1:460533/1‑70 (MQ=255)
                                                                aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                           <  1:3369224/51‑1 (MQ=255)
                                                                aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                           <  1:2177350/51‑1 (MQ=255)
                                                                aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca       <  1:871873/71‑1 (MQ=255)
                                                                aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca       <  1:5361150/71‑1 (MQ=255)
                                                                aCGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcagc        <  1:4004532/70‑1 (MQ=255)
                                                                   aTAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACg                        <  1:3778771/51‑1 (MQ=255)
                                                                   aTAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                  <  1:4950430/57‑1 (MQ=255)
                                                                      gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc  <  1:3250009/70‑1 (MQ=255)
                                                                      |                                                                     
TCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGC  >  minE/1933875‑1934014

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: