Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,252,646	0	T	G	41.9%	65.0 / 13.3	44	V316G (GTA→GGA)	yeiI	putative kinase
Reads supporting (aligned to +/- strand):  ref base T (6/19);  new base G (18/0);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 2.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATG  >  NC_000913/2252502‑2252793                                                                                                                                                  |                                                                                                                                                       aaGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCATAGtgtat                                                                                                                                                    >  1:182795/1‑149 (MQ=255)  aGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑Gc                                                                                                                                                <  1:474945/151‑1 (MQ=255)      gtTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTc                                                                                                                                            <  1:733482/151‑1 (MQ=255)            cAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCTGTGTATA‑ACGTTCtt                                                                                                                                          >  1:422216/1‑147 (MQ=255)             aaCTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCTGGGGTTTA‑GCtttttttt                                                                                                                                        >  1:686953/1‑146 (MQ=255)                 gtttgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCGGGTGTATA‑GGtttttt                                                                                                                                          >  2:389471/1‑142 (MQ=255)                    tgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTCGGCTATA‑GCtttttt                                                                                                                                          >  1:703535/1‑139 (MQ=255)                    tgtttATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑GAttttttt                                                                                                                                         >  2:203238/1‑139 (MQ=255)                       ttATTTGCCCGATGAGTCAGTTTCTTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGACTGGGATATA‑TATTTCTTTTAGGATACAGatt                                                                                                                            >  2:280518/1‑150 (MQ=255)                        tATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTc                                                                                                                           <  1:368779/150‑1 (MQ=255)                                 gATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑ACTTTTTTTTAAGGGACAGTTTTCCCGAc                                                                                                                     >  1:819462/1‑147 (MQ=255)                                      gTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑ACTTTTTTTTAGGATACAGTTTCCGCGACAGCGccc                                                                                                              >  1:703056/1‑149 (MQ=255)                                        cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGGCGGGGAGATA‑AATTTTTTTGAGGAAAAAGTTTTTCCGACAACCCCCCttt                                                                                                          >  2:62933/1‑151 (MQ=255)                                        cAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGGGATA‑TCTTTCTTTTAGGGTACAGTTTTCGCGAGAGCGCCCGgtt                                                                                                          >  1:402141/1‑151 (MQ=255)                                         aGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTTTTGTAGGATACAGTTTTCGCGACAGCGCCCGtttt                                                                                                         >  1:856373/1‑151 (MQ=255)                                              aTTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGt                                                                                                    <  2:23134/151‑1 (MQ=255)                                                tGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTa                                                                                                   <  1:479203/150‑1 (MQ=255)                                                   aGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGATATA‑AATTTTTTTAAGGGTAAAGTTTTCGCGACAGCGCCCGCTTTTCGGTTGcc                                                                                                >  1:777655/1‑150 (MQ=255)                                                     cGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTgc                                                                                             <  2:664814/151‑1 (MQ=255)                                                      gAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATTTA‑GATTTTTTTAAGGATATAATTTTCGCGGCAGCGcccccctttcgggagcc                                                                                                >  2:224633/1‑133 (MQ=255)                                                          aGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGc                                                                                       <  1:61766/151‑1 (MQ=255)                                                             aTGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAAt                                                                                    <  1:767846/151‑1 (MQ=255)                                                               ggAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAAttt                                                                                  <  2:288782/151‑1 (MQ=255)                                                                gagaACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGCGGGGCCTGG‑ATATA‑GGTTTTTTTAAGGATACAATTTTCGCGGCAGCGCCCGTTTTTTGGTAGCCCGCGGGGGCAAttt                                                                                  >  1:581752/1‑151 (MQ=255)                                                                      aaTTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGATATA‑GGTTTTTTTTAGGATACAGGTTCCGCGAGAGAGCGCCTTTTTCGGGGGCCTGC‑GCGGCGAGTTcctgcg                                                                            >  2:133655/1‑145 (MQ=255)                                                                          tttGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAg                                                                        <  2:660393/150‑1 (MQ=255)                                                                              cTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGCCCTCGGATATA‑GCTTTCTTTGATGATACAGGTTTCGCGCCCGCGCCCCCTTTTCGGGAGGCTGC‑GCGGCCAATTcgtgccc                                                                           >  2:349722/1‑141 (MQ=255)                                                                               tGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGTTATT‑ACTTTTTTTTAGGAGACAATTTCCGCGACAGCGCCCCGTTTTCGGGAGCCTGC‑GCGGCAAATTCTCCCGcc                                                                          >  1:524394/1‑143 (MQ=255)                                                                                gACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAg                                                                  <  2:824849/150‑1 (MQ=255)                                                                                  ccGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑Ctt                                                              <  2:163330/151‑1 (MQ=255)                                                                                             cATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑AATTTTTTTGAGGATACAGTTTTCGCGGCAGCGCCCGGTTTTCGGTTGCCTGG‑GGGGCACTTTTTCGCGCCCGCGGGCG‑CCTTAACAaccc                                                     >  2:300681/1‑149 (MQ=255)                                                                                                aCGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGTTATATGATTTTTTG‑AGGGTAAAGTTTTCCCGCCAGCGCCCGGTTTTTGGTAGGCCGG‑GCGGCCATTTTGCCCGGCCGCGGCAGCCTTAAAAAacccaccc                                                 >  1:642727/1‑147 (MQ=255)                                                                                                  gACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑GCTTTCTTTAAGGGTAAAGTTTCCGCGGCCGCGCCCGTGTTTCGGGAGGCCGC‑CCGGCAAATTCCCTCGGCCGCGGCCG‑CCTAaaaaacacccccc                                                >  2:248000/1‑142 (MQ=255)                                                                                                     aGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGGGTATA‑GCATTTTTTTAGGATAAAGTTTTCGCGGCAGCGGCCGCTTTGCGGTAGCCCGC‑GGGGCGAATTTTCGCCCCCGCGGCAG‑CTTAAAAAaccctcccc                                                >  1:239317/1‑146 (MQ=255)                                                                                                         gACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCCGGGATATA‑GCTTTTTTTTAGGGTAAAAGTTTCGCGACAGCGCCCGGTTTTCGGGAGCCTGC‑GCGGCAATTTTGCGCGCCCGCGGCAG‑CTTAAAAAACCCCACCCCCTCTGcc                                        >  1:651482/1‑150 (MQ=255)                                                                                                             gTTTTGGTGCTGACGATGGTTTTACGGC‑GGTCCTCGTATATA‑GCTTTCTTGATGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGCTAGCCTGC‑GCGGCAATTTCACGCGCTAGCGTCAG‑CCTAAACAACGCTACCCTGTCTGCCGATaa                                   <  2:714064/151‑1 (MQ=255)                                                                                                                  ggTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAAcgcg                               <  1:604647/150‑1 (MQ=255)                                                                                                                  ggTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACgcgc                              <  1:203296/151‑1 (MQ=255)                                                                                                                        gACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCAt                        <  2:402880/151‑1 (MQ=255)                                                                                                                         aCGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCAt                        <  2:422083/150‑1 (MQ=255)                                                                                                                                      ggC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGgtg          >  1:519158/1‑151 (MQ=255)                                                                                                                                        c‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTa        <  1:816991/151‑1 (MQ=255)                                                                                                                                            gCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAAc      >  1:649974/1‑150 (MQ=255)                                                                                                                                             ccTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACa     >  2:357267/1‑150 (MQ=255)                                                                                                                                                ggTATATA‑GCTTGCTTGAAGGATACAGTTTCCGGGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATg  <  2:649426/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       AAGGTGTTCAGCAACTGTTTGTTTATTTGCCCGATGAGTCAGTTTATTGCAGCGAAAAGGATGGAGAACAATTTTTGCTGACCGCGCCAGCGCATACGACAGTAGACAGTTTTGGTGCTGACGATGGTTTTATGGC‑GGGCCTGGTATATA‑GCTTTCTTGAAGGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGC‑GCGGCAATTTCGCGCGCCAGCGGCAG‑CTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATG  >  NC_000913/2252502‑2252793

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