Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,064	C→T	A253A (GCG→GCA)	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,064	0	C	T	100.0%	171.1 / NA	50	A253A (GCG→GCA)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/17);  total (33/17)

AATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCACCGCGACC  >  NC_000913/4084933‑4085198                                                                                                                                    |                                                                                                                                       aaTTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGtt                                                                                                                                >  2:755268/1‑140 (MQ=255)      ttCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGaa                                                                                                                           >  1:613129/1‑140 (MQ=255)        cATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAACGAGGATGGATCACAAGCAGCTTTGCGCCGTTGTGAAtt                                                                                                                         >  2:293574/1‑140 (MQ=255)                       aCAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCa                                                                                                               >  2:739279/1‑135 (MQ=255)                       aCAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCa                                                                                                               <  1:739279/135‑1 (MQ=255)                       aCAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGcc                                                                                                          >  1:29076/1‑140 (MQ=255)                       aCAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGcc                                                                                                          >  2:552339/1‑140 (MQ=255)                       aCAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGcc                                                                                                          >  2:241477/1‑140 (MQ=255)                            aaTACGCCTGACAGGAAACCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCg                                                                                                     >  1:1010968/1‑140 (MQ=255)                              tACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGa                                                                                                   >  2:963203/1‑140 (MQ=255)                               aCGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGaa                                                                                                  >  2:719674/1‑140 (MQ=255)                               aCGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGaa                                                                                                  >  1:512924/1‑140 (MQ=255)                               aCGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGaa                                                                                                  >  1:518778/1‑140 (MQ=255)                                   cTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCg                                                                                              >  2:562991/1‑140 (MQ=255)                                             aGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGGCGTTGTGAGTTTGGGCCTCCATCGCGCAGCGGAACCCGACCGGGTGAg                                                                                    >  2:269974/1‑140 (MQ=255)                                                  tGTCAGTACCGGAACGAATAGGGGCATAGTTGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCa                                                                               >  1:830554/1‑140 (MQ=255)                                                  tGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCa                                                                               >  1:175692/1‑140 (MQ=255)                                                    tCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgc                                                                             >  1:648554/1‑140 (MQ=255)                                                     cAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGGATTTTGGCTTCCGTCGCCCAGCGGAGCCCGACCGGGTGAGCTTCAgcg                                                                            >  2:386776/1‑140 (MQ=255)                                                      aGTACCGGAACGAATTGGGGCATAGGAGTGAGCTACCACCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcgg                                                                           <  1:293574/140‑1 (MQ=38)                                                       gTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggc                                                                          >  2:952826/1‑140 (MQ=255)                                                         aCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAGTCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGt                                                                        >  2:637015/1‑140 (MQ=255)                                                           cGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTa                                                                      <  1:562991/140‑1 (MQ=255)                                                             gAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTAcc                                                                    >  2:832935/1‑140 (MQ=255)                                                              aaCGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGAGGACCCCGACCGTGTGAACTTCAGCGGCGTTACCg                                                                   >  2:399543/1‑140 (MQ=255)                                                                     aGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATca                                                            <  2:175692/140‑1 (MQ=255)                                                                     aGGGGCATAGTAGTCACTTCCCCCCGCCGTACGAGTAAACCGAGGATCCATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATca                                                            <  2:456540/140‑1 (MQ=37)                                                                                 gtCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGt                                                >  2:775910/1‑140 (MQ=255)                                                                                 gtCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGt                                                >  1:927595/1‑140 (MQ=255)                                                                                  tCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGtt                                               >  1:847910/1‑140 (MQ=255)                                                                                    aGCTACCGCCGCCGTCCGCGTATAGGGAGCATCGATCCCAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTgg                                             <  1:260256/140‑1 (MQ=37)                                                                                     gCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGc                                            <  1:963203/140‑1 (MQ=255)                                                                                      cTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGGGTTACCGCTCATCGCCACGACGAGGTTGGCg                                           >  2:296315/1‑140 (MQ=255)                                                                                      cTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCg                                           >  2:618117/1‑140 (MQ=255)                                                                                      cTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCg                                           >  2:692332/1‑140 (MQ=255)                                                                                      cTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCg                                           >  2:915634/1‑140 (MQ=255)                                                                                         ccgccgccgTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTc                                        <  2:847910/140‑1 (MQ=255)                                                                                         ccgccgccgTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTc                                        >  1:638383/1‑140 (MQ=255)                                                                                         ccgccgccgTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTc                                        >  1:786374/1‑140 (MQ=255)                                                                                                cgTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGt                                 <  1:755268/140‑1 (MQ=255)                                                                                                cgTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGt                                 <  1:997818/140‑1 (MQ=255)                                                                                                 gTACGCGTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTc                                <  1:915634/140‑1 (MQ=255)                                                                                                     gcgTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAcc                            <  2:660519/140‑1 (MQ=255)                                                                                                      cgTAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAccc                           <  1:241477/140‑1 (MQ=255)                                                                                                      cgTAAAGCGAGGATCGATCACAATCAGATTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAccc                           <  2:951302/140‑1 (MQ=255)                                                                                                        tAAAGCGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAg                         >  1:669301/1‑140 (MQ=255)                                                                                                             cGAGGATCGATCACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGtggtt                    <  2:648554/140‑1 (MQ=255)                                                                                                                     gatcACAATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCg            <  2:613129/140‑1 (MQ=255)                                                                                                                         acaATCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCAcc        <  2:786374/140‑1 (MQ=255)                                                                                                                          caATCAGCTTTGCGCAGTTGTGCATATTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCACcg       <  1:468397/140‑1 (MQ=255)                                                                                                                             tCAGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCACCGCGa    >  2:94878/1‑140 (MQ=255)                                                                                                                               aGCTTTGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCACCGCGAcc  <  1:618117/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       AATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGTCATCGCACCGCGACC  >  NC_000913/4084933‑4085198

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