Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,848,051 T→G F257C (TTT→TGT)  emrD → multidrug efflux MFS transporter EmrD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,848,0510TG96.2% 80.4 / ‑3.4 26F257C (TTT→TGT) emrDmultidrug efflux MFS transporter EmrD
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (21/4);  total (22/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGC  >  NZ_CP009273/3847965‑3848136
                                                                                      |                                                                                     
cTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGtgttt                                                                                    <  1:132040/90‑1 (MQ=255)
  gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                                                  >  2:497914/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:110080/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:400705/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:314090/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:268166/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:223910/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:583542/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:451407/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:383932/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:167500/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:248692/1‑90 (MQ=255)
                                         gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                             >  2:90520/1‑88 (MQ=255)
                                         gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                             <  1:90520/88‑1 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:408224/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  2:29290/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:81026/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:361410/1‑90 (MQ=255)
                                                         tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                           <  1:268166/90‑1 (MQ=255)
                                                         tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                           <  2:451407/90‑1 (MQ=255)
                                                               tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                     >  2:248766/1‑90 (MQ=255)
                                                               tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                     >  1:240560/1‑90 (MQ=255)
                                                                     gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATctgctg               <  1:42593/90‑1 (MQ=255)
                                                                      ccGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCtgctgc              >  2:319041/1‑90 (MQ=255)
                                                                               gCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTgctggct     >  1:116779/1‑90 (MQ=255)
                                                                                 aGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGctggctgg   >  1:310138/1‑90 (MQ=255)
                                                                                  gCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCtggctggc  >  1:101181/1‑90 (MQ=255)
                                                                                      |                                                                                     
CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGC  >  NZ_CP009273/3847965‑3848136

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