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%	152.1 / NA	43	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 (21/22);  total (21/22)

ATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATA  >  NC_000913/4350971‑4351220                                                                                                                           |                                                                                                                                aTAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                                 <  2:369549/139‑1 (MQ=255)     aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                             >  2:606810/1‑139 (MQ=255)               ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGGCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAc                                                                                                                  <  2:726276/124‑1 (MQ=255)               ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTACAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAc                                                                                                                  >  1:726276/1‑124 (MQ=255)                 ccTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCg                                                                                                 >  1:900803/1‑139 (MQ=255)                   ttgCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  2:10482/137‑1 (MQ=255)                   tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  1:606810/139‑1 (MQ=255)                          gCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCAtt                                                                                        >  2:1126556/1‑139 (MQ=255)                           cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAc                                                                                                                  >  2:589370/1‑112 (MQ=255)                           cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAc                                                                                                                  <  1:589370/112‑1 (MQ=255)                                  aGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGtt                                                                                <  2:270256/139‑1 (MQ=255)                                    cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGt                                                                                            <  1:1167189/125‑1 (MQ=255)                                    cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGt                                                                                            >  2:1167189/1‑125 (MQ=255)                                    cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCa                                                                                          <  2:229724/127‑1 (MQ=255)                                    cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCa                                                                                          >  1:229724/1‑127 (MQ=255)                                              tCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGc                                                                    >  2:633060/1‑139 (MQ=255)                                                        ccAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAg                                                          >  1:406423/1‑139 (MQ=255)                                                         cAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAgg                                                         >  1:392509/1‑139 (MQ=255)                                                                aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  1:50787/94‑1 (MQ=255)                                                                aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               >  2:50787/1‑94 (MQ=255)                                                                aaaaGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCg                                                  >  2:949231/1‑139 (MQ=255)                                                                   aGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGcca                                               <  2:1019173/139‑1 (MQ=255)                                                                                 tACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGt                                 <  1:1126556/139‑1 (MQ=255)                                                                                  aCCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGtt                                >  1:864122/1‑139 (MQ=255)                                                                                   ccTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGttt                               <  1:422827/139‑1 (MQ=255)                                                                                           aGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAtttt                       <  2:864122/139‑1 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAaataat                                                                                                                        <  2:1014493/40‑1 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAaataat                                                                                                                        >  1:1014493/1‑40 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAc                                                                                                                      <  1:453002/42‑1 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAc                                                                                                                      >  2:453002/1‑42 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTc                                                                                                      <  1:1030678/58‑1 (MQ=255)                                                                                             aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTc                                                                                                      >  2:1030678/1‑58 (MQ=255)                                                                                              aTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCAt                    >  1:976771/1‑139 (MQ=255)                                                                                                    gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAt                          >  2:135652/1‑127 (MQ=255)                                                                                                    gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGAt                          <  1:135652/127‑1 (MQ=255)                                                                                                     tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt             <  2:406423/139‑1 (MQ=255)                                                                                                         aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                                <  2:180259/52‑1 (MQ=255)                                                                                                         aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGa                                                                                                >  1:180259/1‑52 (MQ=255)                                                                                                              gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTaaa    <  1:214005/139‑1 (MQ=255)                                                                                                               tCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAAt   >  1:373891/1‑139 (MQ=255)                                                                                                                ccGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATa  <  1:633060/139‑1 (MQ=255)                                                                                                                        cATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  2:306950/38‑1 (MQ=255)                                                                                                                        cATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               >  1:306950/1‑38 (MQ=255)                                                                                                                           |                                                                                                                                ATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATA  >  NC_000913/4350971‑4351220

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