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

GTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTG  >  NC_000913/3495130‑3495395                                                                                                                                    |                                                                                                                                       gtggGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                                                                                                                <  2:548998/140‑1 (MQ=255) gtggGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGa                                                                                                                                >  2:723057/1‑140 (MQ=255)     gCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgc                                                                                                                            >  1:661788/1‑140 (MQ=255)      cGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                                                                                                                           >  2:114846/1‑140 (MQ=255)      cGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                                                                                                                           >  2:485852/1‑140 (MQ=255)      cGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                                                                                                                           >  1:425238/1‑140 (MQ=255)      cGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcg                                                                                                                           >  2:575437/1‑140 (MQ=255)        aTACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgcg                                                                                                                         <  1:16134/140‑1 (MQ=255)        aTACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAgcgcg                                                                                                                         >  1:394888/1‑140 (MQ=255)                gACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTg                                                                                                                 >  2:145073/1‑140 (MQ=255)                  cTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCt                                                                                                               <  2:689460/140‑1 (MQ=255)                   tACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTc                                                                                                              <  2:34093/140‑1 (MQ=255)                         aaCCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCtt                                                                                                        <  2:723244/140‑1 (MQ=255)                           ccTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCg                                                                                                      >  1:34752/1‑140 (MQ=255)                             tgctgcAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc                                                                                                    <  2:394888/140‑1 (MQ=255)                             tgctgcAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc                                                                                                    <  1:485852/140‑1 (MQ=255)                             tgctgcAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGAc                                                                                                    <  2:661788/140‑1 (MQ=255)                               ctgcAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGt                                                                                                  >  1:681086/1‑140 (MQ=255)                                   aaCTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCAcc                                                                                              <  2:500733/140‑1 (MQ=255)                                    aCTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCa                                                                                             <  2:561651/140‑1 (MQ=255)                                             gAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATc                                                                                    <  1:114846/140‑1 (MQ=255)                                             gAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATc                                                                                    <  1:575437/140‑1 (MQ=255)                                                     tCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATtgct                                                                            <  2:413930/140‑1 (MQ=255)                                                     tCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATtgct                                                                            <  2:361805/140‑1 (MQ=255)                                                        aaCTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGcc                                                                         <  1:145073/140‑1 (MQ=255)                                                             ccGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATcaa                                                                    <  2:163770/140‑1 (MQ=255)                                                                 aaaaCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACaaa                                                                <  1:42062/140‑1 (MQ=255)                                                                    aCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGc                                                             <  2:425238/140‑1 (MQ=255)                                                                               cTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGt                                                  >  1:196626/1‑140 (MQ=255)                                                                                  aTCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTgc                                               <  1:723057/140‑1 (MQ=255)                                                                                   tCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTgcg                                              <  2:34752/140‑1 (MQ=255)                                                                                        ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCgct                                         >  2:224290/1‑140 (MQ=255)                                                                                            cacaCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCgctgaaa                                     >  2:409831/1‑140 (MQ=255)                                                                                                 tCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCgctgaaagctga                                >  2:590952/1‑140 (MQ=255)                                                                                                  cGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCgctgaaagctgaa                               >  2:323730/1‑140 (MQ=255)                                                                                                     gTAGCGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAAcc                            <  2:306668/140‑1 (MQ=255)                                                                                                         cGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                        >  2:281975/1‑140 (MQ=255)                                                                                                         cGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                        >  1:819294/1‑140 (MQ=255)                                                                                                         cGGCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                        >  2:194046/1‑140 (MQ=255)                                                                                                           gCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCg                      >  2:801127/1‑140 (MQ=255)                                                                                                           gCAAGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCg                      >  1:125088/1‑140 (MQ=255)                                                                                                              aGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                        >  2:404522/1‑135 (MQ=255)                                                                                                              aGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                        <  1:404522/135‑1 (MQ=255)                                                                                                              aGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGt                   >  1:768959/1‑140 (MQ=255)                                                                                                              aGCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGt                   <  1:409831/140‑1 (MQ=255)                                                                                                               gCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTa                  >  1:376057/1‑140 (MQ=255)                                                                                                               gCCGTCTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTa                  <  2:768959/140‑1 (MQ=255)                                                                                                                    cTATCGGTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGc             <  1:263706/140‑1 (MQ=255)                                                                                                                         ggTGTTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCgg        <  1:590952/140‑1 (MQ=255)                                                                                                                            gtTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTgg     <  1:87070/140‑1 (MQ=255)                                                                                                                            gtTGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTgg     <  1:323730/140‑1 (MQ=255)                                                                                                                             ttGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGc    >  2:11062/1‑140 (MQ=255)                                                                                                                             ttGATAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGc    >  1:83565/1‑140 (MQ=255)                                                                                                                               taTAAGCTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTg  >  2:213827/2‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       GTGGGCGATACCAGCGACTACGGTAACCTGCTGCAACTGGTGCTGAACGCTATCGAACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTG  >  NC_000913/3495130‑3495395

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