Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,844,631	A→G	I41T (ATC→ACC)	hycF ←	formate hydrogenlyase subunit HycF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,844,631	0	A	G	100.0%	142.0 / NA	40	I41T (ATC→ACC)	hycF	formate hydrogenlyase subunit HycF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/21);  total (19/21)

GGGCAGACTTCTTCGCAGCGTCCACAGAAGATGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGATGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAGGTAAACATGGCGA  >  NC_000913/2844496‑2844757                                                                                                                                        |                                                                                                                               gggCAGACTTCTTCGCAGCGTCCACAGAAGATGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCa                                                                                                                            <  2:913369/140‑1 (MQ=255)      gACTTCTTCGCAGCGTCCACAGAAGATGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCActgct                                                                                                                       <  2:524753/140‑1 (MQ=255)                               aTGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACg                                                                                              >  2:201373/1‑140 (MQ=255)                               aTGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACg                                                                                              >  1:56333/1‑140 (MQ=255)                                   aGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAg                                                                                          >  2:900387/1‑140 (MQ=255)                                          aaGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTAt                                                                                   >  2:908362/1‑140 (MQ=255)                                          aaGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTAt                                                                                   <  1:596607/140‑1 (MQ=255)                                          aaGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTAt                                                                                   <  2:581866/140‑1 (MQ=255)                                               tAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAAcc                                                                              >  2:1010035/1‑140 (MQ=255)                                                 aaCTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGc                                                                            >  2:369792/1‑140 (MQ=255)                                                    tCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAAt                                                                         <  2:56333/140‑1 (MQ=255)                                                      ccAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                                                                       >  1:154307/1‑140 (MQ=255)                                                      ccAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                                                                       >  2:151631/1‑140 (MQ=255)                                                       cAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCgg                                                                      >  2:532997/1‑140 (MQ=255)                                                       cAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCgg                                                                      >  2:681408/1‑140 (MQ=255)                                                           cAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTcc                                                                  <  1:900387/140‑1 (MQ=255)                                                            aaGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                                                                 <  2:154307/140‑1 (MQ=255)                                                             aGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAg                                                                <  1:191223/140‑1 (MQ=255)                                                                 ctcCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGa                                                            <  1:151631/140‑1 (MQ=255)                                                                    cTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATaa                                                         <  1:532997/140‑1 (MQ=255)                                                                        ggCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                                                                 <  1:391231/128‑1 (MQ=255)                                                                        ggCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                                                                 >  2:391231/1‑128 (MQ=255)                                                                         gCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGa                                                    <  1:369792/140‑1 (MQ=255)                                                                            aGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGt                                                 >  1:635886/1‑140 (MQ=255)                                                                              gTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCg                                               >  2:814382/1‑140 (MQ=255)                                                                              gTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTACACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGa                                                            <  2:446145/127‑1 (MQ=255)                                                                              gTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTACACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGa                                                            >  1:446145/1‑127 (MQ=255)                                                                                      cAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCg                                       >  2:742611/1‑140 (MQ=255)                                                                                       aaCCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCgg                                      <  1:908362/140‑1 (MQ=255)                                                                                            ttAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCATGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTtg                                 >  2:822535/1‑140 (MQ=255)                                                                                                      tGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACttttt                       >  1:19508/1‑140 (MQ=255)                                                                                                         cGGGCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTTTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGa                    <  2:1031534/140‑1 (MQ=255)                                                                                                          gggCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGAt                   <  1:383154/140‑1 (MQ=255)                                                                                                          gggCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGAt                   <  2:208390/140‑1 (MQ=255)                                                                                                          gggCAGGCATTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGAt                   <  2:848540/140‑1 (MQ=255)                                                                                                               ggcaTTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGGTGACTTTTTTGATAAAgg              <  1:1010035/140‑1 (MQ=255)                                                                                                               ggcaTTGACGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAgg              >  2:872335/1‑140 (MQ=255)                                                                                                                      aCGCAGGCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAGGTAAACAt       >  2:587932/1‑140 (MQ=255)                                                                                                                           ggCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAGGTAAACATGGCGa  <  1:742611/140‑1 (MQ=255)                                                                                                                           ggCCGCGCAGCCGGTGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAGGTAAACATGGCGa  <  1:814382/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                               GGGCAGACTTCTTCGCAGCGTCCACAGAAGATGCAGTGCCCAAGATTAAACTCCCAGGCAAGCTCTCCTGTGGCGAGGTCAGTTTCAACCGTTAAGGCGTTTGACGGGCAGGCATTGACGCAGGCCGCGCAGCCGATGCACTGCTGCGGGTTCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCGTGCCGGTTTTGATGACTTTTTTGATAAAGGTAAACATGGCGA  >  NC_000913/2844496‑2844757

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