Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,949 G→A 26.5% A165V (GCG→GTG) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9490GA26.5% 133.7 / 48.5 83A165V (GCG→GTG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base A (11/11);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTC  >  minE/1933889‑1934019
                                                            |                                                                      
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:4792362/70‑1 (MQ=255)
         ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGc                                                      <  1:2095220/70‑1 (MQ=255)
           gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                                   >  1:1441388/1‑71 (MQ=255)
           gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                                   >  1:1526006/1‑71 (MQ=255)
           gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGt                                                   >  1:3136849/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:5118446/71‑1 (MQ=255)
                     tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                         >  1:5022455/1‑71 (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:4588203/1‑71 (MQ=255)
                       gTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCg                                       >  1:2002469/1‑71 (MQ=255)
                        ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                      <  1:3679940/71‑1 (MQ=255)
                          ggAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg                                      >  1:2110155/1‑69 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                 <  1:4376784/71‑1 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                 <  1:1631048/71‑1 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                 <  1:3141770/71‑1 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                 <  1:468999/71‑1 (MQ=255)
                                gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                               >  1:3111569/1‑70 (MQ=255)
                                gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                               >  1:4377217/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:605297/62‑1 (MQ=255)
                                   atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCa                                    >  1:1219607/1‑62 (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:5245815/1‑71 (MQ=255)
                                     cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                         >  1:462623/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:3160493/1‑69 (MQ=255)
                                           aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtt                     >  1:4850361/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:243451/1‑40 (MQ=255)
                                                 gACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                                            >  1:4853913/1‑40 (MQ=255)
                                                 gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag              >  1:1287451/1‑70 (MQ=255)
                                                 gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag              >  1:460533/1‑70 (MQ=255)
                                                 gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag              <  1:408624/70‑1 (MQ=255)
                                                  aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                                <  1:3369224/51‑1 (MQ=255)
                                                  aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                                <  1:2177350/51‑1 (MQ=255)
                                                  aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca            <  1:5361150/71‑1 (MQ=255)
                                                  aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca            <  1:871873/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)
                                                          gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAAc    >  1:3602980/1‑71 (MQ=255)
                                                           cgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGt                 >  1:2058026/1‑57 (MQ=255)
                                                           cgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt   >  1:1744131/1‑71 (MQ=255)
                                                            gcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTc  <  1:3085453/71‑1 (MQ=255)
                                                            |                                                                      
TTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTC  >  minE/1933889‑1934019

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: