Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,093,865	(G)5→6	coding (990/1359 nt)	dgcT →	putative diguanylate cyclase DgcT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,093,860	1	.	G	100.0%	112.1 / NA	32	W329G (TGG→GGG)	dgcT	putative diguanylate cyclase DgcT
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (14/18);  total (14/18)

ACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCT‑GGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAA  >  NC_000913/1093721‑1093994                                                                                                                                             |                                                                                                                                       acaTCGCGATCCCTTAAACAATATATTTAACCGAACTTACTTTTTTAACGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCt                                                                                                                                         <  2:189239/140‑1 (MQ=255)       cgATCCCTTAACCAAAATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCT‑gggggg                                                                                                                                  <  1:296942/140‑2 (MQ=255)               tAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATTGATATCGAGCACGTCAGAAAAGTCAACGACTCCTGGGGGGCATCCGGt                                                                                                                           >  2:225650/1‑140 (MQ=255)                      tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAATAGTCAACGACACCTGGGGGGCATCCGGTTGGCGAt                                                                                                                    >  2:414939/1‑140 (MQ=255)                      tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGAt                                                                                                                    >  1:120204/1‑140 (MQ=255)                      tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGAt                                                                                                                    >  2:698835/1‑140 (MQ=255)                      tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGAt                                                                                                                    >  1:497887/1‑140 (MQ=255)                       atatTTAACAGAAATTATTTTTTTAATGAGCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATc                                                                                                                   <  2:512570/140‑1 (MQ=255)                         atTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAg                                                                                                                 <  1:29170/140‑1 (MQ=255)                         atTTAACAGAAATTAATTTTTTCATTAACGGACAGTTCAATCAGCATCAGCCCAAAAAACGCCATATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAg                                                                                                                 <  1:225650/140‑1 (MQ=255)                                aGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGAt                                                                                                            <  2:628662/138‑1 (MQ=255)                                aGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGAt                                                                                                            >  1:628662/1‑138 (MQ=255)                                aGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATaa                                                                                                          <  1:630430/140‑1 (MQ=255)                                  aaaTTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATaaaa                                                                                                        <  1:414939/140‑1 (MQ=255)                                  aaaTTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATaaaa                                                                                                        >  1:43223/1‑140 (MQ=255)                                     ttattTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACa                                                                                                     <  2:673282/140‑1 (MQ=255)                                        tttttttAATGAACTCACAGTTTAACCAGCATCAGCCTAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAgtg                                                                                                  <  1:187630/140‑1 (MQ=255)                                              tAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAat                                                                                            <  2:497887/140‑1 (MQ=255)                                                  gAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                                        <  1:698835/140‑1 (MQ=255)                                                    aCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCAtt                                                                                      <  2:104461/140‑1 (MQ=255)                                                               aaTCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCAt                                                                           >  1:8378/1‑140 (MQ=255)                                                                  cagcatcagcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACg                                                                        <  1:372472/140‑1 (MQ=255)                                                                    gcatcagcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGAc                                                                      <  2:570512/140‑1 (MQ=255)                                                                        cagcCCAAAAAACGCTTTATTGGGTTATGTTTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGa                                                                  <  1:323123/140‑1 (MQ=255)                                                                          gcCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATg                                                                <  1:410127/140‑1 (MQ=255)                                                                                                        aTGGATATCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGt                                  >  1:473326/1‑140 (MQ=255)                                                                                                              atCGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCg                            >  2:705638/1‑140 (MQ=255)                                                                                                                cGACCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGGCCTTTTAGCGCTCGTCGGCGGGGAAGAGTTTGGCGTc                          >  1:193421/1‑140 (MQ=255)                                                                                                                  aCCACTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCCTCTGCGGCGCAGAGTTTGGCGTCtt                        >  1:230646/1‑140 (MQ=255)                                                                                                                     aCTTCAAAAAAGTCAACGACACCTGGGGGGCTTCCCGTTGGCGATCAGCTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCt                     <  2:193421/140‑1 (MQ=255)                                                                                                                      cTTCAAAAAAGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTg                    >  1:398819/1‑140 (MQ=255)                                                                                                                              aaGTCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACAt            <  2:120204/140‑1 (MQ=255)                                                                                                                                 tCAACGACACCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAGGAGATTTGCGTCTTGCTGACGGACATCGa         >  1:291304/1‑140 (MQ=255)                                                                                                                                        cacCTGGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGaa  <  2:473326/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                       ACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAAAAAAGTCAACGACACCT‑GGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAA  >  NC_000913/1093721‑1093994

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