Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,668	C→T	D567N (GAC→AAC)	fdhF ←	formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,668	0	C	T	100.0%	122.8 / NA	39	D567N (GAC→AAC)	fdhF	formate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

TATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCA  >  NC_000913/4297533‑4297799                                                                                                                                        |                                                                                                                                    tATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                  <  1:879437/139‑1 (MQ=255)           tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                        >  2:542389/1‑139 (MQ=255)           tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                        >  1:389214/1‑139 (MQ=255)                        tcatcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTga                                                                                                           <  1:542389/139‑1 (MQ=255)                                   cgccgcCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGcc                                                                                                >  1:441479/1‑139 (MQ=255)                                             gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                             >  1:967019/1‑132 (MQ=255)                                             gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                             <  2:967019/132‑1 (MQ=255)                                              ccGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                                     <  2:1086318/139‑1 (MQ=255)                                                         ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  2:109989/139‑1 (MQ=255)                                                           ggTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  1:56645/137‑1 (MQ=255)                                                           ggTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          >  2:56645/1‑137 (MQ=255)                                                           ggTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTtc                                                                        >  1:658779/1‑139 (MQ=255)                                                            gTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAACCAGTCGCAGGTGAAGAACTGAGCCAGACCGTTCGGGGTATCAAACTtct                                                                       >  2:450323/1‑139 (MQ=255)                                                              cATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                                        <  2:880711/120‑1 (MQ=255)                                                              cATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                                        >  1:880711/1‑120 (MQ=255)                                                                tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                                   >  2:431376/1‑139 (MQ=255)                                                                tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                                   >  1:457757/1‑139 (MQ=255)                                                                tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                                   >  1:109513/1‑139 (MQ=255)                                                                           gagTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaag                                                        <  2:1059557/139‑1 (MQ=255)                                                                               agtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGt                                                    <  2:389214/139‑1 (MQ=255)                                                                                     cAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAAc                                                                                                      <  2:877619/83‑1 (MQ=255)                                                                                     cAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAAc                                                                                                      >  1:877619/1‑83 (MQ=255)                                                                                                   gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                  >  1:1011864/1‑41 (MQ=255)                                                                                                   gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                  <  2:1011864/41‑1 (MQ=255)                                                                                                   gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAg                                >  1:527980/1‑139 (MQ=255)                                                                                                       aCAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       <  2:1099341/64‑1 (MQ=255)                                                                                                       aCAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       >  1:1099341/1‑64 (MQ=255)                                                                                                         aGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATcgc                          <  2:109513/139‑1 (MQ=255)                                                                                                              cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGc                                                                                                                 <  2:593559/47‑1 (MQ=255)                                                                                                              cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGc                                                                                                                 >  1:593559/1‑47 (MQ=255)                                                                                                                tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCggtatgg                   <  2:658779/139‑6 (MQ=255)                                                                                                                     tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAGCTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  1:627574/79‑1 (MQ=255)                                                                                                                     tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          >  2:627574/1‑79 (MQ=255)                                                                                                                       cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                        >  1:852923/1‑111 (MQ=255)                                                                                                                       cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                        <  2:852923/111‑1 (MQ=255)                                                                                                                       cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtga            >  1:880433/1‑139 (MQ=255)                                                                                                                          gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                                  >  2:501627/1‑82 (MQ=255)                                                                                                                          gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                                  <  1:501627/82‑1 (MQ=255)                                                                                                                                 aGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCa  >  1:119442/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    TATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCA  >  NC_000913/4297533‑4297799

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