Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,852,276	T→C	A47A (GCT→GCC)	hypD →	Fe‑(CN)2CO cofactor assembly scaffold protein HypD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,852,276	0	T	C	100.0%	98.5 / NA	29	A47A (GCT→GCC)	hypD	Fe‑(CN)2CO cofactor assembly scaffold protein HypD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/10);  total (19/10)

GAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGT  >  NC_000913/2852199‑2852353                                                                              |                                                                              gAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATc                                                                            >  2:516313/1‑81 (MQ=255)  aaCGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCt                                                                           >  2:516288/1‑81 (MQ=255)      tGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTaa                                                                       <  2:516255/81‑1 (MQ=255)       gCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTaaa                                                                      >  1:516297/1‑81 (MQ=255)        cTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAAt                                                                     <  1:516256/81‑1 (MQ=255)         ttCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAAtt                                                                    >  1:516293/1‑81 (MQ=255)              aTCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGACGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGcc                                                               >  1:516309/1‑81 (MQ=255)                ctctctTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTc                                                             >  2:516321/1‑81 (MQ=255)                      tACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAg                                                       >  1:516304/1‑81 (MQ=255)                        cacCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGtt                                                     <  1:516257/81‑1 (MQ=255)                         acCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTa                                                    <  2:516258/81‑1 (MQ=255)                               gAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCg                                              >  1:516320/1‑81 (MQ=255)                                aaCGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCg                                              >  2:516326/1‑80 (MQ=255)                                      ctctGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACg                                       >  2:516303/1‑81 (MQ=255)                                            ggATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGt                                 >  2:516336/1‑81 (MQ=255)                                             gATTATGGAAGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGtt                                >  2:516317/1‑81 (MQ=255)                                                     aaGTGTGTGGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCAACg                        <  1:516259/81‑1 (MQ=255)                                                            tgGCGGTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCggg                 >  2:516316/1‑81 (MQ=255)                                                            tgGCGGTCATA‑CCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCggg                 <  1:516260/80‑1 (MQ=255)                                                                 gTCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGcc            >  1:516310/1‑81 (MQ=255)                                                                  tCATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGccc           <  2:516261/81‑1 (MQ=255)                                                                   cATACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCg          >  1:516312/1‑81 (MQ=255)                                                                     tACCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGgt        >  1:516319/1‑81 (MQ=255)                                                                      aCCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGgtg       >  2:516322/1‑81 (MQ=255)                                                                      aCCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGgtg       >  2:516325/1‑81 (MQ=255)                                                                      aCCCACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGgtg       >  1:516338/1‑81 (MQ=255)                                                                       cccACGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGgtgt      <  2:516262/81‑1 (MQ=255)                                                                        ccACGCCATCTTTAAATT‑GGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGc    <  2:516263/81‑1 (MQ=255)                                                                           cGCCATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGt  <  2:516264/81‑1 (MQ=255)                                                                              |                                                                              GAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGATTATGGAAGTGTGTGGCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGT  >  NC_000913/2852199‑2852353

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