Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,482,012 A→G 35.2% T284T (ACT→ACC mdtL ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,482,0120AG35.2% 56.4 / 51.8 88T284T (ACT→ACCmdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (36/21);  new base G (12/19);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 4.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAA  >  minE/2481971‑2482071
                                         |                                                           
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAacac                                   >  1:1594589/1‑68 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:684750/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:382814/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:2483850/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:1956363/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:2049049/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:408172/1‑69 (MQ=255)
aTACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACAcc                                  >  1:1210928/1‑69 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                           >  1:1568999/1‑57 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                           >  1:2320272/1‑57 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                           >  1:840390/1‑57 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                           >  1:3141860/1‑57 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                           >  1:1066143/1‑57 (MQ=255)
   ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATACCCGCCGCCAgg                                           >  1:2273068/1‑57 (MQ=255)
   caaaaCAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCccgc                                                  >  1:665068/3‑50 (MQ=255)
        cAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTg                                <  1:469513/63‑1 (MQ=255)
        cAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTg                                <  1:2946263/63‑1 (MQ=255)
        cAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTg                                <  1:2138426/63‑1 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGAt                                                        >  1:1850567/1‑36 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:394961/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:2114005/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:2999288/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:9460/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:1072154/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:859670/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:1306544/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:731080/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                       >  1:1278551/1‑69 (MQ=255)
           agaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGa                        >  1:308374/1‑68 (MQ=255)
            gaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatc                      <  1:2614429/69‑1 (MQ=255)
            gaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatc                      <  1:3434273/69‑1 (MQ=255)
               aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatcatc                   >  1:564543/1‑69 (MQ=255)
               aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatcatc                   >  1:2496813/1‑69 (MQ=255)
               aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatcatc                   >  1:381238/1‑69 (MQ=255)
               aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatcatc                   >  1:48690/1‑69 (MQ=255)
                   cATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacg               <  1:1882033/69‑1 (MQ=255)
                   cATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacg               <  1:748532/69‑1 (MQ=255)
                   cATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacg               <  1:1401925/69‑1 (MQ=255)
                     tGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgta             <  1:1286640/69‑1 (MQ=255)
                      gggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgta             <  1:722509/68‑1 (MQ=255)
                      gggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgta             <  1:1205902/68‑1 (MQ=255)
                      gggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgta             <  1:1461019/68‑1 (MQ=255)
                       ggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatca                     >  1:2425619/1‑59 (MQ=255)
                       ggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgtacg           >  1:110173/1‑69 (MQ=255)
                       ggAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAacgtac            >  1:3077292/1‑68 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:979822/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2888367/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2833445/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:329237/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:720824/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:648878/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2688990/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2687913/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2458419/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2174230/1‑69 (MQ=255)
                         aaGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2076200/1‑69 (MQ=255)
                         aaGGTGAAACGGAAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGATCATCAACGTACGTg         >  1:2729063/1‑69 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCTCCAGGAATAACACCTGTGAGGTGat                       <  1:1351615/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1542561/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1525164/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1456763/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1328397/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:119141/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:782103/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1830842/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:203581/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:3428771/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1080847/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:1060977/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:959204/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:3090771/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:2763644/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:2685738/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCGCCAGGAATAACACCTGTGAGGTGat                       <  1:2978668/54‑1 (MQ=255)
                          aGGTGAAACGGCAAGGGTGATCCCCGCCACCAGGAATAACACCTGTGAGGTGat                       <  1:1035185/54‑1 (MQ=255)
                          aGGGGAAACGGCAAGGCTGATCCCCGCCGCCAGGAATAACACCTGTGAGGGGat                       <  1:2684112/54‑1 (MQ=25)
                            gTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGc       >  1:1545318/1‑68 (MQ=255)
                            gTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGc       >  1:1778963/1‑68 (MQ=255)
                            gTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGc       >  1:1144929/1‑68 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTTATCATCAACGTACGTGGCTTaaa  <  1:3321248/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:477424/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:1608461/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:412428/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:12370/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:1200600/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:1714474/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:3374453/68‑1 (MQ=255)
                                 aCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTaaa  <  1:3274504/68‑1 (MQ=255)
                                         |                                                           
ATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAA  >  minE/2481971‑2482071

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: