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.2%	57.1 / 21.4	45	L320F (CTT→TTT)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base C (6/20);  new base T (19/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 7.46e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTT  >  NC_000913/2252514‑2252795                                                                                                                                                |                                                                                                                                               aaCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAgg                                                                                                                                         <  2:569320/151‑1 (MQ=255)         gtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑Gtt                                                                                                                                <  2:674818/151‑1 (MQ=255)         gtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑Gt                                                                                                                                 <  2:997777/150‑1 (MQ=255)          tttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATATATATTTTTTTTATGAGACA‑Gttt                                                                                                                               >  1:494408/1‑151 (MQ=255)             aTTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTcc                                                                                                                             <  1:101008/150‑1 (MQ=255)               ttGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCcgc                                                                                                                           <  2:927891/150‑1 (MQ=255)                tGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCcgcg                                                                                                                          <  1:526575/150‑1 (MQ=255)                  cccGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATAGCTTTTTTTTAGGATACA‑ATTTTcgc                                                                                                                           >  2:286835/1‑147 (MQ=255)                   ccGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAg                                                                                                                      <  1:564811/151‑1 (MQ=255)                     gATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAg‑cg                                                                                                                   <  2:422057/151‑1 (MQ=255)                        gAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGTCCCGGGATATAGCTTTCTTTTAGGATAAA‑Attttctcggcacacccc                                                                                                                 >  2:988831/1‑137 (MQ=255)                         aGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATTTATATTTTTTTTAGGAGAAA‑GTTTTcgcg                                                                                                                          >  1:936035/1‑141 (MQ=255)                            cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGTATATAGCTTTTTTTTAGGATAAA‑GTTTTCGCGGCAG‑CG‑ccccgtt                                                                                                           >  2:906002/1‑147 (MQ=255)                               tttATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTCGGATATAACTTTTTTTTTGGGTATA‑ATTTTCGCGACAG‑CG‑CCCGcttt                                                                                                          >  2:555361/1‑149 (MQ=255)                                ttATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGAGTTAGATTTTTTGTAGGATACA‑GATTTCCCGCCAA‑CG‑CCCGcttttt                                                                                                        >  1:591356/1‑148 (MQ=255)                                  aTTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCCGTATATAACTTTCTTTGAGGGTAAA‑GACTCTACGCCCG‑GG‑CCCCTTTTGCggg                                                                                                     >  1:346804/1‑150 (MQ=255)                                   ttGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATAGATTTTTTTAAGGAGAAA‑GTTTTTCCGGCAG‑CG‑ccccctttttg                                                                                                       >  1:978597/1‑140 (MQ=255)                                   ttGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGGTATAGCTTTTTTGTAGGATACA‑GTGTCCGCGACAG‑CG‑CCCGCTTTGCgggg                                                                                                    >  1:63687/1‑149 (MQ=255)                                     gCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑g                                                                                                  <  2:440908/150‑1 (MQ=255)                                      cAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATATAAGTTTTTTTTAGGAAACA‑ATTTTCGCGACAGCCC‑CCCGttttctg                                                                                                       >  1:565014/1‑143 (MQ=255)                                         cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTGCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTgc                                                                                             <  2:65293/151‑1 (MQ=255)                                                gATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCaa                                                                                      <  2:564106/151‑1 (MQ=255)                                                    gagaACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGTATATAGATTTCTTTGCTGAAAAA‑ATTTCCCCTACCA‑CG‑CCCCTTTTGCGGTT‑GGCCCCGCGGCCAttt                                                                                   >  1:429879/1‑150 (MQ=255)                                                    gagaACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATATAGGTTTTTTTTAGGGAAAA‑GGTTTCGCGACAG‑CG‑CCCGGTTTTGGGTT‑GCCCGCGCGGGAAttt                                                                                   >  1:695950/1‑150 (MQ=255)                                                     agaACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGATATAACTTTTTTTTAGGGTACA‑GTTTTCGCGACAG‑CG‑CCCCGTTTTCGGGA‑GCCTGCGCGGCAAttt                                                                                   >  2:262448/1‑149 (MQ=255)                                                        aCAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTcgcgc                                                                              <  2:832207/151‑1 (MQ=255)                                                          aaTTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGGTATAACTTTCTTTTAGGAAACA‑GGTTCCCCCACAG‑CG‑CCCGCTTTTCGGTA‑GCCTGCGCGGCAAttttcc                                                                                >  2:651766/1‑144 (MQ=255)                                                           aTTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGGTATAACTTTTTTGTAGGAGAAA‑GATTCCGCGACAG‑CG‑CCCCGTTTTTGGGA‑GGCTCCGCGGGAAatttccggccc                                                                           >  1:476828/1‑143 (MQ=255)                                                           aTTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGATTTAAATTTTTTTAAGGATACAGGTTTCCCCG‑CAG‑CGCCCCCTTTTTGGTTA‑GCC‑GCGCGGCCAtttttcccgcc                                                                           >  2:1032111/1‑143 (MQ=255)                                                               ttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGGATATAGCTTTTTTTAAGGATACA‑ATTTTCCCGCCAG‑CG‑CCCCGTTTTCGGTA‑GCCCGCGCGGGAATTTCGCGCCCCAGc                                                                        >  2:961196/1‑150 (MQ=255)                                                                tGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCGGGGTATAGCTTTTTTTTAGGAGACA‑GTTTTTCCGGAAG‑CC‑CCCGCTTTTCGGGG‑TCCCGCGCGGCGACTTCGCccgcccgc                                                                        >  1:892598/1‑149 (MQ=255)                                                                   tGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCCCGGATATAACTTTTTTTTAGGGGACA‑GTTTTCGCGGCAA‑CG‑CCCCGTTTTCTGGA‑GCCCGCGCGGGAAATTTGCGCGccccggg                                                                      >  1:70875/1‑143 (MQ=255)                                                                          gcCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAaac                                                            <  1:81686/151‑1 (MQ=255)                                                                             aGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATACGTTTTTTTGGGGATAGC‑GATTTCCCGAAAG‑CG‑ACCGCTTTGCGGGT‑GACTGCGCGGGAAGTTCTCGGGCCCGCCGCGGcgt                                                                >  2:551195/1‑142 (MQ=255)                                                                                gcATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAaccct                                                      <  2:408781/151‑1 (MQ=255)                                                                                  aTACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCTGGGATATAGCTTTTTTGTAGGGTACA‑GATTTCGCGGCAG‑CG‑CCCCTTTTTTGGGAGGCCCGGCGGGCAATTT‑TCGCGCCCGCGGCAGCTTaaaaaa                                                          >  1:711407/1‑145 (MQ=255)                                                                                      gACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCGGGCTATTGGTTTTTTTTAGGGTACA‑AATTTTTGGCCAG‑AG‑CCCCGTTTTTGGTA‑GCCTGCGCGGCAATATTACCCGCCAGCGCCCGCTTTAACAaaccccccccc                                                >  1:165910/1‑144 (MQ=255)                                                                                        cAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTc                                              <  2:306966/151‑1 (MQ=255)                                                                                          gTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTg                                            <  2:163632/151‑1 (MQ=255)                                                                                                         gCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCttt                             <  1:430720/151‑1 (MQ=255)                                                                                                             aCGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCAtt                         <  2:475590/151‑1 (MQ=255)                                                                                                              cGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTa                        <  2:815970/151‑1 (MQ=255)                                                                                                               gATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAg                       <  1:626902/151‑1 (MQ=255)                                                                                                                            cGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtgt           >  1:63719/1‑150 (MQ=255)                                                                                                                                    gTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGtt  >  2:104764/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                               AACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGCCTGGTATATAGCTTTCTTGAAGGATACA‑GTTTCCGCGACAG‑CG‑CCCGTTTTGCGGTA‑GCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTT  >  NC_000913/2252514‑2252795

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