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

GCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGA  >  NC_000913/2252512‑2252806                                                                                                                                                       |                                                                                                                                                     gCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGaa                                                                                                                                                    <  1:443100/151‑1 (MQ=255)   aaCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAgg                                                                                                                                                  <  1:264843/151‑1 (MQ=255)       gtttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGATATA‑GC‑TTTTTTTAAgg                                                                                                                                                  >  2:173037/1‑147 (MQ=255)       gtttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATa                                                                                                                                               <  1:708245/150‑1 (MQ=255)              tATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATA‑GA‑TTTTTTTTAGGAAAC‑AGttt                                                                                                                                        >  2:384814/1‑149 (MQ=255)                      cGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAgc                                                                                                                              <  2:456912/151‑1 (MQ=255)                                   tATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGTGATA‑GC‑TTTTTTTTAGGAGAC‑AGTTTTCgcgcccgcgccccggtttcgg                                                                                                                 >  1:393267/1‑142 (MQ=255)                                        cAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATA‑GT‑TTTTTTTTAGGATAA‑AGATTTCGCGACAACGCCCGCTTTGTgggtgcc                                                                                                            >  2:221930/1‑146 (MQ=255)                                         aGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCt                                                                                                           <  2:673496/151‑1 (MQ=255)                                           cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGGTATA‑TC‑TTTTTTTTAGGGTAA‑AGATTTCGCGCCAGCGCCCCTTTTTCGGttggc                                                                                                            >  1:7622/1‑144 (MQ=255)                                            gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGATATA‑GC‑TTTCTTTTAGGGTAA‑AGTTTCCGCGCCAGCGCCCGCTTTTCgggtgcc                                                                                                            >  2:394583/1‑142 (MQ=255)                                            gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTgc                                                                                                         <  1:496578/150‑1 (MQ=255)                                               aaGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATATA‑AC‑TTTCTTTTAGGATAA‑AGTTTTCGCGACAGCGCCCCCTTTTTGGGAGCCcgcgc‑g                                                                                                     >  2:488715/1‑150 (MQ=255)                                                  gATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCaa                                                                                                 <  2:77499/151‑1 (MQ=255)                                                   aTGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAAt                                                                                                <  2:289279/151‑1 (MQ=255)                                                          aCAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTcgcgc                                                                                         <  1:642664/151‑1 (MQ=255)                                                                 ttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATA‑AC‑TTTTTTTAAGGGTAA‑AGTTTTCCCGCCAGCGCCCGCTTTGCGGTAGGCTGCGC‑GGGAAATTcgcgc                                                                                         >  2:513475/1‑144 (MQ=255)                                                                         cgcgCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGAGATA‑GC‑TTTTTTTAAGGGTAA‑AAGTTTTGCGACCGCGCCCGTTTTTTGGTGGCCTGCGC‑GGCAATTTTTCGCGCCCGCGGCAggct                                                                           >  2:550080/1‑147 (MQ=255)                                                                           cgcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATA‑AA‑TTTTTTTTAGGATAA‑AATTTTCCCGACAGCGCCCGGTTTTCGGGAGGCCGCGCGGGCAATTT‑TCGCGCCCGCGGCCGcctta                                                                         >  2:695538/1‑146 (MQ=255)                                                                             ccAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTATATA‑AC‑TTTTTTTTAGGGTACAAGTTTCC‑CGGCAGCCCCCGTTTTTTGGGAGCCTGCGC‑GGCCAATTTGCGCGCCCGCGGCCGCTTtaa                                                                        >  2:597510/1‑149 (MQ=255)                                                                                  gcATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAaccct                                                                 <  1:706213/151‑1 (MQ=255)                                                                                   cATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATA‑AA‑TTTTTTTTAGGATAC‑AGTTTTCCCGACAACGCCCCGTTTTCGGTAGCCTCCGC‑GGCAAATTTTCGCGCCCGCGGCAGCCTAAACAaccc                                                                  >  2:22743/1‑149 (MQ=255)                                                                                              aGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGGTATA‑GA‑TTTTTTTGAGGGTAA‑AGTTTTCGCGCCAGCGCCCGGTTTTGGGGAGCCCGCGC‑GGCCAATTTGCGCGCCCGCGGCCGCCTTAACAacaccaccccctcttcc                                                     >  1:329926/1‑143 (MQ=255)                                                                                              aGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGACCGGGATATA‑GCTTTTTTTCAGGGAAAA‑AGTTTCC‑CCACAGCGCCCCTTTTTCGGGAGCCCGCGC‑GGCTACTTTTCGCCGCCCCCGCAGCCTAAAAAAACCCACCCTGTc                                                         >  2:348450/1‑147 (MQ=255)                                                                                               gACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTCTATATTC‑TTTTTTTATGG‑TAC‑AGGTTTTCCGCCAGGGCCCGGTTTTCGGGAGCCTGGGC‑GGCAAATTTTTGCGCCCGCCGCCGCCTAaacaaaaccaccc                                                             >  1:321612/1‑134 (MQ=255)                                                                                               gACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGcct                                                    <  2:694516/151‑2 (MQ=255)                                                                                                    ttttGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTTTA‑AC‑TTTTTTTTAGGATAC‑AGTTTTCGCGACAGCGCCCGCTTTTTGGTCGCCTCCGC‑GGCAATTTTTCTCGCCCGCCGCGGCCTCAACAaacctacccctgcttccgcaa                                                 >  1:370764/1‑137 (MQ=255)                                                                                                     tttGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATATA‑GC‑TTTTTTTTAGGAAAC‑AGTTTTCCCGGCAACGCCCGCTTTTTGGGAGGCTGCGC‑GGCCATTTTCCGCGCCCGCCGCAGCTTTAAAAaacccaccccctcctccgctaa                                                >  2:246236/1‑136 (MQ=255)                                                                                                       tGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAAcgc                                             <  2:555056/150‑1 (MQ=255)                                                                                                        ggTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACgcgc                                           <  1:384724/151‑1 (MQ=255)                                                                                                                        ttATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCaa                            <  1:125493/150‑1 (MQ=255)                                                                                                                            ggCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtg                       >  2:363712/1‑151 (MQ=255)                                                                                                                                   cTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAt                >  1:596290/1‑151 (MQ=255)                                                                                                                                   cTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAt                >  1:23374/1‑151 (MQ=255)                                                                                                                                                  ttCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGc    <  2:343239/150‑1 (MQ=255)                                                                                                                                                      tGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGa  <  2:489455/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGA  >  NC_000913/2252512‑2252806

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