Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	503,834	0	T	C	48.1%	52.6 / 24.7	54	D288G (GAT→GGT)	fsr	putative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (23/5);  new base C (0/26);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 9.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGAA  >  NC_000913/503687‑503976                                                                                                                                                      |                                                                                                                                                 aGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑ccc                                                                                                                                              <  2:16029/151‑1 (MQ=255)    aTGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTa                                                                                                                                           <  1:227384/151‑1 (MQ=255)        cAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAgg                                                                                                                                       >  1:697752/1‑151 (MQ=255)         agagAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑c                                                                                                                                     >  2:141355/1‑151 (MQ=255)                  cGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACACCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGa                                                                                                                              >  1:503526/1‑149 (MQ=255)                              aaaTCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGcc                                                                                                                >  2:14106/1‑151 (MQ=255)                              aaaTCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGcc                                                                                                                >  2:241621/1‑151 (MQ=255)                                  caaataacCACCTTAAAACCCCCCGACCAGTGAGGGCGGGGGTAGGGTAAAATCAGCTAAAACGGCGCAACGCCGAGGATAGAGCCCCAACCCCCATTTTCCCCCCCAT‑TTTT‑CC‑CCCCTCCGGG‑CCCGCCGATCACCGTCCCTGCCGc                                                                                                              <  1:659807/147‑1 (MQ=255)                                       aatCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGc                                                                                                        >  2:19297/1‑150 (MQ=255)                                          cACCGTTAAAACCCCCGTCCAGTGCAGGCTGTCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAAc                                                                                                    >  2:743016/1‑151 (MQ=255)                                            ccGTAAAACCCCCCGTCCAGTGCGGGCGGGCGTGGGGAAAAATCAGCTTAACCGGCGCACCGCCGAGGATAGACCCCAAAACCACATATTCCCCCCAAA‑TTTT‑AC‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAg                                                                                                  <  2:412828/151‑1 (MQ=255)                                                tAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCATCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATAACCGTACCTGCCGCAACGGCAAACa                                                                                                   >  2:787641/1‑146 (MQ=255)                                                 aaaaCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAgg                                                                                             >  1:265461/1‑151 (MQ=255)                                                  aaccccccTTCCAGGGCAGGCTGGCTTAGGATAAAACCAGCGAAACGGGCGCACCGCCGAGGATAGACCCCAAAATCACATATTCCCCCCAAA‑TTTT‑AC‑CCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                            <  2:253501/148‑1 (MQ=255)                                                  aaccccccGTCCAGGGCAGGCGGGCGTAGGGAAAAACCAGCGTAACCGGCGAACCGCCGAGGATAGAGCCCAAAACCAAATATTCCCCCCCAA‑TTTA‑AC‑CCCCTACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                            <  2:586639/148‑1 (MQ=255)                                                   aaCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                            >  2:298802/1‑150 (MQ=255)                                                   aaCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                            >  2:540675/1‑150 (MQ=255)                                                    aCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCa                                                                                           >  2:205497/1‑150 (MQ=255)                                                     cccccGCCCAGTGCGGGCTGGCGTAGGGAAAAACCAGCAAAACCGGCGCACCGCCGAGGATAGGCCCCAAAACCACATATTCCCCCCAAA‑TTTT‑AC‑CCCCAACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCaaa                                                                                         <  1:389948/151‑1 (MQ=255)                                                       cccGCCCAGTGCGGGCGGGCGAGGGGAAAAACCAGCGTAACCGCCGCACCGCCGAGGAAAGACCCCCAAACCACAT‑TTTCCCCCCAA‑ATTT‑TC‑CCCCCAAAGGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACa                                                                                       <  1:193650/151‑1 (MQ=255)                                                          gTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAAc                                                                                        >  2:608435/1‑147 (MQ=255)                                                               gTGCGGGCGGGCGTAGGGAAAAACCAGCGAAACCGGCGCACCGCCGAGAATAGACCCCAAAACCACATTTTCCCCCCAAA‑TTTT‑AC‑CCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGaa                                                                                 <  2:138033/149‑1 (MQ=255)                                                                tGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCt                                                                              >  1:539757/1‑151 (MQ=255)                                                                    ggcgggcgTGGGGAAAAACAAGCGAAACCGGCGCACCGCCGAGGATAGCCCCCAAAACCACATTTCCCCGCCAAA‑TTTT‑AC‑CCCCTACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGc                                                                          <  1:677002/151‑1 (MQ=255)                                                                              gggAAAAACCAGCGAAACCGGCCCACCGCCGAGGATAGACCCCAAAACCACATATTCCCGCCAAA‑TTTT‑AC‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTgg                                                                   <  2:867628/148‑1 (MQ=255)                                                                              ggaaaaaaTCAGCGAAACCGGCCCACCGCCGAGGAGAGCCCCCCAACTCACATATTCCCCCCCAT‑TTTT‑AC‑CCCCAACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                                  <  1:205463/145‑1 (MQ=255)                                                                                  aaaaTCAGCAAAACCGGCGCACCCCAGAGAAAAGGCCCCAAAACCACAT‑TTTCCCCCCAATTTTT‑CC‑CCCCAACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGgatag                                                               <  1:73813/148‑1 (MQ=255)                                                                                  aaaaACAGCGAAACGGGCGCACCGCGGAGGATGGACCCCAAAATCAC‑TATTTCCCCCCAATTTTT‑AC‑CCCCTACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGgata                                                                <  1:524198/147‑1 (MQ=255)                                                                                   cgctctgCGAAACCGGCCCAACGCGGAGGAAAGCGCCCCAAACCCAATTTTCCCCCCAAA‑TTTT‑AC‑CCCCGACGGC‑CCCGCCGATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                                  <  1:590447/138‑1 (MQ=255)                                                                                   aaataagcaaaacgggCGCACCGCGGAGAAAAGACCCCAAAACCACATATTCCCGCCCAT‑TTTT‑CC‑CCCCAACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATaa                                                           <  1:375003/137‑1 (MQ=255)                                                                                        aGCGTAACCGGCGCACCCCCGAGAAAAGCCCCCAAAACCACTTTTTCCCCCCCAA‑TTTT‑AC‑CCCCAACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATcc                                                        <  1:183284/149‑1 (MQ=255)                                                                                         gCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                                  >  2:590521/1‑138 (MQ=255)                                                                                         gCGTAAACGGCGCAACCCCGAGGATAGACCCCAAAACCACATATTTCCCCCCAA‑TTTT‑AC‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGaa                                                     <  2:251095/151‑1 (MQ=255)                                                                                             aaaCGGCCCAACGCCAAGGATAAAGCCCCAAAACACATTTCCCCGCCCAA‑TTTT‑AC‑CCCGTACGGG‑CCCGCAGATCACCGTGCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAtt                                                   <  2:505735/149‑1 (MQ=255)                                                                                              aaCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCAGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAttttt                                                >  1:313279/1‑151 (MQ=255)                                                                                                cGGCGAAACGCCGAGGATAGACCCCAAAACCCCAT‑TTTCCCCCCAATTTTT‑AC‑CCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTg                                               <  2:697841/150‑1 (MQ=255)                                                                                                    acAACCCGGAGGATAGCCCCCAAAATACCATATTCCCCCACAT‑TTTA‑AC‑CCACTAAAGC‑CCCGCATATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCAt                                            <  1:141331/148‑1 (MQ=255)                                                                                                     caccGCCGAGGATAGACCCCCAAACCACATTTTCCCCCCCAA‑TTTT‑AC‑CCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCa                                          <  2:308734/147‑1 (MQ=255)                                                                                                          ccGAGGATGGACCCCAAAATCACATATTCCCGCCAAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaa                                     <  2:266401/150‑1 (MQ=255)                                                                                                                 tAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtgt                              >  1:563221/1‑150 (MQ=255)                                                                                                                      ccccAAAGCCCAATTTTCCCGCCAA‑TTTTCAC‑CCCGAAGGGG‑CCCGCGGGTCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATa                          <  1:792702/150‑1 (MQ=255)                                                                                                                      ccccAAAACCAATAATTCCCCCCAA‑TTTT‑TTCCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAg                         <  2:773528/151‑1 (MQ=255)                                                                                                                      cccaAAATCACATATTCCCCCCAAA‑TTTT‑CC‑CCCCAACGGC‑CCCGCCGATCACCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                           <  1:671438/148‑1 (MQ=255)                                                                                                                      cccaAAACCCCATATTCCCCCCAAA‑TTTT‑AC‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtgt                              <  2:242199/145‑1 (MQ=255)                                                                                                                       cccAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgct                       >  2:56402/1‑151 (MQ=255)                                                                                                                         cAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGTAAACAGGAAAGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTAGCATCAGATAAAAGGTGTAATAgctg                      >  1:405025/1‑150 (MQ=255)                                                                                                                         cAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgc                     >  1:87695/1‑151 (MQ=255)                                                                                                                          aacccacaTTTTCCCCCCCAA‑TTTT‑AT‑CCCATAGAGG‑CCCGCCGATCACCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtgt                              <  2:170798/137‑1 (MQ=255)                                                                                                                              caaatatTCCCACCCAT‑TTTT‑AC‑CCCCTACAGG‑CCCGCCGATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGAt                 <  2:503586/147‑1 (MQ=255)                                                                                                                                 atatTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTAACTGCCGCAACGGCAAAAAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTc             >  2:660625/1‑151 (MQ=255)                                                                                                                                       cccccaaa‑TTTT‑AC‑CCCAAACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCata       <  1:278418/145‑1 (MQ=255)                                                                                                                                       ccGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCat        >  2:513759/1‑150 (MQ=255)                                                                                                                                           ccaattttt‑AC‑CCCCTACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAg    <  1:428937/146‑1 (MQ=255)                                                                                                                                             aa‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGaa  >  1:643339/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTT‑AT‑CCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGAA  >  NC_000913/503687‑503976

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