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

TGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCA  >  NC_000913/4297543‑4297799                                                                                                                              |                                                                                                                                    tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                        >  2:1207475/1‑139 (MQ=255) tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                        >  2:1071125/1‑139 (MQ=255)     tAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGt                                                                                                                    <  2:351423/139‑1 (MQ=255)      aGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTc                                                                                                                   >  1:222185/1‑139 (MQ=255)           ggTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATT‑GCGCTACCCAGTCGCAGGt                                                                                                             >  1:585364/1‑139 (MQ=255)                 tcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaaga                                                                                                        >  2:1215800/1‑139 (MQ=255)                  cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       <  2:1419503/139‑1 (MQ=255)                  cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       >  2:1206619/1‑139 (MQ=255)                            cgccAGTGCCGCACAGTGACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                             <  1:578669/139‑1 (MQ=255)                                 gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                                        <  1:1207475/139‑1 (MQ=255)                                   gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                      <  2:9516/139‑1 (MQ=255)                                        acaGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTa                                                                                 <  2:1164523/139‑1 (MQ=255)                                            ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                  >  2:355159/1‑86 (MQ=255)                                            ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                  <  1:355159/86‑1 (MQ=255)                                              aCCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAAc                                                                           >  1:719362/1‑139 (MQ=255)                                               ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  1:1215800/139‑1 (MQ=255)                                               ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  1:996991/139‑1 (MQ=255)                                                     aTCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCtt                                                                    >  2:908715/1‑139 (MQ=255)                                                      tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       <  1:1376824/103‑1 (MQ=255)                                                      tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       >  2:1376824/1‑103 (MQ=255)                                                      tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                                   >  2:645824/1‑139 (MQ=255)                                                       cGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                                  <  2:55513/139‑1 (MQ=255)                                                        gAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATT‑GCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                                <  2:585364/139‑1 (MQ=255)                                                        gAACGGCAAGAGTAGTGACCAACTTAACGCACCGTTGACAGTACCATCGTGTACACGTCGGTGAGTTTGTTGAATGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaa                                                                 <  1:619654/139‑1 (MQ=255)                                                         aaCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                                <  1:1365058/139‑1 (MQ=255)                                                          aCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          >  1:1246168/1‑128 (MQ=255)                                                          aCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  2:1246168/128‑1 (MQ=255)                                                                 gagTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaag                                                        <  2:222185/139‑1 (MQ=255)                                                                      gtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                      <  1:255593/104‑1 (MQ=255)                                                                      gtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                      >  2:255593/1‑104 (MQ=255)                                                                      gtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          >  1:60401/1‑116 (MQ=255)                                                                      gtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  2:60401/116‑1 (MQ=255)                                                                       tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                                <  2:1319784/125‑1 (MQ=255)                                                                       tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                                >  1:1319784/1‑125 (MQ=255)                                                                       tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTcc                                                  >  1:1325830/1‑139 (MQ=255)                                                                         aCCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTcccc                                                >  2:286631/1‑139 (MQ=255)                                                                              cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          <  2:143560/108‑1 (MQ=255)                                                                              cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                          >  1:143560/1‑108 (MQ=255)                                                                                 cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                        <  1:1071125/139‑1 (MQ=255)                                                                                         gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAg                                >  1:1377357/1‑139 (MQ=255)                                                                                                 tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCgg                        >  1:1412659/1‑139 (MQ=255)                                                                                                  aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGc                       >  1:1223442/1‑139 (MQ=255)                                                                                                    cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCaa                     <  2:1377357/139‑1 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       >  2:969289/1‑55 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                       <  1:969289/55‑1 (MQ=255)                                                                                                          gTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCAc               <  2:719362/139‑1 (MQ=255)                                                                                                             cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGAGAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtga            >  2:878164/1‑139 (MQ=255)                                                                                                              tcgtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtgaa           >  2:803726/1‑139 (MQ=255)                                                                                                               cgtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                      >  2:1290737/1‑63 (MQ=255)                                                                                                               cgtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                      <  1:1290737/63‑1 (MQ=255)                                                                                                                gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAg                                                             >  2:107296/1‑87 (MQ=255)                                                                                                                gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAg                                                             <  1:107296/87‑1 (MQ=255)                                                                                                                       aGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCa  <  2:1223442/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                    TGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCA  >  NC_000913/4297543‑4297799

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