Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,282 C→T 10.1% P387P (CCG→CCA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2820CT10.1% 207.3 / 11.1 71P387P (CCG→CCAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (37/25);  new base T (2/5);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGC  >  minE/1933212‑1933349
                                                                      |                                                                   
tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                                                     >  1:3032730/1‑71 (MQ=255)
tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                                                     >  1:597976/1‑71 (MQ=255)
 gATTTCCGGTTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                                                                    >  1:2912796/1‑71 (MQ=255)
  aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCgg                                                                   <  1:694078/71‑1 (MQ=255)
     tCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCt                                                                 >  1:812172/1‑70 (MQ=255)
        gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                                             <  1:548156/71‑1 (MQ=255)
            tAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                                             >  1:2835017/1‑67 (MQ=255)
             aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                                                        >  1:1902277/1‑71 (MQ=255)
               cggcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt                                                      <  1:595333/68‑1 (MQ=255)
                   gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                                                           <  1:114563/62‑1 (MQ=255)
                   gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                                  >  1:3825365/1‑71 (MQ=255)
                   gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                                  <  1:3072780/71‑1 (MQ=255)
                      acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgtg                                                     <  1:960515/65‑1 (MQ=255)
                      acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTg                                               >  1:4096231/1‑71 (MQ=255)
                       cgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTc                                                 >  1:3843282/1‑68 (MQ=255)
                          cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                           >  1:2553201/1‑71 (MQ=255)
                          cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                           >  1:1345709/1‑71 (MQ=255)
                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtg                                                       <  1:4143974/52‑1 (MQ=255)
                                 gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                    <  1:3014788/71‑1 (MQ=255)
                                   ggAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGa                                                         <  1:1030644/48‑1 (MQ=255)
                                    gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                                  >  1:1658654/1‑54 (MQ=255)
                                    gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                                  >  1:265143/1‑54 (MQ=255)
                                    gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCGTCGAAg                                 <  1:536310/71‑1 (MQ=255)
                                      gaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                    >  1:441372/1‑66 (MQ=255)
                                      gaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                    >  1:3434060/1‑66 (MQ=255)
                                        aaCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtc                              <  1:407845/70‑1 (MQ=255)
                                         aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                                                        >  1:1487304/1‑43 (MQ=255)
                                         aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                                                        >  1:1310008/1‑43 (MQ=255)
                                         aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtctc                            >  1:3128126/1‑71 (MQ=255)
                                          cTTCTGCCGGATATTCAGTAATGTAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                    <  1:780667/62‑1 (MQ=255)
                                           ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt                                                <  1:3255240/49‑1 (MQ=38)
                                            tCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTc                                                 >  1:817622/1‑47 (MQ=255)
                                            tCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTg                                               >  1:105183/1‑49 (MQ=255)
                                             cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                           >  1:159308/1‑52 (MQ=255)
                                             cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt                         >  1:2090698/1‑70 (MQ=255)
                                             cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCAAAGATCTCGGt                         >  1:746323/1‑70 (MQ=255)
                                                ccGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGAt                               >  1:3155007/1‑61 (MQ=255)
                                                ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAATAtc                              <  1:371207/62‑1 (MQ=255)
                                                ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAAc                      <  1:860749/70‑1 (MQ=255)
                                                ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAAc                      <  1:858309/70‑1 (MQ=255)
                                                 cGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCt                                                   <  1:680015/40‑1 (MQ=255)
                                                 cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    <  1:2631428/71‑1 (MQ=255)
                                                 cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    <  1:1175725/71‑1 (MQ=255)
                                                 cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    <  1:3512688/71‑1 (MQ=255)
                                                  ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    >  1:2501391/1‑70 (MQ=255)
                                                  ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    >  1:2254113/1‑70 (MQ=255)
                                                  ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGAt                   <  1:358778/71‑1 (MQ=255)
                                                    atatTCCGTAATAAAGGTAGGCTGAATCAGATGTGCTTCTGCCTCTTCATCAAAGATCTCTGTGACAATAc                 <  1:795829/71‑1 (MQ=255)
                                                       ttCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                                  <  1:3335926/35‑1 (MQ=255)
                                                       ttCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGc              >  1:2923334/1‑71 (MQ=255)
                                                         ccgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTcc             >  1:144352/3‑70 (MQ=255)
                                                         cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCa                                            >  1:3785069/1‑39 (MQ=255)
                                                          cgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCa           <  1:1707106/70‑1 (MQ=255)
                                                          aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    <  1:2093972/62‑1 (MQ=255)
                                                          aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                    <  1:2544326/62‑1 (MQ=255)
                                                          aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCa           <  1:3744619/71‑1 (MQ=255)
                                                            tAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGa         >  1:2526580/1‑71 (MQ=255)
                                                              aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                           >  1:1578023/1‑51 (MQ=255)
                                                              aTGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt                         <  1:2412059/53‑1 (MQ=255)
                                                              aTGAAGGTAGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGc              >  1:1472087/1‑64 (MQ=255)
                                                               taaaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAc                      <  1:4032444/53‑1 (MQ=255)
                                                               tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc                 <  1:2278648/60‑1 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGaa                                  >  1:4097567/1‑42 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGaa                                  >  1:2384905/1‑42 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                           >  1:764337/1‑49 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                  >  1:2230538/1‑58 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa                  >  1:849269/1‑58 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc      >  1:3981425/1‑70 (MQ=255)
                                                                gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc      >  1:3735313/1‑70 (MQ=255)
                                                                  aGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt                         <  1:223396/49‑1 (MQ=37)
                                                                    gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc                 >  1:2142712/1‑55 (MQ=255)
                                                                    gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc             >  1:3866125/1‑59 (MQ=255)
                                                                    gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGc  >  1:3615909/1‑70 (MQ=255)
                                                                      |                                                                   
TGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCCAGC  >  minE/1933212‑1933349

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: