Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,495,261	A→G	K417K (AAA→AAG)	nirB →	nitrite reductase catalytic subunit NirB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,495,261	0	A	G	100.0%	92.5 / NA	24	K417K (AAA→AAG)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

CCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAA  >  NC_000913/3495217‑3495306                                             |                                              ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTg                                            <  1:2098308/48‑1 (MQ=255) ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGc                                           >  1:1044539/1‑49 (MQ=255)   aGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCg                                         >  1:1823694/1‑49 (MQ=255)   aGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCg                                         <  1:1809870/49‑1 (MQ=255)    gCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCgg                                        <  1:2228733/49‑1 (MQ=255)    gCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCg                                         >  1:2095172/1‑48 (MQ=255)      acacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGAc                                      <  1:1094532/49‑1 (MQ=255)      acacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGAc                                      <  1:1033370/49‑1 (MQ=255)      acacTCGGGTAGCGGCAAGCCGTCTATCGGTGTAGAAAAGCTGCGGGa                                       >  1:338972/1‑48 (MQ=255)            gggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgc                                 <  1:1992762/48‑1 (MQ=255)            gggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                                  <  1:824047/47‑1 (MQ=255)             ggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgcg                                >  1:2249892/1‑48 (MQ=255)                aGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCaaa                            >  1:1459854/1‑49 (MQ=255)                     cAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGc                       >  1:2216564/1‑49 (MQ=255)                        gCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTcc                    >  1:110196/1‑49 (MQ=255)                             cTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCtt               >  1:1065818/1‑49 (MQ=255)                              tATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTc              >  1:1933076/1‑49 (MQ=255)                              tATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTc              >  1:586007/1‑49 (MQ=255)                                tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCtt               <  1:1145061/46‑1 (MQ=255)                                tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCg             <  1:694874/48‑1 (MQ=255)                                tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGa            >  1:241270/1‑49 (MQ=255)                                    tgtTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGt         >  1:167176/1‑48 (MQ=255)                                          tAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCaaa  <  1:1139492/49‑1 (MQ=255)                                          tAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCaaa  <  1:2283312/49‑1 (MQ=255)                                             |                                              CCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAA  >  NC_000913/3495217‑3495306

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