Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,871,516	C→T	R378C (CGT→TGT)	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,871,516	0	C	T	100.0%	60.1 / NA	18	R378C (CGT→TGT)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/12);  total (6/12)

CGGGATTGAAGTTTCCAGTAAACTGTTTGTTGTGTCTTTTCTGATTTATAACATTTTTCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATCGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAA  >  NC_000913/1871281‑1871722                                                                                                                                                                                                                                            |                                                                                                                                                                                                               cgGGATTGAAGTTTCCAGTAAACTGTTTGTTGTGTCTTTTCTGATTTATAACATTTTTCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATTGGGACATTGTGTGGGTGATAAAg                                                                                                                                                                                                 <  1:229719/251‑1 (MQ=255)           gTTTCCAGTAAACTGTTTGTTGTGTCTTTTCTGATTTATAACATTTTTCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTc                                                                                                                                                                                       <  1:96630/251‑1 (MQ=255)                         gtttgttGTGTCTTTTCTGATTTATAACCTTTTTCAGGTGCTGCAACTCTCCATCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCTGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGTGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGAtatta                                                                                                                                                                         <  1:229340/251‑1 (MQ=255)                                  gtCTTTTCTGATTTATAACATTTTTCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCaaa                                                                                                                                                                <  1:250685/251‑1 (MQ=255)                                                       ttttCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATAttt                                                                                                                                           >  1:81740/1‑251 (MQ=255)                                                          tCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTa                                                                                                                                         <  1:396115/250‑1 (MQ=255)                                                                     aaCTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAg                                                                                                                             <  1:14900/251‑1 (MQ=255)                                                                          tcCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATACTCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTg                                                                                                                         >  1:326435/1‑250 (MQ=255)                                                                           ccAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAATCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGt                                                                                                                          <  2:399082/248‑1 (MQ=255)                                                                           ccAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGNg                                                                                                                         >  2:373330/1‑249 (MQ=255)                                                                                 aaCTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTAAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGCACTTCTGCGTCATGCTGGTTGATATTAATCCGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGAGCTGGTTTCAATTGTAGATATTATACAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGttt                                                                                                                 >  1:193934/1‑251 (MQ=255)                                                                                      gCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCtt                                                                                                            <  1:39954/251‑1 (MQ=255)                                                                                                                                                    attattGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATg                                              <  1:74311/251‑1 (MQ=255)                                                                                                                                                                       gttAAGGGCTTCTGGGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCGCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGCTATTATCCCGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCtt                           >  2:16169/1‑251 (MQ=255)                                                                                                                                                                       gttAAGGACTTCTGTGTCATGCGTGTTGATATTAAGCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTTGGTGCTAAAGTACTGGTTTCAATTGTCAATATTATCCAGCAAAGCATCCGCCCCTATGATATGTTAGCGCGACTGTAGGGTGAGGTGTTTGGCTTTCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGATACGTATTCGTACAAATGTCGATCTCCTGACCGGCt                            <  2:143480/250‑1 (MQ=255)                                                                                                                                                                                           gCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTGTAGTAACAGATATGATGtt        <  1:322701/250‑1 (MQ=255)                                                                                                                                                                                                                     gCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGaa  >  2:182650/1‑230 (MQ=255)                                                                                                                                                                                                                     gCATCAATGCCCAATGGGGACATTGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGaa  <  1:182650/230‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                               CGGGATTGAAGTTTCCAGTAAACTGTTTGTTGTGTCTTTTCTGATTTATAACATTTTTCAGGAGCTGCAACTCTCCAGCAAACTGGCAGTTCATGATGTGCTGACCAATATTTATAATCGGCGCTACTTTTTCAACAGCGTAGAGTCATTATTGTCGCGACCTGTTGTTAAGGACTTCTGTGTCATGCTGGTTGATATTAATCAGTTCAAACGCATCAATGCCCAATGGGGACATCGTGTGGGTGATAAAGTGCTGGTTTCAATTGTCGATATTATCCAGCAAAGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAA  >  NC_000913/1871281‑1871722

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑