Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,302,6580TG44.0% 60.5 / 12.5 50T141P (ACC→CCC) mdtOmembrane translocase (MDR) of MdtNOP efflux pump, PET family
Reads supporting (aligned to +/- strand):  ref base T (11/17);  new base G (22/0);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 2.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAAAAA  >  NC_000913/4302513‑4302805
                                                                                                                                                 |                                                                                                                                                    
gcCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTgtc                                                                                                                                                 >  2:794712/1‑149 (MQ=255)
     aTGCATTTGCCAAATCGCACGACTGGGAATCCACAGCACGCCGATTAACATCATCTGCAAGGTTGGATAGAGGCCAACAACTATACACCACAGCGTTAAGCGCTCTACCACTTCCGTATAGTCCAGCATTGCGGGGAAGGTTTGCCCGTa                                                                                                                                             <  2:505373/150‑1 (MQ=255)
       gCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTaaa                                                                                                                                           <  2:572254/150‑1 (MQ=255)
         aTTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTccccccaata                                                                                                                                         >  2:485031/1‑143 (MQ=255)
           ttGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATa                                                                                                                                         >  1:771873/1‑148 (MQ=255)
             gCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATaac                                                                                                                                       >  2:331040/1‑146 (MQ=255)
              cGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAATAGC‑g                                                                                                                                     >  1:737422/1‑148 (MQ=255)
                aaaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCGAAATAAC‑Gaaaa                                                                                                                                 >  2:335894/1‑150 (MQ=255)
                 aaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATg                                                                                                                               <  1:701567/151‑1 (MQ=255)
                    gtcACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTTCCCGTAAATAGC‑GACAATGGCg                                                                                                                            >  1:251127/3‑151 (MQ=255)
                    ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCg                                                                                                                            <  2:484204/151‑1 (MQ=255)
                       aCGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTTGCCCTAAAAAGG‑GACAAAGgcggcg                                                                                                                         >  1:806093/1‑151 (MQ=255)
                         gACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAAAAGC‑GACAATGGGGGCgg                                                                                                                        >  2:703447/1‑150 (MQ=255)
                          aCTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCg                                                                                                                      >  2:795589/1‑151 (MQ=255)
                          aCTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAATAGC‑GACAAAGGCGACGGCg                                                                                                                      >  1:710332/1‑151 (MQ=255)
                          aCTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGAGCGGGTTTCCCGTTAATAGT‑GGCAAttggg                                                                                                                            >  1:818767/1‑141 (MQ=255)
                           cTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGGAAATAGA‑GACAATGGCGCCGGCg                                                                                                                      >  2:772149/1‑150 (MQ=255)
                              ggAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGaaaa                                                                                                                  <  2:655927/151‑1 (MQ=255)
                                           cGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAATAGG‑GAAAATGGGGAAGGCGAAAAAGAACCGcccc                                                                                                       >  1:603594/1‑149 (MQ=255)
                                             ccGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAAGGTTGCCCGTTAATAAC‑GACAATGGCGACGGCGAAAAAGAACAAGCCCa                                                                                                      >  1:624552/1‑148 (MQ=255)
                                                aTTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGAACAGcccc                                                                                                       >  1:216908/1‑144 (MQ=255)
                                                   aaCGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGGCCGGAAAAAAC‑GACAAAGGCGGCGGCGGAAAAAAACAACCCCCAACGATggg                                                                                             >  2:127254/1‑151 (MQ=255)
                                                         atcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGGCCGTAAATAGC‑GACAAAGGCGACGGCGAAAAAAAACAACCCCAAACGGGGGGTgcgcc                                                                                       >  2:828229/1‑150 (MQ=255)
                                                         atcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTTCCCGTAAAAAGC‑GACAATGGCGACGGCGAAAAAGAACAGCCCCAAACGATGGGGgcg                                                                                         >  2:312676/1‑149 (MQ=255)
                                                         atcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTCAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGCAAATAGG‑GACAATGGTGACCGGGAAAAAGACCAGCCCCCAGCGATGGGGgcgc                                                                                        >  2:42381/1‑150 (MQ=255)
                                                                            ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCAt                                                                    <  1:856773/151‑1 (MQ=255)
                                                                              ccAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTGCCCGTAAATAAC‑GACAAAGGCGGCGGCGAAAAAAACCAGCCCCAAACGAAGGGGGCGCATCAAAAAAAAGCAGCCCATc                                                                   >  1:292073/1‑150 (MQ=255)
                                                                                aacaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGCAAAAAGC‑GACAATGGCGACGGCGGAAAAAACCCGCCCCCAGCGATGGGGGCGCATCAAAAAAAAGCAGCCCAACa                                                                  >  2:89449/1‑149 (MQ=255)
                                                                                  caacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGAt                                                              <  1:224657/151‑1 (MQ=255)
                                                                                   aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAATAAC‑GAAAAAAGCGACGGCGGAAAAAACCCGCCCCCAACGAGGGGTGCGCATAAAAAAAAAATAGCCCATCAgg                                                                >  2:80955/1‑148 (MQ=255)
                                                                                     cGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTGCCCGTAAATAGC‑GACCAAGGGGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCgg                                                           >  1:412213/1‑151 (MQ=255)
                                                                                          caccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACAAGCAGCCCATCAGGATCGGTCCgg                                                      >  2:247879/1‑151 (MQ=255)
                                                                                            ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAATAGC‑GACAAAGGCGACGGCGGAAAAAACCAACCCCCAACGATGGGGGCGCATCAAAAAAAAGCAACCCCTccagaaagggccgg                                                      >  2:879616/1‑136 (MQ=255)
                                                                                             cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCga                                                   <  2:635434/151‑1 (MQ=255)
                                                                                                aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgatga                                                <  2:402238/151‑1 (MQ=255)
                                                                                                 gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgatgat                                               <  2:240572/151‑1 (MQ=255)
                                                                                                   gTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTCAATAGC‑GACAAAGGCGGCGGCGAAAAAGACCAGCCCCAAGCGGAGGGTGCGCCCCAAAAAAATGCCGCCCATCAGGAGCGGTcccggggtg                                                 >  2:142423/1‑140 (MQ=255)
                                                                                                          gcACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGAt                                      <  2:744923/151‑1 (MQ=255)
                                                                                                            aCGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGTCCCTCAATAACGGCCAATGG‑GGCGGCGAAAAAGACCAGCCCCAAGCCAGGGGGGCGCATCAAAAAAATGCCGCCCCTCAGGATCGGTCCGGCGGGCATCAATCgg                                        >  2:455325/1‑147 (MQ=255)
                                                                                                              gACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTCAATAAC‑GACAATGGCGGCGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATCAAAAACATGCAGCCCAACAGGAGCCGTCCGGCGATGATCAAACgggacaa                                   >  2:140879/1‑145 (MQ=255)
                                                                                                              gACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGC‑GACAATGGCGGCGGCGAAAAAGACCAGCCCCAAGCGATGGGTGGGCATCAAAAAAATGCAGCCCATCAGGATCGGTCCGGCGGTGATCAATCGGATaaa                                   >  1:527584/1‑150 (MQ=255)
                                                                                                                      ccGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAAACCAACCCCCAACGATGGGGGGGGATCAAAAACATGCAGCCCAACAGGAACGGGCCGGCGATGATCAAACGGATCAACGGGTTGcc                          >  2:457868/1‑151 (MQ=255)
                                                                                                                                gAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCAt                <  2:491640/151‑1 (MQ=255)
                                                                                                                                gAGCATGGCGGGGAAGGGTTGCCCCGAAAAAGC‑GACAAAGGCGACGGCGGAAAAGACCAACCCCAAGCGATGGGTGCGCATCAAAAAAAAGCAGCCCATCAGGATCGGGCCGGCGGGGGTCAATCGGATCCACGGTTCGCCGCATGACCa                 >  2:298978/1‑150 (MQ=255)
                                                                                                                                  gCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCAttt              <  2:914215/151‑1 (MQ=255)
                                                                                                                                   cATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTg             <  1:330917/151‑1 (MQ=255)
                                                                                                                                   cATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTg             <  1:60332/151‑1 (MQ=255)
                                                                                                                                        cGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGAt        <  1:878392/151‑1 (MQ=255)
                                                                                                                                             aaGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCaaa    >  1:149038/1‑150 (MQ=255)
                                                                                                                                              aGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCaaaaa  <  2:444036/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGC‑GACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAAAAA  >  NC_000913/4302513‑4302805

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