Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,238 C→A 19.4% A69S (GCC→TCC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2380CA19.4% 202.1 / 26.7 77A69S (GCC→TCC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base A (4/11);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTC  >  minE/1934167‑1934304
                                                                       |                                                                  
gCCAACGTCCTGCAGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                                    <  1:1437688/71‑1 (MQ=255)
 ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                                    >  1:3403473/1‑71 (MQ=255)
 ccAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                                    <  1:2043395/71‑1 (MQ=255)
       tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                                              <  1:2322979/71‑1 (MQ=255)
       tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                                              <  1:2141904/71‑1 (MQ=255)
       tCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTc                                                              <  1:2844694/71‑1 (MQ=255)
         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                                             <  1:1274658/70‑1 (MQ=255)
         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                                             <  1:2069430/70‑1 (MQ=255)
         cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa                                                            <  1:2120878/71‑1 (MQ=255)
         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca                                                             <  1:3842180/70‑2 (MQ=255)
         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTca                                                             <  1:1776962/70‑2 (MQ=255)
         cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACAGAGACTTcaa                                                            >  1:1040718/1‑69 (MQ=255)
          tGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACtt                                                               <  1:1167715/67‑1 (MQ=255)
             aGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                        >  1:598419/1‑71 (MQ=255)
             aGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                        >  1:258544/1‑71 (MQ=255)
               ggTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa                                                      <  1:3160277/71‑1 (MQ=255)
                 tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                                             <  1:1337919/62‑1 (MQ=255)
                  aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAgcg                                                   >  1:3533358/1‑71 (MQ=255)
                   aCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                                                           >  1:3665650/1‑62 (MQ=255)
                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                        >  1:696461/1‑64 (MQ=255)
                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa                                                      >  1:3529402/1‑66 (MQ=255)
                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt                                                 >  1:3248566/1‑71 (MQ=255)
                    cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt                                                 <  1:3849325/71‑1 (MQ=255)
                       aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGAc                                                                 >  1:1816200/1‑52 (MQ=255)
                       aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                               >  1:66336/1‑70 (MQ=255)
                         gaggCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAg                                            >  1:2529101/4‑71 (MQ=255)
                            gCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttc                                         <  1:3302059/71‑1 (MQ=255)
                            gCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttc                                         <  1:3534421/71‑1 (MQ=255)
                             cTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttct                                        <  1:3938307/71‑1 (MQ=255)
                               ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                                                           >  1:2104523/1‑50 (MQ=255)
                               ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                                                           >  1:2433148/1‑50 (MQ=255)
                               ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                      >  1:3348784/1‑71 (MQ=255)
                                   ccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt                                       <  1:3335348/66‑1 (MQ=255)
                                    cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                            >  1:2869328/1‑42 (MQ=255)
                                    cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                            >  1:1452122/1‑42 (MQ=255)
                                    cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                                              <  1:2128789/42‑1 (MQ=255)
                                    cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                      <  1:1460137/66‑1 (MQ=255)
                                      tgataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcc                               <  1:2818886/69‑2 (MQ=255)
                                      taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                               <  1:425590/71‑1 (MQ=255)
                                      taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                               <  1:2329692/71‑1 (MQ=255)
                                      taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                               <  1:170121/71‑1 (MQ=255)
                                       aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              >  1:37839/1‑71 (MQ=255)
                                       aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              >  1:3586116/1‑71 (MQ=255)
                                       aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              >  1:3785940/1‑71 (MQ=255)
                                        ataCGACGGGTCATCATTCGGCCAGCAACCGAGCCTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcct                              <  1:860711/70‑3 (MQ=255)
                                        ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa                                                      >  1:372350/1‑46 (MQ=255)
                                         taCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                             <  1:2698314/70‑1 (MQ=255)
                                         taCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              <  1:190178/69‑1 (MQ=255)
                                          acgacgGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttcct                                       <  1:27264/59‑3 (MQ=255)
                                          acgacgGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttcct                                       <  1:461225/59‑3 (MQ=255)
                                            gacgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                              <  1:2382506/50‑1 (MQ=255)
                                             acgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              <  1:760319/65‑1 (MQ=255)
                                             acgGGTCATCATGCGGCCAGCAACGGAGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              <  1:2024128/65‑1 (MQ=255)
                                              cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                             >  1:3508625/1‑65 (MQ=255)
                                               gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                      <  1:3623730/71‑1 (MQ=255)
                                                 gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                    <  1:1904250/71‑1 (MQ=255)
                                                  tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                >  1:509796/1‑58 (MQ=255)
                                                  tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                   >  1:2079312/1‑71 (MQ=255)
                                                   catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGt                                    >  1:3870703/1‑53 (MQ=255)
                                                   catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                   >  1:2449958/1‑70 (MQ=255)
                                                    atcatTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGttc                                         <  1:1635228/47‑1 (MQ=255)
                                                    atcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                               <  1:2672363/57‑1 (MQ=255)
                                                      catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                             >  1:1086959/1‑57 (MQ=255)
                                                      catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                <  1:3895473/70‑1 (MQ=255)
                                                      catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg               >  1:1004369/1‑71 (MQ=255)
                                                       atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc              >  1:3580720/1‑71 (MQ=255)
                                                        tGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                              <  1:475559/38‑1 (MQ=255)
                                                            gCCAGCAACGGCGACTTCTATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa         <  1:3632756/71‑1 (MQ=255)
                                                            gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                    >  1:430752/1‑60 (MQ=255)
                                                             ccAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                   <  1:590094/60‑1 (MQ=255)
                                                              cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt       >  1:797984/1‑71 (MQ=255)
                                                               aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                  >  1:989762/1‑59 (MQ=255)
                                                               aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt      >  1:56651/1‑71 (MQ=255)
                                                               aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt      <  1:2212748/71‑1 (MQ=255)
                                                                gCAACGGCGACTTCGATGTTCAGTGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg     <  1:3438780/71‑1 (MQ=255)
                                                                gCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg     <  1:994345/71‑1 (MQ=255)
                                                                 cAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                    >  1:1360918/1‑55 (MQ=255)
                                                                   aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa         <  1:3577816/64‑1 (MQ=255)
                                                                   aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg     <  1:2617371/68‑1 (MQ=255)
                                                                   aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTc  <  1:2290633/71‑1 (MQ=255)
                                                                       |                                                                  
GCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTC  >  minE/1934167‑1934304

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: