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

ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGA  >  NC_000913/4297538‑4297793                                                                                                                                   |                                                                                                                              aTTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcg                                                                                                                       <  1:136401/139‑1 (MQ=255)      tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                  >  1:525118/1‑139 (MQ=255)                ggTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAgg                                                                                                        >  1:811205/1‑139 (MQ=255)                     atcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGACTAGTGACCAACTTCACGCACCGTCGACAGTACCATCGGGTACTCGTCGGTAAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaag                                                                                                   <  1:533868/139‑1 (MQ=255)                              cgccgcCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGcc                                                                                          >  2:314989/1‑139 (MQ=255)                               gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCa                                                                                         <  2:519152/139‑1 (MQ=255)                                 cgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGa                                                                                       <  2:603063/139‑1 (MQ=255)                                                 ttACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCaa                                                                       >  1:520789/1‑139 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    <  2:811205/139‑1 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcgc                                                                                                                      >  2:69800/1‑82 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcgc                                                                                                                      <  1:69800/82‑1 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 >  1:390043/1‑103 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 <  2:390043/103‑1 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    <  2:339940/132‑1 (MQ=255)                                                           tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    >  1:339940/1‑132 (MQ=255)                                                            cGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTa                                                            <  2:797342/139‑1 (MQ=255)                                                                        gtagtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                <  2:720754/139‑1 (MQ=255)                                                                            tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTgg                                                                                                                         >  2:112328/1‑62 (MQ=255)                                                                            tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTgg                                                                                                                         <  1:112328/62‑1 (MQ=255)                                                                             gACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTCTAAACAGATAAGAAGTccc                                           >  2:532173/1‑139 (MQ=255)                                                                                   cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATc                                     >  2:438247/1‑139 (MQ=255)                                                                                      cACGCACCGTTGACAGTGCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 <  1:43336/76‑1 (MQ=255)                                                                                      cACGCACCGTTGACAGTGCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 >  2:43336/1‑76 (MQ=255)                                                                                      cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                  <  1:438247/139‑1 (MQ=255)                                                                                              gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                          >  2:55248/1‑107 (MQ=255)                                                                                              gTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                          <  1:55248/107‑1 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                               >  2:692659/1‑78 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCgg                                                                               <  1:692659/78‑1 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAg                                                       >  2:703527/1‑102 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAg                                                       <  1:703527/102‑1 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCATTCgg                                                                               >  2:655933/1‑78 (MQ=255)                                                                                                      tACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCATTCgg                                                                               <  1:655933/78‑1 (MQ=255)                                                                                                       aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCggg                                                                              >  2:429479/1‑78 (MQ=255)                                                                                                       aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCggg                                                                              <  1:429479/78‑1 (MQ=255)                                                                                                       aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGc                 >  2:64876/1‑139 (MQ=255)                                                                                                         cATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcgc                                                                                                                      <  2:459431/36‑1 (MQ=255)                                                                                                         cATCGGGTACTCGTCGGTGAGTTTGTTGATTGgcgc                                                                                                                      >  1:459431/1‑36 (MQ=37)                                                                                                         cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGc                                                                                                           <  1:687782/47‑1 (MQ=255)                                                                                                         cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGc                                                                                                           >  2:687782/1‑47 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 >  2:254151/1‑55 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTgaagaa                                                                                                 <  1:254151/55‑1 (MQ=255)                                                                                                           tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg             <  2:525118/139‑1 (MQ=255)                                                                                                              ggTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCa          >  1:721289/1‑139 (MQ=255)                                                                                                                tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    <  2:146799/79‑1 (MQ=255)                                                                                                                tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    >  1:146799/1‑79 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCTCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    <  2:752175/77‑1 (MQ=255)                                                                                                                  cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                    >  1:752175/1‑77 (MQ=255)                                                                                                                      tcgGTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtgaatga  >  1:760227/1‑139 (MQ=255)                                                                                                                        ggTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                >  1:87784/1‑91 (MQ=255)                                                                                                                        ggTGAGTTTGTTGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                <  2:87784/91‑1 (MQ=255)                                                                                                                                tGTTGATTGGCGCTACCCAGTCGCAGGTGCAGAACTGCGCCAGACCGTTCg                                                                                >  2:222092/1‑51 (MQ=255)                                                                                                                                   |                                                                                                                              ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGA  >  NC_000913/4297538‑4297793

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