Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,347,323	A→G	S2P (TCA→CCA)	fumB ←	anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,347,323	0	A	G	100.0%	161.3 / NA	45	S2P (TCA→CCA)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/22);  total (23/22)

CCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCAA  >  NC_000913/4347176‑4347469                                                                                                                                                    |                                                                                                                                                   ccAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCAt                                                                                                                                                 >  1:359428/1‑151 (MQ=255)  cagcagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCAt                                                                                                                                                 >  1:72384/1‑150 (MQ=255)     cagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATa                                                                                                                                                >  2:93278/1‑148 (MQ=255)     cagGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATa                                                                                                                                                <  1:93247/148‑1 (MQ=255)      agGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCtt                                                                                                                                            >  2:667662/1‑151 (MQ=255)         gTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCCTCCa                                                                                                                                         <  2:705969/151‑1 (MQ=255)            aGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGcc                                                                                                                                      >  1:635185/1‑151 (MQ=255)                  tCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTa                                                                                                                                  <  1:726935/149‑1 (MQ=255)                  tCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAAc                                                                                                                                >  2:417095/1‑151 (MQ=255)                     ggTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAActct                                                                                                                             <  1:602470/151‑1 (MQ=255)                       ttCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGt                                                                                                                           >  1:309360/1‑151 (MQ=255)                         ccACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTAt                                                                                                                         <  2:395114/151‑1 (MQ=255)                              ttCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTg                                                                                                                      <  1:688930/149‑1 (MQ=255)                                      ggTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaa                                                                                                             >  2:164969/1‑150 (MQ=255)                                      ggTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaa                                                                                                             <  1:546846/150‑1 (MQ=255)                                       gTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaaa                                                                                                           <  2:361057/151‑1 (MQ=255)                                       gTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaaa                                                                                                           <  2:171543/151‑1 (MQ=255)                                       gTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGaaaa                                                                                                           <  1:213969/151‑1 (MQ=255)                                                  gAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcgc                                                                                                >  1:235400/1‑151 (MQ=255)                                                   aaGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCcgcg                                                                                               <  2:418745/151‑1 (MQ=255)                                                        ggCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGt                                                                                           >  1:565226/1‑150 (MQ=255)                                                                           gAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGt                                                                        >  1:266382/1‑150 (MQ=255)                                                                                   tAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTtat                                                                >  1:593516/1‑150 (MQ=255)                                                                                     gATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAg                                                             <  2:325059/151‑1 (MQ=255)                                                                                      aTAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTtata                                                               >  2:602208/1‑148 (MQ=255)                                                                                        aGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGGTTAAGGGGTTGGTTGGCATAACTTCCAGCCTGTAACTCTGGATTTGTTTTTCGAAAATGGCGGTGCGCGTAGTCACGGGCGGCACGCCAAGTGCAttttttaaga                                                           >  2:324596/1‑143 (MQ=255)                                                                                           aTTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAAccc                                                        <  2:754112/150‑1 (MQ=255)                                                                                           aTTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCg                                                       <  1:487632/151‑1 (MQ=255)                                                                                             tCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGt                                                      <  2:434162/150‑1 (MQ=255)                                                                                              cGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTa                                                     <  2:675541/150‑1 (MQ=255)                                                                                               ggTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTAc                                                    >  2:322849/1‑150 (MQ=255)                                                                                               ggTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTAc                                                    >  1:716834/1‑150 (MQ=255)                                                                                                       cTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCgcggcg                                           >  1:277026/1‑151 (MQ=255)                                                                                                           ccccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGa                                          >  2:250736/1‑148 (MQ=255)                                                                                                             ccATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATc                                     <  1:568753/151‑1 (MQ=255)                                                                                                                 cGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAg                                  <  2:460699/150‑1 (MQ=255)                                                                                                                     aaaGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGaa                             <  2:333511/151‑1 (MQ=255)                                                                                                                     aaaGGTGCCTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGa                              >  2:386224/1‑150 (MQ=255)                                                                                                                            ccTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCAc                      >  2:24/1‑151 (MQ=255)                                                                                                                            ccTGGTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCAc                      <  1:564124/151‑1 (MQ=255)                                                                                                                                gTAGATAAAGGGTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATa                  <  1:562773/151‑1 (MQ=255)                                                                                                                                         gggTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCacac         >  2:724269/1‑151 (MQ=255)                                                                                                                                          ggTTTGTTTGGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACAcc        >  1:318286/1‑151 (MQ=255)                                                                                                                                                tttgGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCa   >  2:90802/1‑150 (MQ=255)                                                                                                                                                 ttgGCATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTCCCGCGTAGTCACGGGCGGCACGCAGAGTGCATTTATAAGAGCCCGTCAATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCaa  <  1:352564/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCAGCAGGGTCAGGGCTTCTGGTTCCACTTTCAGGATGGTTTCGCCGTCGAAGTCGGCAACGCTAACGTAATCGGAAGTGAGTAGATAGTATTCGGTATTGTCTTTCCCCATCGGGAAAGGTGCCTGGTAGATAAAGGGTTTGTTTGACATAGCTTCCAGCCTGTAACTCTGTATTTGTTATTCGAAAATGGCGTGCCGCGTAGTCACGGGCGGCACGCAAAGTGCATTTATAAGAACCCGTACATCGCGGCGAAGATCCAGCCGAAGACGCACGATACGCTCACACCAATCAA  >  NC_000913/4347176‑4347469

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