Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 559,815 A→G 61.8% G354G (GGA→GGG cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE559,8150AG61.8% 53.9 / 59.7 94G354G (GGA→GGGcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base G (15/43);  total (36/58)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGC  >  minE/559748‑559881
                                                                   |                                                                  
cGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGaa                                                                   <  1:1735954/69‑1 (MQ=255)
       ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCa                                                            <  1:424096/69‑1 (MQ=255)
       ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCa                                                            <  1:2635091/69‑1 (MQ=255)
       ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCa                                                            <  1:582350/69‑1 (MQ=255)
       ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCa                                                            <  1:1257728/69‑1 (MQ=255)
       ggTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCa                                                            <  1:1967783/69‑1 (MQ=255)
           cGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATg                                                        >  1:2265642/1‑69 (MQ=255)
           cGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATg                                                        >  1:2411998/1‑69 (MQ=255)
           cGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATg                                                        >  1:1173901/1‑69 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGGAAATGCTGATCttt                                               >  1:2188929/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:783187/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:1306605/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:225920/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:3261027/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:2993431/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCttt                                               >  1:2860506/1‑68 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCtt                                                >  1:2891221/1‑67 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCtt                                                >  1:2879580/1‑67 (MQ=255)
                     tGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCtt                                                >  1:3396550/1‑67 (MQ=255)
                     tGTTTGCAAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCtt                                                >  1:1966056/1‑67 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGt                                      >  1:144186/1‑68 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGt                                      >  1:2264656/1‑68 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:2735010/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:989676/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:3282427/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:2226101/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:500124/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:2108126/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:849286/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:192464/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:2039138/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:1384995/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:1833235/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGTACGGTTTATGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:1645494/1‑69 (MQ=255)
                              gCGATATGAGTAAAGGAACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTa                                     >  1:3096016/1‑69 (MQ=255)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:358454/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:2827355/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:2668597/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:178825/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:584630/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:2461550/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:2442078/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:613982/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:1214162/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:3010343/44‑1 (MQ=38)
                                          aaGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCa                         <  1:2928285/69‑1 (MQ=255)
                                          aaGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCa                         <  1:3408259/69‑1 (MQ=255)
                                          aaGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATCCTGATCTTTACCGTTGGTATTGAAATCAGCa                         <  1:902227/69‑1 (MQ=255)
                                           aGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:1796736/43‑1 (MQ=38)
                                           aGGTACGGTTTCTGCGGCGATGGGGATGCTGCAAATGCTGATc                                                  <  1:1889533/43‑1 (MQ=38)
                                            ggTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCaaa                       <  1:2530259/69‑1 (MQ=255)
                                            ggTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCaaa                       <  1:1344951/69‑1 (MQ=255)
                                             gTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGa                                 >  1:2959808/1‑58 (MQ=255)
                                             gTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGa                                 >  1:2877121/1‑58 (MQ=255)
                                             gTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGa                                 >  1:1097881/1‑58 (MQ=255)
                                             gTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGa                                 >  1:764325/1‑58 (MQ=255)
                                             gTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGa                                 >  1:2201579/1‑58 (MQ=255)
                                              tACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATc                             >  1:1526301/1‑61 (MQ=255)
                                              tACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACa                     >  1:3367192/1‑69 (MQ=255)
                                               aCGGTTTCTGCCGCGATGGGAAGGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACAt                    <  1:2762531/69‑1 (MQ=255)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:424454/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:43417/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:717454/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:724569/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:907064/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:2999204/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:3311067/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:1444577/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:205025/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:2077308/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:2289678/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:1765040/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:2844028/69‑1 (MQ=39)
                                                    ttCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:144238/69‑1 (MQ=39)
                                                     tCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:3253842/68‑1 (MQ=39)
                                                     tCTGCGGCGATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:2146368/68‑1 (MQ=39)
                                                     tCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGGTGGTATTGAAATCAGCAAACAt                    <  1:1639229/63‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTg               <  1:932955/62‑1 (MQ=39)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:1287545/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:1744148/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:988071/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:2711175/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:842780/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:2153979/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:1826653/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:1815064/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:2987115/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:595036/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:2543056/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:1273804/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:558748/69‑1 (MQ=255)
                                                           gcgATGGGGATGCTGCAAATGCTGATCTTTACCGTCGGTATTGAAATCAGCAAACATGCCTGGCTGAAc        <  1:487996/69‑1 (MQ=255)
                                                                 ggAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAAGGGGGGc  <  1:438968/69‑1 (MQ=255)
                                                                 ggAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGc  <  1:1844479/69‑1 (MQ=255)
                                                                 ggAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGc  <  1:1125500/69‑1 (MQ=255)
                                                                   |                                                                  
CGGGACTGGTGCGATTAACCCTGTTTGCCAGCGATATGAGTAAAGGTACGGTTTCTGCCGCGATGGGAATGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGC  >  minE/559748‑559881

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: