Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,153	G→A	S170S (AGC→AGT)	dcuS ←	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,153	0	G	A	100.0%	174.1 / NA	44	S170S (AGC→AGT)	dcuS	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/20);  total (24/20)

GGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTC  >  NC_000913/4351009‑4351293                                                                                                                                                 |                                                                                                                                             ggAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAAc                                                                                                                                        >  2:189864/1‑151 (MQ=255)  gagaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACt                                                                                                                                       >  2:418199/1‑151 (MQ=255)     aTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCaa                                                                                                                                    <  1:721362/151‑1 (MQ=255)      tttCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCa                                                                                                                                     >  2:280756/1‑149 (MQ=255)         cGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGcc                                                                                                                                >  2:471980/1‑151 (MQ=255)         cGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGcc                                                                                                                                >  2:425456/1‑151 (MQ=255)               gTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGc                                                                                                                          <  2:21673/151‑1 (MQ=255)                       ccGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGc                                                                                                                  >  1:62290/1‑151 (MQ=255)                            aaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAAtt                                                                                                             >  2:164479/1‑151 (MQ=255)                              ggATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTg                                                                                                           >  2:727473/1‑151 (MQ=255)                                      ttCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAt                                                                                                   >  2:586065/1‑151 (MQ=255)                                         aGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtttt                                                                                                <  1:686267/151‑1 (MQ=255)                                          gTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTc                                                                                               >  1:727313/1‑151 (MQ=255)                                           tACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTc                                                                                               >  1:546419/1‑150 (MQ=255)                                                tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGt                                                                                          >  2:283174/1‑150 (MQ=255)                                                     agacagCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg                                                                                    <  1:579666/146‑1 (MQ=255)                                                                 ccAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATAcgc                                                                        <  2:89497/151‑1 (MQ=255)                                                                   aaTCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGt                                                                      <  2:860672/151‑1 (MQ=255)                                                                              ccAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTgcgc                                                           <  1:730297/151‑1 (MQ=255)                                                                                       cAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAAcc                                                  <  2:246884/151‑1 (MQ=255)                                                                                        aaaTAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACcg                                                 >  1:355025/1‑151 (MQ=255)                                                                                        aaaTAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACcg                                                 >  1:495735/1‑151 (MQ=255)                                                                                          ataataCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACcgcg                                               <  2:481425/151‑1 (MQ=255)                                                                                            aataCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGa                                              <  2:464259/150‑1 (MQ=255)                                                                                             ataCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGAtt                                            <  2:827229/151‑1 (MQ=255)                                                                                               aCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGa                                          <  2:356940/151‑1 (MQ=255)                                                                                                ccGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGAt                                         >  1:384363/1‑151 (MQ=255)                                                                                                    ccAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCg                                     >  1:636434/1‑151 (MQ=255)                                                                                                      aGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGAc                                   <  2:664734/151‑1 (MQ=255)                                                                                                      aGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGAc                                   <  2:565795/151‑1 (MQ=255)                                                                                                          aaTGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAttt                               >  2:353901/1‑151 (MQ=255)                                                                                                          aaTGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAtt                                <  2:793257/150‑1 (MQ=255)                                                                                                           aTGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAtttt                              >  2:659783/1‑151 (MQ=255)                                                                                                                ccAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTttcttc                          <  2:894682/150‑1 (MQ=255)                                                                                                                 cAGCGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTGTGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTttctt                           <  2:768291/148‑1 (MQ=255)                                                                                                                    cGACTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCAt                     >  2:293159/1‑151 (MQ=255)                                                                                                                      aCTGTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTc                   <  1:717814/151‑1 (MQ=255)                                                                                                                         gTCATTGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAgc                <  1:586115/151‑1 (MQ=255)                                                                                                                              tGATCTGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCttt           >  2:276126/1‑151 (MQ=255)                                                                                                                                   tGTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAgg       >  1:861856/1‑150 (MQ=255)                                                                                                                                    gTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAgg       <  1:750636/149‑1 (MQ=255)                                                                                                                                    gTTGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGAt     >  1:796453/1‑151 (MQ=255)                                                                                                                                     ttGGGTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATg    >  1:228524/1‑151 (MQ=255)                                                                                                                                       gggTCACACGACTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGtc  >  2:653159/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                             GGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTC  >  NC_000913/4351009‑4351293

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