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%	101.4 / NA	32	D320D (GAC→GAT)	fdnG	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/12);  total (20/12)

ATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACC  >  NC_000913/1548241‑1548487                                                                                                                        |                                                                                                                                aTTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTacg                                                                                                              >  1:226600/1‑139 (MQ=255)           tGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGaaaa                                                                                                    >  1:27976/1‑139 (MQ=255)               tGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                >  1:175922/1‑139 (MQ=255)               tGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATT‑‑‑‑‑‑ACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCa                                                                                                                       <  2:251145/110‑1 (MQ=255)               tGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATT‑‑‑‑‑‑ACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCa                                                                                                                       >  1:251145/1‑110 (MQ=255)                 gCGTTTTGCGCTACCTGATCGAAAACAACGAAATCAACGCCGAATTCGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCa                                                                                              >  2:642418/1‑139 (MQ=255)                       tGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGa                                                                                        >  1:615366/1‑139 (MQ=255)                           cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                    >  1:897364/1‑139 (MQ=255)                                tGATCGAAAACTACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGATGAAAAACGTCAATACGATAAATCGTc                                                                               <  1:720064/139‑1 (MQ=255)                                tGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                               <  1:642418/139‑1 (MQ=255)                                        aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACt                                                                       >  2:19268/1‑139 (MQ=255)                                            acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTac                                                                                                               >  2:822127/1‑95 (MQ=255)                                            acaaAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTac                                                                                                               <  1:822127/95‑1 (MQ=255)                                                aaTCAACGCCGAATACGTTAAGCATTACACCAACGTCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTc                                                               >  1:222026/1‑139 (MQ=255)                                                aaTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTc                                                               >  1:1089534/1‑139 (MQ=255)                                                                 ttAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGGAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaa                                                         <  2:797107/128‑1 (MQ=255)                                                                 ttAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaa                                                         >  1:797107/1‑128 (MQ=255)                                                                   aaGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGa                                            >  2:689254/1‑139 (MQ=255)                                                                       aTTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                                        >  2:1090936/1‑139 (MQ=255)                                                                         tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                                      >  1:680755/1‑139 (MQ=255)                                                                         tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                                      >  2:853920/1‑139 (MQ=255)                                                                         tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcgcg                                      >  2:299020/1‑139 (MQ=255)                                                                          acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCTGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                                        <  2:796265/136‑1 (MQ=255)                                                                          acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAAcg                                        >  1:796265/1‑136 (MQ=255)                                                                                           ctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATc                    >  1:94476/1‑139 (MQ=255)                                                                                               tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                <  2:27976/139‑1 (MQ=255)                                                                                               tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                <  1:853920/139‑1 (MQ=255)                                                                                                 cGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                >  2:797672/1‑57 (MQ=255)                                                                                                 cGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGt                                                                                                <  1:797672/57‑1 (MQ=255)                                                                                                 cGTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCt              <  1:19268/139‑1 (MQ=255)                                                                                                  gTGATGATTTTGCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCtg             <  2:175922/139‑1 (MQ=255)                                                                                                             gCTTTCGAAGATGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAAcc  <  2:897364/139‑1 (MQ=255)                                                                                                                        |                                                                                                                                ATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTTCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACC  >  NC_000913/1548241‑1548487

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