Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	65.0%	3.9 / 13.9	20	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base A (13/0);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-04

GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTG  >  NC_000913/2386725‑2386991                                                                                                                                          |                                                                                                                                    gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGGAAAAAAGACATGCAAatt                                                                                                                        >  2:76502/1‑151 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt                                                                                                                   <  1:185610/151‑1 (MQ=255)           aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAAAGAACaatttatttttc                                                                                                               >  2:167405/1‑148 (MQ=255)                aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                                         <  2:18972/151‑1 (MQ=255)                     cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAAGAAAAGCACATTTTATTTTTCCTCCTCTtttt                                                                                                    >  1:256965/1‑151 (MQ=255)                       cGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACATTTTATTTTTCCTCATCTGttttt                                                                                                  >  1:228200/1‑150 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                               <  1:256768/151‑1 (MQ=255)                              aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTATTTTTCCTCATCTGTTTTTAAA‑TGCg                                                                                          >  1:260787/1‑151 (MQ=255)                                                                  gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGGAAAAAAGACACAGAAATTTTTTTTTTTCTCATCCGGTTTTAGATTGGGAAAAAGAAAAACAAAGGGG‑GGAAAAAAAAATGGGa                                                      >  2:249726/1‑151 (MQ=255)                                                                      tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAt                                                 <  1:55688/151‑1 (MQ=255)                                                                         aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAAA‑TGCGGAACAAGAAATTTAACGGA‑GGAAAAAAAAAATGGGATGATTgg                                              >  1:134218/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTTTTTTTAAA‑AGCGGAAAAAGAAAATTAAAGGG‑GGAAAAATAAAATGGGATGATTTGCtt                                           >  2:215030/1‑150 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAA‑TGCGGAACATTAAATTTAAAGGA‑GGGAAAATAAATGGGAATGATTTGCtt                                           >  1:142978/1‑150 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGAACATTTTTTTTTTCCTCATCTGTTTTTAGA‑AGCGAAACATGAAAATCAACGGA‑GGGAAAATAAATGGGGATGATTGGCTa                                           >  2:267263/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGA‑TGCGAAAAAAGAAAATTAACGGG‑GGGAAAAAAAATTGGGATGATTTGCTa                                           >  1:75668/1‑151 (MQ=255)                                                                                     cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCCAAAAAAAAGCATGCACCTTTTATTTTTCCTCTTCTTTTTTTAGA‑TTTGGAAAATACAATTAAAAGGAGGGAAAAAAAATTGGGAAAAATTGG‑TAtttttgaca                                  >  2:232280/1‑146 (MQ=255)                                                                                            ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAACGACAAGCACATTTTTTTTTTTCTCATCTGTTTTTAGA‑AGCGAAAAATGAAAATTAAAGGA‑GGAAAAAAAAAAGGGGATGATTGGGTATTTTGAAGAAAATGaa                           >  1:125640/1‑150 (MQ=255)                                                                                                             aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa          <  1:208349/151‑1 (MQ=255)                                                                                                                     ctcGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg  <  1:311738/151‑1 (MQ=255)                                                                                                                     ctcGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg  <  1:168706/151‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGA‑TGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTG  >  NC_000913/2386725‑2386991

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