Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,093	G→A	G190G (GGC→GGT)	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,093	0	G	A	100.0%	97.6 / NA	29	G190G (GGC→GGT)	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 (13/16);  total (13/16)

ACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGG  >  NC_000913/4350975‑4351213                                                                                                                       |                                                                                                                         aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                      >  1:339648/1‑139 (MQ=255)  cTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAcc                                                                                                              <  2:346419/130‑1 (MQ=255)  cTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAcc                                                                                                              >  1:346419/1‑130 (MQ=255)   tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                          <  1:405528/133‑1 (MQ=255)   tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                          >  2:405528/1‑133 (MQ=255)     cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                          >  1:692393/1‑131 (MQ=255)     cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                          <  2:692393/131‑1 (MQ=255)     cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATg                                                                                                  >  1:272230/1‑139 (MQ=255)             ccTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCg                                                                                          >  2:190268/1‑139 (MQ=255)               tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                        <  1:190268/139‑1 (MQ=255)               tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                        <  2:339648/139‑1 (MQ=255)               tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                        <  2:272230/139‑1 (MQ=255)                    ttGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCa                                                                                   <  1:473339/139‑1 (MQ=255)                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAataata                                                                                                                >  1:234274/1‑107 (MQ=255)                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAataata                                                                                                                <  2:234274/107‑1 (MQ=255)                       cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTg                                                                                <  2:146830/139‑1 (MQ=255)                         cAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGAt                                                                              <  1:685884/139‑1 (MQ=255)                            aCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAaataa                                                                                                                  >  2:12501/1‑100 (MQ=255)                            aCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAaataa                                                                                                                  <  1:12501/100‑1 (MQ=255)                                cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTgg                                                                       >  1:261659/1‑139 (MQ=255)                                                    ccAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAg                                                   >  2:745735/1‑139 (MQ=255)                                                       ggCCGTAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                            >  1:9282/1‑79 (MQ=255)                                                       ggCCGTAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                            <  2:9282/79‑1 (MQ=255)                                                               aGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGcca                                        <  1:490817/139‑1 (MQ=255)                                                                               ccTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGttt                        <  2:93327/139‑1 (MQ=255)                                                                                    aCCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAt                   >  2:604848/1‑139 (MQ=255)                                                                                                     aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgggg  <  2:140746/139‑1 (MQ=255)                                                                                                          gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                        <  1:710083/48‑1 (MQ=255)                                                                                                          gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                        >  2:710083/1‑48 (MQ=255)                                                                                                                       |                                                                                                                         ACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGG  >  NC_000913/4350975‑4351213

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