Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,293,781	2 bp→AT	coding (239‑240/597 nt)	nrfG →	putative formate‑dependent nitrite reductase complex subunit NrfG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,293,781	0	G	A	85.0%	45.2 / 0.8	20	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (8/9);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
*	NC_000913	4,293,782	0	C	T	85.0%	46.5 / 2.2	20	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (8/9);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACGC  >  NC_000913/4293544‑4294024                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                   tgaaacgacacAAAATACCGGTGAAAATGCTGAAAACGCGCACGATCTTGATGGTATTTCTGTGTGTCGGAAGTGAGCTGTTAAGTCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGGCAGCGCGATCCGCTACATCAGTTGGCCAGTCAGCAAACCACGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCTTTACTGGATGAGTACTATCt                                                                                                                                                                                                                                         <  2:2559/240‑1 (MQ=255) ggAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCtg                                                                                                                                                                                                                                        <  2:110217/250‑1 (MQ=255)  gAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCtg                                                                                                                                                                                                                                        <  1:342298/250‑1 (MQ=255)                   cGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCAACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGc                                                                                                                                                                                                                      >  1:11084/1‑251 (MQ=255)                                  caacaCTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGTCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGAGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCg                                                                                                                                                                                                       >  1:52585/1‑251 (MQ=255)                                      acTCACGATATGGATGGTATTTCGCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTAcc                                                                                                                                                                                                   <  2:22401/251‑1 (MQ=255)                                                        ttttCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGCCAGCGCGATCCGGTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTa                                                                                                                                                                                  <  2:95406/249‑1 (MQ=255)                                                                      gCAGTTATCCGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAATCGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACTTGGTTAGAACGCt                                                                                                                                                                   <  1:285602/251‑1 (MQ=255)                                                                                                      gTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTtta                                                                                                                                   <  1:51786/251‑1 (MQ=255)                                                                                                        aCGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTTCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTTGCGACGGTGCTttat                                                                                                                                  <  1:158882/250‑1 (MQ=255)                                                                                                        aCGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                                 >  1:171359/1‑251 (MQ=255)                                                                                                                aGTAGCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGTTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAggtc                                                                                                                         <  2:118170/251‑3 (MQ=255)                                                                                                                 gTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgcc                                                                                                                         >  1:34328/1‑250 (MQ=255)                                                                                                                   aTCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagc                                                                                                                      >  1:45097/1‑251 (MQ=255)                                                                                                                     aagcgtaagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTACAACTACGTGGTGAGAACGCTGCACTGGATGCGGCGCTGGCGACGGGGATTTATTACCAGgccagac                                                                                                                     >  2:260171/2‑248 (MQ=255)                                                                                                                                                  tGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCCGCATATGACCGCCCAGACTCGCGCAATg                                                                                        <  1:248248/250‑1 (MQ=255)                                                                                                                                                                ccccGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGCCAAAGCCCTc                                                                         >  2:8991/1‑251 (MQ=255)                                                                                                                                                                      aaGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCGGGCGACCGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCtaccgc                                                                     >  2:52914/1‑244 (MQ=255)                                                                                                                                                                                                                       cAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAACACGATTACAGCACTTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGttt                   >  1:19175/1‑250 (MQ=255)                                                                                                                                                                                                                                       ttACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAAGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTAcgc  >  2:201698/1‑251 (MQ=255)                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                   TGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACGC  >  NC_000913/4293544‑4294024

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