Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,026,449	(C)5→6	coding (1569/1695 nt)	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,444	1	.	C	100.0%	88.2 / NA	25	V525G (GTA→GGA)	dgcQ	putative diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/11);  total (14/11)

TTCTCCGCCCGTTAAGCGTTATCGCTCGCGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTA‑CCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAA  >  NC_000913/2026312‑2026567                                                                                                                                      |                                                                                                                            ttCTCCGCCCGTTAAACGTTATCGCTCGGGAATACCCGATTGCGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTa                                                                                                                              <  2:249001/133‑1 (MQ=255)        cccGTTAAGCGTTATCGCTCGCGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTA‑cccccca                                                                                                                      >  1:10327/1‑131 (MQ=255)              aaGCGTTATCGCTCGCGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAgg                                                                                                                 <  2:463982/133‑1 (MQ=255)                     aTCGCTCGCGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGAt                                                                                                          >  1:627553/1‑133 (MQ=255)                         ctcgcGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGt                                                                                                      >  2:406693/1‑133 (MQ=255)                         ctcgcGAATACCCGATTACGCCCACCCTGTTTCGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGt                                                                                                      <  2:203031/133‑1 (MQ=255)                            gcgAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCg                                                                                                       <  1:1230/129‑1 (MQ=255)                            gcgAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCg                                                                                                       >  2:1230/1‑129 (MQ=255)                                            gCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGa                                                                                   >  2:468705/1‑133 (MQ=255)                                                   cTGTTTAGCGAGATATAGACGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACAt                                                                            <  1:171251/133‑1 (MQ=255)                                                      tttAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAAGCACTGATGCGTGTCGACGTACTCTGGGCGATCAACAtctc                                                                         >  2:177005/1‑133 (MQ=255)                                                            gagaTAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTc                                                                   >  2:441710/1‑133 (MQ=255)                                                              gaTAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCAt                                                                 <  1:454439/133‑1 (MQ=255)                                                                tAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCAttt                                                               >  1:574386/1‑133 (MQ=255)                                                                   aGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAAt                                                            <  1:406693/133‑1 (MQ=255)                                                                     ccGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAAttt                                                          <  1:468705/133‑1 (MQ=255)                                                                     ccGACGGTCGCCCAGCGACTGGAGTTGTCCAACATCATAATGACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAAttt                                                          <  2:176761/133‑1 (MQ=255)                                                                         cGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGt                                                      <  2:574386/133‑1 (MQ=255)                                                                         cGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGACGAGCAAGATCTCTTTGTCATATAATTTCAGt                                                      >  2:206774/1‑133 (MQ=255)                                                                                                ttCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGc                               >  2:31760/1‑133 (MQ=255)                                                                                                ttCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGc                               >  2:629473/1‑133 (MQ=255)                                                                                                 tCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCa                              <  1:629473/133‑1 (MQ=255)                                                                                                  cAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAg                             >  1:386733/1‑133 (MQ=255)                                                                                                                 cGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTcg              >  1:510095/1‑133 (MQ=255)                                                                                                                 cGGTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTACTGCGACTTCCGCATCCTCCGACAGACTcg              >  1:186767/1‑133 (MQ=255)                                                                                                                   gTTTCCTCGCTGCTGCTTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTcgcg            <  2:331117/133‑1 (MQ=255)                                                                                                                             tgctgctTACCCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGaa  >  2:389902/1‑133 (MQ=255)                                                                                                                                      |                                                                                                                            TTCTCCGCCCGTTAAGCGTTATCGCTCGCGAATACCCGATTACGCCCAGCCTGTTTAGCGAGATAAAGCCGACGGTCGGCCAGTGACTGGAGTTGTTCAAAATCATAATCACCGGTTTCCTCGCTGCTGCTTA‑CCCCCAGCGAGGCACTGATGCGTATCGTCGTACTCTTGGCGATCAACATCTCTTTTTCATTTAATTTCAGTCGAATACGTTCTGCGACCTCCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAA  >  NC_000913/2026312‑2026567

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