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%	89.8 / NA	23	K417K (AAA→AAG)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/11);  total (12/11)

CCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACC  >  NC_000913/3495217‑3495303                                             |                                           ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTg                                         <  1:1918061/48‑1 (MQ=255)  cAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGcc                                       >  1:2033269/1‑49 (MQ=255)   aGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCg                                      >  1:503405/1‑49 (MQ=255)      acacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                    >  1:665538/1‑48 (MQ=255)       cacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                    >  1:1207310/1‑47 (MQ=255)       cacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACa                                  <  1:445859/49‑1 (MQ=255)        acTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAg                                 <  1:315980/49‑1 (MQ=255)         cTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgc                                <  1:660895/49‑1 (MQ=255)         cTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgc                                >  1:24835/1‑49 (MQ=255)          tCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                               >  1:47948/1‑49 (MQ=255)            gggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgc                              <  1:1427947/48‑1 (MQ=255)               tAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCaa                          <  1:1898343/49‑1 (MQ=255)                  cGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATc                       >  1:11193/1‑49 (MQ=255)                   ggCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCt                      <  1:747766/49‑1 (MQ=255)                    gCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTg                     <  1:487559/49‑1 (MQ=255)                       aGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTc                  >  1:640777/1‑49 (MQ=255)                         ccGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCt                <  1:1853212/49‑1 (MQ=255)                          cGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTg               >  1:931799/1‑49 (MQ=255)                            tCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCt             >  1:651987/1‑49 (MQ=255)                                tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCg          <  1:707334/48‑1 (MQ=255)                                     gtTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCa    >  1:238935/1‑49 (MQ=255)                                      ttGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCa    <  1:1249091/48‑1 (MQ=255)                                        gATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAcc  >  1:1196992/1‑48 (MQ=255)                                             |                                           CCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACC  >  NC_000913/3495217‑3495303

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