Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,153,569	0	A	G	44.0%	56.9 / 13.3	50	R90R (CGA→CGG)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (22/0);  total (35/15)
Fisher's exact test for biased strand distribution p-value = 2.07e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAA  >  NC_000913/1153419‑1153719                                                                                                                                                       |                                                                                                                                                         ggCGCATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGa                                                                                                                                                          <  1:661636/151‑1 (MQ=255)        atCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAgggg                                                                                                                                                      >  1:528188/1‑147 (MQ=255)               gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAggggggaggacaa                                                                                                                                             >  2:683287/1‑140 (MQ=255)               gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGGAC‑Att                                                                                                                                           >  1:84324/1‑151 (MQ=255)                      ggATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGc                                                                                                                                    <  2:801136/151‑1 (MQ=255)                               tCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑ATTTAACa                                                                                                                                      >  1:703891/1‑140 (MQ=255)                                cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGcac‑aattaaaaacgg                                                                                                                                  >  1:386812/1‑132 (MQ=255)                                 aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGCAC‑ATTTaaaaaggggcc                                                                                                                               >  1:222030/1‑136 (MQ=255)                                      ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAAAAGCGGACCGGCCACACgg                                                                                                                     >  1:863990/1‑150 (MQ=255)                                           gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAAAAGCGGACCGGCAACGCGGAGtct                                                                                                                >  1:408145/1‑150 (MQ=255)                                           gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGCCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATtct                                                                                                                <  1:477904/150‑1 (MQ=255)                                                ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGt                                                                                                          <  2:559856/151‑1 (MQ=255)                                                  tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑AATGAACAGAGGACCGGCCACGCGGGGTTTCTcc                                                                                                            >  1:349202/1‑147 (MQ=255)                                                    tgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGCAC‑AATGAACAGCGGGCACGGAACGCGGGGtctctc                                                                                                             >  2:215250/1‑144 (MQ=255)                                                     gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACg                                                                                                      <  2:163092/150‑1 (MQ=255)                                                     gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCAC‑ATTTAACAACGGGCCGGGGAAGCGGGTTCTCTCCCTTTCgg                                                                                                     >  1:678304/1‑151 (MQ=255)                                                      aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑ATTTAAAAACGGACCCGGGACGCGGGGTCTCTCCCTTTCg                                                                                                      >  2:519477/1‑149 (MQ=255)                                                      aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAACAGCGGACCGGCAACGCGGGGTCTCTCCGTTTCgg                                                                                                     >  2:565608/1‑150 (MQ=255)                                                               cTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCAA‑AATTAAAAGCGGGCCGGGCACGCGGGGTTTCTCCGTTTCgggctattc                                                                                              >  2:711102/1‑141 (MQ=255)                                                                 cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGCCC‑ATTTAACAGCGGAGCGGGAACACGGGGTTTCTCCCCTTCTGGTTATCCTg                                                                                            >  1:594447/1‑148 (MQ=255)                                                                      tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCCC‑AATGAAAAGGGGACCGGGGACCCGGAGTCTCTCCGTTACGGCTTATCCTGCACAtt                                                                                      >  2:674334/1‑149 (MQ=255)                                                                        aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGGAC‑AATGGAAAGCGGACCGGGGACGCGGGGTCTCTTCCTTTCCGCTTAttatccacaa                                                                                       >  2:216926/1‑137 (MQ=255)                                                                        aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa                                                                                  <  2:691932/151‑1 (MQ=255)                                                                            agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑AATGTACAACGGACCGGGGAAGCGGGGTCTCTCCCTTTTGGGTTTTTCTGCACATTACg                                                                                   >  2:607029/1‑146 (MQ=255)                                                                                aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACAGCGGACCCGCAACGCGGGTTTTCTCCCTTTCGGCCTTTTCTTCCCATTACGAcccgtt                                                                            >  1:836484/1‑145 (MQ=255)                                                                                      aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                    <  2:252538/151‑1 (MQ=255)                                                                                       cGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGGACAATTGAA‑AACGGGCCGGCAACGCGGGGTTTCTCCGTTTCGGCTTATTCTGCACATTAAGACCGTTTGCg                                                                         >  2:181455/1‑145 (MQ=255)                                                                                        gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                  <  2:261736/151‑1 (MQ=255)                                                                                              agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGGAA‑ACTGAAAAGCGGACCGGCCACGCGGGGTCTCTCCCTTTCGGCCTTTCCTGCACACTACGCCCGCTTGCGTAAAGAgggggg                                                             >  2:190262/1‑148 (MQ=255)                                                                                                 agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAGCAC‑ATTGAAAAACGGACCGGCAACCCGGGGTCTCTCCGCTATGGCTTATTCTCCACAATACGACCGCTTGCGTGAAGAGGGGGTTACg                                                         >  2:739910/1‑151 (MQ=255)                                                                                                  gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt                                                        <  1:508103/151‑1 (MQ=255)                                                                                                   aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGAAC‑ATTGAAAAGCGGACCGGCAACGCGGGGTTTCTCCCTTACGGCTTATCCTGCACATTAAGACCGGTTGCGTAAAGAGGGGGTTACGtt                                                       >  1:499824/1‑151 (MQ=255)                                                                                                      aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgc                                                    <  1:216949/151‑1 (MQ=255)                                                                                                             tGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCAC‑AATTAAAAACGGGCCGGCAACGCGGGTTCTCTCCGTTTCGGCTTATTCTGCACATTACGACCGCTTGCGTGAAGAGGGGGGTACGTTTGCGC‑Caa                                              >  2:582587/1‑149 (MQ=255)                                                                                                              ggTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGc                                            <  2:354261/150‑1 (MQ=255)                                                                                                                tgtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCAC‑AATTTACAGCGGGCCGGCCACCCGGGTTCTCTCCCCTTCCGGTTTTCCCTCCCAATACCACCCTTTTCGTAACGGGGGGGTTAAGtttggggcaaacc                                            >  1:756868/1‑138 (MQ=255)                                                                                                                 gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACAGCGGGGCGGCAACGCGGGTTCTCTCCGTTACGGCTTATTCTGCACATTACGACCGGTTGCGTAACGAGGGGGGTTCGTTTGCGC‑TAAGCCCgg                                        >  1:330117/1‑151 (MQ=255)                                                                                                                 gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCg                                         <  1:120122/150‑1 (MQ=255)                                                                                                                        aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAAAAGCGGACCGGCAACGCGGGTTTTCTCCCTTTCCGCTTTTTCTGCACATTACGACCGTTTTCGGAACGAGGGGGTTACGTTGGCGC‑TAAACCCggggggg                                   >  1:66803/1‑144 (MQ=255)                                                                                                                           tCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGAC‑ATTTAAAAACGGACCCGGAACGCGGGGTCTCTCCGTTTCGGCTTATCCTGCACCTTAAGACCGGTTGCGTGACGAGGGGGGTACGTTGGCGC‑TGAGACCCGGGCGGCgggg                              >  2:110138/1‑151 (MQ=255)                                                                                                                             gTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGgc                            <  1:611026/151‑1 (MQ=255)                                                                                                                             gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGATAGGAC‑ATTGTACAGCGGGCCGGCAACACGGGGTCTCTCCCCTAATGGTTTTCCTGGACCATACGACCGTTTGCGTTACGAGGGGGGTACGTTTGGGCGTAGGCCCgggcgggttgggg                             >  1:672844/1‑140 (MQ=255)                                                                                                                                aTCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGCAC‑ATTTAAAAACGGACCGGCAAAACGGGGTCTCTCCCCTTCGGCTTTTTCTGCACACTACGACCGCTTTCGTAACGAGGGGGGTATGGTTGCGC‑GCAACCCCGGGCGGCGgggg                             >  2:566777/1‑147 (MQ=255)                                                                                                                                aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCa                         <  2:597076/151‑1 (MQ=255)                                                                                                                                  cAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAACAACGGGACGGCAACGCGGGTTCTCTCCGGTACGGCCTATCCTGCACATTACGACCGTTTGTGTAACGAGGGGGTTACGTTGGCGC‑TAAGCCCCGGGGGGCTgggggg                           >  1:83204/1‑145 (MQ=255)                                                                                                                                   aGTCGCGGTAGTGGCGGGCGGGGGGACAACAC‑ATTGAAAAGCGGGCCGGCCACGCGGGGTTTCTCCCCTTTGGCCTATTCTTCCCAATACGACCGGTTTCGGAAAGAGGGGGGTTCGTTGGCGC‑CAAGCCCCGGGCGGGGGGGgggcga                        >  1:659573/1‑147 (MQ=255)                                                                                                                                      cgcgGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTc                   <  2:607135/151‑1 (MQ=255)                                                                                                                                              tGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGc            >  2:277104/1‑150 (MQ=255)                                                                                                                                               ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGc            >  2:357622/1‑149 (MQ=255)                                                                                                                                               ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAg          <  1:441169/151‑1 (MQ=255)                                                                                                                                                      gAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa    >  1:598403/1‑150 (MQ=255)                                                                                                                                                       aGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaaa  <  1:399579/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         GGCGCATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAA  >  NC_000913/1153419‑1153719

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