Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,252,660	0	G	T	48.5%	41.3 / 19.3	33	E321* (GAA→TAA)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base G (3/14);  new base T (16/0);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCG  >  NC_000913/2252518‑2252805                                                                                                                                                   |                                                                                                                                                      gtttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATa                                                                                                                                                <  2:403172/150‑1 (MQ=255)  tttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGGGT‑AT‑AGA‑TTTTTTTTAGGATaa                                                                                                                                               >  1:685127/1‑149 (MQ=255)         aTTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCTGG‑AT‑AT‑AGG‑TTTTTTTTAGGATACAGatt                                                                                                                                          >  2:520462/1‑148 (MQ=255)                cGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGGT‑AT‑ATC‑TTTTTTTTAGGGTATAATTTTCGCGAc                                                                                                                                   >  1:69623/1‑148 (MQ=255)                     aGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATAAT‑AGT‑TTTTTTTATGGAT‑CAGTTTTCTCGACAGCGcc                                                                                                                             >  2:608234/1‑149 (MQ=255)                        cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑ttt                                                                                                                       <  1:635298/151‑1 (MQ=255)                           tttATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGTT‑AT‑AGC‑TTTTTTTTAGGGTACAGTTTTCGCGACAGCGCCCG‑gttt                                                                                                                      >  2:206263/1‑149 (MQ=255)                             tATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGTAG‑TT‑AAA‑TTTTTTGTAGGAGAAAGTTTTTTCGACAGCGCCCG‑cttttggg                                                                                                                  >  1:263910/1‑147 (MQ=255)                                     cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTgc                                                                                                          <  1:409842/151‑1 (MQ=255)                                     cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGAT‑AT‑AGC‑TTTTTTTGAGGATACAGGTTCCGCGCCAACGCCCC‑TTTTTCGGGAGcc                                                                                                             >  1:191989/1‑148 (MQ=255)                                      gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGGAT‑AT‑AAA‑TTTTTTTGAAGATACAGTTTTTTCGACCACGCCCG‑GTTTGCGGGAGGCcgcg                                                                                                         >  2:473958/1‑151 (MQ=255)                                       aaaaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGGT‑AT‑TGT‑TTTTTTTAAGGATAAAGATTTTGCCACAGCGCCCGCTTTTCCGGtgac                                                                                                              >  2:624863/1‑143 (MQ=255)                                       aaaaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGGT‑AT‑AAC‑TTTTTTTTAAGGTAAAGTTTTCGCGCCAGCGGCCC‑CTTTTCGGTAGcc                                                                                                             >  1:406063/1‑146 (MQ=255)                                            gATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGAT‑TT‑AAC‑TTTTTTTTAGGATAAAGTTTTCGCGGCAGCGCCCG‑CTTTTCGGGGGCCTCCGC‑gg                                                                                                     >  2:665518/1‑148 (MQ=255)                                                          tttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCGGGTGT‑AT‑AAC‑TTTTTTTTAGGGGACAGGTTTCGCGACCGCGGCCG‑GTTTTCGGTTGCCTGCGC‑GGCAAatttccccgcc                                                                                       >  2:633175/1‑140 (MQ=255)                                                                gACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAg                                                                              <  2:434810/150‑1 (MQ=255)                                                                gACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGGT‑AT‑AGCTTTTTTTTAGGGGAAAAGTTTCC‑CGCCAGCGCCCG‑CTTTTCGGGGGCCTGCGC‑GGGCATTTTTCCCGCCCG‑Cgg                                                                                 >  1:233065/1‑147 (MQ=255)                                                                 aCCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGTAT‑ATAAGT‑TTTTTTTAGGGAAAAATTTTCCCCG‑CAACGCCCC‑GTTTTCGGGAGGCCGCGCGGGCAATTTC‑CGCGCCCG‑ccg                                                                                 >  2:478553/1‑144 (MQ=255)                                                                  ccGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑Ctt                                                                          <  2:531990/151‑1 (MQ=255)                                                                              aTACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCTGG‑AT‑TT‑AGA‑TTTTTTTTAGGAGACAGTTTTTGCGGCAGCGCCCC‑GTTTTTGGGGGCCTGCGC‑GGCCATTTTTCGCGCCCG‑CGGCGGCCTTAAAAAacc                                                                  >  1:480531/1‑148 (MQ=255)                                                                                      gTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCCGGGAT‑AT‑AAC‑TTTTTTTTAGGGTAAAGTTTTCGCGGCCGCGCCCG‑CTTTGCGGTAGACTTCGC‑GGCGATTTTTCGCCGCCG‑CGGCCG‑CCTATACAaacccaccccctc                                                        >  1:547548/1‑145 (MQ=255)                                                                                           cAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGAt                                                 <  2:268786/151‑1 (MQ=255)                                                                                            aGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATa                                                <  1:193427/151‑1 (MQ=255)                                                                                                ttGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAAcg                                             <  2:112740/150‑1 (MQ=255)                                                                                                       tGACGATGGTTTTATGGC‑GGGGCCGGGGT‑TT‑AAA‑TTTTTTTTAGGAGACAATTTTCGCGGCCGCGCCCC‑CTTTTCGGTAGGCTGCGC‑GGCCATTTTTCGCGCCCGCCGGCGC‑CTTAAAAAACCC‑ACCCCGTCTGCCGATTAAGCGCctt                                       >  1:477549/1‑149 (MQ=255)                                                                                                        gACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCAt                                    <  1:504891/151‑1 (MQ=255)                                                                                                                       gC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtgt                     <  2:50710/151‑1 (MQ=255)                                                                                                                           ggCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAAc                  <  2:173359/151‑1 (MQ=255)                                                                                                                                 gTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGtt            <  2:687329/151‑1 (MQ=255)                                                                                                                                  tat‑at‑aGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTg           >  1:316520/1‑151 (MQ=255)                                                                                                                                   at‑at‑aGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTg           <  2:225679/150‑1 (MQ=255)                                                                                                                                         aGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCgg       <  2:245309/150‑1 (MQ=255)                                                                                                                                              ttCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCg  >  1:553153/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTAT‑AT‑AGC‑TTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCG‑TTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAG‑CGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCG  >  NC_000913/2252518‑2252805

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