Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,276 G→A 9.7% F389F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2760GA 9.7% 216.5 / 10.8 72F389F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/5);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCC  >  minE/1933206‑1933345
                                                                      |                                                                     
ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                                                       <  1:3276287/71‑1 (MQ=255)
 gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                                                      <  1:804056/71‑1 (MQ=255)
    tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg                                                                   >  1:946646/1‑71 (MQ=255)
      tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                                                    <  1:1782671/68‑1 (MQ=255)
      tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                                                 >  1:3032730/1‑71 (MQ=255)
      tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                                                 >  1:597976/1‑71 (MQ=255)
      tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGt                                                                  <  1:1036006/70‑1 (MQ=255)
       gATTTCCGGTTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                                                                >  1:2912796/1‑71 (MQ=255)
        aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                                                  >  1:441733/1‑68 (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)
                       tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                                                  >  1:2519386/1‑53 (MQ=255)
                         gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGaa                                                                     >  1:2960271/1‑48 (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)
                              gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                                                    <  1:626717/44‑1 (MQ=255)
                              gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                                                    <  1:2994512/44‑1 (MQ=255)
                               acgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                                                  <  1:1586248/45‑1 (MQ=255)
                                cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                       >  1:1345709/1‑71 (MQ=255)
                                cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                       >  1:2553201/1‑71 (MQ=255)
                                     cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                                                      >  1:229937/1‑35 (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:265143/1‑54 (MQ=255)
                                          gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                                              >  1:1658654/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:1310008/1‑43 (MQ=255)
                                               aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                                                    >  1:1487304/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:1175725/71‑1 (MQ=255)
                                                       cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                <  1:3512688/71‑1 (MQ=255)
                                                       cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                <  1:2631428/71‑1 (MQ=255)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                >  1:2254113/1‑70 (MQ=255)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                >  1:2501391/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)
                                                                     tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATAc             <  1:2278648/60‑1 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGaa                              >  1:2384905/1‑42 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGaa                              >  1:4097567/1‑42 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                       >  1:764337/1‑49 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa              >  1:849269/1‑58 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa              >  1:2230538/1‑58 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc  >  1:3981425/1‑70 (MQ=255)
                                                                      gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc  >  1:3735313/1‑70 (MQ=255)
                                                                      |                                                                     
GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCC  >  minE/1933206‑1933345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: