Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,686	0	G	T	39.6%	42.9 / 14.3	53	intergenic (+40/+25)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (21/11);  new base T (6/15);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCT  >  NC_000913/2151542‑2151831                                                                                                                                                     |                                                                                                                                                  tataCGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTATTCTTTTTTTGTTACCTATCCAAAAAAAAAAACCGACGAGCTAAGAGCAAAAATAATTTG‑CTAATCAGCCGGTTTA‑TTTTTTTTTTTAAAT‑TGC‑GGGTGGAttttttttt                                                                                                                                             >  1:1332372/1‑147 (MQ=255)  ataCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCt                                                                                                                                             <  1:434750/150‑1 (MQ=255)    aCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCa                                                                                                                                          <  2:1009106/151‑1 (MQ=255)     cgcgcgGGGATAAAATTTACGGCGGGTTTACCATTTATCCCTTTTTTGTCACCTATACAAAAAACAAAGCAGACGCGCTAAGAGCAAAAATAATTCT‑CTCATCAACCGCTATATTTTTTTTTTTAAACA‑GGG‑GGTTGAAtttttttttt                                                                                                                                            >  2:90154/1‑144 (MQ=255)      gcgcgGGGATAAAATTTACGGTGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCaa                                                                                                                                         >  2:1030212/1‑150 (MQ=255)           gggATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCa                                                                                                                                   <  1:158624/151‑1 (MQ=255)                                 taCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCcac                                                                                                              >  1:707586/1‑150 (MQ=255)                                    cATTGATCTTTTATTGGCCCCCTACCAAAAAAAAGAAGCCGACGAGAAAAGAAAAAAATTAGTCCG‑CTAACCAGCCTTTATA‑TTTATTTTTAAAACG‑GGC‑GGTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCcaca                                                                                                             <  1:334207/148‑1 (MQ=255)                                       tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc                                                                                                        >  1:1312786/1‑150 (MQ=255)                                           ccTTTATGGGCCCCCTATCAAAAAAAAGAACCCGACGAGAAAGGAACAAAATTATTTCG‑CAAATCGGCCTTTTTT‑TTTATTTTTGAAGCG‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                       <  1:1164015/147‑1 (MQ=255)                                             tttATTGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCTTTATA‑TTTATTTTTTAAACG‑GGC‑GTTTGAATTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                       <  1:1034720/145‑1 (MQ=255)                                             tttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTTAACCCCCCTGAGAAAATCACTCTTCCCACAGCGCCAGGGcc                                                                                                  >  1:306474/1‑150 (MQ=255)                                                aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTTAAACGCCATGAGAAAATAACTCTGCCCACCGCGCCCGTGCCAACg                                                                                              >  1:198859/1‑151 (MQ=255)                                                  tgtcaacttATCAAAAAAAAAAAGCCGACGAAAAAAGAGAAAAATTTTTCCC‑CTAATCGGCCTTTATT‑TTTTTTTTTTAAACG‑GGC‑GTTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                                                                                                <  1:734678/139‑1 (MQ=255)                                                    gcccccgaccaaaaaGAAGAGCCCGAGGAGCAAAAAGAAAAATTGGTTCG‑AAAATCGGCTTTTATT‑TTTTTTGTTAAACCT‑GGC‑GGTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtttt                                                                                          <  1:1047700/140‑1 (MQ=255)                                                     tCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAAATGTTGTCTTCAACAGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGtttt                                                                                          >  2:901693/1‑150 (MQ=255)                                                      ccccTATCAAAAAAAAGAAGCCGACGAGTAAAAAGAAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTTTTTTTTAAACT‑GGC‑GGTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                                                                        <  1:864875/149‑1 (MQ=255)                                                              cAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCCACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATcc                                                                                >  2:892224/1‑151 (MQ=255)                                                                  ataagaagCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTa                                                                             >  1:818013/3‑150 (MQ=255)                                                                    aagaagCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGGATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGCGCCAACGTTTTATCATCATCCCTAc                                                                            >  1:319167/1‑149 (MQ=255)                                                                     agaagCCGACGAGCAAAGAAAAAAATTATTTCG‑CTAATCAGCCTTTATA‑TTTATTTTTTAAACG‑GGC‑GGTTGATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCt                                                                         <  2:1013333/151‑1 (MQ=255)                                                                             aaaaagtaaaaataaaaATAATTCGCCTAATCAGCTATTATT‑TTTATTATTAAATCT‑GGC‑GGTTAATTTTTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTa                                                                   <  2:927046/137‑1 (MQ=255)                                                                                 gaaaaGAGAAAATTTATTTCG‑CAAATCAGCCTTTATT‑TTTTTTTTTTAACCG‑GGG‑GTTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTa                                                                   <  1:666474/142‑1 (MQ=255)                                                                                   tAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAAAGTTTTTTCATCATCCGTACGCTCGTTTACCGCCgg                                                            >  2:235142/1‑150 (MQ=255)                                                                                   tAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCggg                                                           >  1:907747/1‑151 (MQ=255)                                                                                     agagCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGGTTACCGCCGGGc                                                          >  1:799054/1‑150 (MQ=255)                                                                                       agaaaaaTTAGTTCG‑CTAATCAGCCTTTATT‑TTTATTTTTTAACCG‑GGC‑GTTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTg                                                        <  1:117639/147‑1 (MQ=255)                                                                                         aaaaaTTATTTCT‑CAAATCAGCCTTTATT‑TTTTTTTTTTAAACG‑GGC‑GGTTGATTTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTgg                                                       <  2:264105/148‑1 (MQ=255)                                                                                          aaaaTTATTGGG‑TTAATCAGCCTTTTTT‑TTAATTTTTTAACCT‑GGC‑GTTTGAATTTTTTTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGaaa                                                    <  1:149297/151‑3 (MQ=255)                                                                                          aaaaTTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGGTGGAc                                                     >  2:694611/1‑150 (MQ=255)                                                                                           aaaTTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCTTCCGTACGCTCGTTTACCGCCGGGCTGGAc                                                     >  1:797039/1‑149 (MQ=255)                                                                                            aaattattcctttAATCAGCCTTTATT‑TTTTTTTTTAAAGCT‑GGC‑GTTTAATTTGTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaa                                                   <  2:1331625/139‑1 (MQ=255)                                                                                            aaTTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTTAATTTTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTACCGCCGGGCTGGACaa                                                   >  2:1295967/1‑150 (MQ=255)                                                                                              ttAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGTATTTTTTTCTTCAAACCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGGCAAAGTTTTATCATTATCCGCACGCCCGTTTACCGCCCGGGTGGaaaaaa                                                 >  2:169843/1‑150 (MQ=255)                                                                                                  ttCG‑CTAATCAGCCTTTATA‑TTTTTTTTTTAACCG‑GGG‑GGTTAAATTGTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaaa                                                <  1:1133123/147‑1 (MQ=255)                                                                                                       cTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAAATGCttt                                         >  2:746775/1‑150 (MQ=255)                                                                                                         aaTCAGCCTTTATT‑TTTTTTTTTTAACCG‑GGC‑GGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCttt                                         <  1:100370/148‑1 (MQ=255)                                                                                                              gCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTTACTCTGCCCCCAGCGCCAGTGCCCAAGTTTTATTATCATCCCGACCCTCGTTTAACGCCCGGCTGGACAAAAAATGCTTTAACAAtt                                  >  2:54601/1‑150 (MQ=255)                                                                                                               ccGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCACTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCCGGACAAAAACTGCTTTAACAAtt                                  >  1:462822/1‑149 (MQ=255)                                                                                                                 gTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTTTTGTCTTTAACCCCCCTGAGAAACTTACTCTGCCCACAGCGCCCGTGCCAACGTTTTTTCATCATCCGGACGCTCGGTTAACGCCGGGGTGGACAAAAAATGCTTTAAAAATTCAgg                              >  1:666812/1‑151 (MQ=255)                                                                                                                  ttATA‑TTTTTTGTTAAAATG‑GGC‑GTTTGAATTGTTGTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                              <  1:753519/150‑1 (MQ=255)                                                                                                                      a‑tttattGTTTAAACT‑GGC‑GGTTGAATTTTTGTCTTCAACCGCCATTAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATTATCATCCGTACGCTCGTTTACCGCCGGGCGGGAAAAAAACTTCTTTAACAATTCAGGgaa                           >  2:496977/1‑149 (MQ=255)                                                                                                                         ttttttgttAAACTGGGC‑GTTTGATTTTTTTTTTACACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaata                         <  2:136386/143‑1 (MQ=255)                                                                                                                         ttattGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaataa                        >  1:761430/1‑150 (MQ=255)                                                                                                                          ttttttttAACCT‑GGC‑GTTAAAATGTTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATaaa                       <  2:629439/145‑1 (MQ=255)                                                                                                                          tattGTTTAAACT‑GGC‑GGTTGAATTGTTTTCTTCAACCGCCCTGAGAAACTCACTCTGGCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATaaa                       >  2:482430/1‑150 (MQ=255)                                                                                                                            ttGTTAAAACG‑GGG‑GGTTAATTTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                      <  2:1315640/149‑1 (MQ=255)                                                                                                                              tttaaaaCC‑GGCGGGTTGAATTTTTCTCTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCg                    <  2:1130506/146‑1 (MQ=255)                                                                                                                                 tAAACT‑GGC‑GGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAggt                >  2:838570/1‑148 (MQ=255)                                                                                                                                  aaaCG‑GGC‑GTTTAATTTTTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttg               <  1:345484/150‑1 (MQ=255)                                                                                                                                  aaaCG‑GGC‑GGTTGATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttg               <  1:1276866/150‑1 (MQ=255)                                                                                                                                     agcggggggTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt             <  1:237946/144‑1 (MQ=255)                                                                                                                                      t‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCt          >  1:600801/1‑151 (MQ=255)                                                                                                                                                taaTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct  >  1:832127/2‑151 (MQ=255)                                                                                                                                                gAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct  >  1:452601/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCG‑CTAATCAGCCGTTATA‑TTTATTGTTTAAACT‑GGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCT  >  NC_000913/2151542‑2151831

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