Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,381,504	T→C	E272G (GAA→GGA)	frdA ←	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,504	0	T	C	100.0%	190.3 / NA	53	E272G (GAA→GGA)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/22);  total (31/22)

GCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGC  >  NC_000913/4381356‑4381647                                                                                                                                                     |                                                                                                                                                gCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTccc                                                                                                                                               >  1:1175753/1‑151 (MQ=255) gCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTccc                                                                                                                                               >  1:214240/1‑151 (MQ=255)     aGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGc                                                                                                                                           >  1:86724/1‑151 (MQ=255)          aCGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAt                                                                                                                                       >  2:959494/1‑150 (MQ=255)              aaGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGt                                                                                                                                  <  2:899161/151‑1 (MQ=255)                         aaCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGa                                                                                                                       <  2:308619/151‑1 (MQ=255)                          aCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGAt                                                                                                                      <  2:63654/151‑1 (MQ=255)                           ccACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGAt                                                                                                                      >  2:789750/1‑150 (MQ=255)                           ccACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATa                                                                                                                     <  1:672810/151‑1 (MQ=255)                            cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATaa                                                                                                                    <  2:888198/151‑1 (MQ=255)                               aTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACgg                                                                                                                 <  2:991020/151‑1 (MQ=255)                                 cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTa                                                                                                               <  2:898455/151‑1 (MQ=255)                                 cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTa                                                                                                               <  2:524803/151‑1 (MQ=255)                                   gcACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGc                                                                                                             >  1:662424/2‑151 (MQ=255)                                   ccACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGc                                                                                                             >  1:527317/1‑151 (MQ=255)                                   ccACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGAAGc                                                                                                             >  1:819546/1‑151 (MQ=255)                                     aCGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGcc                                                                                                            >  2:711767/1‑150 (MQ=255)                                          gcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                      <  2:1097332/151‑1 (MQ=255)                                          gcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                      <  2:650797/151‑1 (MQ=255)                                             tGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgtt                                                                                                   >  1:414417/1‑151 (MQ=255)                                              ggAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgttg                                                                                                  <  1:467926/151‑1 (MQ=255)                                              ggAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgttg                                                                                                  >  2:350513/1‑151 (MQ=255)                                                          gCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATAcc                                                                                       <  2:669803/150‑1 (MQ=255)                                                                aCGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTc                                                                                 <  2:731034/150‑1 (MQ=255)                                                                      ttCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATCTATTTGTTTTTCGGCTCGCCCAGCGGCGCTCCCGGGCCCATGCCGTAATCTTGCAGCTAACGGTAGCCATTTTTGGTGACCAGAATAACGGCTTCACcgc                                                                            >  1:1099942/1‑149 (MQ=255)                                                                       tCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGc                                                                         >  2:351363/1‑151 (MQ=255)                                                                               gAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCg                                                                 >  1:915683/1‑151 (MQ=255)                                                                                aaGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTGGGCCCGCCCAGCGGAGTTCCCGGCCCCATGCCGTAATCTTGCAGATAGCGGTAGCCATTTTTGTTGCCCAGAATACCGCCTTCTCCGCGGCAACCTTCgg                                                                <  1:1009678/151‑1 (MQ=255)                                                                                aaGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCgg                                                                <  2:333290/151‑1 (MQ=255)                                                                                   gCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtc                                                              >  1:881116/1‑150 (MQ=255)                                                                                      tGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcatc                                                           >  2:56984/1‑150 (MQ=255)                                                                                        agagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAgg                                                        >  1:1066000/1‑151 (MQ=255)                                                                                        agagaCTTTGGCGCGTGGACCCAGATCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAAGACGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAgg                                                        >  1:108693/1‑151 (MQ=255)                                                                                            aCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATAc                                                    <  2:32353/151‑1 (MQ=255)                                                                                                tGTCGCGTGGACCCAGTTCCATATTTTTGTTTTGCGGCTCGCCCAGCGGATTCCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCg                                                  <  1:1009155/149‑1 (MQ=255)                                                                                                     cgTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTg                                           <  2:733303/151‑1 (MQ=255)                                                                                                        ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCa                                        <  2:1178642/151‑1 (MQ=255)                                                                                                         gACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGc                                         >  1:355736/1‑149 (MQ=255)                                                                                                          aCCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGc                                         >  2:273901/1‑148 (MQ=255)                                                                                                                 tCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGtt                                <  1:711601/150‑1 (MQ=255)                                                                                                                   cATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTggg                             >  1:1092133/1‑151 (MQ=255)                                                                                                                       tatTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGa                          <  2:1176014/150‑1 (MQ=255)                                                                                                                        atTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGAt                         >  1:743224/1‑150 (MQ=255)                                                                                                                            gTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTg                     >  2:548549/1‑150 (MQ=255)                                                                                                                               tttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACg                 >  1:251925/1‑151 (MQ=255)                                                                                                                                     cTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTcc           >  1:1048758/1‑151 (MQ=255)                                                                                                                                       cGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCAt         >  1:400436/1‑151 (MQ=255)                                                                                                                                        gCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg        >  1:1026819/1‑151 (MQ=255)                                                                                                                                          ccAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc      >  2:1089605/1‑151 (MQ=255)                                                                                                                                          ccAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc      >  2:327563/1‑151 (MQ=255)                                                                                                                                             gCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCAc    <  1:505379/150‑1 (MQ=255)                                                                                                                                             gCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACg   >  2:1003138/1‑151 (MQ=255)                                                                                                                                              cGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc  >  1:492585/1‑151 (MQ=255)                                                                                                                                               ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc  >  2:533698/1‑150 (MQ=255)                                                                                                                                                gAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc  >  1:609077/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGC  >  NC_000913/4381356‑4381647

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