Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,456	G→T	Q123K (CAA→AAA)	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,456	0	G	T	100.0%	141.9 / NA	41	Q123K (CAA→AAA)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/17);  total (24/17)

CGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATA  >  NC_000913/4085319‑4085592                                                                                                                                          |                                                                                                                                         cgcgGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTttct                                                                                                                                        >  2:80307/1‑140 (MQ=255) cgcgGAAGCACACAGCATCCCGGTAGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTttct                                                                                                                                        >  1:403525/1‑140 (MQ=255)                   cccGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTg                                                                                                                      >  1:537123/1‑140 (MQ=255)                        tGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCa                                                                                                                 <  1:355526/140‑1 (MQ=255)                          gagaGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACg                                                                                                               <  2:221957/140‑1 (MQ=255)                             ggCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTa                                                                                                            >  1:459587/2‑140 (MQ=255)                                       tAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAAc                                                                                                  >  2:733912/1‑140 (MQ=255)                                            gTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAg                                                                                             <  2:424586/140‑1 (MQ=255)                                               aCGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACg                                                                                          >  2:746495/1‑140 (MQ=255)                                               aCGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACg                                                                                          >  1:671730/1‑140 (MQ=255)                                                cGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACgg                                                                                         >  2:505822/1‑140 (MQ=255)                                                     tCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTc                                                                                    >  2:244986/1‑140 (MQ=255)                                                        gCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGa                                                                                 <  1:25455/140‑1 (MQ=255)                                                        gCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGa                                                                                 <  2:473567/140‑1 (MQ=255)                                                             ttGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGa                                                                            >  1:230223/1‑140 (MQ=255)                                                                        aGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATccacc                                                                 <  1:195499/140‑1 (MQ=255)                                                                          ttAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAg                                                               >  1:388677/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  2:607374/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  1:555929/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  2:12741/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  2:225392/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  2:231135/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  1:347862/1‑140 (MQ=255)                                                                               aTCGCGGTCTTCTTTCATCAGTTTAGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAg                                                          >  1:475586/1‑140 (MQ=255)                                                                                   cgGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCTGGCCAGCGcc                                                      >  2:510565/1‑140 (MQ=255)                                                                                   cgGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGcc                                                      >  2:601400/1‑140 (MQ=255)                                                                                    ggTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCt                                                     <  1:284810/140‑1 (MQ=255)                                                                                         tcttTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCgg                                                <  1:496064/140‑1 (MQ=255)                                                                                         tcttTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCgg                                                <  2:25138/140‑1 (MQ=255)                                                                                                tcaGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCTGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAgg                                         >  2:343257/1‑139 (MQ=255)                                                                                                   gTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCa                                      <  1:231135/140‑1 (MQ=255)                                                                                                     ttGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCAcc                                    >  1:198294/1‑140 (MQ=255)                                                                                                       ggCGATGCGATCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACcgc                                  <  2:658813/140‑1 (MQ=255)                                                                                                         cgatgcgatCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCgg                                >  1:572334/1‑140 (MQ=255)                                                                                                                 tCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCgg                        <  2:403525/140‑1 (MQ=255)                                                                                                                 tCAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCgg                        <  2:572334/140‑1 (MQ=255)                                                                                                                  cAAACGCCTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCggg                       <  1:80307/140‑1 (MQ=255)                                                                                                                        ccTCTTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATcc                 >  2:648157/1‑140 (MQ=255)                                                                                                                           ctTCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGa              >  1:266942/1‑140 (MQ=255)                                                                                                                            ttCCCAACTGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGAt             <  1:554087/140‑1 (MQ=255)                                                                                                                                    tGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGTTTTCGGGGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCg     <  2:293963/140‑1 (MQ=255)                                                                                                                                    tGATTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCg     <  1:192553/140‑1 (MQ=255)                                                                                                                                      aTTTTCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGat   >  1:365236/1‑140 (MQ=255)                                                                                                                                       ttttCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGata  <  2:216400/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGACAAAGTGCACCGCGGTTGACCGGGTGATCCGGATCGCCTTCGATA  >  NC_000913/4085319‑4085592

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