Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,346,086	A→G	I414T (ATC→ACC)	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,346,086	0	A	G	100.0%	232.4 / NA	64	I414T (ATC→ACC)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/32);  total (32/32)

CTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGAT  >  NC_000913/4345957‑4346219                                                                                                                                  |                                                                                                                                      ctgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCg                                                                                                                              >  1:1227759/1‑139 (MQ=255)  tgcCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCgg                                                                                                                             <  2:1522647/139‑1 (MQ=255)    ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGaa                                                                                                                           >  1:535058/1‑139 (MQ=255)      gccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGt                                                                                                                         >  1:457089/1‑139 (MQ=255)          tGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCtt                                                                                                                     >  2:1338365/1‑139 (MQ=255)                      agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCaa                                                                                                         <  2:901935/139‑1 (MQ=255)                       gATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAAt                                                                                                        >  2:1397196/1‑139 (MQ=255)                        aTCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATc                                                                                                       >  2:144939/1‑139 (MQ=255)                                  gAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                                                               >  2:678606/1‑137 (MQ=255)                                  gAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                                                               <  1:678606/137‑1 (MQ=255)                                    gTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGc                                                                                           <  2:1250391/139‑1 (MQ=255)                                    gTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGc                                                                                           <  2:250532/139‑1 (MQ=255)                                     tCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCt                                                                                          >  1:279902/1‑139 (MQ=255)                                         tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTg                                                                                                                                     >  2:487692/1‑92 (MQ=255)                                         tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                                        >  1:149215/1‑137 (MQ=255)                                         tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                                        <  2:149215/137‑1 (MQ=255)                                         tACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                                      <  1:1397196/139‑1 (MQ=255)                                           cGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgt                                                                                    <  2:1227759/139‑1 (MQ=255)                                                  cGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAat                                                                             >  2:1345012/1‑139 (MQ=255)                                                            gCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGTGTGATCTTTGGAGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGccc                                                                   >  1:958263/1‑139 (MQ=255)                                                            gCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGccc                                                                   <  1:1338365/139‑1 (MQ=255)                                                               aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    <  1:1556525/135‑1 (MQ=255)                                                               aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    >  2:1556525/1‑135 (MQ=255)                                                               aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACg                                                                >  1:664584/1‑139 (MQ=255)                                                                     aCCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATg                                                          >  1:1156128/1‑139 (MQ=255)                                                                      ccTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATgg                                                         >  1:619293/1‑139 (MQ=255)                                                                      ccTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATgg                                                         >  1:956453/1‑139 (MQ=255)                                                                         gATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGc                                                      <  2:99121/139‑1 (MQ=255)                                                                          aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATc                                                                                                       <  2:1604679/89‑1 (MQ=255)                                                                          aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATc                                                                                                       >  1:1604679/1‑89 (MQ=255)                                                                          aTGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGcc                                                     >  2:244931/1‑139 (MQ=255)                                                                               tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTg                                                                                                                                     >  1:126327/1‑54 (MQ=255)                                                                               tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                          >  1:1218907/1‑65 (MQ=255)                                                                               tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAg                                                                                                                          <  2:1218907/65‑1 (MQ=255)                                                                               tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCtt                                                                                         >  1:839224/1‑98 (MQ=255)                                                                               tAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCtt                                                                                         <  2:839224/98‑1 (MQ=255)                                                                                   cGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGa                                            <  2:619293/139‑1 (MQ=255)                                                                                    ggCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    >  2:622306/1‑114 (MQ=255)                                                                                    ggCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    <  1:622306/114‑1 (MQ=255)                                                                                           gTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAg                                    <  2:535058/139‑1 (MQ=255)                                                                                           gTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAg                                    <  2:457089/139‑1 (MQ=255)                                                                                            ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTAc                                                                                                                                  >  2:1506187/1‑44 (MQ=255)                                                                                            ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTAc                                                                                                                                  <  1:1506187/44‑1 (MQ=255)                                                                                            ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   <  2:279902/139‑1 (MQ=255)                                                                                                     gACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                   >  1:282030/1‑82 (MQ=255)                                                                                                     gACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCgtg                                                                                   <  2:282030/82‑1 (MQ=255)                                                                                                     gACCCGCGTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCgg                          <  1:1345012/139‑1 (MQ=255)                                                                                                          gcgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGt                                                  <  1:1584140/110‑1 (MQ=255)                                                                                                          gcgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGt                                                  >  2:1584140/1‑110 (MQ=255)                                                                                                           cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   >  1:421907/1‑124 (MQ=255)                                                                                                           cgTAGTAGATCGGGTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   <  2:421907/124‑1 (MQ=255)                                                                                                                      gggTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTggg         <  2:958263/139‑1 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAg                                    <  1:219259/110‑1 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAg                                    >  2:219259/1‑110 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   >  1:1301413/1‑111 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   <  2:1301413/111‑1 (MQ=255)                                                                                                                        gTGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCg       >  2:796036/1‑139 (MQ=255)                                                                                                                         tGATCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGa      <  2:1156128/139‑1 (MQ=255)                                                                                                                           aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    >  1:452989/1‑75 (MQ=255)                                                                                                                           aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGcc                                                                    <  2:452989/75‑1 (MQ=255)                                                                                                                           aTCTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAgg    <  1:244931/139‑1 (MQ=255)                                                                                                                             cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   <  2:1411716/106‑1 (MQ=255)                                                                                                                             cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGt                                   >  1:1411716/1‑106 (MQ=255)                                                                                                                             cTTTGGTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGAt  <  2:664584/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      CTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGAT  >  NC_000913/4345957‑4346219

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