Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,252,646	0	T	G	44.7%	55.7 / 13.5	47	V316G (GTA→GGA)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base T (8/18);  new base G (21/0);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 3.76e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTA  >  NC_000913/2252500‑2252787                                                                                                                                                   |                                                                                                                                                   gCAAGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGtatat                                                                                                                                                <  1:608269/151‑1 (MQ=255)      gtgtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGGTATAGGtt                                                                                                                                           >  1:831126/1‑151 (MQ=255)       tgtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGGTTTAGCttt                                                                                                                                          >  1:586639/1‑151 (MQ=255)          tCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGGTATAACttttt                                                                                                                                        >  2:925484/1‑148 (MQ=255)          tCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTTGtatatatgtttttt                                                                                                                                       >  2:538940/1‑143 (MQ=255)           cAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGgatatttttttttt                                                                                                                                       >  2:575334/1‑141 (MQ=255)           cAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCtt                                                                                                                                       <  2:879633/150‑1 (MQ=255)            aGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑a                                                                                                                                    <  1:18553/151‑1 (MQ=255)               aaCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATAACttttttttagg                                                                                                                                  >  2:835312/1‑146 (MQ=255)                 cTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTGTATAGCttttttttagg                                                                                                                                  >  2:167952/1‑144 (MQ=255)                     ttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACa                                                                                                                            <  2:878992/150‑1 (MQ=255)                          tATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTc                                                                                                                       <  2:102401/150‑1 (MQ=255)                            tttGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTTTATAACTTTTTTT‑TAGGAAATAGttttc                                                                                                                      >  2:876103/1‑147 (MQ=255)                              tGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCCGGATATAGCTTTTTTT‑TAGGATACAGatt                                                                                                                        >  2:337150/1‑145 (MQ=255)                                  ccattgatccgTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGa                                                                                                                  <  1:131811/137‑1 (MQ=255)                                  cGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAgc                                                                                                              <  2:586834/151‑1 (MQ=255)                                   gATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAgcg                                                                                                             <  1:92189/151‑1 (MQ=255)                                   gATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATAGGTTTCTTT‑TAGGAGACAGttttctccgcagc                                                                                                              >  1:800579/1‑140 (MQ=255)                                          cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGttt                                                                                                     <  1:115906/151‑1 (MQ=255)                                                    cAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGca                                                                                           <  1:838248/151‑2 (MQ=255)                                                        gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATATCTTTTTTT‑TAGGGGACAGTTTTCGCGACAGCG‑CCCGCTTTGCGGGAGcc                                                                                           >  2:909495/1‑147 (MQ=255)                                                        gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCGGGGATAAAGATTTTTTT‑TAGGGTACAGATTCCCCGAAACAG‑CCCGGTTTGCg                                                                                                 >  2:619238/1‑141 (MQ=255)                                                         aaaaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTTTTT‑TAAGATAAAGTTTTCTCCACAAAG‑CCCGttttttgg                                                                                                >  2:149830/1‑137 (MQ=255)                                                           aaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATAACTTTTTTT‑GAGGATACAGTTTTCCCGACAGAG‑CCCGCTTTTCGGTAGcc                                                                                           >  2:820349/1‑144 (MQ=255)                                                            aGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGc                                                                                  <  1:113683/151‑1 (MQ=255)                                                            aGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTGTATAGCTTTCTTT‑GAGGATACAAATTTCGCGACAGGG‑CCCGGTTTGCGGTTGTCTCCGC‑gg                                                                                   >  1:315741/1‑150 (MQ=255)                                                              gATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGTTATAGCTTTTTTG‑TAGGATAAAGTTTTCGCGGCAGCG‑CCCGCTTTGCGGGGGCCTGCGC‑gg                                                                                   >  1:268981/1‑148 (MQ=255)                                                               aTGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAAt                                                                               <  2:580994/151‑1 (MQ=255)                                                                         aTTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGGTTTATCTTTTTTT‑TAGGATAAAGTTTTCCCGGCCGCC‑CCCCCTTTTCGGTTGCCCTCGC‑GGCGAtttttcccccc                                                                     >  1:834507/1‑143 (MQ=255)                                                                          tttttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGATATAACTTTTTTT‑AAGGGTAAAGTTTTCCCGGCAGCG‑CCCGGTTTTCGGGAGCCTGCGC‑GGCAAtttttcccgcc                                                                     >  2:196258/1‑142 (MQ=255)                                                                            tttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGGTATAGAGTTTTTT‑TAGGGTAAAGCTTTTGCGGCCCCG‑CCCGCTTTTTGGTTGCCCTCGC‑GGCCATTTccccgg                                                                       >  1:362895/1‑141 (MQ=255)                                                                             ttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGc                                                                  <  2:982456/150‑1 (MQ=255)                                                                                               cATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTCGGATATAACTTTCTTT‑TAGGGTACAGTTTTCCCGGCAGCG‑CCCGTTTTTCGGGAGCCTGCGCGGGCAATTT‑TCGCGCCCGCGGCAG‑CGTTAaaaaac                                                 >  2:986683/1‑148 (MQ=255)                                                                                                aTACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAacccta                                             <  2:1010421/150‑1 (MQ=255)                                                                                                  aCGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGATATAGCTTTCTTT‑TAGGAGAAAGCTTTTCCGACAACG‑CCCGTTTTTCTGGATACTTCGC‑GGCAATTTCTTGCCCCCGCCGCCG‑C‑CTAAAAAaaccccccc                                          >  1:964574/1‑146 (MQ=255)                                                                                                   cGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAaccctaccct                                         <  2:383649/151‑1 (MQ=255)                                                                                                    gACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAACTTTTTTT‑TAGGGTAAAGTTTTCCCGGCAGCG‑CCCGGTTTTCGGTAGCCTGCGC‑GGCAATTTTGCGCGCCCCCGGCGG‑C‑TTAAAAAaacccaccc                                          >  2:754640/1‑149 (MQ=255)                                                                                                    gACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATAGCTTTTTTT‑TAGGATACAAATTTCGCGGCAGCG‑CCCGCTTTTCGGGGGCCTGCGC‑GGGAATTTTTCGCGCCCGCGGCCG‑G‑TTAAAAAaaccccccc                                          >  1:961428/1‑144 (MQ=255)                                                                                                      cAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGTAAGGA‑AAAGTTTTCGCGACAACG‑CCCGCTTTTCGGGAGCCCGCGC‑GGCAAATTTTAGCGCCAGCGGGGG‑C‑TTAAAAAaacccacc                                           >  2:592730/1‑146 (MQ=255)                                                                                                       aGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTTGGATATAGCTTTCTTT‑AAGGATACAGTTTTCGCGCCAGCG‑CCCGTTTTTCGGTTGCCTGCGC‑GGGGAATTTGTGCGCCAGCGGCAG‑C‑TTTaaaaaaaccacccc                                         >  1:432567/1‑138 (MQ=255)                                                                                                       aGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTATATAACTTTTTTT‑AAGGAGACAGTTTTCGCGACAGCG‑CCCGTTTTGCGGGAGCCCGCGC‑GGCAAATTTTCGCGCCCGCCGCAG‑C‑TTAAACAaccccaccc                                          >  2:42375/1‑146 (MQ=255)                                                                                                       aGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGTTTTAGCTTTTTTT‑TAGGTTAAAGGTTCCCCGAAAGCC‑CCCGGTTTTTGGTAGCCT‑CGCGGGCAATTTTCCGCGCCCGGGGCGC‑C‑GTTAAAAaccctcccc                                          >  2:429936/1‑146 (MQ=255)                                                                                                           gACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCg                                 <  2:706816/151‑1 (MQ=255)                                                                                                            aCAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGa                                <  2:883249/151‑1 (MQ=255)                                                                                                                  ttGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCCTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGATAAc                            <  1:63931/149‑1 (MQ=255)                                                                                                                   tGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATAAGTTTCTTT‑TAGGGTACAGATTTCGCGACAACG‑CCCGTTTTTCGGTAGCCCGCGC‑GGCACTTTTGCGCGCCCGCGGCCG‑C‑TTTAACAAACCCCCCCCGTCTGCCGCTaa                             >  2:398555/1‑147 (MQ=255)                                                                                                                         tGACGATGGTTTTATGGCGGGCCTGGGATATAGCTTTTTTT‑TAGGGTAAAGGTTTCGCGGCAGCG‑CCCGCTTTTCGGGAGCCTTCGC‑GGGAATTTTGCGCGCCCGCGGCCGCC‑TTAAAAAACCC‑ACCCCGTCTTCCGGTAAAGCGCctt                     >  2:590915/1‑149 (MQ=255)                                                                                                                           aCGATGGTTTTATGGCGGGCCCGGGATATAAATTTTTTT‑TAGGGTAAAGTTTTCTCGGCAGCGCCCCGTTTTTGGGGGGCC‑GCGC‑GGCGATTTTTCGCGCCCGCGGGGG‑CGTTAAAAAAACC‑ACCCCCCCTGCCGCTAACGCGCCTTCttt                 >  1:611111/1‑151 (MQ=255)                                                                                                                                ggTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTg             <  1:228943/150‑1 (MQ=255)                                                                                                                                           gggCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTa  >  1:359943/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCAAGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTG‑AAGGATACAGTTTCCGCGACAGCG‑CCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAG‑C‑TTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTA  >  NC_000913/2252500‑2252787

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