Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,249 A→G 12.7% 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,2490AG12.7% 148.1 / 12.9 55S398S (TCT→TCClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (3/4);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAG  >  minE/1933181‑1933318
                                                                    |                                                                     
cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                                     <  1:2515502/71‑1 (MQ=255)
  accaccGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGgggaca                                                                   >  1:3168460/1‑69 (MQ=255)
   ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                                  <  1:2176779/71‑1 (MQ=255)
    caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                 <  1:3398231/71‑1 (MQ=255)
    caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                 <  1:188691/71‑1 (MQ=255)
     accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                                >  1:1417476/1‑71 (MQ=255)
        aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTg                                                             >  1:3490247/1‑71 (MQ=255)
          tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                 <  1:517299/65‑1 (MQ=255)
            aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTc                                                               >  1:821439/1‑65 (MQ=255)
            aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                                          <  1:3500476/70‑1 (MQ=255)
             agaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                        <  1:900427/71‑1 (MQ=255)
                 cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                        <  1:3607278/67‑1 (MQ=255)
                    aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGata                                                      <  1:2608841/66‑1 (MQ=255)
                    aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAg                                                 >  1:592188/1‑71 (MQ=255)
                     aaacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGt                                                <  1:2930847/68‑1 (MQ=255)
                         ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                            <  1:3276287/71‑1 (MQ=255)
                          gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                           <  1:804056/71‑1 (MQ=255)
                             tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                                     <  1:4118789/42‑1 (MQ=255)
                             tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                        <  1:1018364/55‑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)
                               tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGacac                                                                  <  1:1363723/43‑1 (MQ=16)
                               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)
                                 aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACtt                                                                >  1:1314136/1‑43 (MQ=16)
                                  tttCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                                    <  1:2840527/54‑1 (MQ=255)
                                    tCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                                     <  1:2955099/35‑1 (MQ=255)
                                    tCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCt                                  >  1:812172/1‑70 (MQ=255)
                                       gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                              <  1:548156/71‑1 (MQ=255)
                                           tAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                              >  1:2835017/1‑67 (MQ=255)
                                           tAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGat                                                       >  1:3854578/1‑42 (MQ=16)
                                            aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCa                                                  >  1:2515140/1‑46 (MQ=255)
                                            aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt                         >  1:1902277/1‑71 (MQ=255)
                                              cggcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt                       <  1:595333/68‑1 (MQ=255)
                                                tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                        <  1:1600906/36‑1 (MQ=255)
                                                tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGt                                                <  1:1739025/44‑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)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTcc                                                  >  1:3453221/2‑39 (MQ=19)
                                                     acgacgCGCCAGCGGAGAAACTTCTGCCGGATATTCa                                                  >  1:194979/1‑37 (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:2994512/44‑1 (MQ=255)
                                                       gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                                         <  1:626717/44‑1 (MQ=255)
                                                        acgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                       <  1:1586248/45‑1 (MQ=255)
                                                         cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc            >  1:2553201/1‑71 (MQ=255)
                                                         cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc            >  1:1345709/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)
                                                                    |                                                                     
CGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAG  >  minE/1933181‑1933318

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: