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%	114.3 / NA	32	Y586Y (TAC→TAT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/14);  total (18/14)

ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTC  >  NC_000913/4297538‑4297686                                                                        |                                                                              aTTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGAcc                                                                       <  2:881103/80‑1 (MQ=255)  tttGTGCGTAGCCAGGTTCATCAGCCAGCGCCGTCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCaa                                                                     >  1:881160/1‑81 (MQ=255)      tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACtt                                                                  >  1:881150/1‑80 (MQ=255)        cGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCAc                                                               >  2:881149/1‑81 (MQ=255)          tAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCACCGAACGGCAAGAATAGTGACCAACTTCACg                                                              >  2:881152/1‑80 (MQ=255)            gCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCAc                                                           <  2:881104/81‑1 (MQ=255)             ccAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCAcc                                                          >  1:881144/1‑81 (MQ=255)              cAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCg                                                         >  1:881148/1‑81 (MQ=255)                 gTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGtt                                                       >  1:881169/1‑80 (MQ=255)                    catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGAc                                                    <  2:881105/80‑1 (MQ=255)                        agccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTa                                                <  1:881106/80‑1 (MQ=255)                           cagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCAt                                            >  2:881166/1‑81 (MQ=255)                             gcgcCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCg                                          <  1:881107/81‑1 (MQ=255)                              cgccgcCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCgg                                         >  1:881159/1‑81 (MQ=255)                               gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCggg                                        >  2:881153/1‑81 (MQ=255)                                   ccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACt                                    <  2:881109/81‑1 (MQ=255)                                   ccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACt                                    <  2:881108/81‑1 (MQ=255)                                    cAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtc                                   >  1:881155/1‑81 (MQ=255)                                      gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgt                                 <  2:881110/81‑1 (MQ=255)                                      gTGCCGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgt                                 >  2:881206/1‑81 (MQ=255)                                         ccGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgtcg                               <  2:881111/80‑1 (MQ=255)                                          cGCACAGTTACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGt                             <  2:881112/81‑1 (MQ=255)                                                  tACCGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGt                     >  1:881158/1‑81 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTc                    <  1:881115/80‑1 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTc                    <  1:881114/80‑1 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCg                   <  1:881113/81‑1 (MQ=255)                                                      ggTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAt                 >  1:881170/1‑81 (MQ=255)                                                      ggTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAt                 <  2:881116/81‑1 (MQ=255)                                                       gTCATCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAtt                >  2:881172/1‑81 (MQ=255)                                                           tCGAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgcg            >  1:881167/1‑81 (MQ=255)                                                             gAACGGCAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgcgc           >  1:881161/1‑80 (MQ=255)                                                                   cAAGAATAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTTATCGGTGAGTTTGTCGATTGGCGCTACCCAg    >  2:881171/1‑81 (MQ=255)                                                                     agaaTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTc  <  1:881117/81‑1 (MQ=255)                                                                        |                                                                              ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTC  >  NC_000913/4297538‑4297686

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