Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,083,951	A→G	G624G (GGT→GGC)	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,083,951	0	A	G	100.0%	189.8 / NA	55	G624G (GGT→GGC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/26);  total (29/26)

CCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACG  >  NC_000913/4083820‑4084086                                                                                                                                    |                                                                                                                                        cccTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTtctc                                                                                                                                 <  1:661135/140‑1 (MQ=255)           cATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTGCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAg                                                                                                                       >  2:726590/1‑140 (MQ=255)              aCATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAg                                                                                                                    >  2:18149/1‑140 (MQ=255)               cATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAgg                                                                                                                   >  1:61390/1‑140 (MQ=255)                aTCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGc                                                                                                                      <  2:452881/136‑1 (MQ=255)                aTCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGc                                                                                                                      >  1:452881/1‑136 (MQ=255)                aTCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGt                                                                                                                  <  1:699610/140‑1 (MQ=255)                aTCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGt                                                                                                                  >  1:410976/1‑140 (MQ=255)                     gcgcAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATg                                                                                                             >  1:678880/1‑140 (MQ=255)                     gcgcAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATg                                                                                                             >  1:619698/1‑140 (MQ=255)                     gcgcAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATg                                                                                                             >  1:45841/1‑140 (MQ=255)                               aGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGaa                                                                                                   <  2:538432/140‑1 (MQ=255)                               aGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGaa                                                                                                   <  2:704719/140‑1 (MQ=255)                                    gATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACtt                                                                                              <  1:726590/140‑1 (MQ=255)                                      tACCGGAGAGGATCTCGCCGTCAGTAAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACAAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCg                                                                                            <  2:240888/140‑1 (MQ=255)                                      tACCGGAGAGGATCTCGCCGTCAGTAAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACAAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCg                                                                                            >  1:240888/1‑140 (MQ=255)                                           gagaGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTg                                                                                       >  2:105337/1‑140 (MQ=255)                                           gagaGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTg                                                                                       <  1:333748/140‑1 (MQ=255)                                           gagaGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTg                                                                                       >  1:652738/1‑140 (MQ=255)                                           gagaGGATCTCGCCGTCAGTCAACGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTg                                                                                       >  1:804201/1‑140 (MQ=255)                                              agGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGaa                                                                                                   <  2:144522/125‑1 (MQ=255)                                              agGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGaa                                                                                                   >  1:144522/1‑125 (MQ=255)                                              agGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGAt                                                                                    >  1:351146/1‑140 (MQ=255)                                                gATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCt                                                                                  >  2:694173/1‑140 (MQ=255)                                                  tctcGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTc                                                                                >  2:543703/1‑140 (MQ=255)                                                  tctcGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTc                                                                                <  1:105337/140‑1 (MQ=255)                                                   ctcGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTc                                                                                         <  1:326659/130‑1 (MQ=255)                                                   ctcGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTc                                                                                         >  2:326659/1‑130 (MQ=255)                                                          tcagtcagCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTc                                                                        >  2:196375/1‑140 (MQ=255)                                                            agtcagCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCaa                                                                      <  2:619698/140‑1 (MQ=255)                                                            agtcagCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCaa                                                                      >  2:408804/1‑140 (MQ=255)                                                              tcagCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACt                                                                    <  2:410976/140‑1 (MQ=255)                                                                   gcAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGtt                                                               <  1:196375/140‑1 (MQ=255)                                                                        ccccGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCt                                                          <  2:678880/140‑1 (MQ=255)                                                                        ccccGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCt                                                          <  1:199118/140‑1 (MQ=255)                                                                        ccccGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCt                                                          <  2:45841/140‑1 (MQ=255)                                                                                 gTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGt                                                 <  2:351146/140‑1 (MQ=255)                                                                                  tCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGtt                                                >  1:792835/1‑140 (MQ=255)                                                                                   ccGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGtt                                                               <  1:432684/124‑1 (MQ=255)                                                                                   ccGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGtt                                                               >  2:432684/1‑124 (MQ=255)                                                                                     gCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTg                                             <  1:543703/140‑1 (MQ=255)                                                                                        ccTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCa                                          >  2:817150/1‑140 (MQ=255)                                                                                        ccTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATCGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCa                                          <  1:49302/140‑1 (MQ=255)                                                                                           ttCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGaa                                       >  1:563813/1‑140 (MQ=255)                                                                                            tCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAg                                      >  2:325361/1‑140 (MQ=255)                                                                                              cAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGtt                                    <  2:151596/140‑1 (MQ=255)                                                                                                   ccACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAgg                               >  2:381472/1‑140 (MQ=255)                                                                                                   ccACTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAgg                               >  1:829732/1‑140 (MQ=255)                                                                                                      cTGTAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGttt                            >  2:740044/1‑140 (MQ=255)                                                                                                         tAACCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTttcag                         <  1:740044/140‑1 (MQ=255)                                                                                                           aCCAGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTttcagtg                       >  2:246491/1‑140 (MQ=255)                                                                                                              aGCGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTttcagtgttc                    >  2:462672/1‑140 (MQ=255)                                                                                                                cGACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTttcagtgttcag                  <  1:243890/140‑1 (MQ=255)                                                                                                                 gACCTGAGTTAACGATAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTttcagtgttcagt                 <  1:694173/140‑1 (MQ=255)                                                                                                                                tAGAGCCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACg  <  2:652738/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        CCCTGTTCGGCATACATCTTGCGCAAGCGCAGGAAGATACCGGAGAGGATCTCGCCGTCAGTCAGCGCAATCCCCGGGGCGTCCGCACCTTTCCAGTGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACG  >  NC_000913/4083820‑4084086

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