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

GGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTA  >  NC_000913/4568369‑4568637
                                                                                                                                       |                                                                                                                                      
gggCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaaa                                                                                                                                    <  1:394085/140‑1 (MQ=255)
gggCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa                                                                                                                                    >  2:385946/1‑140 (MQ=255)
     aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaag                                                                                                                               >  1:313812/1‑138 (MQ=255)
       aCTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTg                                                                                                                             <  2:427645/140‑1 (MQ=255)
          gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCa                                                                                                                          <  2:645231/140‑1 (MQ=255)
          gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCa                                                                                                                          <  1:716448/140‑1 (MQ=255)
                ggCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAACGGCATaaaaa                                                                                                                    >  2:792132/1‑140 (MQ=255)
                ggCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGCAAAAAAACCTGGCATaacaa                                                                                                                    >  1:148216/1‑140 (MQ=255)
                 gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                               <  2:42646/128‑1 (MQ=255)
                 gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg                                                                                                                               >  1:42646/1‑128 (MQ=255)
                 gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGaaaaaaaa                                                                                                                   >  2:588511/1‑140 (MQ=255)
                        aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTcc                                                                                                            <  1:22450/140‑1 (MQ=255)
                            aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACCCCtttt                                                                                                        >  1:820461/1‑138 (MQ=255)
                               tACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACAGGGGCTAACAAACAAATCCTTGT‑ccc                                                                                                    >  1:801691/1‑139 (MQ=255)
                                     aTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAAc                                                                                              <  1:792132/140‑1 (MQ=255)
                                            aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCTTGT‑CCCGAGAACCCCAACg                                                                                       >  1:174646/1‑140 (MQ=255)
                                              cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCCACggg                                                                                     >  1:416395/1‑140 (MQ=255)
                                               ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACCAGGCATAACAAACAAATCCCTGT‑CCAGCGGACCCCACCtggg                                                                                    >  2:60858/1‑140 (MQ=255)
                                               ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAAATCCTTCT‑TCGGGGAACCCCAACgggg                                                                                    >  1:359082/1‑140 (MQ=255)
                                               ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAAGCACAAAAAAAAAACCCTTTT‑CCAGGGGACCCCAACgggg                                                                                    >  2:423671/1‑140 (MQ=255)
                                               ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAACAAAAAACTCCTTGT‑CCGGGGAACCCCaaacggg                                                                                    >  1:636141/1‑135 (MQ=255)
                                               ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCACAAAAAAAAAACCCCTGTCCCGGGGAACCCCCGCggg                                                                                     >  2:603452/1‑140 (MQ=255)
                                                     tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGTAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACg                                                                              <  2:569627/140‑1 (MQ=255)
                                                     tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACg                                                                              <  2:148986/140‑1 (MQ=255)
                                                     tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACg                                                                              <  1:116762/140‑1 (MQ=255)
                                                      ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGa                                                                             <  2:347428/140‑1 (MQ=255)
                                                        aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGaaa                                                                           >  1:1168/1‑140 (MQ=255)
                                                        aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGaaa                                                                           >  1:503924/1‑140 (MQ=255)
                                                        aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAACAAAAAACTCCTTGT‑CCGGAGAACCCCCACGGGGAAAAAGaaa                                                                           >  2:791086/1‑140 (MQ=255)
                                                            gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCa                                                                       >  1:344581/1‑140 (MQ=255)
                                                            gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCa                                                                       <  2:344581/140‑1 (MQ=255)
                                                             acaggTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCaa                                                                      <  1:467300/136‑1 (MQ=255)
                                                             aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCaa                                                                      >  1:657740/1‑140 (MQ=255)
                                                             aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGCATAAAAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCaa                                                                      >  1:711516/1‑140 (MQ=255)
                                                             aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGAAAAAAAAAAGGGCATAAAAAAAAAATCCTTTT‑CCGGAGAACCCCAACGGGGAAAACGaaaataaa                                                                      >  2:695831/1‑135 (MQ=255)
                                                                 gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATAcc                                                                  <  2:657740/140‑1 (MQ=255)
                                                                   cGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAAAAACTCCTTGT‑CCGGAGAACCCCAACGGGGGAAAAGAAATTCAATAcccc                                                                >  2:280502/1‑138 (MQ=255)
                                                                        gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                           <  1:634370/140‑1 (MQ=255)
                                                                         aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGGAAAAGAAAATCAAAACCGTCGGGaa                                                          >  1:761209/1‑140 (MQ=255)
                                                                            cGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAACAAAAAACTCCTTTT‑CCGGAGAACCCCAAAGGGGGAAAAGAAAATCAAAACCGTCgggaccgt                                                       >  2:467300/1‑134 (MQ=255)
                                                                                 caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                  <  1:695831/140‑1 (MQ=255)
                                                                                 caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                  <  1:423671/140‑1 (MQ=255)
                                                                                   gaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAgg                                                <  1:723563/139‑1 (MQ=255)
                                                                                   aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAAATCCTTGT‑CCGGGGAACCCCAAAGGGGGAAAAGAAAATCAATAACCTCGGGAATGACGaaaag                                                >  2:198646/1‑138 (MQ=255)
                                                                                    aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACTCCTTGT‑CCGGGGAAACCCCACGGGGGAAACGAAAATCAATACCGTCGGTAATGACGGAAggg                                               >  2:634370/1‑139 (MQ=255)
                                                                                      tCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAg                                             <  1:198646/140‑1 (MQ=255)
                                                                                         tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAAAGGGGGAAAAGAAAATTAATACCGGCGGGAATGACGAAAAGCAGGTg                                          >  2:783029/1‑140 (MQ=255)
                                                                                           cAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgag                                        <  2:1168/140‑1 (MQ=255)
                                                                                            aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                       <  2:333974/140‑1 (MQ=255)
                                                                                            aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                       <  1:588511/140‑1 (MQ=255)
                                                                                               aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGCTAAAAAAAAAATCCCTGT‑CCGGGGAAACCCCACGGGGGAAAAGAAAATTAATAACCTCGGGAATGACGAAAAGGAGGGGAGAAgg                                    >  2:723563/1‑140 (MQ=255)
                                                                                                tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACCCCCTGT‑CCGGAGAACCCCAAAGGGGGAAAAGAAAATTAATAACCGCGGGAATGACGAAAAGCAGGGGGGAAagg                                   >  2:497876/1‑140 (MQ=255)
                                                                                                   aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                <  1:603452/140‑1 (MQ=255)
                                                                                                      ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                             <  1:497876/140‑1 (MQ=255)
                                                                                                         tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                          <  2:636141/140‑1 (MQ=255)
                                                                                                                    aGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc               <  1:385946/140‑1 (MQ=255)
                                                                                                                      cGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAACAACTCCTTGT‑CCGTAGAACCCCCACGGGGGAAACGAAAATCAATACCGTCGGGAATGACGAAAGGGAGGGGGGAAGGGGAACGAGACTGCAAATGTACTg             >  2:261761/1‑140 (MQ=255)
                                                                                                                            ggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGTAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTc       <  1:261761/140‑1 (MQ=255)
                                                                                                                                 ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa  <  1:791086/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGT‑CCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTA  >  NC_000913/4568369‑4568637

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: ‑