Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,153,573	0	T	G	43.1%	73.0 / 30.3	51	Y92D (TAC→GAC)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (1/28);  new base G (22/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.47e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA  >  NC_000913/1153423‑1153721                                                                                                                                                       |                                                                                                                                                        cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                         <  1:453700/151‑1 (MQ=255)    atCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACa                                                                                                                                                      <  1:494572/151‑1 (MQ=255)          cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑At                                                                                                                                               <  1:313572/151‑1 (MQ=255)                   gATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACa                                                                                                                                         >  1:760377/1‑148 (MQ=255)                      gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGAAC‑Att                                                                                                                                              >  2:344073/1‑140 (MQ=255)                      gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAACAACgg                                                                                                                                     >  1:292607/1‑149 (MQ=255)                      gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGa                                                                                                                                    <  2:491514/150‑1 (MQ=255)                         tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGAcc                                                                                                                                  <  1:780359/149‑1 (MQ=255)                         tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCg                                                                                                                                 <  1:511323/150‑1 (MQ=255)                          gTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGc                                                                                                                               <  2:526575/151‑1 (MQ=255)                                ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAAcgcg                                                                                                                         <  2:283639/151‑1 (MQ=255)                                   tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGAAC‑AATGAAAAGAGGACCGGGAACGCgg                                                                                                                        >  2:604182/1‑149 (MQ=255)                                   tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACAATTGAA‑AACGGACCGGCCAAGCgg                                                                                                                        >  1:655054/1‑149 (MQ=255)                                    gatgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATTAACAGCGGGCCCGCAACGCGGGtt                                                                                                                     >  1:91468/1‑151 (MQ=255)                                          ttCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATtctctc                                                                                                                <  2:480880/150‑1 (MQ=255)                                           tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGAAC‑ATTGAACAAGGGAACGGCCACGCGGGTTTTCTcc                                                                                                               >  2:265525/1‑150 (MQ=255)                                              tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑AATGTAAAGCGGGCCGGCAAAGCGGGTTCTCTcccct                                                                                                            >  2:294573/1‑147 (MQ=255)                                               gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGtt                                                                                                            <  2:584206/149‑1 (MQ=255)                                               gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTAc                                                                                                          <  2:887182/151‑1 (MQ=255)                                                      cTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCAC‑ATTTTAAAAAGGACCGGCAACGCGGGGTTTCTCCGTTTCgg                                                                                                        >  1:793909/1‑146 (MQ=255)                                                        ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGACC‑ACTGAACAGCGGACCAGCAACGCGGGGTCTCTTCGCTACGGccttt                                                                                                   >  1:754191/1‑145 (MQ=255)                                                        ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCACGGC‑ATTTAACAGCGGACCGGCCAAGCGGGTGCTCCCCCTTTCTGGTTAt                                                                                                   >  1:118552/1‑149 (MQ=255)                                                            tCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGGCAGGGC‑AATCAAAAGGGGGCCGGCGACGCGGGTTTTCTCTTTTTCGGGTTAt                                                                                                   >  1:885424/1‑145 (MQ=255)                                                                   gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACg                                                                                      <  2:676228/151‑1 (MQ=255)                                                                    aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGCAC‑AATTAAAAGCGGGCGGGCAAAGCGGGGTCTCTCCGTTACGGCGTATCCTGCACAATACg                                                                                      >  2:270780/1‑150 (MQ=255)                                                                       aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAACAGCGGGACGGCAACGCGGGTTTTCTCCCTTTCGGCTTATTCTGCACATTACG‑Acc                                                                                  >  1:712648/1‑150 (MQ=255)                                                                            gtGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGc                                                                            <  2:154717/150‑1 (MQ=255)                                                                                       ggCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAAAAACGGGCCGGCCACGCGGGGTCTCTCCCTTACGGCTTTTCCTGCACATTACG‑ACCCCTTTCGTGACGAgggg                                                                 >  2:95215/1‑151 (MQ=255)                                                                                         cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑CTTTAACAACGGACCGGCAACCCGGGTTCTCTCCCTTTCGGCCTTTTCTGCCCCCTACC‑ACCCGTTTCGTAACGAggggg                                                                >  1:80967/1‑149 (MQ=255)                                                                                         cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGAt                                                               <  1:282166/151‑1 (MQ=255)                                                                                           gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑AATTAACAACGGACCCGCCACCCGGGTTTTCTCCCCTTCCGCTTTTCCCCCACAATAAC‑ACCGCTTTCGTAAAAAgggg                                                                 >  2:59314/1‑147 (MQ=255)                                                                                           gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCTGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTa                                                             <  2:628119/151‑1 (MQ=255)                                                                                           gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTa                                                             <  1:574559/151‑1 (MQ=255)                                                                                                  aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTgg                                                      <  2:457787/150‑1 (MQ=255)                                                                                                   gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGgcg                                                    <  2:592282/151‑1 (MQ=255)                                                                                                     aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCt                                                  <  1:245823/151‑1 (MQ=255)                                                                                                           gtgtGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGccc                                            <  1:780916/151‑1 (MQ=255)                                                                                                           gtgtGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGccc                                            <  1:193547/151‑1 (MQ=255)                                                                                                                   aaaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGTAAAACGGACCGGCAACGCGGGTTCTCTCCGTTACGGCTTATCCTCCACACTACG‑ACCCTTTGCGTAACGAGGGGGTTTCG‑TTTGGGCTAAGCCCggggggg                                     >  1:550961/1‑145 (MQ=255)                                                                                                                   aaaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGc                                    <  1:939630/151‑1 (MQ=255)                                                                                                                       tCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGAAGCCC‑ATTTGACAGGGGAGCGGCCACGCGGCTTCTCTCCCTTACGGCTTATTCTGCCCATTACG‑ACCGTTTGCGTCACGAGGGGGTTACG‑TTGGGGCCAAGACCCGTGCGGCTggg                                >  1:810518/1‑151 (MQ=255)                                                                                                                         gTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGgc                              <  1:527485/151‑1 (MQ=255)                                                                                                                           gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGgcgc                            <  2:623891/151‑1 (MQ=255)                                                                                                                            aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCa                           <  1:752538/151‑1 (MQ=255)                                                                                                                             tCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAAAAACGGACCGGCGAAGCGGGATTTTTCCGTTTCGGGCTTTACACCCCATTAAGAACCCTTTG‑GGAACCAGGGGGGTACG‑TTTGGGCTAAACCCCGGGCGGCTgggggg                             >  2:423737/1‑146 (MQ=255)                                                                                                                               aGTCGCGGTAGTGGCGGGCGGGGGGACAACAC‑AATGTAAAGCGGGGCGGCCACGCGGGGTCTCTCCCTTTCGGCCTTTCCTGCACCTTACG‑ACCCGTTTCGTGACGAGGGGGGTAAGTTTGGGCCTAAGCCCGG‑GCGGGTGGGgggcga                          >  2:922844/1‑147 (MQ=255)                                                                                                                                gTCGCGGTAGTGGCGGGCGAGGGGACAGCCC‑ATTGAAAAACGGGCCGGCAACGCGGGTTCTCTCCCCTTCGGCTTTTCCCGCCCATTACG‑ACCGTTTTCGGAACGAGGGGGTTTCG‑TTGGCGCTAAGCCCCGGGCGGCGGGGgggccaacc                       >  2:134185/1‑146 (MQ=255)                                                                                                                                        aGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTg               <  1:359573/151‑1 (MQ=255)                                                                                                                                           ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCa             <  2:72595/150‑1 (MQ=255)                                                                                                                                           ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCa             <  1:92237/150‑1 (MQ=255)                                                                                                                                                  gAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa      >  1:294931/1‑150 (MQ=255)                                                                                                                                                     ggTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACa  <  1:42160/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACG‑ACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA  >  NC_000913/1153423‑1153721

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