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.4%	65.9 / 11.5	46	V316G (GTA→GGA)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base T (4/21);  new base G (20/0);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.95e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTG  >  NC_000913/2252504‑2252796                                                                                                                                                |                                                                                                                                                          ggTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCt                                                                                                                                                   <  2:44329/151‑1 (MQ=255)  gtgtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCGGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCtt                                                                                                                                                  <  2:437985/151‑1 (MQ=255)   tgtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGG‑CTGGGATATAGCttt                                                                                                                                                 >  2:140050/1‑151 (MQ=255)    gtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTc                                                                                                                                                <  2:212648/151‑1 (MQ=255)      tCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCt                                                                                                                                               <  1:348518/150‑1 (MQ=255)            aCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGG‑CTGGGATATAACttttttttagg                                                                                                                                         >  2:716688/1‑145 (MQ=255)              tgtttgttTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCTGGGTATAGCtttttt                                                                                                                                              >  1:326628/1‑143 (MQ=255)              tgtttgttTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGggtataaatttttt                                                                                                                                              >  2:132826/1‑135 (MQ=255)               gtttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGAt                                                                                                                                       <  2:177257/149‑1 (MQ=255)                  tgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGGATATAGCTTTTTTGTAGGATACAG‑t                                                                                                                                 >  2:322399/1‑151 (MQ=255)                   gtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGTTATAGGTTTTTTTTAGGGTACAA‑tt                                                                                                                                >  2:195105/1‑151 (MQ=255)                    tttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTGTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑tg                                                                                                                                <  1:253283/150‑2 (MQ=255)                          tGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGAAAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGCCGCTGGAAT‑TAAATTTCTTGGAGGTTAAAG‑GTTCcgtgc                                                                                                                         >  1:649313/1‑148 (MQ=255)                               gATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAgcg                                                                                                                    <  2:282336/151‑1 (MQ=255)                                 tGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGGGTATAAATTTTTTTTAGGGTACAG‑ATTTCCCGCCAGCGcc                                                                                                                  >  1:56654/1‑151 (MQ=255)                                  gAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGccc                                                                                                                 <  1:449316/151‑1 (MQ=255)                                      cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGG‑CCGGGAAATAACTTTTTTTAAGGAAAAAG‑GTTCCGCCCCAGCGCCCG‑gtt                                                                                                            >  1:6460/1‑151 (MQ=255)                                          ttATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGG‑CCGGGATATAGCTTTCTTTTTGGGTACAG‑TTTTCTCGGCAGCGcccctcgt                                                                                                            >  1:177332/1‑143 (MQ=255)                                          ttATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGG‑CCGGGATATAACTTTTTTTTAAGATAAAG‑GTTTCCCCACAGCGCCCC‑tttt                                                                                                           >  2:545535/1‑148 (MQ=255)                                                    gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGGTATTACTTTCTTTTTGGATACAG‑GTTTCCCGGCAGCcccccctttttgggagcc                                                                                                   >  1:393687/1‑133 (MQ=255)                                                                        tttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGGATATAACTTTCTTGAAGGAGAAAG‑ATTCCGCGGCAACGCCCG‑TTTTGCGGTAGCCTGCGCGGCAAATTCGCGCGcc                                                                             >  2:809804/1‑148 (MQ=255)                                                                          tGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGGATATAGCTTTTTTTTAGGGTACAGTTTTCCCCGCCA‑CCCCCC‑CTTTGCCGGGGCCCCCCCGGCAATTTCGCGCGCCacc                                                                          >  1:552642/1‑147 (MQ=255)                                                                              gACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGGATATAGCTTTTTTTTAGGATAAAA‑TTTTCCCGACAGCGCCCC‑TTTTTCGGTAGCCTGCGCGGGCATTTTGCGCGCCCGCgg                                                                        >  2:444369/1‑147 (MQ=255)                                                                                    gcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAaac                                                             <  1:162297/151‑1 (MQ=255)                                                                                            aTACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAaccctac                                                     <  1:353285/151‑1 (MQ=255)                                                                                              aCGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAaccctaccc                                                   <  2:599515/151‑1 (MQ=255)                                                                                               cGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGATATAGCTTTTTTTGAGGGTAAAA‑TTTTCGCGGCAGCGCCCG‑CTTTTCGGGAGCCTGCGCGGCGATTTCGCGCGCCCGCGGCCG‑CTTAAAACaccctcacc                                                   >  1:176583/1‑146 (MQ=255)                                                                                                   aGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCt                                              <  2:294403/151‑1 (MQ=255)                                                                                                      aGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGc                                            <  1:709318/150‑1 (MQ=255)                                                                                                       gACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGATATAACTTTTTTTTAGGATAAAG‑TTTTCCCCGCAACGCCCC‑GTTTTCGGGAGCCTGCGCGGCCATTTTTCGCGCCCCCGGGAGCCTTAAAAAACCC‑ACCCCCTCTGcc                                           >  2:537604/1‑150 (MQ=255)                                                                                                        aCAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGATTTAGCTTTTTTTAAGGGTAAAA‑TTTTCGCGGCCACGCCCC‑CTTTTCGGGAGCCCGCGCGGCCATTTTGCCCGCCCGCGGCCG‑CTTAAACAACCCCCCCCCGTCt                                              >  2:174760/1‑146 (MQ=255)                                                                                                               tGGTGCTGACGATGGTTTTATGGCGGGC‑CCGGGATATAGCTTTTTTTTAGGGTACAG‑TTTCCCCCACAGCGCCCG‑CTTTTCGGTAGCCTGCGCGGCAATTTTGCGCGCCCGCGGCAG‑CTTTAACAACCCCACCCCGTCTGCCGATaa                                      >  2:853306/1‑147 (MQ=255)                                                                                                                ggTGCTGACGATGGTTTTATGGCGGGC‑CTGGGATATAGCTTTCTTTAAGGATAAAG‑GTTTCGCGACCACGCCCC‑CTTTTTGGGAGCCTGCGCGGCAAATTTTCGCGCCCCCGGGCG‑CTTAAACAAACCCACCCCCTCTTCCGGTaa                                      >  2:101006/1‑146 (MQ=255)                                                                                                                 gTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGGTTTTTTTTAGGGTAAAG‑TTTTCGCCGCAGCGCCCG‑TTTTGCGGTAGCCTGCGGGGCAATTTTTCGCGCCCGCGGCAGCCTTAAAAAACCC‑ACCCCGTCTGcccctaaaac                                   >  2:86852/1‑140 (MQ=255)                                                                                                                    cTGACGATGGTTTTATGGCGGGG‑CCGGGGTTTAACTTTTTTTTAGGAGACAG‑GTTTCCCGGCAGCGCCCGCTTTTGGGGT‑GCCCGCGCGGCCAATTTTCGCGCCCGCGGGCG‑CTTAAAAAAACCCACCCCTTCTGCCGCTaaaacgcct                               >  1:282137/1‑142 (MQ=255)                                                                                                                      gACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCa                            <  2:176616/150‑1 (MQ=255)                                                                                                                          aTGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGt                       <  1:217621/151‑1 (MQ=255)                                                                                                                            ggTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGc                     <  2:393727/151‑1 (MQ=255)                                                                                                                             gTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGcc                    <  2:385431/151‑1 (MQ=255)                                                                                                                              ttttATGGCGTGC‑CTGGTATAGAGTTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGTCCG‑TTTTGCGTTATCCTTCGCGGCAATTTCGCGCGTCACCGGCAG‑CTTAAACAACCCTTCCCTGTTTGCCGATAACGCGTTTTCATTAgc                       <  1:319867/147‑2 (MQ=255)                                                                                                                                  aTGGCGGGC‑CTGGGATATAGCTTTCTTTAAGGAAACAG‑TTTCCGCGACAGCGCCCG‑GTTTTCGGTAGCCTGCGCGGCAATTTTGCGCGCCAGCGGCAG‑CTTTAACAAACCCACCCCGTCTGCCGATAAAGCGCCTTCCTTAATGCCCATgg               >  1:211410/1‑151 (MQ=255)                                                                                                                                    ggCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtg             >  1:821083/1‑151 (MQ=255)                                                                                                                                      cGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtgt            >  1:297690/1‑150 (MQ=255)                                                                                                                                        tgC‑CTGGTATATAGCTTTATTGAAGGATACAG‑GTTCCGCGACAGCGACCG‑TTTTGCGGTATCCTGCGCGGCAATTGCGCGCGCCTGAGGCAG‑TTTAAAAAACCCTACACTGTCTGCCGATAACGCGCTTTCATTAGAGCTAATGGGGTaa          <  1:500776/149‑1 (MQ=255)                                                                                                                                         gC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACa        <  1:4662/151‑1 (MQ=255)                                                                                                                                          c‑cTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACaa       <  1:444311/151‑1 (MQ=255)                                                                                                                                             tGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATg     >  2:125014/1‑151 (MQ=255)                                                                                                                                                tatataGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTg  >  1:32666/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                          GGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGC‑CTGGTATATAGCTTTCTTGAAGGATACAG‑TTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTG  >  NC_000913/2252504‑2252796

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