Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,568,5040TG33.9% 88.5 / 11.4 57T6P (ACC→CCC) mdtMmultidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (19/0);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 5.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAT  >  NC_000913/4568374‑4568639
                                                                                                                                  |                                                                                                                                        
aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGTTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                                 <  1:927310/140‑1 (MQ=255)
aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                                 >  1:751452/1‑140 (MQ=255)
aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                                 <  2:102238/140‑1 (MQ=255)
            gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaa                                                                                                                     >  1:48356/1‑140 (MQ=255)
            gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATaaaaaa                                                                                                                     >  2:941975/1‑140 (MQ=255)
             cATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATaaaaaaa                                                                                                                    >  2:56494/1‑139 (MQ=255)
              aTTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAATGGCATaaaaaaaa                                                                                                                   >  2:328076/1‑138 (MQ=255)
                   aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTcc                                                                                                             <  2:300307/140‑1 (MQ=255)
                   aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTcc                                                                                                             <  2:12677/140‑1 (MQ=255)
                       aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGt                                                                                                         <  1:550811/140‑1 (MQ=255)
                       aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACATGGCGTAAAAAAAAACTCCTTGt                                                                                                         >  2:300690/1‑140 (MQ=255)
                             cacaTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGgag                                                                                                   <  2:343116/140‑1 (MQ=255)
                               caTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGGGaa                                                                                                 >  2:174580/1‑140 (MQ=255)
                                aTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCCTTTCCGGAGaaa                                                                                                >  1:501637/1‑139 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAAAAACAACTCCCTGTCCGGGGAACCCCAACg                                                                                         >  1:409719/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAAAAAACAACTCCTTGTCCGGAGAAACCCAACg                                                                                         >  1:218544/1‑140 (MQ=255)
                                       aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGATAAAAAAAAACTCCTTGTTCGGAGAACCCCCACg                                                                                         >  2:22638/1‑140 (MQ=255)
                                         cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGCGGCATAAAAAACAACACCTTCTCCGGAGGACCCCCACggg                                                                                       >  2:220168/1‑140 (MQ=255)
                                         cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGCGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg                                                                                       >  1:122584/1‑140 (MQ=255)
                                         cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGCAAAAAAACATGGCATAACAAACAACTCCTTTTTCGGAGAAGCCCCACggg                                                                                       >  1:854763/1‑140 (MQ=255)
                                          ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTCAAAAAAAAGGGAATAAAAAAACACCCCTTGTTCATGGAAACCCCACgggg                                                                                      >  2:445311/1‑140 (MQ=255)
                                          ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACAAGGCCTAAAAAAAAAATCCCTTTGCGGGGGAACCCCAAgggg                                                                                       >  2:943570/1‑140 (MQ=255)
                                             ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCAGAAAAAAAAAATCCTTTTCCCGAGAAACCCCACGGGGgaa                                                                                   >  2:889383/1‑140 (MQ=255)
                                                tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                <  1:22638/140‑1 (MQ=255)
                                                tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                <  1:445311/140‑1 (MQ=255)
                                                tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                <  2:48356/140‑1 (MQ=255)
                                                  gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTATCCGGAGAACCCCAACGGGGAAAAAGaa                                                                              >  1:649815/1‑140 (MQ=255)
                                                   aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGTGCATAAAAAACAACTCCTTGTCCGGAGAAACCCCACGGGGGAAAAGaaa                                                                             >  2:27017/1‑140 (MQ=255)
                                                      aGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAAAAAATCCCTGTCCGGAGAAACCCAACGGGGGAAACGAAAATc                                                                          >  2:48206/1‑140 (MQ=255)
                                                       gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCa                                                                         <  1:416491/140‑1 (MQ=255)
                                                            gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATAcc                                                                    <  1:114764/140‑1 (MQ=255)
                                                             tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTTTCCGGGGAACCCCCACGGGGGAAAAGAAAATCaaaaacc                                                                   >  2:770129/1‑135 (MQ=255)
                                                             tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAACAAAAAACTCCTTCTCCGGAGAAACCCAACGGGGGAAACAAAAATCAATAccc                                                                   >  2:114764/1‑139 (MQ=255)
                                                             tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAATTCAATAccc                                                                   >  2:731236/1‑139 (MQ=255)
                                                                         cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAAATCCCTTTTCGGGGAAACCCCACAGGGGAAAAGAAAATCAAAACCCTCGGGAAAGACg                                                       >  2:799379/1‑140 (MQ=255)
                                                                            caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAAATCAATACCGTCGGGAATGAAGaaa                                                    >  1:90459/1‑140 (MQ=255)
                                                                              aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAAAAAAAAACACCCTTTCCGGAGAACCCCCACGGGGGAAAAGAAAATCAATACCGTCGGGAATGACGaaaaa                                                  >  2:196894/1‑138 (MQ=255)
                                                                               aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc                                                 <  1:300690/140‑1 (MQ=255)
                                                                                   aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGt                                             <  1:32508/140‑1 (MQ=255)
                                                                                   aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGt                                             <  2:751452/140‑1 (MQ=255)
                                                                                    tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCAGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTg                                            >  2:955401/1‑140 (MQ=255)
                                                                                      cAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgag                                          <  1:220168/140‑1 (MQ=255)
                                                                                       aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                         <  2:804270/140‑1 (MQ=255)
                                                                                       aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                         <  2:187194/140‑1 (MQ=255)
                                                                                       aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                         <  1:196894/140‑1 (MQ=255)
                                                                                           tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAggg                                     >  1:51577/1‑140 (MQ=255)
                                                                                           tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGGGAACCCCAAAGGGGGAAAAGAAAATCAATAACGTCGGGAAAGACGAAAAGCAGGGGAGAAggg                                     >  2:32508/1‑140 (MQ=255)
                                                                                              gagacagcGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                  <  1:65916/133‑1 (MQ=255)
                                                                                               cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAAc                                 <  1:889383/140‑1 (MQ=255)
                                                                                                 ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTTGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                               <  1:27017/140‑1 (MQ=255)
                                                                                                 ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                               <  1:943570/140‑1 (MQ=255)
                                                                                                  cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag                              <  1:792133/140‑1 (MQ=255)
                                                                                                     cGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGGATAAAAAAAAAACCCTTGTTCCGGGAACCCCCACGGGGGAAAAGAAAATCAATAACCTCCGGAATGACGAAAAGCAGGGGGGAAGGGCAACGAGAcc                           >  2:597380/1‑139 (MQ=255)
                                                                                                           aaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAAt                     <  2:187067/140‑1 (MQ=255)
                                                                                                                gCGTGGCGGCATGGCGGGGAAAAAAAAGCGGCCTAAAAAAAAAATCCCTCTCCGGAGAAACCCAACCGGGGAAAAGAAAATCAATACCCTCGGGAAAGAAGAAAAGGAGGTGAGAAAGGCAACGAGAAAGCAACTGtctc                >  2:65916/1‑137 (MQ=255)
                                                                                                                 cGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAACTCCTTGTCCGGGGAACCCCCACGGGGAAAAAGAAAATCAAAAACGTCGGGAATGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCAAATGGACTg               >  2:31335/1‑140 (MQ=255)
                                                                                                                    ggcggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAAt            <  2:501637/140‑1 (MQ=255)
                                                                                                                            ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa    <  2:879087/140‑1 (MQ=255)
                                                                                                                              cGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAt  <  1:31335/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
AGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAT  >  NC_000913/4568374‑4568639

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