Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,870,885	A→G	Q167Q (CAA→CAG)	cdgI →	putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,870,885	0	A	G	100.0%	101.7 / NA	29	Q167Q (CAA→CAG)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/12);  total (17/12)

TCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAACTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCA  >  NC_000913/1870749‑1871017                                                                                                                                         |                                                                                                                                     tCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTc                                                                                                                                   <  2:563492/140‑1 (MQ=255)        aaGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTa                                                                                                                            <  2:949609/140‑1 (MQ=255)        aaGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTa                                                                                                                            <  1:700410/140‑1 (MQ=255)               gCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGttt                                                                                                                     >  2:441659/1‑140 (MQ=255)               gCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGttt                                                                                                                     <  1:441659/140‑1 (MQ=255)                          tGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTgg                                                                                                          <  2:377679/140‑1 (MQ=255)                           ggCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGgc                                                                                                         >  1:487015/1‑140 (MQ=255)                                       ttAGTTGAAACCATTGTAATTATCCATAAACCACTTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTa                                                                                             >  2:395202/1‑140 (MQ=255)                                       ttAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTa                                                                                             >  1:670257/1‑140 (MQ=255)                                           ttGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGa                                                                                         >  1:416746/1‑140 (MQ=255)                                                aCCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGa                                                                                    >  1:483703/1‑140 (MQ=255)                                                aCCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGa                                                                                    >  2:832527/1‑140 (MQ=255)                                                            aTCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTGTGTTATGGAAAAGACAACATACTTGa                                                                        >  2:858969/1‑140 (MQ=255)                                                               cATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGacaa                                                                     >  1:756111/1‑140 (MQ=255)                                                                  aaaCCAATTAACGGCAGTACACTGATCCAGACAAATTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGAcaacaa                                                                  <  2:670257/140‑1 (MQ=255)                                                                   aaCCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAAt                                                                 <  2:793292/140‑1 (MQ=255)                                                                             cGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGaaaaaaa                                                       <  2:416746/140‑1 (MQ=255)                                                                                          aTCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGtt                                          >  2:713706/1‑140 (MQ=255)                                                                                           tCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTg                                         >  2:505428/1‑140 (MQ=255)                                                                                             cAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCt                                       >  1:164622/1‑140 (MQ=255)                                                                                               gACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTgg                                     >  1:211889/1‑140 (MQ=255)                                                                                                  aaaGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGgcgc                                  >  1:201242/1‑140 (MQ=255)                                                                                                    aGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTg                                >  1:24284/1‑140 (MQ=255)                                                                                                                tttctattttctatATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTa                    >  1:834328/1‑140 (MQ=255)                                                                                                                  tctattttctatATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGt                  <  1:713706/140‑1 (MQ=255)                                                                                                                        ttctatATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTcc            <  2:834328/140‑1 (MQ=255)                                                                                                                          ctatATTTTCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTcccc          <  2:767703/140‑1 (MQ=255)                                                                                                                                tttCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTgg    >  1:375955/1‑140 (MQ=255)                                                                                                                                  tCCGCCAGCTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGca  <  2:483703/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     TCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAGACAAAGTCGAATGATGTTTCTATTTTCTATATTTTCCGCCAACTCAGTTTTATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCA  >  NC_000913/1870749‑1871017

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: ‑