Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,548,360	C→T	D320D (GAC→GAT)	fdnG →	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,548,360	0	C	T	100.0%	89.3 / NA	29	D320D (GAC→GAT)	fdnG	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)

GTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCT  >  NC_000913/1548252‑1548475                                                                                                             |                                                                                                                    gTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaaaa                                                                                       >  1:1199519/1‑139 (MQ=255)         ttttGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATa                                                                               <  2:1199519/139‑1 (MQ=255)                           aaaaCAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                   >  2:98065/1‑133 (MQ=255)                           aaaaCAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                   <  1:98065/133‑1 (MQ=255)                             aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                    <  2:630427/114‑1 (MQ=255)                             aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                    >  1:630427/1‑114 (MQ=255)                                      aTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                   <  1:575411/122‑1 (MQ=255)                                      aTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                   >  2:575411/1‑122 (MQ=255)                                          aCGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGa                                              >  2:545408/1‑139 (MQ=255)                                            gCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                 <  1:749159/102‑1 (MQ=255)                                            gCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATAc                                                                              <  1:564430/105‑1 (MQ=255)                                            gCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATAc                                                                              >  2:564430/1‑105 (MQ=255)                                            gCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaaa                                            <  2:1222674/139‑1 (MQ=255)                                            gCCGAATACGTTAAGCATAACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCaa                                                                                 >  2:749159/1‑102 (MQ=255)                                                            aTTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                            <  2:233029/139‑1 (MQ=255)                                                              tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                          >  1:292570/1‑139 (MQ=255)                                                              tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                          >  2:439820/1‑139 (MQ=255)                                                               acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGa                         >  1:1318114/1‑139 (MQ=255)                                                                    aaCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGaaa                    >  2:118381/1‑139 (MQ=255)                                                                     aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacg                                                                                                  <  1:643361/60‑1 (MQ=255)                                                                     aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacg                                                                                                  >  2:643361/1‑60 (MQ=255)                                                                     aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAac                   >  1:52805/1‑139 (MQ=255)                                                                     aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAac                   >  2:969275/1‑139 (MQ=255)                                                                           gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGAc             >  1:244366/1‑139 (MQ=255)                                                                               gctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTg                                                                <  2:422777/84‑1 (MQ=255)                                                                               gctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTg                                                                >  1:422777/1‑84 (MQ=255)                                                                                ctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATc        <  1:439820/139‑1 (MQ=255)                                                                                    ttcGTGATTATTTTGGTTTCGAAGTTGGTCTGTTCCGCGGCTACGACGGTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCCTCcgcg    <  2:606067/137‑1 (MQ=255)                                                                                      cGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCt  <  2:244366/139‑1 (MQ=255)                                                                                                             |                                                                                                                    GTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCT  >  NC_000913/1548252‑1548475

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