Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,473,900	C→T	G52R (GGA→AGA)	bdcA ←	c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,473,900	0	C	T	100.0%	68.0 / NA	20	G52R (GGA→AGA)	bdcA	c‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/12);  total (8/12)

TCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAA  >  NC_000913/4473654‑4474101                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          tCACGGAGCCGATGATTAAGATGCGCCCGCCTGCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACATGCATTTACCACCAGGATATCCAATGCGACGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGtc                                                                                                                                                                                                       >  1:31522/1‑251 (MQ=255)                       gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCg                                                                                                                                                                                  <  1:182430/250‑1 (MQ=255)                                          tGCCGGGCGGCTTCAACAGAGGCATGATCTGGAGCATGACTATTGATTTTGAAAACGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGACATCCAATGCGCCGCTCTTACGAACGACATCACTGACAGCGTCTCTGTCAGCCCTCCCTGTGAATCCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                               >  2:122063/1‑250 (MQ=255)                                           gCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAg                                                                                                                                                             <  1:220954/251‑1 (MQ=255)                                                         aCAGAGGCATGATAAGGAGCATAAATATTGATTTTGCAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCAAAAGACGCCAATTCCTGCATTTACCAACAGGATATCCAATGCGCTGCTATTACGAAAGACCTCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATAATGCTGTCGCTCTAGTCTCTTGTGACAGGCGTTTAGCAGCATCTTTCGACCCCGCATAGGTGAATCGTACAtt                                                                                                                                               >  2:164031/1‑251 (MQ=255)                                                              ggCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTGTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGccc                                                                                                                                          <  2:46988/251‑1 (MQ=255)                                                                             aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                                                                                                           >  1:5018/1‑251 (MQ=255)                                                                                 atatTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgacg                                                                                                                       <  2:198663/251‑1 (MQ=255)                                                                                      gATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCCGGGCCTCGCCACAGACGCCAATACATGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACCGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATCGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACg                                                                                                                   >  1:231298/1‑250 (MQ=255)                                                                                                                           tCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATga                                                                              <  2:200583/250‑1 (MQ=255)                                                                                                                                      gCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAc                                                                  <  1:233438/251‑1 (MQ=255)                                                                                                                                        cAAAGACGCCAATACCTGCATTTACCACCAGCATATCCAATGAGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                >  2:195223/1‑251 (MQ=255)                                                                                                                                         aaaGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGGACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                               <  2:220526/251‑1 (MQ=255)                                                                                                                                              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaaa                                                          <  2:171529/251‑1 (MQ=255)                                                                                                                                                                   ccaGGAGATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCTATTCTCTTGTGCCAGGCGTTGATCCGCAGCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAGACGACGTACGATAGCGGCATCGATACCACGACCGTCACCGAGGATGAGAACGGTCTTAGCTGGAAAAACGCCCaattggagttcttga                                     <  2:236900/251‑15 (MQ=255)                                                                                                                                                                    caGGATATCCAATGCGACGCTCTTTCGAACGACCTCCATGACAGCGTCTCTGTCAGCAATATCTGTGAATACTGATGTCGCTCTAGTCTCTTGTGCCAGGAGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCACATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAAGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt                                    >  2:229851/1‑251 (MQ=255)                                                                                                                                                                       gATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATtta                                 <  1:118435/251‑1 (MQ=255)                                                                                                                                                                                 gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACa                       <  2:187469/251‑1 (MQ=255)                                                                                                                                                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCCCCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACACAACTg        >  2:195954/1‑251 (MQ=255)                                                                                                                                                                                                       cAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTGGATTTATTATGTAACATGCATTACAAAACTGTTTTaa  <  1:32360/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                          TCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAA  >  NC_000913/4473654‑4474101

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