Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	637,742	G→C	intergenic (‑124/+85)	ybdO ← / ← dsbG	putative DNA‑binding transcriptional regulator/thiol:disulfide interchange protein, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	637,742	0	G	C	100.0%	186.4 / NA	48	intergenic (‑124/+85)	ybdO/dsbG	putative DNA‑binding transcriptional regulator/thiol:disulfide interchange protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (23/25);  total (23/25)

GATTCAGCAGCTTTGCTAATGCTTAAATGCTGGTAAATACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACAGCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTATTTTCCTTACTCATGTAATAGATAGCCGGCGTGACATT  >  NC_000913/637521‑637931                                                                                                                                                                                                                              |                                                                                                                                                                                              gATTCAGCAGCTTTGCTAATGCTTAAATGCTGGTAAATACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCagag                                                                                                                                                                                <  1:367413/237‑1 (MQ=255)  aTTCAGCAGCTTTGCTAATGCTTAAATGCTGGTAAATACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAgagcga                                                                                                                                                                             >  2:111372/1‑239 (MQ=255)            tttGCTAATGCTTAAATGCTGGTAAATACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATg                                                                                                                                                                   <  1:111372/239‑1 (MQ=255)                                     aTACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAg                                                                                                                                          <  1:243871/239‑1 (MQ=255)                                          cTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGata                                                                                                                                     <  1:772677/239‑1 (MQ=255)                                           tCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGatat                                                                                                                                    >  1:290924/1‑239 (MQ=255)                                           tCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTGCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGatat                                                                                                                                    >  2:532067/1‑239 (MQ=255)                                                          agaTTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAccatccat                                                                                                                     >  2:433952/1‑239 (MQ=255)                                                                                      agagATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATaa                                                                                          <  1:433192/238‑1 (MQ=255)                                                                                            tGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCaa                                                                                                                             >  2:600773/1‑197 (MQ=255)                                                                                            tGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCaa                                                                                                                             <  1:600773/197‑1 (MQ=255)                                                                                             ggCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTa                                                                                  >  1:702953/1‑239 (MQ=255)                                                                                                  tACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATg                                                                                                                                                                   <  2:569496/153‑1 (MQ=255)                                                                                                  tACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATg                                                                                                                                                                   >  1:569496/1‑153 (MQ=255)                                                                                                         ctCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGaa                                                                                                                                                <  1:232238/165‑1 (MQ=255)                                                                                                         ctCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGaa                                                                                                                                                >  2:232238/1‑165 (MQ=255)                                                                                                                 atatTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCt                                                                                                            >  2:465056/1‑193 (MQ=255)                                                                                                                 atatTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCt                                                                                                            <  1:465056/193‑1 (MQ=255)                                                                                                                      tttCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCa                                                         <  1:433952/239‑1 (MQ=255)                                                                                                                          ttcCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGc                                                     <  2:702953/239‑1 (MQ=255)                                                                                                                                           ttAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGCTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTAttt                                     >  2:279929/1‑238 (MQ=255)                                                                                                                                            tAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGc                                                                                                                                            <  1:46520/134‑1 (MQ=255)                                                                                                                                            tAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGc                                                                                                                                            >  2:46520/1‑134 (MQ=255)                                                                                                                                            tAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTAttt                                     >  1:438953/1‑237 (MQ=255)                                                                                                                                              aTCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCt                                                                                                            <  2:557810/164‑1 (MQ=255)                                                                                                                                              aTCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCt                                                                                                            >  1:557810/1‑164 (MQ=255)                                                                                                                                                        tttGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGc                                                                <  2:512144/198‑1 (MQ=255)                                                                                                                                                        tttGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGc                                                                >  1:512144/1‑198 (MQ=255)                                                                                                                                                           gCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTATTTTCCTTACTCATGTAata                    <  2:290924/239‑1 (MQ=255)                                                                                                                                                                         aTTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCTACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGtt                                                                                                                                >  2:860387/1‑117 (MQ=255)                                                                                                                                                                         aTTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGtt                                                                                                                                <  1:860387/117‑1 (MQ=255)                                                                                                                                                                         aTTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAcc                                                                                                                           <  2:542450/122‑1 (MQ=255)                                                                                                                                                                         aTTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAcc                                                                                                                           >  1:542450/1‑122 (MQ=255)                                                                                                                                                                             aTTAGGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTATTTTCCTTACTCATGTAATAGATAGCCGGCGTGACAtt  <  2:59213/239‑1 (MQ=255)                                                                                                                                                                                                     gATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAcc                                                                                                                           <  2:2972/94‑1 (MQ=255)                                                                                                                                                                                                     gATTAAGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAcc                                                                                                                           >  1:2972/1‑94 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGc                                                                                                                                            >  2:389423/1‑72 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGc                                                                                                                                            <  1:389423/72‑1 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCaa                                                                                                                             >  1:829823/1‑87 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCaa                                                                                                                             <  2:829823/87‑1 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAccatcc                                                                                                                       <  1:648588/93‑1 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAccatcc                                                                                                                       >  2:648588/1‑93 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAAtt                                                                                                                >  1:105330/1‑100 (MQ=255)                                                                                                                                                                                                          aGTTTTTCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAAtt                                                                                                                <  2:105330/100‑1 (MQ=255)                                                                                                                                                                                                               ttCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAccat                                                                                                                         >  2:552720/1‑86 (MQ=255)                                                                                                                                                                                                               ttCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAAccat                                                                                                                         <  1:552720/86‑1 (MQ=255)                                                                                                                                                                                                               ttCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCggg                                                                 >  1:336956/1‑142 (MQ=255)                                                                                                                                                                                                               ttCTAAAATAGCACACCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCggg                                                                 <  2:336956/142‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                              GATTCAGCAGCTTTGCTAATGCTTAAATGCTGGTAAATACACTCAAATATGACAAGAAGATTTAAGTCGAACTTTTTCAAGTCGTAGAGATTGGCCATACACTACTCCGGTTATATTTTTCTTCCATGATATATTCACTTAATCAATGTTTTTTGCGTATATATTTTTATTGATTATGTTTTTTTGCTTAACTCAGGATTAAGTTTTTCTAAAATAGCACAGCCACTTAGGGCAGAGCGAGCAATTCATGCTTTTTTAGGATGGTTGAAAAGCAGGGATATAGTTCAACCATCCATGAATTACCTTTTATTTATTCCCCATAATGATATTAAGCGTTTTCTGATCGGGCAACCCCACGGCCTGTTGTAGCGTATTTTCCTTACTCATGTAATAGATAGCCGGCGTGACATT  >  NC_000913/637521‑637931

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