Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,494,086	T→G	S26A (TCT→GCT)	nirB →	nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,494,086	0	T	G	100.0%	203.3 / NA	53	S26A (TCT→GCT)	nirB	nitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/31);  total (22/31)

ATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACGGCAT  >  NC_000913/3493942‑3494231                                                                                                                                                 |                                                                                                                                                  aTTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATg                                                                                                                                             >  1:390938/1‑151 (MQ=255)            tAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACttt                                                                                                                                  >  1:178421/1‑151 (MQ=255)             aaGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACttt                                                                                                                                  <  1:493820/150‑1 (MQ=255)              aGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTg                                                                                                                                 >  1:125589/1‑150 (MQ=255)                  gggTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATAtt                                                                                                                            <  1:666119/151‑1 (MQ=255)                  gggTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATAtt                                                                                                                            <  2:76257/151‑1 (MQ=255)                   ggTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATAtt                                                                                                                            >  1:107236/1‑150 (MQ=255)                       gCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCg                                                                                                                        <  2:54073/150‑1 (MQ=255)                         tGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGttt                                                                                                                     >  1:225783/1‑151 (MQ=255)                           aaTCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTc                                                                                                                   >  1:607450/1‑151 (MQ=255)                              cGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCtgt                                                                                                                <  2:367257/151‑1 (MQ=255)                                ttAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCtgtg                                                                                                               >  2:212140/1‑150 (MQ=255)                                   aGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTgaaga                                                                                                           <  1:778217/151‑1 (MQ=255)                                     gTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAAc                                                                                                         <  2:254452/151‑1 (MQ=255)                                       aGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACcg                                                                                                       <  1:910069/151‑1 (MQ=255)                                          cGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCgcgc                                                                                                    <  2:73978/151‑1 (MQ=255)                                          cGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCgcgc                                                                                                    <  2:530886/151‑1 (MQ=255)                                          cGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCgcgc                                                                                                    <  1:381158/151‑1 (MQ=255)                                              aaTAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCg                                                                                                >  2:457902/1‑151 (MQ=255)                                              aaTAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCg                                                                                                >  1:754937/1‑151 (MQ=255)                                                 agaaaagaaaTCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCtt                                                                                             >  2:700566/1‑151 (MQ=255)                                                 agaaaagaaaTCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCt                                                                                              <  2:104367/150‑1 (MQ=255)                                                       gaaaTCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGAcc                                                                                       <  2:536806/151‑1 (MQ=255)                                                          aTCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCACTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACcgcg                                                                                    <  1:83413/151‑1 (MQ=255)                                                            cGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGt                                                                                   <  1:674885/150‑1 (MQ=255)                                                                           cAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTAc                                                                   <  1:149997/151‑1 (MQ=255)                                                                              aGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTtc                                                                >  1:702814/1‑151 (MQ=255)                                                                                  aGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTtctctc                                                            <  1:419905/151‑1 (MQ=255)                                                                                    aCTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTcac                                                          <  2:287893/151‑1 (MQ=255)                                                                                         cAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCacac                                                     <  2:65542/151‑1 (MQ=255)                                                                                            ttATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACAccgc                                                  >  2:496421/1‑151 (MQ=255)                                                                                             tATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACAccgc                                                  <  1:421320/150‑1 (MQ=255)                                                                                              aTCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACAccgccg                                                >  1:17267/1‑151 (MQ=255)                                                                                              aTCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACAccgcc                                                 <  1:894925/150‑1 (MQ=255)                                                                                               tCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACAccgccg                                                >  1:335010/1‑150 (MQ=255)                                                                                                       ggTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGt                                       <  1:910461/151‑1 (MQ=255)                                                                                                            ggTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCt                                   <  1:407200/150‑1 (MQ=255)                                                                                                                ggCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTgc                              >  1:731175/1‑151 (MQ=255)                                                                                                                 gCCATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTgcg                             <  1:711709/151‑1 (MQ=255)                                                                                                                  ccATCGCTTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTgcgc                            <  2:97281/151‑1 (MQ=255)                                                                                                                      cGCCTTATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGa                          >  2:266336/1‑149 (MQ=255)                                                                                                                         tttATCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCt                     >  2:544530/1‑151 (MQ=255)                                                                                                                            aTCGAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTc                   <  2:178418/150‑1 (MQ=255)                                                                                                                               gAAGATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACg               >  2:357343/1‑151 (MQ=255)                                                                                                                                  gATCTTCTTGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACg               >  2:548941/1‑148 (MQ=255)                                                                                                                                         ttGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACgg     >  1:909096/1‑151 (MQ=255)                                                                                                                                         ttGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACgg     >  2:653301/1‑151 (MQ=255)                                                                                                                                         ttGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCAGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACgg     >  2:235511/1‑151 (MQ=255)                                                                                                                                          tGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAacac       <  2:867306/148‑1 (MQ=255)                                                                                                                                          tGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACgg     <  2:107235/150‑1 (MQ=255)                                                                                                                                          tGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACgg     <  2:38982/150‑1 (MQ=255)                                                                                                                                          tGATAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACg      <  2:277182/149‑1 (MQ=255)                                                                                                                                            aTAAAGCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACGGCAt  <  1:484743/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGAACCGCGCATCGCTTATGACCGCGTACACCTCTCGTCTTACTTCTCTCACCACACCGCCGAAGAGCTGTCGCTGGTGCGCGAAGGCTTCTACGAGAAACACGGCAT  >  NC_000913/3493942‑3494231

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