Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,986	G→T	intergenic (‑154/‑87)	comR ← / → bhsA	DNA‑binding transcriptional repressor ComR/DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,986	0	G	T	100.0%	68.4 / NA	22	intergenic (‑154/‑87)	comR/bhsA	DNA‑binding transcriptional repressor ComR/DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/13);  total (9/13)

CATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTG  >  NC_000913/1168746‑1169222                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             cgtGGCCTTTGCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTGACACATTGTGTTGAGTAAGTTGATATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCTATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATcc                                                                                                                                                                                                                                    <  2:132659/249‑1 (MQ=255)               gcgcGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATc                                                                                                                                                                                                                      <  1:216331/251‑1 (MQ=255)                                              gcCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttg                                                                                                                                                                                       >  1:93494/1‑251 (MQ=255)                                                                 caTTGTGTTGAGTCAGTTGCCATGAGGCCTTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCTAGATCCCATAAAAATAATTGCTTTATATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                                                                                                                                                                    <  2:112038/251‑1 (MQ=255)                                                                              cAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACtatt                                                                                                                                                        <  1:228368/250‑1 (MQ=255)                                                                                gTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATg                                                                                                                                                     <  1:405383/251‑1 (MQ=255)                                                                                        ggtgCCGTTGTACCTGGTGACTTTGAATGAAAGGTTATTATAAAAATAATCACCTCAGTTCACCAGTCCAGATCCCATAAAAATGATTGCTTTCTATTTACCTGAAATTTAAAGATTTTTAAATTACTTAATGCTTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGTTCGTTAAGTAATTGCTTGCGACGTCATTCTTCTGCATAAGGACACTATTATGAAATACg                                                                                                                                              >  2:121672/3‑250 (MQ=255)                                                                                                    ccTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAATTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCAt                                                                                                                                 <  2:197072/251‑1 (MQ=255)                                                                                                       ggTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTatcg                                                                                                                                                                                                          <  1:344564/175‑1 (MQ=255)                                                                                                       ggTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTatcg                                                                                                                                                                                                          >  2:344564/1‑175 (MQ=255)                                                                                                          gACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAACCGTAAAAACCCTCATCGCTg                                                                                                                            <  2:325385/250‑1 (MQ=255)                                                                                                                              tataAAAATAATCACCTACGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGt                                                                                                       >  1:385702/1‑251 (MQ=255)                                                                                                                               ataAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCGTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTc                                                                                                      <  2:317748/251‑1 (MQ=255)                                                                                                                                                           ccAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACta                                                                                                                                                          <  1:1607/171‑1 (MQ=255)                                                                                                                                                           ccAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACta                                                                                                                                                          >  2:1607/1‑171 (MQ=255)                                                                                                                                                                    aTAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAAcc                                                                 >  1:251176/1‑250 (MQ=255)                                                                                                                                                                                           aaCTGAAATTTAAAGATTTTTAAACTAAATAATGATGGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCAGTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCACCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCTATGACATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCg                                           <  1:372846/250‑1 (MQ=255)                                                                                                                                                                                                    ttAAAGATTTTTCAATTAATTAATGATTGTTATAAAAAATATCTTTTATGTTATCCAGATCACATCTATCATTTAGATATCGATCGTTATGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGGAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGACAACACCAGAAGGCcaacacacaggcggaacaaatac                                 >  2:397144/1‑233 (MQ=255)                                                                                                                                                                                                                aaaTTAATTAATGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGCTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCggg                     <  1:57004/251‑1 (MQ=255)                                                                                                                                                                                                                 aattaattaattATTGTTATAAAAAATATCTTTTATGTGAATCAGATTACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCgggg                    <  1:44628/251‑1 (MQ=255)                                                                                                                                                                                                                          atGATTGTTATAAAAAATATCTTTTATGTGATCCAGATCACATCTATCATTTAGATATCGCTCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATCAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGCTTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACACAAAGTCGGTACAATCATTGCTAACGCGGGGACAAATCTg           >  2:184374/1‑251 (MQ=255)                                                                                                                                                                                                                                   ataAAAAATATCTTTTCTGTGATCCAGATCACATCTATCATTTCGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGACAACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGCGCTAACGCGGGGACACATCTGGGCTCGCTg  >  1:314523/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             CATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTG  >  NC_000913/1168746‑1169222

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