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

CCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAC  >  NC_000913/4568372‑4568643
                                                                                                                                     |                                                                                                                                           
ccAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaaaaa                                                                                                                                      <  2:261741/140‑1 (MQ=255)
         cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAATAAACGTGGCATa                                                                                                                             <  2:713394/140‑1 (MQ=255)
         cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATa                                                                                                                             >  1:347348/1‑140 (MQ=255)
             ggCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCGGCCTaaaaa                                                                                                                         >  2:220553/1‑140 (MQ=255)
                    aaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTc                                                                                                                  <  2:407420/140‑1 (MQ=255)
                    aaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTc                                                                                                                  <  1:664740/140‑1 (MQ=255)
                     aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTcc                                                                                                                 <  1:47285/140‑1 (MQ=255)
                        cACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTg                                                                                                              <  1:311054/140‑1 (MQ=255)
                              ccACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGga                                                                                                        >  1:534351/1‑140 (MQ=255)
                                       tGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAACAACTCCTTGTTCGGAGAAccccca                                                                                               >  2:509206/1‑138 (MQ=255)
                                         aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTTTCCCGAGAACCCCAACg                                                                                             >  1:275053/1‑140 (MQ=255)
                                         aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCCTAAAAAACAAAACCCTGTCCGGAGAAACCCCAcc                                                                                             >  2:330480/1‑139 (MQ=255)
                                         aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGGGGCGGCATGGCGGGGAAAAAAAAAGGGGGTAAAAAAAAAACCCCTTTCCCGGGGAACCCCaag                                                                                              >  1:576811/1‑139 (MQ=255)
                                           cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg                                                                                           >  1:96289/1‑140 (MQ=255)
                                           cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAAAAAAAAAACCCTTGTCCGGGGAACCCCCACggg                                                                                           >  2:341799/1‑140 (MQ=255)
                                            ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTTTCCGGAGAACCCCAACgggg                                                                                          >  1:391984/1‑140 (MQ=255)
                                            ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAACAAACAACTCCTTTTCCGGAGAACCCCCACgggg                                                                                          >  2:365170/1‑140 (MQ=255)
                                              aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCCTAAAAAAAAAAACCTTTTCCGGAGGACCCCCACcgggga                                                                                        >  2:738458/1‑138 (MQ=255)
                                              aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGGAAAAAAGAAGGCATAAAAAAAAACTCCTTGTCCCGAGAACCCCAACGGGGaa                                                                                        >  1:59259/1‑140 (MQ=255)
                                              aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACAGGGATTACAAACAAATCCTTCTCCGGAGAAACCCAACggggga                                                                                        >  2:175652/1‑138 (MQ=255)
                                                gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTTTCCGGAGAACCCCAACGGGGgaaa                                                                                      >  1:293552/1‑140 (MQ=255)
                                                gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAAACCTTGTCCGGAGAACCCCAACGGGGaaaa                                                                                      >  2:393230/1‑140 (MQ=255)
                                                gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCGGCATAACAAAAAAAACCTTGTCCGGGGAACCCCCACCGGGaaaa                                                                                      >  2:17494/1‑140 (MQ=255)
                                                  tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                    <  1:222539/140‑1 (MQ=255)
                                                  tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                    <  2:534351/140‑1 (MQ=255)
                                                  gagATCAGATCCGTCGACAGAGACGCAGCAAAGTCATACAAAATCAACGCCATCGGCAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACACCTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                    <  2:362055/138‑1 (MQ=255)
                                                    gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAAAAAAAAAAACCTTTTTCGGAGAAACCCCACGGGGGAAAAGaa                                                                                  >  2:766278/1‑140 (MQ=255)
                                                     aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                    >  2:481788/1‑137 (MQ=255)
                                                     aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                    <  1:481788/137‑1 (MQ=255)
                                                     aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCGGCATAAAAAAAAAATCCCTGTCCGGAGAACCCCCACGGGGGAAAAGaaa                                                                                 >  2:37467/1‑140 (MQ=255)
                                                           tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAAt                                                                           <  1:365170/140‑1 (MQ=255)
                                                               tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTTTCCGGAGAACCCCAACGGGGGAAAAGAAAATCAATaacc                                                                       >  2:552951/1‑137 (MQ=255)
                                                                      aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                                <  2:561167/140‑1 (MQ=255)
                                                                      aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCCACGGGGAAAACGAAAATCAAAACCGTCGGTaa                                                                >  2:589769/1‑140 (MQ=255)
                                                                              caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                        <  2:452395/140‑1 (MQ=255)
                                                                                aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAgg                                                      <  2:278952/140‑1 (MQ=255)
                                                                                aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCTTGTCCGGAGAACCCCCACGGGGGAAAAGAAAATCAATACCCTCGGTAATGACGaaaag                                                      >  2:457652/1‑138 (MQ=255)
                                                                                 aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAAATCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAAATCAATAACGTCGGGAATGACGaaaagg                                                     >  2:401213/1‑137 (MQ=255)
                                                                                     aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGt                                                 <  1:552951/140‑1 (MQ=255)
                                                                                         aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                             <  2:754045/140‑1 (MQ=255)
                                                                                         aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                             <  1:37467/140‑1 (MQ=255)
                                                                                             tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGCGGCAAAAAAAAAAACTCCTTTTCCGGGGAAGCCCCACAGGGAAAACGAAAATCAAAACCCTCGGGAATGACGAAAAGCAGGGGAGAAagg                                         >  2:289965/1‑140 (MQ=255)
                                                                                                 cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGAAATTCAATACCGTCGGGAATGACGAAAGGCAGGTGAGAAGGGCAAc                                     >  1:449729/1‑140 (MQ=255)
                                                                                                 cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAAc                                     <  1:738458/140‑1 (MQ=255)
                                                                                                           aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                           <  2:576811/140‑1 (MQ=255)
                                                                                                             aaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAAt                         <  1:766278/140‑1 (MQ=255)
                                                                                                                 aGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                     <  2:96289/140‑1 (MQ=255)
                                                                                                                    gTGGCGGCATGGCGGGTAAGAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGa                  <  2:279708/140‑1 (MQ=255)
                                                                                                                    gTGGCGGCATGGCGGGTAAGAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGa                  >  1:279708/1‑140 (MQ=255)
                                                                                                                      ggcggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAAt                <  1:289965/140‑1 (MQ=255)
                                                                                                                        cggcATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTGTCCGGAGAACCCCAACGGGGAAAAAGAAAATCAATACCGTCGGTAATGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCCAAAGTACTGAAttt              >  2:493359/1‑140 (MQ=255)
                                                                                                                        cggcATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCCTTTCCGGAGAAACCCAACGGGGGAAACGAAAATCAAAACCCTCGGTAAAGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCAAATGTACTGAAttt              >  2:244661/1‑140 (MQ=255)
                                                                                                                            aTGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATg          <  1:509206/140‑1 (MQ=255)
                                                                                                                             tGGCGGGTAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa        <  2:660096/140‑1 (MQ=255)
                                                                                                                              ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa        <  1:31044/140‑1 (MQ=255)
                                                                                                                                cGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAt      <  1:401213/140‑1 (MQ=255)
                                                                                                                                   ggTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCa    <  2:216309/140‑1 (MQ=255)
                                                                                                                                     tAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAc  <  2:116058/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                           
CCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAC  >  NC_000913/4568372‑4568643

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