Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,873,384	A→G	K448K (AAA→AAG)	dgcJ →	putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,873,384	0	A	G	100.0%	99.9 / NA	28	K448K (AAA→AAG)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/15);  total (13/15)

CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873248‑1873516                                                                                                                                         |                                                                                                                                     cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                   <  2:319092/140‑1 (MQ=255)   ttattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCt                                                                                                                                 >  2:255798/1‑140 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                 <  2:130151/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                 <  1:255798/140‑1 (MQ=255)                             gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAgaga                                                                                                       >  1:352436/1‑140 (MQ=255)                                          gATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTc                                                                                          <  2:320803/140‑1 (MQ=255)                                                   aTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTa                                                                                 >  1:318126/1‑140 (MQ=255)                                                   aTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTa                                                                                 >  1:4865/1‑140 (MQ=255)                                                      cGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt                                                                                                <  2:101466/122‑1 (MQ=255)                                                      cGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCtt                                                                                                >  1:101466/1‑122 (MQ=255)                                                         cTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGAt                                                                           >  2:145027/1‑140 (MQ=255)                                                                   aTGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAAc                                                                 >  2:112265/1‑140 (MQ=255)                                                                   aTGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCGTCCGGTATTTACGCGATGAAAGAAAAc                                                                 <  2:215443/140‑1 (MQ=255)                                                                      aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                              <  1:112265/140‑1 (MQ=255)                                                                           tGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTGCGCGATGAAAGAAAACGATACGtt                                                         >  1:161488/1‑140 (MQ=255)                                                                                 atTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatga                                                   <  2:318126/140‑1 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                  <  2:352436/140‑1 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                  <  1:228694/140‑1 (MQ=255)                                                                                  ttCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                  <  2:334000/140‑1 (MQ=255)                                                                                          gATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATaa                                          <  2:4865/140‑1 (MQ=255)                                                                                            ttCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAg                                        >  2:318847/1‑140 (MQ=255)                                                                                             tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                  >  2:91213/1‑129 (MQ=255)                                                                                             tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatgat                                                  <  1:91213/129‑1 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTGCGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                <  2:161488/140‑1 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                >  1:75041/1‑140 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGAt                                >  1:13622/1‑140 (MQ=255)                                                                                                               gcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtata                     <  2:75041/140‑1 (MQ=255)                                                                                                                cacaACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTatat                    <  2:13622/140‑1 (MQ=255)                                                                                                                                  tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaaa  >  2:348621/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873248‑1873516

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