Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 465,858 T→C I511T (ATC→ACC)  mdlA → SmdA family multidrug ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273465,8580TC75.7% 54.5 / 21.1 37I511T (ATC→ACC) mdlASmdA family multidrug ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base C (12/16);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGC  >  NZ_CP009273/465725‑465995
                                                                                                                                     |                                                                                                                                         
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGTCGCACTGAGCACCAGATCCTg                                                                                                                                       >  1:223184/1‑138 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTg                                                                                                                                       >  1:178817/1‑138 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTg                                                                                                                                       >  2:816345/1‑138 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTg                                                                                                                                       >  2:178683/1‑138 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTg                                                                                                                                       >  1:815039/1‑138 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGc                                                                                                                                      >  2:218523/1‑139 (MQ=255)
ggCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGCGGACGGACCCCCTGAGCACCAGGCCCTg                                                                                                                                       >  2:499817/1‑138 (MQ=255)
               aTGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGt                                                                                                                       <  1:645440/139‑1 (MQ=255)
               aTGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       <  1:218523/139‑1 (MQ=255)
               aTGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTAGTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       <  1:499817/139‑1 (MQ=255)
                     tCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTggggg                                                                                                                  >  1:866355/1‑138 (MQ=255)
                     tCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTggggg                                                                                                                  <  1:816345/138‑1 (MQ=255)
                       cggcggGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCa                                                                                                                <  1:865372/138‑1 (MQ=255)
                       cggcggGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGGCAGTGGGGGCAg                                                                                                               >  2:240265/1‑139 (MQ=255)
                                aaaaCAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACg                                                                                                      <  2:384134/139‑1 (MQ=255)
                                        gTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAatcat                                                                                              <  1:95763/139‑1 (MQ=255)
                                                       gTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAg                                                                                           >  2:109679/1‑127 (MQ=255)
                                                       gTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAg                                                                                           <  1:109679/127‑1 (MQ=255)
                                                         gCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCttt                                                                             >  1:129603/1‑139 (MQ=255)
                                                               ttaGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCAc                                                                       <  2:129603/139‑1 (MQ=255)
                                                               ttaGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCAc                                                                       <  2:178817/139‑1 (MQ=255)
                                                                       cgcgGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACgg                                                                 >  1:802076/1‑137 (MQ=255)
                                                                       cgcgGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       >  2:348646/1‑83 (MQ=255)
                                                                       cgcgGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       <  1:348646/83‑1 (MQ=255)
                                                                                    aTCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAatta                                                  >  2:731796/1‑139 (MQ=255)
                                                                                           atgatgCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGAt                                           <  2:802076/139‑1 (MQ=255)
                                                                                            tgatgCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATg                                          <  1:623076/139‑1 (MQ=255)
                                                                                            tgatgCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATg                                          <  2:8123/139‑1 (MQ=255)
                                                                                              atgCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATg                                          >  1:635059/1‑137 (MQ=255)
                                                                                               tgCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATgcag                                       <  1:670079/139‑1 (MQ=255)
                                                                                                gcgcTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATtg                                               >  1:762961/1‑130 (MQ=255)
                                                                                                gcgcTTTCGGCGGTGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATtg                                               <  2:762961/130‑1 (MQ=255)
                                                                                                             tGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       <  2:95370/45‑1 (MQ=255)
                                                                                                             tGGACGGACGCACTGAGCACCAGACCCTGCATAACCTGCGTCAGt                                                                                                                       >  1:95370/1‑45 (MQ=255)
                                                                                                                         cTGAGCACCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCaa             <  1:630219/139‑1 (MQ=255)
                                                                                                                               aCCAGACCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCatga       <  2:815039/139‑1 (MQ=255)
                                                                                                                                    aTCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGc  <  1:731796/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
GGCGAGCGCGGTGTGATGCTTTCCGGCGGGCAAAAACAGCGTATCTCCATTGCTCGTGCGTTATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGC  >  NZ_CP009273/465725‑465995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑