Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,084,906	+G	coding (917/3051 nt)	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,084,904	1	.	G	100.0%	172.1 / NA	49	N307H (AAC→CAC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/27);  total (22/27)

GTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑GGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGA  >  NC_000913/4084824‑4084975                                                                                  |                                                                        gTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGtt                                                                          >  1:825008/1‑81 (MQ=255)  tACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑g                                                                        >  2:825023/1‑81 (MQ=255)  tACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑g                                                                        >  2:825009/1‑81 (MQ=255)        cGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggtat                                                                   <  1:824944/80‑5 (MQ=255)         gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAg                                                                  >  1:825020/1‑81 (MQ=255)           ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGc                                                                >  1:825007/1‑81 (MQ=255)                 gcgTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGgt                                                          <  1:824945/81‑1 (MQ=255)                   gtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGgtgt                                                        <  2:824946/81‑1 (MQ=255)                   gtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGgtgt                                                        <  1:824947/81‑1 (MQ=255)                    tcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTa                                                       <  2:824948/81‑1 (MQ=255)                         gCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTcg                                                   <  1:824949/80‑1 (MQ=255)                            ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGt                                               >  1:825011/1‑81 (MQ=255)                             gTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGtt                                              <  1:824950/81‑1 (MQ=255)                             gTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGtt                                              >  1:825016/1‑81 (MQ=255)                             gTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGtt                                              >  2:825029/1‑81 (MQ=255)                               gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTg                                             >  2:825037/1‑80 (MQ=255)                               gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTg                                             >  2:825021/1‑80 (MQ=255)                                aaCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGa                                            <  2:824951/80‑1 (MQ=255)                                 aCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAt                                          >  1:825006/1‑81 (MQ=255)                                  cAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAtt                                         <  2:824952/81‑1 (MQ=255)                                  cAGGCCATCCTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAtt                                         >  1:825032/1‑81 (MQ=255)                                     gCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAttttt                                      >  2:825019/1‑81 (MQ=255)                                      ccATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTc                                     <  2:824953/81‑1 (MQ=255)                                        aTCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCAt                                   >  2:825038/1‑81 (MQ=255)                                             cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTg                                 <  1:824955/78‑1 (MQ=255)                                             cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGtt                               >  2:824996/1‑80 (MQ=255)                                             cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGt                                <  1:824954/79‑1 (MQ=255)                                              aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGtt                               >  1:825026/1‑79 (MQ=255)                                              aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTc                              <  2:824956/80‑1 (MQ=255)                                               aaGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcag                            <  2:824957/81‑1 (MQ=255)                                                   cGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAgg                        <  1:824959/81‑1 (MQ=255)                                                   cGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAgg                        <  1:824958/81‑1 (MQ=255)                                                    gTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGt                       <  2:824960/81‑1 (MQ=255)                                                    gTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGt                       <  2:824961/81‑1 (MQ=255)                                                       aTCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACa                    <  1:824962/81‑1 (MQ=255)                                                       aTCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACa                    <  2:824963/81‑1 (MQ=255)                                                        tCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAg                   <  1:824964/81‑1 (MQ=255)                                                         ccTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAg                   >  2:825031/1‑80 (MQ=255)                                                          cTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGc                  <  2:824965/80‑1 (MQ=255)                                                          cTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGc                  <  2:824967/80‑1 (MQ=255)                                                          cTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCa                 <  1:824966/81‑1 (MQ=255)                                                              cgcacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATAc             <  1:824968/81‑1 (MQ=255)                                                               gcacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACg            >  2:825054/1‑81 (MQ=255)                                                                 acgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGcc          >  1:825057/1‑81 (MQ=255)                                                                 acgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGcc          <  2:824969/81‑1 (MQ=255)                                                                  cgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt         <  2:824970/81‑1 (MQ=255)                                                                   gATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt         >  1:825044/1‑80 (MQ=255)                                                                   gATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTg        <  2:824971/81‑1 (MQ=255)                                                                    aTCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGa       >  2:825034/1‑81 (MQ=255)                                                                      cAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGa       >  1:825012/1‑79 (MQ=255)                                                                       aGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAg    >  1:825027/1‑81 (MQ=255)                                                                       aGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAg    >  1:825025/1‑81 (MQ=255)                                                                           tGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGa  >  2:825039/1‑79 (MQ=255)                                                                                  |                                                                        GTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑GGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGA  >  NC_000913/4084824‑4084975

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