Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,548,145	A→C	N249H (AAC→CAC)	fdnG →	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,548,145	0	A	C	100.0%	199.3 / NA	57	N249H (AAC→CAC)	fdnG	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/29);  total (28/29)

AAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAAT  >  NC_000913/1548006‑1548290                                                                                                                                            |                                                                                                                                                  aaGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAacgac                                                                                                                                        >  2:206043/1‑151 (MQ=255)    tCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGc                                                                                                                                      >  1:186150/1‑150 (MQ=255)         cTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCtt                                                                                                                                >  2:513518/1‑151 (MQ=255)             aaCATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATc                                                                                                                            >  2:7696/1‑151 (MQ=255)             aaCATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATc                                                                                                                            >  1:493701/1‑151 (MQ=255)                    ggTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCg                                                                                                                     >  2:545945/1‑151 (MQ=255)                     gTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGa                                                                                                                    <  2:76886/151‑1 (MQ=255)                      tCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGAt                                                                                                                   >  2:421631/1‑151 (MQ=255)                          ggTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATccc                                                                                                                <  1:480599/150‑1 (MQ=255)                            tGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGt                                                                                                             >  2:51747/1‑151 (MQ=255)                             gCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGtt                                                                                                            >  1:520942/1‑151 (MQ=255)                              cGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGttt                                                                                                           >  1:716511/1‑151 (MQ=255)                               gATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGtttt                                                                                                          <  2:182496/151‑1 (MQ=255)                                 tGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTAc                                                                                                        >  1:269260/1‑151 (MQ=255)                                     caaccaCTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGt                                                                                                    <  2:247992/151‑1 (MQ=255)                                     caaccaCTGGGTGGATATCAAAAACGCTAACGGCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGt                                                                                                    <  1:668776/151‑1 (MQ=255)                                          aCTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGc                                                                                               >  2:808522/1‑151 (MQ=255)                                                             cGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTAc                                                                            <  2:265653/151‑1 (MQ=255)                                                                tAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTAc                                                                            >  1:158992/1‑148 (MQ=255)                                                                 aaCGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTAcgc                                                                          <  1:539168/149‑1 (MQ=255)                                                                      cgtGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTAtt                                                                   >  2:154957/1‑151 (MQ=255)                                                                       gtgatggtgatgGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTc                                                                  <  1:638360/151‑1 (MQ=255)                                                                         gatggtgatggGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGt                                                                <  1:855581/151‑1 (MQ=255)                                                                            ggtgatggGCGGTAACGCGGCTCAAGCGCATCCCGTCGGTTTCCGCTGGGCCATGGAAGCGAAACACAACAACCACGCAACCTTGATCGTTGTCGAGCCCCGTTTTACGTGTCCCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTcc                                                             <  1:1430/151‑1 (MQ=255)                                                                                  gggCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTAc                                                        <  2:868848/150‑1 (MQ=255)                                                                                  gggCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACg                                                       >  1:575206/1‑151 (MQ=255)                                                                                          aCGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTAc                                               <  1:585779/151‑1 (MQ=255)                                                                                           cGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTAc                                               <  2:89909/150‑1 (MQ=255)                                                                                                  gAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTg                                        <  2:8878/150‑1 (MQ=255)                                                                                                  gAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTg                                        >  1:478963/1‑150 (MQ=255)                                                                                                  gAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGt                                       >  1:355523/1‑151 (MQ=255)                                                                                                   aaGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTc                                      <  1:531517/151‑1 (MQ=255)                                                                                                   aaGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTc                                      >  1:123738/1‑151 (MQ=255)                                                                                                   aaGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTc                                      >  2:545868/1‑151 (MQ=255)                                                                                                      cgcATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTgg                                   >  2:90084/1‑151 (MQ=255)                                                                                                       gcATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGc                                  <  2:747433/151‑1 (MQ=255)                                                                                                       gcATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGc                                  <  2:729355/151‑1 (MQ=255)                                                                                                        cATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCg                                 <  1:729202/151‑1 (MQ=255)                                                                                                        cATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCg                                 >  2:718734/1‑151 (MQ=255)                                                                                                        cATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCg                                 <  2:50109/151‑1 (MQ=255)                                                                                                         aTCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGt                                >  2:621718/1‑151 (MQ=255)                                                                                                           cccGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGtt                               <  2:397180/150‑1 (MQ=255)                                                                                                             cGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTg                            <  1:455025/151‑1 (MQ=255)                                                                                                             cGTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGAACCGCTTCTGCGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTAAGTTCCTGTCTGGCGttt                              >  1:247023/1‑149 (MQ=255)                                                                                                              gTCGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTgc                           <  2:635444/151‑1 (MQ=255)                                                                                                                cGGTTTCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTgcgc                         >  2:456186/1‑151 (MQ=255)                                                                                                                   tttCCGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTAc                      >  1:318889/1‑151 (MQ=255)                                                                                                                      ccGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTg                   <  2:855558/151‑1 (MQ=255)                                                                                                                       cGCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTg                   >  2:859348/1‑150 (MQ=255)                                                                                                                        gCTGGGCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGAt                 <  1:704174/151‑1 (MQ=255)                                                                                                                           gggCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATc                <  2:676956/149‑1 (MQ=255)                                                                                                                           gggCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCg               >  1:303333/1‑150 (MQ=255)                                                                                                                           gggCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGa              <  1:585196/151‑1 (MQ=255)                                                                                                                            ggCGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGaa             <  1:850695/151‑1 (MQ=255)                                                                                                                              cGATGGAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGaaaa           <  2:918596/151‑1 (MQ=255)                                                                                                                                   gAAGCGAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAacaaca      <  1:409269/151‑1 (MQ=255)                                                                                                                                        gAAACACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAAt  >  1:674731/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                                  AAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGCGTACCGCTTCTGTGGCGGATATTTACGCGCCTATTCGTTCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAAT  >  NC_000913/1548006‑1548290

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