Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,301,182	0	G	T	42.0%	56.1 / 25.0	50	Q633K (CAA→AAA)	mdtO	membrane translocase (MDR) of MdtNOP efflux pump, PET family
Reads supporting (aligned to +/- strand):  ref base G (7/22);  new base T (21/0);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 2.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAGAAGATCGCG  >  NC_000913/4301036‑4301326                                                                                                                                                   |                                                                                                                                                 ggCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTgagg                                                                                                                                               <  2:905541/150‑1 (MQ=255)  gCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGGTCTTTACACGTTTTTTtagggg                                                                                                                                             >  2:298624/1‑149 (MQ=255)                 tCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCg                                                                                                                               <  1:438031/150‑1 (MQ=255)                         gggcgTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTTTCTTCAAGCGTTTTTTTAGGAGAACGGCTGTGGGCGGCTGcccg                                                                                                                       >  1:124664/1‑148 (MQ=255)                            cgTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGGTATTTGAGGAGAACGGCTGAGTGCGGGTGCCa                                                                                                                        >  1:144130/1‑146 (MQ=255)                             gTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCg                                                                                                                   <  1:247918/150‑1 (MQ=255)                             gTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCg                                                                                                                   <  2:876011/150‑1 (MQ=255)                                 aTGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGc                                                                                                              <  1:778320/151‑1 (MQ=255)                                 aTGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGc                                                                                                              <  1:601482/151‑1 (MQ=255)                                 aTGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTCAAGCGGTTTTTTAGGAGAAAGGCTGGGTGGGGGTGGCCGACCgg                                                                                                                  >  2:892440/1‑147 (MQ=255)                                  tGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTATTTTAGGGGAACGGCCGAGTGCGGGTGCCAGACCGGCAGCg                                                                                                             >  1:52490/1‑151 (MQ=255)                                      ggggCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTAtt                                                                                                         <  1:433425/151‑1 (MQ=255)                                         gCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTtct                                                                                                      <  1:560937/151‑1 (MQ=255)                                            gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTGTCTTCCAACGTTATTTTTGGAGAAAGGCTGGGGGCGGTTTCCAGACCGGGAGCGGATTTTTcc                                                                                                    >  1:788152/1‑150 (MQ=255)                                            gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTtctc                                                                                                     <  1:283357/149‑1 (MQ=255)                                              ccAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAgg                                                                                                 <  2:208058/151‑1 (MQ=255)                                              ccAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAgg                                                                                                 <  1:678587/151‑1 (MQ=255)                                                    tGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTTAAGCGTTATTTTAGGAGAACGGCTGAGTGCGGTTGCCAGACCCGCAGCGGATTTTTCCCggggg                                                                                              >  1:779868/1‑146 (MQ=255)                                                    tGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGGCTCTTTAAACGGTATTTGAGGAGAACGGCTGAGGGCGGGTGTCCGACCGGCAGCGGATTTTTCCCGGGCGTCg                                                                                           >  1:111049/1‑151 (MQ=255)                                                            gggCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACt                                                                                    <  2:558945/150‑1 (MQ=255)                                                                gAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTTAGGAGAAAGGCTGAGGGCGGTTGCCAGACCGGCAGCGGATTTTTCCAGGGcggcggggaacgg                                                                                   >  1:341412/1‑140 (MQ=255)                                                                gAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCgg                                                                                          <  2:341404/140‑1 (MQ=255)                                                                gAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGGTATTTTAGGAGAACGGCTGAGTGGGGGTGCCAGACCGGCAGCGTATTTCTCCCGGGGGTCGGCAAACTGAcc                                                                                >  1:682168/1‑150 (MQ=255)                                                                 aaaaGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGGCTCTTCAAGCGTTATTTTAGGGGAAAGGCTGAGTGTGGGTTTCCGACCGGACGCGGATTTTTCCAGGGGGTCGGCAAAATGAc                                                                                 >  1:657440/1‑148 (MQ=255)                                                                  aaaGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGGTTTTTTAGGAGAAAGGCTGGGGGCGGGTGCCAGACCGGCAGCGTATTTCTCCCGGGGGGCGGCCAAATGAcc                                                                                >  1:524708/1‑148 (MQ=255)                                                                        ggCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTg                                                                         <  1:837530/149‑1 (MQ=255)                                                                        ggCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTTTTTAAGCGTTATTTTAGGGGAACGGCTGAGTGCGGTTGCCAGACCGGCCGCGGATTTTTCCAGGGCGTCGGCAAAATGACCCGCTCTGt                                                                        >  2:804642/1‑150 (MQ=255)                                                                          cAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTTAAGCGGTTTTTTAGGGGAAAGGGTGAGGGCGGGTGCCCGACCGGCAGCGGATTTTTCCAGGGGGGCGGCAAAATGACCTGCTCTGTcc                                                                      >  1:312801/1‑150 (MQ=255)                                                                            gACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAACGGTATTTTAGGAGGACGGCTGAGGGCGGTTGCCCGACCGGCAGCGGATTTTTCCCGGGCGTCGGCCAACTGACCCGCTCTGTcc                                                                      >  2:591029/1‑148 (MQ=255)                                                                              cGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTg                                                                 <  1:551164/151‑1 (MQ=255)                                                                              cGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTATTTTAGGAGAAAGGGTGGGTGGGGGTGGCAGACCGGCAGCGGATTTTTTCCGGGGGTCGGCAAAATGACCCGCCCTGTCCGGGTg                                                                 >  2:787656/1‑151 (MQ=255)                                                                                  ttGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTTAAGCGTTTTTTTAGGGGAACGGCTGAGGGGGGGTGCCCGACCGGGAGCGTATTTTTCCCGGGCGTCGGCAAAATGACCCGCTCTGTCCgg                                                                    >  1:785601/1‑144 (MQ=255)                                                                                   tGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGGTATTTTAGGGGGACGGCTGAGTGGGGGTGCCAAACCGGCAGCGGATTTTTTCCGGGGGGCGGCCAAACGACCCGCCCCGTCCCGCCGTAggg                                                            >  2:65059/1‑149 (MQ=255)                                                                                     cTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGGCTCTTCAAGCGTTTTTTTAGGAGAAAGGCTGGGTGCGGGTGGCAGACCGGCAGCGTATTTTTCCAGGGCGGCGGGCAACTGACCTGCTCTGTCCGGCTGTAGTGCg                                                          >  1:81893/1‑151 (MQ=255)                                                                                     cTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGCGTTTCTTTAAACGTTTTTTTGGGGGAACGGCTGGGGGGGGGTGCCAGACCGGCAGCGGATTTTTTCAGGGGGGCGGCCAACCGACCTGCCCTTTCCGGCTGGAGTGCg                                                          >  2:488615/1‑151 (MQ=255)                                                                                               aataaGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAgcg                                                <  1:106189/151‑1 (MQ=255)                                                                                               aataaGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCCTCAAACGGTATTTGAGGAGAACGGCTGGGTGGGGGTGGCCGACCCGCAGCGGATTTCTCCAGGGCGGCGGGCAACCGACCTGCCCTGTCCGGCTGTAGTGCGGTTTcc                                                    >  2:718190/1‑147 (MQ=255)                                                                                                 taaGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCCAACGGTTTTTTAGGAGAACGGGTGGGGGGGGGTGGCAGACCGGCAGCGGATTTTTCCAAGGCGGCGGGAAACCGACCTCCTCTGTCCCGGTGTTGTGCGGTTTTCCGCGcc                                              >  2:643603/1‑151 (MQ=255)                                                                                                 taaGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGGTTTTTGAGGAGAACGGCTGAGGGCGGTTGCCAGACCGGCAGCGGATTTCTCCCGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGcc                                              >  1:594967/1‑151 (MQ=255)                                                                                                        ggCAGGATGGCCTGAGAGGCGGGTGTCTCTTTAAGCGTTATTTTAGGAGAACGGCTGAGTGGGGTTTCCCGACCGGCAGCGGATTTCTCCCGGGCGGCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGGTATCCAGCGCCCGCGcc                                        >  1:855093/1‑150 (MQ=255)                                                                                                           aGGATGGCCTGAGAGGCGGGGGGCTCTTCAAGCGGTATTTTAGGAGAACGGGTGGGGGGGGGTGGCAGACCGGCAGCGTATTTTTCCCGGGCGGCGGCAAACTGACCTGCTCTGTCCGGCTGTGGTGCGTTTTTCCGCGCCTGCGCCg                                       >  2:665232/1‑148 (MQ=255)                                                                                                                  ccTGAGAGGCGGGGGTCTCCTCAAGCGTTTTTTTAGGGGAAAGGCTGGGTGCGGGTGCCAGACCCGCCGCGGATTTTTCCCGGGCGTCGGCAAAATGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTccccgg                             >  2:535045/1‑147 (MQ=255)                                                                                                                        agGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATc                       <  1:883517/151‑1 (MQ=255)                                                                                                                           cGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATccc                     >  1:603985/1‑150 (MQ=255)                                                                                                                             ggTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAgg                  <  2:536506/151‑1 (MQ=255)                                                                                                                                gtCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCg                >  2:117428/1‑150 (MQ=255)                                                                                                                                  ctctTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTg              >  2:615787/1‑150 (MQ=255)                                                                                                                                    ctTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTg              <  1:573065/148‑1 (MQ=255)                                                                                                                                      tCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTgaag           <  2:273158/149‑1 (MQ=255)                                                                                                                                       cAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGaagaag        <  1:821889/151‑1 (MQ=255)                                                                                                                                            gTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAGAAGATcgc   <  1:226464/151‑1 (MQ=255)                                                                                                                                              tATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAGAAGATcgcg  <  1:57815/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGACCTGCTCTGTCCGGCTGTAGTGCGTTATCCAGCGCCTGCGCCGAGTTCCAGGTGGCATCCCAGGCGTGAAGAAGATCGCG  >  NC_000913/4301036‑4301326

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