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%	94.7 / NA	28	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 (15/13);  total (15/13)

CACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA  >  NC_000913/4473655‑4474136                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             cACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGtct                                                                                                                                                                                                                                         >  1:103442/1‑251 (MQ=255)        ccGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGc                                                                                                                                                                                                                                  <  1:226130/251‑1 (MQ=255)                  aGATGCGCCCGCCTTCTGACATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTTGCCAAAGACGCCAATACCTGCATTTACCTCCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCTCTATCTGTGAATCCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTa                                                                                                                                                                                                                        <  1:161084/251‑1 (MQ=255)                      gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCgg                                                                                                                                                                                                                    <  1:73005/251‑1 (MQ=255)                                      aTCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGc                                                                                                                                                                                                    >  2:75300/1‑251 (MQ=255)                                              ggcggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTg                                                                                                                                                                                             >  1:269794/1‑250 (MQ=255)                                                         cAGAGGCATGATAAGGAGCGTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAAATCCAGGGCCCCGCCATAGACGCCAATACCTGCATTTACCACCAGGATATCCGATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACGATCGGTGAATTCTGCTGTCGCTCTAGTCTCTAGTGCCAGGCGTATAGCGGCTTCTTTCTACCCCGGATAGGGGAATCGTTCATTg                                                                                                                                                                                 <  2:290024/251‑1 (MQ=255)                                                                      aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAg                                                                                                                                                                                                                       <  1:54867/200‑1 (MQ=255)                                                                      aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCAAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTCTCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAg                                                                                                                                                                                                                       >  2:54867/1‑200 (MQ=255)                                                                                              aaGGCGATCAATATCGTCGGCCTTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                                                                                            <  2:315794/251‑1 (MQ=255)                                                                                                            cgtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGGGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGCTAGCGGCACCGATACCACGACTg                                                                                                                              >  1:226620/1‑251 (MQ=255)                                                                                                                     ttaacgccAGTGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACCTCAAGGACAGCTTCTCTGTCAGCACGATCTTTGAAATCTGCTGTCGCTCTAGTCGCTTGTTCCAGGCGTTTAGCGGCATCTTTCGGCCCCGCATAGGTGAATCGTATATTGGCCCCATCCGTGACGAAACGACGTACGATATCGGCCCCGATACCACGACTGCCACCGAg                                                                                                                      <  2:122736/244‑1 (MQ=255)                                                                                                                           ccAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                  >  1:179958/1‑248 (MQ=255)                                                                                                                                 ccTCTCCAAAGACGCCAATTCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACAGCAATGACAGCGTCTCTGTCATCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACGGCCACCGAGGATGAGAACTGTc                                                                                                         <  1:245856/251‑1 (MQ=255)                                                                                                                                        aaaGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                                                                  >  1:31562/1‑251 (MQ=255)                                                                                                                                                  aTACCTGCATTTACAACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGcc                                                                                        <  1:158481/251‑1 (MQ=255)                                                                                                                                                    aCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGAACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                                                                                       >  2:197618/1‑250 (MQ=255)                                                                                                                                                       tGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                                                                                    >  1:221013/1‑250 (MQ=255)                                                                                                                                                                   caGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACATCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTATCGGCATCTTTCAACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAATCGCCCATAAGTACTCCTTGAt                                                                       <  1:125500/251‑1 (MQ=255)                                                                                                                                                                                   ccgcTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCCCTCTCTGTGAATACTGCTGTCGCTCTAGTCTCTTGAGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGCCTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTacaatgc                                                       >  2:158034/1‑245 (MQ=255)                                                                                                                                                                                   ccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTCTCTGAGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACAGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGc                                                       >  1:96700/1‑251 (MQ=255)                                                                                                                                                                                       tctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTAGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATACGGACTCCTTGATTTATTATGTAACATGCAtt                                                    >  1:259693/1‑250 (MQ=255)                                                                                                                                                                                                  aCATCAATGTCTGCGTCTATGGCAGCGCTATCTGTGAATACTGCTGTCGCTCTGGTCTCTTGTTCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAACCGTACATTGGCCCCATCGGTGACGAATCGACGTCCGATAGCGGCACCGATACCACGACTGCTACCGAGGATGAGAACTTTCTTACCTTTAAAAGCGCCCATAATGACTACTTGATTTATTATGTAACATGCATTATAAAACTGttt                                        <  2:103442/251‑1 (MQ=255)                                                                                                                                                                                                                            aCTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTa               >  1:144196/1‑250 (MQ=255)                                                                                                                                                                                                                               aTCTGTGAATACTGCTGTCGCTCTAGTCTCTTGGTCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAAt             <  2:221013/249‑1 (MQ=255)                                                                                                                                                                                                                                  tgtgAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGCCGTACGATAGCGGCACCGATAcc                                                                                                                                     >  1:189065/1‑126 (MQ=255)                                                                                                                                                                                                                                         aCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaa  >  1:132450/1‑250 (MQ=255)                                                                                                                                                                                                                                         aCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaaaa  <  1:133915/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             CACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAA  >  NC_000913/4473655‑4474136

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