Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,948 C→G 26.5% A165A (GCG→GCC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9480CG26.5% 131.3 / 52.7 83A165A (GCG→GCC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (24/37);  new base G (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 = 9.98e-01

TTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACT  >  minE/1933889‑1934018
                                                           |                                                                      
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:397563/71‑1 (MQ=255)
             gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                                <  1:2935649/71‑1 (MQ=255)
             gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                                                <  1:4292694/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:962843/71‑1 (MQ=255)
                     tCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATc                                        >  1:1704494/1‑71 (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:468999/71‑1 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                <  1:1631048/71‑1 (MQ=255)
                             gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                <  1:3141770/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: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:2749457/1‑71 (MQ=255)
                                     cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                        >  1:462623/1‑71 (MQ=255)
                                     cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                        >  1:5245815/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)
                                                           |                                                                      
TTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACT  >  minE/1933889‑1934018

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: