Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,026,558	C→T	G487D (GGC→GAC)	dgcQ ←	putative diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,026,558	0	C	T	100.0%	140.0 / NA	40	G487D (GGC→GAC)	dgcQ	putative diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

TCCTCGCTGCTGCTTACCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATC  >  NC_000913/2026429‑2026692                                                                                                                                  |                                                                                                                                       tCCTCGCTGCTGCTTACCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAAt                                                                                                                              >  1:424139/1‑140 (MQ=255)   cTCGCTGCTGCTTACCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATca                                                                                                                            >  2:344152/1‑140 (MQ=255)                       gCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCg                                                                                                        >  2:705212/1‑140 (MQ=255)                       gCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCg                                                                                                        >  2:446061/1‑140 (MQ=255)                       gCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCg                                                                                                        >  1:609887/1‑140 (MQ=255)                         gAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGAc                                                                                                      <  1:696196/140‑1 (MQ=255)                         gAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGAc                                                                                                      <  1:705212/140‑1 (MQ=255)                           ggCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGAccc                                                                                                    >  1:72536/1‑140 (MQ=255)                            gCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTcacc                                                                                                            >  2:820007/1‑131 (MQ=255)                            gCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTcacc                                                                                                            <  1:820007/131‑1 (MQ=255)                              aCTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGAcccgcc                                                                                                 <  2:172625/140‑1 (MQ=255)                                      gTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACg                                                                                         >  1:609476/1‑140 (MQ=255)                                        aTCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTc                                                                                       <  2:609887/140‑1 (MQ=255)                                         tcgtcgTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTcc                                                                                      >  1:297883/1‑140 (MQ=255)                                         tcgtcgTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTcc                                                                                      >  2:388287/1‑140 (MQ=255)                                            tcgtACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTgc                                                                                   >  1:284832/1‑140 (MQ=255)                                              gtACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTgcgc                                                                                 <  1:772560/140‑1 (MQ=255)                                                      ggCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAgg                                                                         <  1:936997/140‑1 (MQ=255)                                                      ggCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAgg                                                                         >  2:936997/1‑140 (MQ=255)                                                                    ctTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCtaattaat                                                           <  2:253056/140‑1 (MQ=255)                                                                        ttCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGTCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCgg                                                       <  2:609476/140‑1 (MQ=255)                                                                        ttCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCgg                                                       <  2:297883/140‑1 (MQ=255)                                                                               aaTTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATcc                                                         >  2:272751/1‑131 (MQ=255)                                                                               aaTTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATcc                                                         <  1:272751/131‑1 (MQ=255)                                                                               aaTTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATg                                                <  1:525095/140‑1 (MQ=255)                                                                                              cGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTc                                 <  2:72536/140‑1 (MQ=255)                                                                                                  cTGCGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTcgccc                             <  2:488814/140‑1 (MQ=255)                                                                                                     cGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTcgcccgcc                          >  2:111767/1‑140 (MQ=255)                                                                                                     cGACCTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTcgcccgcc                          >  2:805265/1‑140 (MQ=255)                                                                                                        ccTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCtga                       >  2:578845/1‑140 (MQ=255)                                                                                                         cTCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCtgat                      >  1:188397/1‑140 (MQ=255)                                                                                                          tCCGCAGCCTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCtgatg                     >  2:414493/1‑140 (MQ=255)                                                                                                                 ccTCCGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAg              <  2:188397/140‑1 (MQ=255)                                                                                                                     cGTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGt          <  1:578845/140‑1 (MQ=255)                                                                                                                      gTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAg              <  1:240640/135‑1 (MQ=255)                                                                                                                      gTCAGACTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAg              >  2:240640/1‑135 (MQ=255)                                                                                                                           aCTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTaa    >  2:687604/1‑140 (MQ=255)                                                                                                                           aCTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTaa    <  1:388287/140‑1 (MQ=255)                                                                                                                           aCTCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTaa    <  2:947205/140‑1 (MQ=255)                                                                                                                             tCGCGTCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATc  <  2:424139/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                       TCCTCGCTGCTGCTTACCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATC  >  NC_000913/2026429‑2026692

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