Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,249 A→G 33.3% S398S (TCT→TCC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2490AG33.3% 56.2 / 46.1 54S398S (TCT→TCClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (8/28);  new base G (7/11);  total (15/39)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGA  >  minE/1933180‑1933316
                                                                     |                                                                   
aCGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGga                                                                     <  1:2053879/70‑1 (MQ=255)
  gACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGggga                                                                   >  1:3336015/1‑70 (MQ=255)
  gACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                                   <  1:1762637/70‑1 (MQ=255)
  gACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                                   <  1:2000826/70‑1 (MQ=255)
  gACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgag                                                                    >  1:2691720/1‑69 (MQ=255)
   accaccAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGaaa                                                                 <  1:489697/71‑1 (MQ=255)
     caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                                <  1:2885904/70‑1 (MQ=255)
     caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                               <  1:986997/71‑1 (MQ=255)
        cgaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGgggac                                                                  <  1:2261431/63‑2 (MQ=255)
        caATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                                                            <  1:1541563/71‑1 (MQ=255)
         gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                                                           <  1:2143574/70‑1 (MQ=255)
         gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                                                           <  1:990147/70‑1 (MQ=255)
            gaagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                                        >  1:767057/1‑71 (MQ=255)
               gaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                      >  1:3235564/1‑70 (MQ=255)
                 attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGatat                                                   <  1:1320466/69‑1 (MQ=255)
                  cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                                  <  1:3411527/71‑1 (MQ=255)
                    cAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCgg                                                       <  1:722270/64‑1 (MQ=255)
                    cAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATAtt                                                  <  1:704346/69‑1 (MQ=255)
                     aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTc                                                 <  1:1405310/69‑1 (MQ=255)
                        gCGGTCTGTGATTTCCGGGTTAACGGCGTTACGACGCGCCAGCGgaga                                                                   <  1:2889304/48‑1 (MQ=255)
                         cggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGata                                                    >  1:245831/1‑62 (MQ=255)
                         cGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGga                                                                     <  1:3150566/45‑1 (MQ=255)
                          ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                          <  1:147228/71‑1 (MQ=255)
                           gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                         <  1:1878731/71‑1 (MQ=255)
                           gTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAAtaa                                         <  1:3803418/71‑3 (MQ=255)
                            tCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                              <  1:2256603/49‑1 (MQ=255)
                             cggtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                       >  1:2340843/3‑71 (MQ=255)
                              tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                      <  1:1264421/55‑1 (MQ=255)
                               gcgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTaa                                            <  1:3846345/62‑1 (MQ=255)
                                 gATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGga                                                                     >  1:2496929/1‑37 (MQ=255)
                                  aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTgg                                  >  1:3352935/1‑71 (MQ=255)
                                      ccGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTg                               >  1:2317851/1‑70 (MQ=255)
                                         ggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                     <  1:892255/61‑1 (MQ=255)
                                         ggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                            >  1:3920323/1‑70 (MQ=255)
                                          gTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                            >  1:1595899/1‑69 (MQ=255)
                                           ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                                   <  1:3165632/61‑1 (MQ=255)
                                           ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                         <  1:867695/71‑1 (MQ=255)
                                           ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                         <  1:1366373/71‑1 (MQ=255)
                                               cgtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt                     <  1:971606/71‑1 (MQ=255)
                                                gtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgtg                    >  1:3734586/1‑71 (MQ=255)
                                                 tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGAATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGc                   <  1:2590603/71‑1 (MQ=255)
                                                   gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGc                   >  1:1879522/1‑69 (MQ=255)
                                                    ttACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGAtg                      <  1:544849/65‑1 (MQ=255)
                                                     tACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt               >  1:1157883/1‑71 (MQ=255)
                                                      acggcgTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt               <  1:1208896/70‑1 (MQ=255)
                                                        gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGcc            <  1:2321248/71‑1 (MQ=255)
                                                         acgCGCCATCGGAGAAACTTCTGCCGGATATTCAg                                               <  1:405096/35‑1 (MQ=255)
                                                           gtgcCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCAct         <  1:3935828/69‑1 (MQ=255)
                                                               cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt     <  1:1288964/71‑1 (MQ=255)
                                                               cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt     >  1:3056262/1‑71 (MQ=255)
                                                               cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt     >  1:1575206/1‑71 (MQ=255)
                                                               cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt     <  1:2490657/71‑1 (MQ=255)
                                                                aGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACttcttc    <  1:2505579/71‑1 (MQ=255)
                                                                 gCGGGGACACTTCTGCCGTATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc  <  1:1798835/71‑2 (MQ=255)
                                                                 gCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc  <  1:1455906/71‑2 (MQ=255)
                                                                 gCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc    >  1:2304842/1‑70 (MQ=255)
                                                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                         <  1:3688637/49‑1 (MQ=255)
                                                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACttcttc    >  1:3209953/1‑70 (MQ=255)
                                                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg   <  1:2743481/71‑1 (MQ=255)
                                                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg   <  1:1282024/71‑1 (MQ=255)
                                                                  cGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCActt        <  1:1325994/65‑1 (MQ=255)
                                                                  cGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGa  <  1:3332325/71‑1 (MQ=255)
                                                                   ggAGAAACTTCTGCCGTATATTCTTTAATGAAGGTGGGCTGAATCAGAtg                      <  1:1280261/50‑1 (MQ=255)
                                                                     |                                                                   
ACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGA  >  minE/1933180‑1933316

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: