Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,641	A→G	M61T (ATG→ACG)	fdoG ←	formate dehydrogenase O subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,085,641	0	A	G	100.0%	203.4 / NA	57	M61T (ATG→ACG)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/32);  total (25/32)

GCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCC  >  NC_000913/4085505‑4085771                                                                                                                                         |                                                                                                                                   gCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCa                                                                                                                                 <  1:1012556/140‑1 (MQ=255)  cTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCaa                                                                                                                                <  1:424090/140‑1 (MQ=255)  cTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCaa                                                                                                                                <  2:254955/140‑1 (MQ=255)         aGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCca                                                                                                                         <  1:1032780/140‑1 (MQ=255)         aGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCca                                                                                                                         <  2:337850/140‑1 (MQ=255)         aGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCca                                                                                                                         <  2:736855/140‑1 (MQ=255)            ggATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAg                                                                                                                      >  1:474849/1‑140 (MQ=255)             gATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGc                                                                                                                     >  1:199428/1‑140 (MQ=255)              aTGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGcc                                                                                                                    >  2:271384/1‑140 (MQ=255)                gAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTa                                                                                                                  <  2:474849/140‑1 (MQ=255)                gAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTa                                                                                                                  <  1:911505/140‑1 (MQ=255)                   aTCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACgg                                                                                                               >  1:559561/1‑140 (MQ=255)                       accaGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACa                                                                                                           <  2:595872/140‑1 (MQ=255)                       accaGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACa                                                                                                           >  2:751304/1‑140 (MQ=255)                        ccaGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGANNCTGTACGTCAACAGCCCACAGCCTACGGAACaa                                                                                                          <  1:298453/140‑1 (MQ=255)                        ccaGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACaa                                                                                                          <  2:525139/140‑1 (MQ=255)                         caGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAAt                                                                                                         <  2:992380/140‑1 (MQ=255)                         caGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGGCACCGAGGCTGTACGTCAACAGCCTACAGCCTACGGAACAAt                                                                                                         <  2:665590/140‑1 (MQ=255)                           ggCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCca                                                                                                                         <  1:691772/122‑1 (MQ=255)                           ggCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCca                                                                                                                         >  2:691772/1‑122 (MQ=255)                            gCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAgg                                                                                                      >  1:684394/1‑140 (MQ=255)                              cAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTg                                                                                                    <  2:684394/140‑1 (MQ=255)                                gcgcCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCa                                                                                                  <  1:271384/140‑1 (MQ=255)                                      ttCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTAtt                                                                                            >  2:532598/1‑140 (MQ=255)                                       tCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTa                                                                                           <  1:1107932/140‑1 (MQ=255)                                         ggACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACg                                                                                         <  1:532598/140‑1 (MQ=255)                                               aGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTc                                                                                   >  2:1089503/1‑140 (MQ=255)                                               aGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTc                                                                                   >  2:111094/1‑140 (MQ=255)                                                 tGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTcgc                                                                                 >  2:766442/1‑140 (MQ=255)                                                        cgGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTCTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACg                                                                                         <  2:146513/125‑1 (MQ=255)                                                        cgGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTCTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACg                                                                                         >  1:146513/1‑125 (MQ=255)                                                          gTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgca                                                                        <  2:644915/140‑1 (MQ=255)                                                           ttGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgcag                                                                       <  1:1089503/140‑1 (MQ=255)                                                           ttGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgcag                                                                       >  1:870206/1‑140 (MQ=255)                                                           ttGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgca                                                                        <  1:80194/139‑1 (MQ=255)                                                           ttGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCgca                                                                        >  2:80194/1‑139 (MQ=255)                                                                     gATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATAc                                                             >  2:1133308/1‑140 (MQ=255)                                                                             tCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGt                                                     >  1:450220/1‑140 (MQ=255)                                                                                ccTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTc                                                  >  2:1139936/1‑140 (MQ=255)                                                                                ccTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTc                                                  >  1:205235/1‑140 (MQ=255)                                                                                ccTTCGATATGGAAGATAGATGCTTTGGCG‑TTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTcc                                                 >  2:777211/1‑140 (MQ=255)                                                                                 cTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTcc                                                 >  1:1040408/1‑140 (MQ=255)                                                                                 cTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAACAGTTTATACTGCCGGGTTTcc                                                 >  1:643971/1‑140 (MQ=255)                                                                                   tCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCcgc                                               <  2:199428/140‑1 (MQ=255)                                                                                    cGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCcgcg                                              <  2:1040408/140‑1 (MQ=255)                                                                                        atGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTg                                          <  2:205235/140‑1 (MQ=255)                                                                                            aaGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTAc                                      <  1:1139936/140‑1 (MQ=255)                                                                                             agatagatGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACg                                     <  1:111094/140‑1 (MQ=255)                                                                                                 agatGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTgg                                 <  2:559561/140‑1 (MQ=255)                                                                                                   atGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGt                               >  2:261439/1‑140 (MQ=255)                                                                                                     gCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTg                                  <  1:1139142/135‑1 (MQ=255)                                                                                                     gCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTg                                  >  2:1139142/1‑135 (MQ=255)                                                                                                      cTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAATAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCa                            >  1:613448/1‑140 (MQ=255)                                                                                                      cTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACg                                     <  1:32703/131‑1 (MQ=255)                                                                                                      cTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACg                                     >  2:32703/1‑131 (MQ=255)                                                                                                                     cTCCGTCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCt             <  2:343066/140‑1 (MQ=255)                                                                                                                          tCACCGAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCgt        >  1:974268/1‑140 (MQ=255)                                                                                                                               gAGGCTGTACGTCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGtgc   <  1:261439/140‑1 (MQ=255)                                                                                                                                gggCTGTACGTCAACAGCCCACACCCTCCGGAACAATAGGCGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTCCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTACCGCTGCCGTGGtgcc  <  2:178875/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   GCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATATGGAAGATAGATGCTTTGGCGTTTTTTGCTCCGTCACCGAGGCTGTACATCAACAGCCCACAGCCTACGGAACAATAGGTGCAGGTATTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCC  >  NC_000913/4085505‑4085771

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