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	96.2%	84.1 / NA	26	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/16);  minor base C (1/0);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	4,293,782	0	C	T	96.2%	82.3 / NA	26	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/15);  minor base G (0/1);  total (10/16)
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.97e-01

GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGG  >  NC_000913/4293542‑4294014                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                         gATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATc                                                                                                                                                                                                                                <  1:54529/251‑1 (MQ=255)           cccAAAATACCGGTGAAAATGCTGACAACACTCACGAAATTGATGGTATTTCTCTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCCGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGCACTATCTGTGGCaca                                                                                                                                                                                                                       >  1:184404/1‑248 (MQ=255)            ccAAAATACGTGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGGGTCG‑GCATTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAAc                                                                                                                                                                                                                     <  2:102923/251‑1 (MQ=255)                                     aacaCTCACGATATTGTTTGTATTTCCCTGTGTCGTGC‑GTTATCTGTTAAGCCCAACATGGCAGGCGTTACGTGCGGAGTAACATCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGt                                                                                                                                                                                            <  2:52522/251‑1 (MQ=255)                                                          aTTTCTCTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTAc                                                                                                                                                                       >  2:198996/1‑251 (MQ=255)                                                             tctcTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtg                                                                                                                                                                    >  1:66126/1‑251 (MQ=255)                                                              ctctGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtgg                                                                                                                                                                   >  2:4500/1‑251 (MQ=255)                                                                      cG‑GCAGTGATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTTGTGAGAAc                                                                                                                                                            <  1:62023/250‑1 (MQ=255)                                                                               aTCTGTGAAGCCCAAAATGGCAGGCGGTACGTGCGGAGGATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCTCGGAAGCGCAGCTGCATGCATTGCAGGATAAAATCCGCGTGAATCCACAAAACAGCGAACAGTGGGCGTTACTGGAGGAGTACTATCTGTGGCAAAACGAGTACGGCAATTCGTTGCTGGCGTCCCGGGAGGCGTTGCAGCTACGTGGTGAGTTCGCTGAACt                                                                                                                                                    <  2:117312/250‑1 (MQ=255)                                                                                      aaGCACAAAATGGCAGGAGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCAGTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTGAGGCGTTGAAACTACGTGGTGAGAACGCTGAACTGTATgcgg                                                                                                                                            <  1:203912/251‑1 (MQ=255)                                                                                                tGGCAGGCGGTACGTGCGGAGTATCAGCGTCCGCGCGATCCGCTACATCAGTTTGCCCGCCAGCACACCCCGGAAGCGCAGCTTCAGCCATTGCAGGATCAACTCCGCGCTAATCCACAAAACAGCGCACAGTGGGCGTTACTGGCTGAATACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGCGAGAACGCTGAACTGAATGCGGCGCTGGCGAc                                                                                                                                  >  1:116428/1‑251 (MQ=255)                                                                                                           aCGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTttatt                                                                                                                       <  1:243097/251‑1 (MQ=255)                                                                                                                  gAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgc                                                                                                                <  1:54313/251‑1 (MQ=255)                                                                                                                        cagcgtcagcgCGATCCGCTACATCAGGTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACTATTACAGCAATTCGTTGCTGGCGTACCTTCAGGCGTTGCAACTACGTGGTGAGAACGCGGAACTGTAGGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagcca                                                                                                          <  1:27891/251‑1 (MQ=255)                                                                                                                           cgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATACGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGAAGGTGCTTTATTACCAGGCcagccagc                                                                                                        >  2:231540/1‑250 (MQ=255)                                                                                                                                                       ccagacagCCAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGTTCTAATCCACAAAACAGCGAACAGGTGGCGTTACTTTATGAGTACTATCCGTGGCAAAACGATTACAGCAATGCGTTGCTGGCGTACCGGCATGCGGTGCAACTACGTGGTGAGAAGGCTGAACTGTATGCGTCGCTGGCGACGTTGCTTTAGTCCTAGGGCAGCGAGAATATGATCGCCCAGACTctgtcaatgatc                                                                           <  2:127126/251‑12 (MQ=255)                                                                                                                                                        cagccagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATGAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCg                                                                          <  1:22064/251‑1 (MQ=255)                                                                                                                                                            cagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACa                                                                       >  1:70881/1‑250 (MQ=255)                                                                                                                                                              gcAAACCACGGAAGCGCAGCTTCCGGCCTTGCAGGCTAAAATCCGCGCTAATCCACCAAACAGCGAACACTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACaa                                                                      >  2:17835/1‑249 (MQ=255)                                                                                                                                                                    cccGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAACGCCCTcg                                                              >  2:144223/1‑251 (MQ=255)                                                                                                                                                                                              aGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGGGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGTCAGTAATGAAATCACCGcc                                    <  1:77047/251‑1 (MQ=255)                                                                                                                                                                                                      aTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgct                            >  2:99346/1‑251 (MQ=255)                                                                                                                                                                                                        ccGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgctg                           <  1:4582/250‑1 (MQ=255)                                                                                                                                                                                                         cgcgCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATgctgct                         <  2:73116/251‑1 (MQ=255)                                                                                                                                                                                                               aaTCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCtt                    <  1:123916/250‑1 (MQ=255)                                                                                                                                                                                                                                aCAGTGGGCGTTACTGGATGAGTACGATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGGGGCGCTGGCGACGGTGCTTTATTGCCAGGCCAGCCAGCATATGACCGGCCAGGCTCGCGCAATGATGGACAAAGCCCGCGCGGTGGACAGTAATGAAAGCACCGCCCTGATGCTGCTGGCTTCCGATGCGTGTTTGCAgg  <  2:87291/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                         GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCG‑GCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGG  >  NC_000913/4293542‑4294014

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