Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,609	G→A	Y586Y (TAC→TAT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,609	0	G	A	100.0%	128.5 / NA	36	Y586Y (TAC→TAT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/24);  total (12/24)

GTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTC  >  NC_000913/4297532‑4297686                                                                              |                                                                              gTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGaatag                                                                            <  1:813089/81‑1 (MQ=255) gTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGaatag                                                                            >  1:813143/1‑81 (MQ=255)         ttGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAAc                                                                    >  1:813142/1‑81 (MQ=255)          tgtgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAAc                                                                    <  1:813090/80‑1 (MQ=255)            tgtgTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACtt                                                                  >  2:813146/4‑80 (MQ=255)            tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTc                                                                 <  1:813091/81‑1 (MQ=255)              cGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCAc                                                               <  2:813093/81‑1 (MQ=255)              cGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCAc                                                               <  2:813092/81‑1 (MQ=255)              cGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCAc                                                               <  2:813091/81‑1 (MQ=255)                    cAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCg                                                         >  2:813144/1‑81 (MQ=255)                     aGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGt                                                        <  1:813095/81‑1 (MQ=255)                          catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGAc                                                    <  1:813096/80‑1 (MQ=255)                           atcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACa                                                   <  1:813097/80‑1 (MQ=255)                              agccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTAc                                               <  1:813098/81‑1 (MQ=255)                                ccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCa                                             >  2:813147/1‑81 (MQ=255)                                  agcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATc                                           <  2:813098/81‑1 (MQ=255)                                     gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCggg                                        >  1:813160/1‑81 (MQ=255)                                      ccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGt                                       <  1:813100/81‑1 (MQ=255)                                           aGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcg                                  >  2:813151/1‑81 (MQ=255)                                            gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgt                                 <  1:813101/81‑1 (MQ=255)                                             tGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgtc                                <  1:813102/81‑1 (MQ=255)                                               ccGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCgg                              <  2:813102/81‑1 (MQ=255)                                                cGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGt                             <  1:813104/81‑1 (MQ=255)                                                 gCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGt                             >  1:813158/1‑80 (MQ=255)                                                 gCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                            <  1:813105/81‑1 (MQ=255)                                                  cacaGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGa                           <  1:813106/81‑1 (MQ=255)                                                  cacaGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGa                           <  2:813106/81‑1 (MQ=255)                                                   acaGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAg                          >  1:813151/1‑81 (MQ=255)                                                   acaGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAg                          >  1:813159/1‑81 (MQ=255)                                                    caGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGt                         <  2:813107/81‑1 (MQ=255)                                                         aCCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTc                    >  1:813165/1‑81 (MQ=255)                                                          ccGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTc                    <  1:813109/80‑1 (MQ=255)                                                           cGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGa                  >  2:813148/1‑81 (MQ=255)                                                            ggTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAt                 >  2:813154/1‑81 (MQ=255)                                                             gTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTGCTCGTCGGTGAGTTTGTCGAtt                <  1:813111/81‑1 (MQ=255)                                                             gTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAtt                <  2:813109/81‑1 (MQ=255)                                                                    aaCGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTa         <  2:813111/81‑1 (MQ=255)                                                                          aaGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGt   <  1:813113/81‑1 (MQ=255)                                                                           agaaTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTc  <  1:813114/81‑1 (MQ=255)                                                                              |                                                                              GTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTC  >  NC_000913/4297532‑4297686

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