Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 387,871 G→T 8.3% S175I (AGT→ATT)  mepA → MATE family multi‑antimicrobial extrusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730387,8710GT 8.3% 172.3 / 5.0 48S175I (AGT→ATT) mepAMATE family multi‑antimicrobial extrusion protein
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base T (2/2);  total (25/23)
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.10e-01

GGTGAAACGTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTGA  >  CP000730/387732‑388009
                                                                                                                                           |                                                                                                                                          
ggTGAAACGTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAgt                                                                                                                                           >  2:306805/1‑141 (MQ=255)
      aCGTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGc                                                                                                                                     <  2:115858/141‑1 (MQ=255)
      aCGTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGc                                                                                                                                     <  2:249820/141‑1 (MQ=255)
        gTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTATTGTAGGCtt                                                                                                                                   >  1:330924/1‑141 (MQ=255)
        gTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTATTGTAGGCtt                                                                                                                                   >  1:330832/1‑141 (MQ=255)
          tAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTaa                                                                                                                                 <  1:442896/141‑1 (MQ=255)
               ttAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAtatt                                                                                                                            >  1:189581/1‑141 (MQ=255)
               ttAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAtatt                                                                                                                            >  2:457422/1‑141 (MQ=255)
               ttAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAtatt                                                                                                                            >  2:350436/1‑141 (MQ=255)
               ttAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAtatt                                                                                                                            >  1:424919/1‑141 (MQ=255)
               ttAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAtatt                                                                                                                            >  1:423833/1‑141 (MQ=255)
                  aCAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATattatt                                                                                                                         >  2:199685/1‑141 (MQ=255)
                  aCAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATattatt                                                                                                                         >  2:254189/1‑141 (MQ=255)
                      gTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAg                                                                                                                     >  2:146387/1‑141 (MQ=255)
                       gaaTTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGa                                                                                                                    >  1:436423/2‑141 (MQ=255)
                                gAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCaatttt                                                                                                           <  2:379705/141‑1 (MQ=255)
                                gAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCaatttt                                                                                                           <  2:165670/141‑1 (MQ=255)
                                gAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCCTCATATTAGAACAATTTGCGCGTGCAATTGGGGCACCGATGGGTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAAGATTATGTTAGATCCaattct                                                                                                           <  2:15352/141‑3 (MQ=255)
                                  aaGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCaattttaa                                                                                                         <  1:328421/141‑1 (MQ=255)
                                     tAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCaattttaattt                                                                                                      <  2:81351/141‑1 (MQ=255)
                                     tAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCaattttaattt                                                                                                      <  2:423903/141‑1 (MQ=255)
                                            tttAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGttt                                                                                               <  1:430294/141‑1 (MQ=255)
                                                    cACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTaa                                                                                       <  2:316976/141‑1 (MQ=255)
                                                          ttGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACgttg                                                                                 >  1:399622/1‑141 (MQ=255)
                                                          ttGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACgttg                                                                                 >  2:17031/1‑141 (MQ=255)
                                                           tGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACgttgt                                                                                >  1:186734/1‑141 (MQ=255)
                                                           tGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTA‑ATTCAATTTTTATTTTTGGGTTTTGTTTTAACgttttt                                                                               >  2:15967/1‑141 (MQ=255)
                                                            gTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTATTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACgttgtt                                                                               <  2:330924/141‑1 (MQ=255)
                                                            gTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTATTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACgttgtt                                                                               <  2:330832/141‑1 (MQ=255)
                                                              aaTTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTgg                                                                             <  2:375948/141‑1 (MQ=255)
                                                                               tAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTg                                                            >  1:247030/1‑141 (MQ=255)
                                                                                  aaCAATTTGCACGTGCAATTGGGGCACCAACGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCaa                                                         >  2:414075/1‑141 (MQ=255)
                                                                                            aCGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGt                                               <  2:293171/141‑1 (MQ=255)
                                                                                                          cACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGt                                 <  2:281068/141‑1 (MQ=255)
                                                                                                                 ggTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTAt                          >  1:176800/1‑141 (MQ=255)
                                                                                                                     tCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCatt                      >  2:435861/1‑141 (MQ=255)
                                                                                                                     tCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCatt                      >  1:87711/1‑141 (MQ=255)
                                                                                                                            ttGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTTGTTTTGTTTTTAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCAtttattttt               >  2:55445/1‑140 (MQ=255)
                                                                                                                              ggTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATg             <  1:17031/141‑1 (MQ=255)
                                                                                                                                  tGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGaaaa         <  1:199685/141‑1 (MQ=255)
                                                                                                                                  tGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGaaaa         <  1:254189/141‑1 (MQ=255)
                                                                                                                                    ttAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAAt       >  2:456237/1‑141 (MQ=255)
                                                                                                                                    ttAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGATTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAAt       >  1:135101/1‑141 (MQ=255)
                                                                                                                                      aGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAg     <  1:403686/141‑1 (MQ=255)
                                                                                                                                       gCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGt    <  2:399622/141‑1 (MQ=255)
                                                                                                                                        cTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTg   <  1:69506/141‑1 (MQ=255)
                                                                                                                                        cTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTg   <  1:86854/141‑1 (MQ=255)
                                                                                                                                         tAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTGa  >  1:300772/1‑141 (MQ=255)
                                                                                                                                         tAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTGa  >  1:193793/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
GGTGAAACGTTAGCTTTAACAAGTAATTATTTGAAAGTAATGTTTTTAAGTGCACCTTTTGTAATTTTGTTCTTCATATTAGAACAATTTGCACGTGCAATTGGGGCACCAATGGTTTCTATGATTGGTATGTTAGCTAGTGTAGGCTTAAATATTATTTTAGATCCAATTTTAATTTTTGGTTTTGATTTAAACGTTGTTGGTGCAGCTTTGGGTACTGCAATCAGTAATGTTGCTGCTGCTCTGTTCTTTATCATTTATTTTATGAAAAATAGTGA  >  CP000730/387732‑388009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑