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%	117.5 / NA	34	N307H (AAC→CAC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (17/17);  total (17/17)

TACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑GGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAG  >  NC_000913/4084825‑4084978                                                                                 |                                                                           tACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑g                                                                           <  1:757403/81‑1 (MQ=255)   cTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gg                                                                          <  1:757404/80‑1 (MQ=255)      gCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggta                                                                       <  2:757405/80‑4 (MQ=255)      gCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggt                                                                        <  2:757404/79‑3 (MQ=255)       cGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggtat                                                                      <  2:757406/80‑5 (MQ=255)       cGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggta                                                                       >  2:757450/1‑76 (MQ=255)       cGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑gggta                                                                       >  1:757457/1‑76 (MQ=255)          ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGc                                                                   >  1:757461/1‑81 (MQ=255)             tCCGCGTCATAGCCGGTGAGCAGGCCATCTACAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCtt                                                                 <  2:757408/80‑1 (MQ=255)              ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTc                                                                <  1:757410/80‑1 (MQ=255)              ccGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTc                                                                <  1:757411/80‑1 (MQ=255)                 cgtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGgtg                                                            <  2:757411/81‑1 (MQ=255)                     gtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTAtt                                                        <  2:757412/81‑1 (MQ=255)                     gtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTAtt                                                        >  1:757468/1‑81 (MQ=255)                      tAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTc                                                       <  2:757413/81‑1 (MQ=255)                         ccGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTcgcg                                                    <  1:757415/81‑1 (MQ=255)                           ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGt                                                  >  2:757465/1‑81 (MQ=255)                               aaCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGaa                                              >  2:757453/1‑81 (MQ=255)                                  aGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAttt                                           <  1:757416/81‑1 (MQ=255)                                  aGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAttt                                           >  1:757456/1‑81 (MQ=255)                                     ccATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAAttttt                                         >  2:757459/1‑80 (MQ=255)                                     ccATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTc                                        <  2:757416/81‑1 (MQ=255)                                            cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGtt                                  <  2:757418/80‑1 (MQ=255)                                            cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGt                                   <  1:757418/79‑1 (MQ=255)                                                        ccTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAcgg                      >  1:757465/1‑78 (MQ=255)                                                        ccTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAg                      >  2:757462/1‑80 (MQ=255)                                                         cTCACGCACGATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGc                     <  2:757419/80‑1 (MQ=255)                                                           cacgcacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGcta                   >  2:757457/1‑78 (MQ=255)                                                            acgcacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATa                 >  2:757474/1‑81 (MQ=255)                                                              gcacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACg               >  1:757459/1‑81 (MQ=255)                                                               cacgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGc              >  2:757477/1‑81 (MQ=255)                                                                acgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGcc             <  2:757421/81‑1 (MQ=255)                                                                 cgATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt            <  2:757422/81‑1 (MQ=255)                                                                  gATCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTg           <  2:757423/81‑1 (MQ=255)                                                                   aTCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGa          >  2:757446/1‑81 (MQ=255)                                                                   aTCAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGa          >  2:757478/1‑81 (MQ=255)                                                                     cAGGCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGa          >  1:757473/1‑79 (MQ=255)                                                                        gCTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAgg      >  1:757467/1‑80 (MQ=255)                                                                         cTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGaa    <  2:757424/81‑1 (MQ=255)                                                                         cTGGCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGa     >  1:757470/1‑80 (MQ=255)                                                                            gCGTTGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAg  <  1:757426/80‑1 (MQ=255)                                                                                 |                                                                           TACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTT‑GGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAG  >  NC_000913/4084825‑4084978

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