Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,482,022 G→A 8.9% A281V (GCG→GTG)  mdtL ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,482,0220GA 8.9% 129.2 / 3.9 45A281V (GCG→GTG) mdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACCGGCGCAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCAAATG  >  minE/2481952‑2482088
                                                                      |                                                                  
aCCGGCGCAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCccg                                                                    >  1:128614/1‑71 (MQ=255)
 ccGGCGCAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCccgc                                                                   >  1:54908/1‑71 (MQ=255)
    gcgcAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCcgccgc                                                                >  1:1691522/1‑71 (MQ=255)
     cgcAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCgccgcc                                                               <  1:1490279/71‑1 (MQ=255)
         aaTCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                                            <  1:1088288/70‑1 (MQ=255)
         aaTCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAgg                                                            >  1:1171600/1‑70 (MQ=255)
           tCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAAt                                                         <  1:1851434/71‑1 (MQ=255)
             aGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGa                                                           >  1:193104/1‑67 (MQ=255)
             aGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAAt                                                         <  1:2054980/69‑1 (MQ=255)
              gCGTGATACCAAACAGAGAAACCGCATGGGAAGGTTAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAac                                                      <  1:561633/71‑1 (MQ=255)
                gTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAacac                                                    <  1:994597/71‑1 (MQ=255)
                gTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAacac                                                    <  1:441787/71‑1 (MQ=255)
                  gATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCACCGCCAGGAATaa                                                       >  1:1363891/1‑66 (MQ=255)
                    tACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGc                                                <  1:277846/71‑1 (MQ=255)
                     aCCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGaa                                                          >  1:1105783/1‑60 (MQ=255)
                      ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCg                                               <  1:1197677/70‑1 (MQ=255)
                      ccAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGa                                              <  1:484331/71‑1 (MQ=255)
                        aaaCAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAAt                                                         <  1:1388743/58‑1 (MQ=255)
                        aaaCAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAg                                             >  1:243865/1‑70 (MQ=255)
                         aaCAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGt                                           >  1:1177861/1‑71 (MQ=255)
                            agagAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCTAGGTGat                                        <  1:370610/71‑1 (MQ=255)
                            agagAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCACCGCCAGGAATAACACCTGCGAGGTGat                                        <  1:735030/71‑1 (MQ=255)
                             gagaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCACCGCCAGGAATAACACCTGCGAGGTGat                                        <  1:886916/70‑1 (MQ=255)
                               gaAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGatca                                      <  1:1945832/70‑1 (MQ=255)
                                  aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCaa                                  <  1:877213/71‑1 (MQ=255)
                                  aCCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCaa                                  <  1:967994/71‑1 (MQ=255)
                                     gCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGATTAACACCTGCGAGGTGATCATCAccgt                               <  1:2006328/71‑1 (MQ=255)
                                             aGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCt                       >  1:2117952/1‑71 (MQ=255)
                                              ggTGAAACGGCAAGAGTTATCCCCGCCGCCAGGAAt                                                         <  1:1240401/36‑1 (MQ=255)
                                                  aaaCGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGat                                        <  1:986165/49‑1 (MQ=255)
                                                    aCGGCAAGAGTGATCCCCACCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAAt                >  1:311566/1‑71 (MQ=255)
                                                       gCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTcc             >  1:114525/1‑71 (MQ=255)
                                                       gCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTcc             >  1:893996/1‑71 (MQ=255)
                                                       gCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTcc             >  1:1565349/1‑71 (MQ=255)
                                                          agagTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAg          >  1:2002625/1‑71 (MQ=255)
                                                           gagTGATCCCCGCCGCCAGTAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAgc         >  1:386020/1‑71 (MQ=255)
                                                            agTGATCCCCGCCGCCAGGAATAACACCTGCGAgg                                            >  1:427927/1‑35 (MQ=255)
                                                            agTGATCCCCGCCGCCAGGAATAACACCTGCGAgg                                            >  1:266233/1‑35 (MQ=255)
                                                            agTGATCCCCGCCGCCAGGAATAACACCTGCGAgg                                            >  1:804922/1‑35 (MQ=255)
                                                            agTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAgcg        <  1:589970/71‑1 (MQ=255)
                                                             gTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAgcg        >  1:1378462/1‑70 (MQ=255)
                                                             gTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAgcg        >  1:1924668/1‑70 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaaa    >  1:1213567/1‑71 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaaa    >  1:689687/1‑71 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaaa    <  1:2199174/71‑1 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaaa    <  1:1584309/71‑1 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaa     <  1:278120/70‑1 (MQ=255)
                                                                aTCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCaa     <  1:241557/70‑1 (MQ=255)
                                                                  ccccGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCATCGCAAATg  >  1:1531965/1‑71 (MQ=255)
                                                                    ccgccgCCAGGAATAACACCTCCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCa      <  1:1300793/65‑1 (MQ=255)
                                                                      |                                                                  
ACCGGCGCAAATCAGCGTGATACCAAACAGAGAAACCGCATGGGAAGGTGAAACGGCAAGAGTGATCCCCGCCGCCAGGAATAACACCTGCGAGGTGATCATCAACGTACGTGGCTTAAAAATTCCCAGCGCAAATG  >  minE/2481952‑2482088

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: