Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	62.5%	0.0 / 13.2	17	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/6);  new base A (10/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-03

GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAG  >  NC_000913/2386725‑2386993                                                                                                                                          |                                                                                                                                     gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGCTCCTAAAAAAACGACaaaccaa                                                                                                                           >  1:47739/1‑143 (MQ=255)      aaaCAATACCACTTGCCCGGTGATGCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACattttatt                                                                                                                   <  2:386496/151‑1 (MQ=255)      aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACattttatt                                                                                                                   <  1:378869/151‑1 (MQ=255)         cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATTGCCAattttttttt                                                                                                                 >  2:260042/1‑151 (MQ=255)           aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAAAAACA‑AGAACaatttatttttc                                                                                                               >  1:130390/1‑148 (MQ=255)              ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAACACA‑CGCACCTTTTTTTTTTCCttt                                                                                                           >  1:209290/1‑149 (MQ=255)                          tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCCCATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                               <  2:269666/151‑1 (MQ=255)                              aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGAGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACCGAAAAAAAAAAA‑AGCACAATTTTTTTTTCCTCCTCTGTTTTTAGaagag                                                                                           >  2:373577/1‑147 (MQ=255)                               ccGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTGAGAAAAACAGACA‑TGCACAATTTTTTTTTTCTCATCTGTTTTTTGATACGa                                                                                          >  1:91314/1‑151 (MQ=255)                                                           ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGAAC‑AGCACATTTTTTTTTTTCTTCTCCGTTTTTATATTGGGAAAATGAAAATTAAAGGAGGGAaaaaaa                                                              >  1:16151/1‑151 (MQ=255)                                                                 ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTTTCATAAAAAACGACA‑TGCAACTTTTATTTTTCCTTCTCTGGTTTTAAATTCGAAAAATGAAAATCAAACGCGGTAAAATAAATggtg                                                        >  2:138053/1‑149 (MQ=255)                                                                     gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAA‑AGCACATTTTTTTTTTCCTCATCTGTTTTTAGAAGCGAAAAATGAAAATCAAAGGAGGAAAAATAAATTGGAATGa                                                    >  2:327343/1‑151 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAA‑AAAACAAATTATTTTTTCTCCACCGTTTTTAGAAGAGAAAAAAGAAAATTAAACGATGAAAAAAAAATAGGAAAGagttggtt                                             >  1:46508/1‑144 (MQ=255)                                                                            gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACA‑AGCACAATTTATTTTTTCTCAAATGTTTTTAAATGCGGAAAATGAAATTTAACGTAGGAAAAAAAAATAGGTAAGATTTGCtt                                             >  2:141881/1‑150 (MQ=255)                                                                                          ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt                                <  2:130420/150‑1 (MQ=255)                                                                                                                       cGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg  <  1:194384/151‑1 (MQ=255)                                                                                                                       cGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg  <  1:115456/151‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACA‑TGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAG  >  NC_000913/2386725‑2386993

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