Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,509,782 G→T 12.2% V454L (GTA→TTA)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,509,7820G→T12.2% 2.8 ‑4.4 ‑0.1331360acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (8/21):  total (12/21)
Fisher's exact test strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.11e-01

TGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGC  >  REL606/2509747‑2509816
                                   |                                  
tGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTg                                    <  1:6585062/36‑1 (MQ=255)
       gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                             <  1:766295/36‑1 (MQ=255)
       gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                             <  1:6537706/36‑1 (MQ=255)
       gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                             <  1:6128479/36‑1 (MQ=255)
        aTTGCGATGGTTCTTTCGGCGGTATTTGTACCAATg                            <  1:810727/36‑1 (MQ=255)
        aTTGCGATGGTTCTTTCGGCGGTATTTGTACCAATg                            <  1:6978742/36‑1 (MQ=255)
         ttGCGATGGTTCTTTCGGCGGTATTTTTTCCAAttg                           >  1:6176863/1‑34 (MQ=21)
         ttGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                           <  1:2906457/36‑1 (MQ=255)
         ttGCGATGGTTCTTTCGGCGGTAGTTGTACCaagtt                           >  1:3758325/1‑33 (MQ=25)
           gCGATGGTTCTTTCGGCGGTATTTGTACCAATGGcc                         <  1:7470044/36‑1 (MQ=255)
           gCGATGGTTCTTTCGGCGGTATTTGTACCAATGGcc                         <  1:2629623/36‑1 (MQ=255)
           gCGATGGTACTGTCGGCGGTATTTGTATCCATgtgc                         >  1:2803219/1‑33 (MQ=9)
             gATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                       <  1:2206140/36‑1 (MQ=255)
               tGGTTCTTTCGGCGGTATTTGTACCAATGGCcttct                     <  1:6297066/36‑1 (MQ=255)
                ggTTCTTTCGGCGGTATTTGTACCAATGGCcttctt                    <  1:2735566/36‑1 (MQ=255)
                 gTTCTTTCGGCGGTATTTGTAGCAATGGCCttcttc                   <  1:1921926/36‑1 (MQ=255)
                 gTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                   <  1:2509742/36‑1 (MQ=255)
                   tCTTTCGGCGGTATTTTTTCCCATTGCCttcttctg                 >  1:603394/1‑34 (MQ=16)
                   tCTTTCGGCGGTATTTTTACCAATTGCcttcttttt                 >  1:5364942/1‑33 (MQ=18)
                   tCTTTCGGCGGTATTTTTACCAATTGCcttcttttg                 >  1:4754381/1‑33 (MQ=18)
                   tCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcgg                 <  1:3317904/36‑1 (MQ=255)
                         ggcggTATTTGTACCAATGGCCTTCTTTGGCGGcac           >  1:2528428/1‑36 (MQ=255)
                          gcggTATTTGTACCAATGGCCTTCTTCGGCGGcacc          <  1:239754/36‑1 (MQ=255)
                          gcggTATTTGTACCAATGGCCTTCTTCGGCGGcacc          >  1:2722514/1‑36 (MQ=255)
                          gcggTATTTGTACCAATGGCCTTCTTCGGCGGcacc          <  1:1498862/36‑1 (MQ=255)
                          gcggTATTTGTACCAATGGCCTTCTTCGGCGGcacc          >  1:1422292/1‑36 (MQ=255)
                            ggTATTTGTACCAATGGCCTTCTTCGGCGGcaccac        <  1:581249/36‑1 (MQ=255)
                            ggTATTTGTACCAATGGCCTTCTTCGGCGGcaccac        >  1:6637248/1‑36 (MQ=255)
                              tATTTGTACCAATGGCCTTCTTCGGCGGCACCACCg      <  1:1529581/36‑1 (MQ=255)
                              tATTTGTACCAATGGCCTTCTTCGGCGGCACCACCg      <  1:1168740/36‑1 (MQ=255)
                               aTTTGTACCAATGGCCTTCTTCGGCGGCACCACCgg     >  1:698201/1‑36 (MQ=255)
                                tttGTACCAATGGCCTTCTTCGGCGGCACCACCGGt    <  1:6608929/36‑1 (MQ=255)
                                tttGTACCAATGCCCTTCTTCGGCGGCACCACCGGt    >  1:1382233/1‑36 (MQ=255)
                                  tGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGc  <  1:2803140/36‑1 (MQ=255)
                                   |                                  
TGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGC  >  REL606/2509747‑2509816

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG