Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,236 G→C 16.7% A69A (GCC→GCG) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2360GC16.7% 183.3 / 26.1 80A69A (GCC→GCG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base C (3/10);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTC  >  minE/1934166‑1934304
                                                                      |                                                                    
cGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                                                                      <  1:3731146/71‑1 (MQ=255)
 gcaaaCGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACgg                                                                     >  1:1476048/4‑71 (MQ=255)
 gCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAAAgg                                                                     >  1:1291757/1‑71 (MQ=255)
 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)
                   aaCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACg                                                                      <  1:78977/52‑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:3849325/71‑1 (MQ=255)
                     cGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCt                                                 >  1:3248566/1‑71 (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:2433148/1‑50 (MQ=255)
                                ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGAt                                                           >  1:2104523/1‑50 (MQ=255)
                                ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                      >  1:3348784/1‑71 (MQ=255)
                                    ccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt                                       <  1:3335348/66‑1 (MQ=255)
                                     cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                            >  1:1452122/1‑42 (MQ=255)
                                     cATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                            >  1:2869328/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:3785940/1‑71 (MQ=255)
                                        aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                              >  1:3586116/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)
                                                                      |                                                                    
CGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTC  >  minE/1934166‑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: