Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,037,140	0	A	G	64.3%	3.7 / 11.3	14	D131G (GAC→GGC)	hyaF	hydrogenase‑1 protein nickel incorporation factor
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (9/0);  total (11/3)
Fisher's exact test for biased strand distribution p-value = 2.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03

ATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAGTTACTGGCACATGTACGTAACCC  >  NC_000913/1037004‑1037290                                                                                                                                         |                                                                                                                                                       aTCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCCGACGGCGCGGCGGAAAAATGTTGGGGGGCAAAAaggaggggg                                                                                                                                          >  1:114087/1‑148 (MQ=255)                 gggaaggggaagTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCCGACGGCGCGGCGGAAAAATGTTGGGGGGGAAACTGGAGGGGGGGGGCGcgccccggg                                                                                                                           >  1:135722/1‑145 (MQ=255)                       gggaagTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGGGGACGGCCCGGCGGGAAAGTGCTGGGGGGAAAACCGGGGGGTGGCCGGGCGCCCCCGGCGGtgtg                                                                                                                    >  2:241453/1‑151 (MQ=255)                            gtgtCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCGGACGGCGCGGGGGAAAAATGGTGGGGGGAAAAATGGAGGGGGGGCGCGCGCCGCCGGCGGTgggggggg                                                                                                               >  1:166696/1‑147 (MQ=255)                            gtgtCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCGGAAGGCGCGGGGGAAAAATGCTGGGGGGGAAAATGGAGGGTGGGTGCGCGCCGCCGGCGGTGTgggggg                                                                                                               >  1:96999/1‑147 (MQ=255)                                cAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCCGAAGGCGCGGCGGAAAAATGCTGGGGGGCAAAAAGGAGGGTGGCTGCGCGCCCCCGGCGGTGTGGGAggcggg                                                                                                           >  1:223492/1‑150 (MQ=255)                                         ttCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGGCGCGGGGGAAAAGTGGTGGGGGGAAAACTGGAGGGTGGCTGGGCGCCGCCGGGGGTGTGGGAggggggcaagcacaa                                                                                                  >  1:115881/1‑141 (MQ=255)                                                                   gaaaTACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGGCGAAGGCGCGGCGGAAAAATGCTGGGGGGGAAAATGGGGGGTGGCTGCGCGCCGCCGGCGGTGTGGGAGGGGGGGACGCAAAAAATCTTTCCGCCAGATTTGCttttt                                                                         >  1:23740/1‑149 (MQ=255)                                                                   gaaaTACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGGGCCGACGGCGCGGCGGAAAGATGCTGGGGGGCAAAATGGAGGCTGGCCGCGCGCCGCCGGCGGTGGGGGGGGGGGGCACGCCAAAAATCTTTCCCGCAGATTCTCtttttt                                                                        >  1:261120/1‑149 (MQ=255)                                                                                                aTGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATgg                                            >  1:171393/1‑150 (MQ=255)                                                                                                aTGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATgg                                            >  1:200275/1‑150 (MQ=255)                                                                                                               aCGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGgcg                               <  2:217092/148‑1 (MQ=255)                                                                                                                  tCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAg                         <  1:241455/151‑1 (MQ=255)                                                                                                                                tGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAGTTACTGGCACATGt           <  1:291701/151‑1 (MQ=255)                                                                                                                                        gaCAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAGTTACTGGCACATGTACGTAAc    <  1:22414/150‑1 (MQ=255)                                                                                                                                         aCAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAGTTACTGGCACATGTACGTAAccc  <  2:145127/151‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                       ATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCCTACCGTTGGCGCATGAGTTACTGGCACATGTACGTAACCC  >  NC_000913/1037004‑1037290

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