Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,236 G→C 7.8% 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,2360GC 7.8% 269.0 / 10.1 90A69A (GCC→GCG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base C (4/3);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

CCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCA  >  minE/1934168‑1934305
                                                                    |                                                                     
ccAACGTCATGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGc                                                                     >  1:4095849/1‑71 (MQ=255)
     gTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACtt                                                                >  1:1098166/1‑71 (MQ=255)
      tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACttg                                                               <  1:1051536/71‑2 (MQ=255)
      tCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTc                                                               <  1:2166119/71‑1 (MQ=255)
      tCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTc                                                               <  1:437585/71‑1 (MQ=255)
      tCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTc                                                               <  1:1555851/71‑1 (MQ=255)
        cTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa                                                             >  1:1566708/1‑71 (MQ=255)
        cTGCAGGGTAACGAAAGACGCTTTACCCATAAGACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa                                                             >  1:3726255/1‑71 (MQ=255)
        cTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCGGCCAGCAACCGAGACTTcaa                                                             >  1:636344/1‑69 (MQ=255)
               gTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCa                                                       >  1:1977437/1‑70 (MQ=255)
                tAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATg                                                           <  1:857792/65‑1 (MQ=255)
                  aCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa                                                             >  1:3775650/1‑61 (MQ=255)
                  aCGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                         <  1:281911/65‑1 (MQ=255)
                     aaaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                                <  1:1108702/71‑1 (MQ=255)
                      aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACtt                                                                >  1:1983487/1‑54 (MQ=255)
                      aaGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                               <  1:2172686/71‑1 (MQ=255)
                       aGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                                >  1:1663461/1‑69 (MQ=255)
                        gACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAg                                             <  1:1960104/71‑1 (MQ=255)
                         aCGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAg                                                      >  1:840469/1‑61 (MQ=255)
                         aCGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGt                                            >  1:600392/1‑71 (MQ=255)
                          cGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGtt                                           <  1:3470108/71‑1 (MQ=255)
                          cGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGtt                                           <  1:2684230/71‑1 (MQ=255)
                            cTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATg                                                           <  1:779022/53‑1 (MQ=255)
                             tttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                         <  1:100945/54‑1 (MQ=255)
                             tttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttctt                                        >  1:1023094/1‑71 (MQ=255)
                              ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGa                                                             >  1:2187415/1‑49 (MQ=255)
                              ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGtt                                                         >  1:2649902/1‑53 (MQ=255)
                              ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                                <  1:2620888/62‑1 (MQ=255)
                              ttACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttcttc                                       >  1:2543666/1‑71 (MQ=255)
                                aCCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCg                                      <  1:1130741/70‑1 (MQ=255)
                                 cccATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCg                                      <  1:908885/69‑1 (MQ=255)
                                   cATGATACGACGGGTCATCATTCGGCCAGCAACCTAGACTTCAATGTTTAAGGATTCCAGTTCCTGGttat                                  >  1:3574618/1‑69 (MQ=255)
                                   cATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtct                                  <  1:3419247/71‑1 (MQ=255)
                                    ataataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAgcg                                                    >  1:2465860/1‑52 (MQ=255)
                                    ataataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                 >  1:2254108/1‑71 (MQ=255)
                                     tgataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTc                                                               <  1:3355678/38‑1 (MQ=255)
                                     taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATg                                                           <  1:1301943/44‑1 (MQ=255)
                                      aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                               >  1:3156672/1‑71 (MQ=255)
                                       ataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                              >  1:515733/1‑71 (MQ=255)
                                       ataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                              >  1:1575533/1‑71 (MQ=255)
                                       ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCg                                                              >  1:3043812/1‑39 (MQ=255)
                                       ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                              >  1:3368729/1‑71 (MQ=255)
                                       ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                              >  1:2479681/1‑71 (MQ=255)
                                       ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                              <  1:1308797/71‑1 (MQ=255)
                                       ataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTT‑‑AGTTCTTCGTTCTCTTTGc                            >  1:1943096/1‑71 (MQ=255)
                                           gacgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCg                          >  1:1101812/1‑71 (MQ=255)
                                            acgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGt                         <  1:1299553/71‑1 (MQ=255)
                                             cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGtt                                    <  1:1606635/59‑1 (MQ=255)
                                             cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTc                        <  1:1364761/71‑1 (MQ=255)
                                              gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                       <  1:2595525/71‑1 (MQ=255)
                                               gctcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                                <  1:1176970/43‑1 (MQ=255)
                                                gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                     <  1:3162616/71‑1 (MQ=255)
                                                  catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                   >  1:3651078/1‑71 (MQ=255)
                                                   atcatGCGGCCAGCAACGGCGACTTCGATGTTCAg                                                      >  1:2355558/1‑35 (MQ=255)
                                                    tcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGtt                                    <  1:2365464/52‑1 (MQ=255)
                                                    tcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGtt                                    <  1:718834/52‑1 (MQ=255)
                                                    tcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTc                  >  1:288390/1‑70 (MQ=255)
                                                    tcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                 <  1:3803520/71‑1 (MQ=255)
                                                    tcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                 <  1:3993051/71‑1 (MQ=255)
                                                     catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                     >  1:2687599/1‑66 (MQ=255)
                                                     catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                >  1:1345016/1‑71 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTcc                                               >  1:4091039/1‑39 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCa                                              >  1:1421734/1‑40 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtct                                  >  1:1736851/1‑52 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                       <  1:2176049/63‑1 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGa                      >  1:2482577/1‑64 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGa                      >  1:2994118/1‑64 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc               >  1:1994772/1‑71 (MQ=255)
                                                      atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc               >  1:15203/1‑71 (MQ=255)
                                                        tcGGCCAGCAACGGCGACTTCGATGTTCAGCGCATCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                <  1:712372/67‑1 (MQ=255)
                                                        gCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTc                                                >  1:1256312/1‑36 (MQ=255)
                                                        gCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt             >  1:38936/1‑71 (MQ=255)
                                                         cGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAg                                             >  1:365492/1‑38 (MQ=255)
                                                         cGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTg            >  1:2315813/1‑71 (MQ=255)
                                                          ggCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGc           <  1:3742653/71‑1 (MQ=255)
                                                          ggCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGc           <  1:313067/71‑1 (MQ=255)
                                                           gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGttc                                          <  1:1138750/39‑1 (MQ=255)
                                                           gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGc           <  1:3659099/70‑1 (MQ=255)
                                                            ccAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCaa         >  1:3955698/1‑71 (MQ=255)
                                                             cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                 <  1:2559555/46‑1 (MQ=255)
                                                             cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCaa         <  1:182546/70‑1 (MQ=255)
                                                             cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt        <  1:2295899/71‑1 (MQ=255)
                                                              aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                 <  1:1178894/45‑1 (MQ=255)
                                                              aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTc                  >  1:848756/1‑60 (MQ=255)
                                                                  aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                 <  1:1142122/41‑1 (MQ=255)
                                                                  aCGGCGACTTCGATGTACAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                >  1:2168962/1‑58 (MQ=255)
                                                                   cGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                 <  1:1212803/40‑1 (MQ=255)
                                                                   cGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa  <  1:3160453/71‑1 (MQ=255)
                                                                   cGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa  <  1:81045/71‑1 (MQ=255)
                                                                   cGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa  <  1:1081731/71‑1 (MQ=255)
                                                                    ggCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                   >  1:3320817/1‑53 (MQ=255)
                                                                    ggCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                   >  1:3316134/1‑53 (MQ=255)
                                                                    |                                                                     
CCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCA  >  minE/1934168‑1934305

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: