Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,225 C→A 31.1% A216A (GCG→GCT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,2250CA31.1% 90.7 / 44.8 77A216A (GCG→GCTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (45/8);  new base A (11/13);  total (56/21)
Fisher's exact test for biased strand distribution p-value = 7.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-03

TGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGT  >  minE/312163‑312292
                                                              |                                                                   
tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCa                                                               >  1:1139815/1‑69 (MQ=255)
tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCa                                                               >  1:2364843/1‑69 (MQ=255)
   gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc                                                            >  1:372151/1‑69 (MQ=255)
   gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc                                                            >  1:2975600/1‑69 (MQ=255)
   gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc                                                            >  1:725117/1‑69 (MQ=255)
   gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAAcc                                                            >  1:2028324/1‑69 (MQ=255)
         ataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCaa                                                              >  1:2658836/1‑61 (MQ=255)
            gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc                                                   >  1:430561/1‑69 (MQ=255)
            gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc                                                   >  1:3264492/1‑69 (MQ=255)
            gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTc                                                   >  1:2786224/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCt                                                  >  1:2250333/1‑68 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:2259199/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:2164412/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:1431834/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:3295229/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:2181951/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:835399/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:2696083/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:484426/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:1320212/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:2878777/1‑69 (MQ=255)
              ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTg                                                 >  1:1499852/1‑69 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTTACCCGCCGCAACCTGGGCGTTCTGCGCttt                                           <  1:809402/69‑1 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCACCCTGGGCGTTCTGCGCttt                                           >  1:1705023/1‑69 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt                                           <  1:2277039/69‑1 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt                                           >  1:1302374/1‑69 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt                                           <  1:2255684/69‑1 (MQ=255)
                    aaGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCttt                                           >  1:2232589/1‑69 (MQ=255)
                        tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg                                       >  1:3118144/1‑69 (MQ=255)
                        tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg                                       >  1:1897063/1‑69 (MQ=255)
                        tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg                                       >  1:600790/1‑69 (MQ=255)
                        tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATg                                       >  1:3252444/1‑69 (MQ=255)
                             gCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  >  1:3324381/1‑69 (MQ=255)
                             gCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  >  1:2686746/1‑69 (MQ=255)
                                tttCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAAt                               >  1:2842782/1‑69 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTc                                                   >  1:2855617/1‑45 (MQ=38)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:2105950/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:889310/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:676711/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:1694619/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:696110/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:2452437/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:1816651/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:1844992/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:1248518/1‑68 (MQ=255)
                                    aCCGGCGGCGTACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTAATGAc                            >  1:3344012/1‑68 (MQ=255)
                                       ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa                        >  1:1422177/1‑69 (MQ=255)
                                       ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCa                        >  1:2179174/1‑69 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:3185091/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:1478262/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:603127/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:101998/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:2941744/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:1410470/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:2375443/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:2178183/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:2112361/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:145919/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:1957006/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:1604366/53‑1 (MQ=255)
                                              tACCACCGAGCTGACCAGCTGCAACCTGAGCGTTCTGTGCTTTGATGGCGGTa                                 <  1:863649/53‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  <  1:3017226/50‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  <  1:2763983/50‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  <  1:1736255/50‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  <  1:3425509/50‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGt                                  <  1:1647967/50‑1 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTaa                                >  1:3361060/1‑52 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:89342/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:366514/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:3324596/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:3138870/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:2531617/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:216625/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:1992864/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGtt               >  1:1114657/1‑69 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGt                >  1:1848539/1‑68 (MQ=255)
                                                ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGt                >  1:1507303/1‑68 (MQ=255)
                                                             ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt  <  1:455265/69‑1 (MQ=255)
                                                             ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt  <  1:3401244/69‑1 (MQ=255)
                                                             ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt  <  1:1071213/69‑1 (MQ=255)
                                                              |                                                                   
TGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGT  >  minE/312163‑312292

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: