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%	127.7 / NA	37	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 (17/20);  total (17/20)

TCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAAT  >  NC_000913/4350965‑4351219                                                                                                                                 |                                                                                                                               tcttTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAataata                                                                                                                      <  1:655636/139‑1 (MQ=255)     tGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                                  <  2:748506/139‑1 (MQ=255)      gATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAc                                                                                                                 <  1:501842/139‑1 (MQ=255)         aGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                                  <  1:128107/135‑1 (MQ=255)         aGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGa                                                                                                                  >  2:128107/1‑135 (MQ=255)           aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                            >  2:149759/1‑139 (MQ=255)           aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                            >  1:666650/1‑139 (MQ=255)             tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATaa                                                                                                          >  2:694689/1‑139 (MQ=255)                    tGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCa                                                                                                   >  2:42039/1‑139 (MQ=255)                     ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAg                                                                                                  >  2:147773/1‑139 (MQ=255)                        cTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                               <  1:526442/139‑1 (MQ=255)                         tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                              <  2:666650/139‑1 (MQ=255)                         tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                              <  1:17385/139‑1 (MQ=255)                         tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                              <  1:149759/139‑1 (MQ=255)                         tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                              <  1:694689/139‑1 (MQ=255)                         tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGG‑CTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACGGACCAGATAATGCTCCAGCGACt                                                                                             >  1:400711/1‑139 (MQ=255)                                     aaCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCt                                                                                  >  2:613429/1‑139 (MQ=255)                                          cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTgg                                                                             >  1:146365/1‑139 (MQ=255)                                           gTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTggg                                                                            >  2:739969/1‑139 (MQ=255)                                                gaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                       <  2:164892/139‑1 (MQ=255)                                                     cGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                                  <  2:59031/139‑1 (MQ=255)                                                         gggTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCt                                                                                                      <  1:427764/99‑1 (MQ=255)                                                         gggTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCt                                                                                                      >  2:427764/1‑99 (MQ=255)                                                          ggTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACt                                                             >  2:186935/1‑139 (MQ=255)                                                                      aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                               >  1:776941/1‑93 (MQ=255)                                                                      aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                               <  2:776941/93‑1 (MQ=255)                                                                                       tACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGt                                <  1:739969/139‑1 (MQ=255)                                                                                            tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATg                           >  1:412801/1‑139 (MQ=255)                                                                                               ccAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtt                        <  1:492261/139‑1 (MQ=255)                                                                                                   aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                       >  2:648233/1‑88 (MQ=255)                                                                                                   aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                       <  1:648233/88‑1 (MQ=255)                                                                                                     tGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                                <  2:299190/45‑1 (MQ=255)                                                                                                     tGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                                >  1:299190/1‑45 (MQ=255)                                                                                                          gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATTCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGa             <  2:478602/139‑1 (MQ=255)                                                                                                           tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt            <  1:42039/139‑1 (MQ=255)                                                                                                           tGCCAATCAGTCCGACCAGCATACCAAATAAGACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt            <  2:740726/139‑1 (MQ=255)                                                                                                                     tCCGACCACCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAAt  >  2:522114/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                               TCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAAT  >  NC_000913/4350965‑4351219

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