Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,645,851 A→G 35.4% G58G (GGT→GGC yhgN ← predicted antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,645,8510AG35.4% 60.8 / 57.1 96G58G (GGT→GGCyhgNpredicted antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base G (12/22);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

ATGCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATG  >  minE/2645792‑2645908
                                                           |                                                         
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  <  1:2668077/69‑1 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:756534/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  <  1:1108127/69‑1 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:1110938/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  <  1:3385173/69‑1 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:1183369/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:120244/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:598817/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:2906636/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:1390350/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  <  1:1423432/69‑1 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:2726315/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:1445360/1‑69 (MQ=255)
atgCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAAc                                                  >  1:464730/1‑69 (MQ=255)
   ccgccAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               <  1:1327397/69‑1 (MQ=255)
      ccAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAac                                            <  1:1103612/69‑1 (MQ=255)
           aaTGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:596374/1‑61 (MQ=255)
           aaTGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2713737/1‑61 (MQ=255)
           aaTGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:446383/1‑61 (MQ=255)
           aaTGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2323675/1‑61 (MQ=255)
             tGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGc                                      >  1:2804111/1‑68 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:481715/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:476270/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2616089/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:731575/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:230365/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2661911/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:645919/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2080244/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2838308/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:1522572/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAgg                                               >  1:2004068/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACACGCAAACAgg                                               >  1:595282/1‑57 (MQ=255)
               gagaCGGTTTCTGCTCGTAGGCTAAGAAATGCAGAATTTTCTCACCCGCAAACAgg                                               >  1:652623/1‑56 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGTCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:2498373/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAgg                                               >  1:769971/1‑52 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGa                                              >  1:1744944/1‑53 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:1784227/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:278358/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:2381918/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:589539/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:3290474/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:3321459/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:2731781/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:509656/1‑66 (MQ=255)
                    ggTTTCTGCCCGCAGGCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCAca                                >  1:2772214/1‑66 (MQ=255)
                                   gCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:286491/39‑1 (MQ=255)
                                   gCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:2796386/39‑1 (MQ=255)
                                   gCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:2390920/39‑1 (MQ=255)
                                   gCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:1757102/39‑1 (MQ=255)
                                   gCTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:1431632/39‑1 (MQ=255)
                                    cTAAGAAATGCCAGAATTTTCTCGCCCGCAAACAGGaa                                             <  1:2191594/38‑1 (MQ=255)
                                    cTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAAc              <  1:2870715/69‑1 (MQ=255)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:3023577/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:3059658/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:375713/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:369667/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:347122/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:521450/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:525392/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:2417073/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:2255677/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:2151246/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:1803429/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:1763779/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:1560081/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:867879/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:1117783/51‑1 (MQ=38)
                                      aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                              <  1:1003002/51‑1 (MQ=38)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCTCAc      <  1:1588184/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:973679/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      >  1:231812/1‑69 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      >  1:1053818/1‑69 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:2121346/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:2125378/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:2329313/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:2536556/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:3398187/69‑1 (MQ=255)
                                            tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc      <  1:2722648/69‑1 (MQ=255)
                                                ggAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:2058965/68‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:623404/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:872645/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:868751/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:1214504/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:1632783/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:1764084/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:1897031/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:3323271/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:606147/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:3406497/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:3101180/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:244963/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:2964101/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:2627273/69‑1 (MQ=255)
                                                aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:2648695/69‑1 (MQ=255)
                                                 gAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg  <  1:3195410/68‑1 (MQ=255)
                                                           |                                                         
ATGCCGCCAGAAATGGAGACGGTTTCTGCTCGTAGGCTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATG  >  minE/2645792‑2645908

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