Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,548,360	C→T	D320D (GAC→GAT)	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,360	0	C	T	100.0%	216.8 / NA	68	D320D (GAC→GAT)	fdnG	formate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/38);  total (30/38)

TCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACCTGCT  >  NC_000913/1548229‑1548491                                                                                                                                    |                                                                                                                                     tCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGt                                                                                                                               <  2:1082941/139‑1 (MQ=255)   cGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑c                                                                                                                            >  2:344246/1‑139 (MQ=255)     gTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAg                                                                                                                          >  2:1324309/1‑139 (MQ=255)           aCATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGAATTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTa                                                                                                                    <  2:1196691/139‑1 (MQ=255)                    tCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAACCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGa                                                                                                           >  1:1081269/1‑139 (MQ=255)                     ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGtt                                                                                                                              <  2:1338555/120‑1 (MQ=255)                     ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGtt                                                                                                                              >  1:1338555/1‑120 (MQ=255)                     ccTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGaa                                                                                                          >  2:408210/1‑139 (MQ=255)                         tCTGGCGTTTTGCGCTACCTGATCTAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAAc                                                                                                      >  2:979362/1‑139 (MQ=255)                              cGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCaa                                                                                                 >  2:998966/1‑139 (MQ=255)                                  ttGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTacg                                                                                                                  >  1:1280661/1‑118 (MQ=255)                                  ttGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTacg                                                                                                                  <  2:1280661/118‑1 (MQ=255)                                   tGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACAAACGCCAGACTGCTGGTGCGTGTTAATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGa                                                                                            >  2:931242/1‑139 (MQ=255)                                       cTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                        <  1:998966/139‑1 (MQ=255)                                        tACCTGATCTAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                       <  1:979362/139‑1 (MQ=255)                                        tACCTGATCGAAAACAACACAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGCTTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                       <  2:1081269/139‑1 (MQ=255)                                          ccTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCg                                                                                     >  2:358856/1‑139 (MQ=255)                                              aTCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCt                                                                                 >  2:1178452/1‑139 (MQ=255)                                                cGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGTCTGCTGGTGCGTGATGATTTTGTTTTCGAAGATGGTCTGGT‑CAGCGGCTATGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                               <  1:1053695/139‑1 (MQ=255)                                                    aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                               <  2:720442/136‑1 (MQ=255)                                                    aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTTCCAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTgg                                                                               >  1:720442/1‑136 (MQ=255)                                                    aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGaaa                                                                                                         <  1:674812/109‑1 (MQ=255)                                                    aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGaaa                                                                                                         >  2:674812/1‑109 (MQ=255)                                                    aacaacAAAATCAACGCCGAATACGTTAAGCATTACACCAAAGCCAGCCTGATGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACt                                                                           >  2:1022743/1‑139 (MQ=255)                                                             aTCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCg                                                                  <  1:1022743/139‑1 (MQ=255)                                                                  cGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaa                                                             <  1:1178452/139‑1 (MQ=255)                                                                  cGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGaa                                                             <  1:678250/139‑1 (MQ=255)                                                                              tAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAGACGTCAATACGATaaa                                                                                        >  2:633502/1‑100 (MQ=255)                                                                              tAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATaaa                                                                                        <  1:633502/100‑1 (MQ=255)                                                                                 gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGa                                                                 <  1:414561/120‑1 (MQ=255)                                                                                 gCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGa                                                                 >  2:414561/1‑120 (MQ=255)                                                                                   aTTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGTTCGATGAAAACGGCTATGCGAAAcg                                            <  2:275382/139‑1 (MQ=255)                                                                                     tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAAc                                                                                                      <  1:980942/79‑1 (MQ=255)                                                                                     tACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAAc                                                                                                      >  2:980942/1‑79 (MQ=255)                                                                                      acacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAACCGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGa                                         >  2:1184019/1‑139 (MQ=255)                                                                                       cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCaa                                                                                                 <  1:961083/82‑1 (MQ=255)                                                                                       cacCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCaa                                                                                                 >  2:961083/1‑82 (MQ=255)                                                                                            aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCa                                                                                                  >  1:14362/1‑76 (MQ=255)                                                                                            aCGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCa                                                                                                  <  2:14362/76‑1 (MQ=255)                                                                                             cGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   <  1:1645007/90‑1 (MQ=255)                                                                                             cGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   >  2:1645007/1‑90 (MQ=255)                                                                                               ccAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTa                                                                                                                    <  2:668582/55‑1 (MQ=255)                                                                                               ccAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTa                                                                                                                    >  1:668582/1‑55 (MQ=255)                                                                                               ccAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACt                                <  2:177650/139‑1 (MQ=255)                                                                                                 aGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGa                              <  1:408210/139‑1 (MQ=255)                                                                                                  gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   <  1:1271301/85‑1 (MQ=255)                                                                                                  gCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   >  2:1271301/1‑85 (MQ=255)                                                                                                    ctgctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTc                           >  2:551030/1‑139 (MQ=255)                                                                                                       ctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCCGAAAAACGTCa                                                                                                  >  1:1166169/1‑65 (MQ=255)                                                                                                       ctgGTGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCCGAAAAACGTCa                                                                                                  <  2:1166169/65‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  1:344246/139‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  1:358856/139‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  1:705663/139‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  1:739421/139‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  2:1244766/139‑1 (MQ=255)                                                                                                           tGCGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCcgcg                    <  2:220763/139‑1 (MQ=255)                                                                                                             cGTGATGATTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCt                  <  1:1324309/139‑1 (MQ=255)                                                                                                               tgatgaTTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGAt                                        <  2:1425351/115‑1 (MQ=255)                                                                                                               tgatgaTTTTGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGAt                                        >  1:1425351/1‑115 (MQ=255)                                                                                                                     ttttGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTgg          <  1:1058607/139‑1 (MQ=255)                                                                                                                       ttGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGaa        <  1:551030/139‑1 (MQ=255)                                                                                                                        tGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   <  2:720184/63‑1 (MQ=255)                                                                                                                        tGCTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTc                                                                                   >  1:720184/1‑63 (MQ=255)                                                                                                                          cTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAACCACTGACTCATCCGCGCTGTGTGTGGa         <  2:986763/135‑1 (MQ=255)                                                                                                                          cTTTCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTAGGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGa         >  1:986763/1‑135 (MQ=255)                                                                                                                             tCGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACctgct  >  1:1168710/1‑139 (MQ=255)                                                                                                                              cGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                       >  1:1596810/1‑53 (MQ=255)                                                                                                                              cGAAGATGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAAt                                                                                       <  2:1596810/53‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     TCCGGTACGGACATTACGTTCCTGTCTGGCGTTTTGCGCTACCTGATCGAAAACAACAAAATCAACGCCGAATACGTTAAGCATTACACCAACGCCAGCCTGCTGGTGCGTGATGATTTTGCTTTCGAAGACGGTCTGTT‑CAGCGGCTACGACGCTGAAAAACGTCAATACGATAAATCGTCCTGGAACTATCAGCTCGATGAAAACGGCTATGCGAAACGCGATGAAACACTGACTCATCCGCGCTGTGTGTGGAACCTGCT  >  NC_000913/1548229‑1548491

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