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%	174.1 / NA	54	D567N (GAC→AAC)	fdhF	formate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/28);  total (26/28)

TGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATG  >  NC_000913/4297541‑4297792                                                                                                                                |                                                                                                                             tgtgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTa                                                                                                                   >  1:838987/1‑139 (MQ=255)                 catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaa                                                                                                   <  1:1078271/139‑1 (MQ=255)                    cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                <  2:838987/139‑1 (MQ=255)                    cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                <  2:1020848/139‑1 (MQ=255)                    cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGGACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                <  2:1298664/139‑1 (MQ=255)                             ccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAg                                                                                       <  1:1278703/139‑1 (MQ=255)                                     gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                               <  1:879621/139‑1 (MQ=255)                                              ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcg                                                                                                                      <  2:1384710/91‑1 (MQ=255)                                              ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcg                                                                                                                      >  1:1384710/1‑91 (MQ=255)                                              ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACGCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCaa                                                                      >  2:997692/1‑139 (MQ=255)                                              ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCaa                                                                      >  1:50590/1‑139 (MQ=255)                                                aCCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAAc                                                                    <  2:1507643/139‑1 (MQ=255)                                                 ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                   <  2:638335/139‑1 (MQ=255)                                                  cGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACtt                                                                  <  2:496117/139‑1 (MQ=255)                                                         cGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                           <  2:1406927/139‑1 (MQ=255)                                                                        gtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTc                                            >  2:119717/1‑139 (MQ=255)                                                                              aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                           >  1:366512/1‑54 (MQ=255)                                                                              aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                           <  2:366512/54‑1 (MQ=255)                                                                              aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGa                                      >  2:1462692/1‑139 (MQ=255)                                                                               aCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                           <  2:920440/53‑1 (MQ=255)                                                                               aCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                           >  1:920440/1‑53 (MQ=255)                                                                               aCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGAt                                     >  2:687276/1‑139 (MQ=255)                                                                                cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                      >  1:1404519/1‑89 (MQ=255)                                                                                cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                      <  2:1404519/89‑1 (MQ=255)                                                                                 ttCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCg                                   <  2:50590/139‑1 (MQ=255)                                                                                          cGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                           >  2:394830/1‑106 (MQ=255)                                                                                          cGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                           <  1:394830/106‑1 (MQ=255)                                                                                          cGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGaa                          >  2:1271216/1‑139 (MQ=255)                                                                                          cGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGaa                          >  2:1315710/1‑139 (MQ=255)                                                                                          cGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGaa                          >  2:718572/1‑139 (MQ=255)                                                                                                cAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              <  1:369882/81‑1 (MQ=255)                                                                                                cAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              >  2:369882/1‑81 (MQ=255)                                                                                                  gTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              <  1:499276/79‑1 (MQ=255)                                                                                                  gTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              >  2:499276/1‑79 (MQ=255)                                                                                                    aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                               <  1:1257896/76‑1 (MQ=255)                                                                                                    aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                               >  2:1257896/1‑76 (MQ=255)                                                                                                    aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                            >  1:401555/1‑95 (MQ=255)                                                                                                    aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                            <  2:401555/95‑1 (MQ=255)                                                                                                      cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgggg                                                                            >  1:1396046/1‑77 (MQ=255)                                                                                                      cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgggg                                                                            <  2:1396046/77‑1 (MQ=255)                                                                                                        tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                         <  1:81375/94‑1 (MQ=255)                                                                                                        tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                         >  2:81375/1‑94 (MQ=255)                                                                                                        tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg            >  1:110629/1‑139 (MQ=255)                                                                                                             tACTCGTCGGTTAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCTAGGCCACt       <  1:997692/139‑1 (MQ=255)                                                                                                             tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                               >  1:1481061/1‑67 (MQ=255)                                                                                                             tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                               <  2:1481061/67‑1 (MQ=255)                                                                                                             tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                   <  2:414315/79‑1 (MQ=255)                                                                                                             tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                   >  1:414315/1‑79 (MQ=255)                                                                                                               cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtga     >  2:659986/1‑139 (MQ=255)                                                                                                                  gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                >  2:1068081/1‑45 (MQ=255)                                                                                                                  gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                <  1:1068081/45‑1 (MQ=255)                                                                                                                  gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtgaatg  <  2:110629/139‑1 (MQ=255)                                                                                                                          gTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              <  1:665106/55‑1 (MQ=255)                                                                                                                          gTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                              >  2:665106/1‑55 (MQ=255)                                                                                                                                |                                                                                                                             TGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATG  >  NC_000913/4297541‑4297792

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: ‑