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%	150.3 / NA	44	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/23);  total (21/23)

CTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATA  >  NC_000913/4350966‑4351220                                                                                                                                 |                                                                                                                                cttTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAc                                                                                                                      <  2:207613/139‑1 (MQ=255) cttTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATAc                                                                                                                      <  1:398045/139‑1 (MQ=255)          aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                             >  2:697875/1‑139 (MQ=255)          aCTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGAt                                                                                                             >  1:145170/1‑139 (MQ=255)            tGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATaa                                                                                                           >  2:642751/1‑139 (MQ=255)            ggCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTCGGGTTCCAGGCCGAAAAGGA‑TTTTTTGCAGTACCTTAACCAGAATGCATGTACCAATCAGTCCGACCAACATACCAAATAATACCGACCAGATac                                                                                                           <  2:732086/138‑2 (MQ=255)              cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAAt                                                                                                          >  1:82401/1‑139 (MQ=255)              cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAAtt                                                                                                         >  2:686860/1‑138 (MQ=255)                    ggCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAg                                                                                                   >  2:704241/1‑139 (MQ=255)                        tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  2:668689/139‑1 (MQ=255)                        tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  1:697875/139‑1 (MQ=255)                        tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  1:704241/139‑1 (MQ=255)                        tGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  2:145170/139‑1 (MQ=255)                         ggCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACt                                                                                              >  2:308688/1‑139 (MQ=255)                                   aaaCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATc                                                                                    <  2:736444/139‑1 (MQ=255)                                         cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTgg                                                                              >  1:251202/1‑139 (MQ=255)                                          gTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTggg                                                                             >  1:260714/1‑139 (MQ=255)                                               gaTTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                        <  1:421087/139‑1 (MQ=255)                                               gaTTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                        <  1:308688/139‑1 (MQ=255)                                                     gTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCt                                                                   <  2:82401/139‑1 (MQ=255)                                                        gggTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAAc                                                               >  2:781854/1‑139 (MQ=255)                                                           ttCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCa                                                            >  1:695174/1‑139 (MQ=255)                                                              cAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAgg                                                         <  2:251202/139‑1 (MQ=255)                                                                ggCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGcc                                                       <  1:251109/139‑1 (MQ=255)                                                                     aaaaGGA‑TTTTTTTCCGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                        >  2:538815/1‑117 (MQ=255)                                                                     aaaaGGA‑TTTTTTTCAGTTCCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGGCCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTcaca                                                                        <  1:538815/117‑1 (MQ=255)                                                                        aGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGcca                                               <  2:404201/139‑1 (MQ=255)                                                                                ttttCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAGCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGcc                                        >  1:574799/1‑139 (MQ=255)                                                                                  ttCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCaa                                      <  1:781854/139‑1 (MQ=255)                                                                                            tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               <  2:431569/72‑1 (MQ=255)                                                                                            tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGAc                                                                                               >  1:431569/1‑72 (MQ=255)                                                                                            tAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATg                            >  1:745713/1‑139 (MQ=255)                                                                                                   aaTGCAGGTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCa                     >  1:70548/1‑139 (MQ=255)                                                                                                          gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                                 >  1:752084/1‑40 (MQ=255)                                                                                                          gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGAcc                                                                                                                 <  2:752084/40‑1 (MQ=255)                                                                                                          gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGa              >  2:180069/1‑139 (MQ=255)                                                                                                          gTGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGa              >  2:773507/1‑139 (MQ=255)                                                                                                           tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt             <  2:695174/139‑1 (MQ=255)                                                                                                           tGCCAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGAt             <  1:773507/139‑1 (MQ=255)                                                                                                             ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg           >  2:29453/1‑139 (MQ=255)                                                                                                             ccAATCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgg           >  1:659076/1‑139 (MQ=255)                                                                                                               aaTCAGTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATgggg         <  1:180069/139‑1 (MQ=255)                                                                                                                    gTCCGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTaaa    <  2:659076/139‑1 (MQ=255)                                                                                                                      ccGACCAGCATACCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATa  <  2:745713/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                CTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGA‑TTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATA  >  NC_000913/4350966‑4351220

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