Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,143,684	(G)6→7	coding (419/3318 nt)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,143,684	1	.	G	100.0%	73.7 / NA	22	G140? (GGT→GNT)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/14);  total (8/14)

TTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTTCATTATTGTGCTGTT  >  NC_000913/2143436‑2143933                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          ttATCGCCTGGCTGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGGTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCAGTGTTCCGCCTCTGGTggggggg                                                                                                                                                                                                                                                           <  2:343638/250‑2 (MQ=255)                             ggAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTGCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTTACGCCCGgag                                                                                                                                                                                                                              <  1:138938/251‑1 (MQ=255)                                        gAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCctc                                                                                                                                                                                                                     <  2:73253/249‑1 (MQ=255)                                        gAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCtctc                                                                                                                                                                                                                   <  2:185764/251‑1 (MQ=255)                                                           atCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATggg                                                                                                                                                                                                >  1:259018/1‑251 (MQ=255)                                                              gccgcATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGCAGGGGCAGAGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCCACATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACCCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                                                                                                                                                                                             >  1:171729/1‑251 (MQ=255)                                                               ccgcATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACt                                                                                                                                                                                            >  1:144165/1‑251 (MQ=255)                                                                   atccatccTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCa                                                                                                                                                                                        <  2:55612/251‑1 (MQ=255)                                                                                     caccaGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTg                                                                                                                                                                      >  1:204715/1‑251 (MQ=255)                                                                                          gCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACt                                                                                                                                                                 <  1:72014/251‑1 (MQ=255)                                                                                                                gacCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAttt                                                                                                                                           <  1:21936/251‑1 (MQ=255)                                                                                                                gacCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAttt                                                                                                                                           >  2:315067/1‑251 (MQ=255)                                                                                                                            tgttgAATCCGTGGTCGGGGCAGTGCTGCTACTTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGGCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGGCGCCCGGAGACGAGCCTCTCAGGGCATTTTTGATAGGGGTACTGTCAGAATCCATCGGCGCTCTGTCACTGGTGCCGCTGGGATTGTTAGTTAAACCACACTa                                                                                                                                <  2:341330/250‑1 (MQ=255)                                                                                                                                  aGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATctgctg                                                                                                                         >  2:109681/1‑251 (MQ=255)                                                                                                                                                       gctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTg                                                                                                    <  2:282124/251‑1 (MQ=255)                                                                                                                                                           cGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCttt                                                                                                <  1:248900/251‑1 (MQ=255)                                                                                                                                                                   gctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACGGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACACACGGTTGCTTTTTGAGTc                                                                                         <  2:264449/250‑1 (MQ=255)                                                                                                                                                                      gcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAAGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgctg                                                                                     <  2:259018/251‑1 (MQ=255)                                                                                                                                                                      gcCGGGGTGTAACCCCTTGCCAAATCGGGCTGACTGGCTGCGTCTGGCACTCGGCGGCGCTATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCgctg                                                                                     <  2:79644/251‑1 (MQ=255)                                                                                                                                                                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTc                                                                                  >  1:248794/1‑250 (MQ=255)                                                                                                                                                                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTc                                                                                  >  1:38280/1‑250 (MQ=255)                                                                                                                                                                                                 ggctgccgggCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGt                                                            <  2:204715/241‑1 (MQ=255)                                                                                                                                                                                                                                         ccGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACtt                  <  1:358545/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ggggggg‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTGTTGAGTCGCTGCTCACGTTAGCCATCACACTGACTTTAAGCTGGCTTTCGAGGCTGGATCTGCCGTGGCCTTTTACTGTCATTATtgt        <  2:228648/249‑1 (MQ=255)                                                                                                                                                                                                                                                      gggg‑TGTTCTGGTTGGCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTTCATTATTGTGCTg    <  2:266563/251‑1 (MQ=255)                                                                                                                                                                                                                                                        gg‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTTCATTATTGTGCTGtt  <  2:156795/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          TTATCGCCATGCCGGGCGCATGTGGCCGGGAATTGCGCTGAGCTGTTCGCTGGGAAATATCGCCGCATCCATCCTGCTTTTTTCCACCAGCTCGCTGAACATGACCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACCACACTATCTGCTGCGCCATCGCAACCCACGGTTGCTTTTTGAGTCGCTGCTCACGTTAGCCATCACACTGACGTTAAGCTGGCTTTCGATGCTGTATCTGCCGTGGCCTTTTACTTTCATTATTGTGCTGTT  >  NC_000913/2143436‑2143933

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