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,0510TG91.8% 144.1 / 3.4 49F257C (TTT→TGT) emrDmultidrug efflux MFS transporter EmrD
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (32/13);  total (35/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCT  >  NZ_CP009273/3847972‑3848137
                                                                               |                                                                                      
ggCGTGCTGATGGGCGCGGTGTTAGGGCTGGGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGGGCa                                                                              >  1:1200399/1‑90 (MQ=255)
  cGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATg                                                                            <  2:1468793/90‑1 (MQ=255)
   gTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATgg                                                                           >  1:453536/1‑90 (MQ=255)
   gTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGgcgc                                                                               >  2:835858/1‑86 (MQ=255)
   gTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGgcgc                                                                               <  1:835858/86‑1 (MQ=255)
           gggCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCgg                                                                   <  1:938743/90‑1 (MQ=255)
               gcgGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGt                                                               >  1:725001/1‑90 (MQ=255)
                cgGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTc                                                              <  2:941886/90‑1 (MQ=255)
                cgGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTc                                                              >  1:320271/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:219298/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1244722/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1547419/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  <  2:1492318/90‑1 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1280263/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:520997/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1459536/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:669174/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:1578122/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:1492318/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  1:157882/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:871413/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:219408/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:192120/1‑90 (MQ=255)
                            tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                  >  2:930680/1‑90 (MQ=255)
                                  gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCa                                            >  2:855143/1‑90 (MQ=255)
                                  gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCa                                            >  2:658445/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:640710/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:181338/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:1228276/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:1243128/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  2:691081/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:390991/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:71141/1‑90 (MQ=255)
                                            cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                  >  1:704903/1‑90 (MQ=255)
                                                  tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  2:1265444/90‑1 (MQ=255)
                                                  tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                            <  1:520997/90‑1 (MQ=255)
                                                   ttGGTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCg                           >  1:808234/1‑90 (MQ=255)
                                                      tttATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGtt                        <  2:157882/90‑1 (MQ=255)
                                                        tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTc                                                              <  1:977064/50‑1 (MQ=255)
                                                        tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTc                                                              >  2:977064/1‑50 (MQ=255)
                                                        tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                      >  2:781837/1‑90 (MQ=255)
                                                           tCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATctg                   <  2:640710/90‑1 (MQ=255)
                                                              gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATctgctg                <  1:930680/90‑1 (MQ=255)
                                                               ccGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCtgctgc               >  2:1221229/1‑90 (MQ=255)
                                                                          aGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGctggctgg    <  2:390991/90‑1 (MQ=255)
                                                                           gCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCtggctggc   >  2:1508739/1‑90 (MQ=255)
                                                                            cGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                              <  2:1545103/46‑1 (MQ=255)
                                                                            cGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                              >  1:1545103/1‑46 (MQ=255)
                                                                            cGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTggctggct  <  1:938048/90‑1 (MQ=255)
                                                                               |                                                                                      
GGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCT  >  NZ_CP009273/3847972‑3848137

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