Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,252,657	0	C	T	42.9%	40.5 / 20.1	35	L320F (CTT→TTT)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base T (15/0);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 4.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGC  >  NC_000913/2252512‑2252804                                                                                                                                                    |                                                                                                                                                      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)                                                                                                                                                    |                                                                                                                                                      GCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATA‑GC‑TTTCTTGAAGGATAC‑AGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGC‑GGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGC  >  NC_000913/2252512‑2252804

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