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

GCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAC  >  NC_000913/3495216‑3495302                                              |                                          gCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTg                                        >  1:749806/1‑49 (MQ=255)   cAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGcc                                      <  1:908832/49‑1 (MQ=255)   cAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGcc                                      >  1:2492107/1‑49 (MQ=255)    aGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCg                                     >  1:1831048/1‑49 (MQ=255)      cacaCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                   >  1:184579/1‑49 (MQ=255)        cacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                   >  1:2811974/1‑47 (MQ=255)        cacTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                   >  1:2567751/1‑47 (MQ=255)          cTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgc                               <  1:2118704/49‑1 (MQ=255)             gggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                              <  1:1875231/47‑1 (MQ=255)              ggTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgcg                            >  1:590639/1‑48 (MQ=255)                 aGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCaaa                        >  1:845469/1‑49 (MQ=255)                 aGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCaaa                        >  1:1509004/1‑49 (MQ=255)                   cGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATc                      >  1:1467494/1‑49 (MQ=255)                    ggCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCt                     >  1:498973/1‑49 (MQ=255)                       aaGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCt                  >  1:1222006/1‑49 (MQ=255)                         gCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTcc                <  1:2399670/49‑1 (MQ=255)                          ccGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCt               >  1:2575507/1‑49 (MQ=255)                           cGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTg              >  1:419773/1‑49 (MQ=255)                            gTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGc             <  1:2187754/49‑1 (MQ=255)                            gTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGc             <  1:1310917/49‑1 (MQ=255)                              cTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCtt           >  1:498698/1‑49 (MQ=255)                                aTCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTc          >  1:2082636/1‑48 (MQ=255)                                 tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCtt           <  1:46498/46‑1 (MQ=255)                                 tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTc          <  1:2110318/47‑1 (MQ=255)                                 tCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGa        >  1:1049718/1‑49 (MQ=255)                                  cGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc       >  1:15367/1‑49 (MQ=255)                                   ggTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACg      <  1:2114999/49‑1 (MQ=255)                                        tGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAc  >  1:2311550/1‑48 (MQ=255)                                              |                                          GCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAC  >  NC_000913/3495216‑3495302

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: ‑