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

ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGG  >  NC_000913/4297538‑4297782                                                                                                                                   |                                                                                                                   aTTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcg                                                                                                            <  1:619847/139‑1 (MQ=255)      tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                       >  2:1310915/1‑139 (MQ=255)           aGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTc                                                                                                  >  1:729822/1‑139 (MQ=255)             ccAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGc                                                                                                <  1:663403/139‑1 (MQ=255)                ggTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAgg                                                                                             <  2:20377/139‑1 (MQ=255)                       cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  1:335467/139‑1 (MQ=255)                              cgccgcCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGcc                                                                               >  2:28682/1‑139 (MQ=255)                               gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCa                                                                              <  2:571141/139‑1 (MQ=255)                                      gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                       <  2:301247/139‑1 (MQ=255)                                      gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                       <  1:761452/139‑1 (MQ=255)                                        gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                     <  2:11712/139‑1 (MQ=255)                                         ccGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  2:807226/121‑1 (MQ=255)                                         ccGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      >  1:807226/1‑121 (MQ=255)                                         ccGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                    <  1:28682/139‑1 (MQ=255)                                           gCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgggg                                                                  <  1:724054/139‑1 (MQ=255)                                                 ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCaa                                                            >  2:337081/1‑139 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                         >  1:1071302/1‑139 (MQ=255)                                                            cGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                 <  1:516022/139‑1 (MQ=255)                                                              aaCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                               <  1:1057559/139‑1 (MQ=255)                                                              aaCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                               <  1:1222293/139‑1 (MQ=255)                                                                cGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACa                                             <  2:116625/139‑1 (MQ=255)                                                                        gtagtGACCAACTTCACGCACCGTTGACATTACCATCGGGTACTCGTCGGTGAGTTTGTTGAGTGGCGCTACCCAGTAGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                     <  1:882565/139‑1 (MQ=255)                                                                        gtagtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                     <  1:337081/139‑1 (MQ=255)                                                                        gtagtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                     <  1:1310915/139‑1 (MQ=255)                                                                             gACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                     <  2:456750/102‑1 (MQ=255)                                                                             gACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                     >  1:456750/1‑102 (MQ=255)                                                                              aCCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTcccc                               >  1:32878/1‑139 (MQ=255)                                                                              aCCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTcccc                               >  1:591330/1‑139 (MQ=255)                                                                                 aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGa                            >  2:1256463/1‑139 (MQ=255)                                                                                   cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  1:78786/79‑1 (MQ=255)                                                                                   cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      >  2:78786/1‑79 (MQ=255)                                                                                        cGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                       <  2:1065252/89‑1 (MQ=255)                                                                                        cGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt                                                                       >  1:1065252/1‑89 (MQ=255)                                                                                        cGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCAt                     <  1:1256463/139‑1 (MQ=255)                                                                                                       aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGc      >  1:975941/1‑139 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      >  1:1217028/1‑55 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  2:1217028/55‑1 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg  <  1:48193/139‑1 (MQ=255)                                                                                                                tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                               >  2:1226939/1‑89 (MQ=255)                                                                                                                tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                               <  1:1226939/89‑1 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  2:661336/48‑1 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      >  1:661336/1‑48 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                         >  1:1327306/1‑77 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTAGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                         <  2:1327306/77‑1 (MQ=255)                                                                                                                     gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTGTCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                       >  2:142539/1‑108 (MQ=255)                                                                                                                     gtcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTGTCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                       <  1:142539/108‑1 (MQ=255)                                                                                                                        ggTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      <  1:700883/42‑1 (MQ=255)                                                                                                                        ggTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                      >  2:700883/1‑42 (MQ=255)                                                                                                                         gTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAg                                                                             <  2:1180101/50‑1 (MQ=255)                                                                                                                         gTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAg                                                                             >  1:1180101/1‑50 (MQ=255)                                                                                                                         gTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                     <  1:527841/58‑1 (MQ=255)                                                                                                                         gTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                     >  2:527841/1‑58 (MQ=255)                                                                                                                            aGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATcgc         >  2:1370006/1‑115 (MQ=255)                                                                                                                            aGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATcgc         <  1:1370006/115‑1 (MQ=255)                                                                                                                                   |                                                                                                                   ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGG  >  NC_000913/4297538‑4297782

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