Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,494,795	(G)5→6	coding (785/2544 nt)	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,494,790	1	.	G	100.0%	57.0 / NA	18	R260R (CGT→CGG)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (7/11);  total (7/11)

AGGTGTTGAAGCGCGTAAAACCATGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGT‑GGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCA  >  NC_000913/3494658‑3494912                                                                                                                                      |                                                                                                                           aGGTGTTGAAGCGCGTAAGACCGTGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCAGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGt                                                                                                                             <  1:252425/133‑1 (MQ=255)   gtgtTGAAGCGCGTAAAACCATGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGT‑gg                                                                                                                          <  2:283360/133‑1 (MQ=255)    tgtTGAAGCGCGTAAAACCATGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGT‑ggg                                                                                                                         >  1:414776/1‑133 (MQ=255)              cgTAAAACCATGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTc                                                                                                                <  2:414776/133‑1 (MQ=255)                       aTGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGAt                                                                                                       <  2:274322/133‑1 (MQ=255)                                       gCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCg                                                                                       <  2:479218/133‑1 (MQ=255)                                       gCAGCGAACTCGAAGTCGACTTTATCGTTTTCCCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCg                                                                                       <  2:130906/133‑1 (MQ=255)                                            gAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCg                                                                                  <  1:262925/133‑1 (MQ=255)                                             aaCTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCgg                                                                                 <  1:266952/133‑1 (MQ=255)                                                          cTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCAt                                                                    <  2:442062/133‑1 (MQ=255)                                                             tATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCgg                                                                 >  2:77381/1‑133 (MQ=255)                                                                    ttctCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTATGCCATCGGTGAATgc                                                          >  2:258581/1‑133 (MQ=255)                                                                         tACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATgcg                                                         <  2:118875/129‑1 (MQ=255)                                                                         tACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATgcg                                                         >  1:118875/1‑129 (MQ=255)                                                                             ggTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTgg                                                 <  2:304326/133‑1 (MQ=255)                                                                              gTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGa                                                >  1:370388/1‑133 (MQ=255)                                                                                                         cccAGTGTGGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTATGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCAc                     <  1:258581/133‑1 (MQ=255)                                                                                                               gtgGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCt               >  1:68028/1‑133 (MQ=255)                                                                                                                tgGTCTGGACGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTa              <  1:293315/133‑1 (MQ=255)                                                                                                                        aCGTTGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATgg      >  1:17914/1‑133 (MQ=255)                                                                                                                            tGCTCCGCGTGGGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAAACGTGTATTTGGTTTGGAAGCACCTGGCTACCAAATGGCGCa  >  2:147509/1‑133 (MQ=255)                                                                                                                                      |                                                                                                                           AGGTGTTGAAGCGCGTAAAACCATGCGTTTTGCCGACGGCAGCGAACTGGAAGTCGACTTTATCGTCTTCTCTACCGGTATCCGTCCGCGCGATAAGCTGGCAACCCAGTGTGGTCTGGACGTTGCTCCGCGT‑GGGGGTATTGTCATTAATGATTCCTGCCAGACTTCCGATCCGGATATCTACGCCATCGGTGAATGCGCAAGCTGGAACAACCGTGTATTTGGTCTGGTAGCACCTGGCTACAAAATGGCGCA  >  NC_000913/3494658‑3494912

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