Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,272,083	C→T	S453F (TCC→TTC)	pdeN →	putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,272,083	0	C	T	100.0%	76.7 / NA	23	S453F (TCC→TTC)	pdeN	putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/9);  total (14/9)

AAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTG  >  NC_000913/2271839‑2272329                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       aaaCTGCTCACTTCCCTGCCCGCACACCATTTCCAGATTGTGCTGGTAATTACCGAGCGCGATATGTTGAACGCGCAAGAAGCCACACAACTCTTCGCCTGGCTGAACTCGGTAGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAATTCGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGt                                                                                                                                                                                                                                                  >  1:264471/1‑251 (MQ=255)                                tCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCCAGAAGCGACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTAGGCACAGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGc                                                                                                                                                                                                                   >  1:32747/1‑251 (MQ=255)                                                      cGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTg                                                                                                                                                                                             >  1:58193/1‑251 (MQ=255)                                                            cgATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAgggg                                                                                                                                                                                       <  1:78810/251‑1 (MQ=255)                                                                atGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCg                                                                                                                                                                                    >  1:133603/1‑250 (MQ=255)                                                                             aagaagCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAAc                                                                                                                                                                        >  1:21649/1‑249 (MQ=255)                                                                               gaagCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATAGGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGgc                                                                                                                                                                    <  1:240699/251‑1 (MQ=255)                                                                                                                       tAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGc                                                                                                                            >  1:208096/1‑251 (MQ=255)                                                                                                                                           ttCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGtt                                                                                                        >  2:135942/1‑251 (MQ=255)                                                                                                                                                  ccGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTg                                                                                                 <  1:135942/251‑1 (MQ=255)                                                                                                                                                     ggCACAGCGCGCTTACCTATCGTGAGCGGTTTACGCTCGATTATCTGAAAATTGGCCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTTACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAATCGATCGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTGTCCGCTGGAc                                                                                              <  2:37047/251‑1 (MQ=255)                                                                                                                                                            cgcgcTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTg                                                                                       >  1:198114/1‑251 (MQ=255)                                                                                                                                                                 ttATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGCTGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGc                                                                                  >  2:156133/1‑251 (MQ=255)                                                                                                                                                                  tatctatcTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGTCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCGCGTTACCGCTGGATGATTTTGTTCGc                                                                                  <  2:267894/250‑1 (MQ=255)                                                                                                                                                                   atctatctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTg                                                                                >  1:41111/1‑251 (MQ=255)                                                                                                                                                                              gCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTGCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAAc                                                                     <  1:179882/251‑1 (MQ=255)                                                                                                                                                                                     aCGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATAc                                                              <  2:190559/251‑1 (MQ=255)                                                                                                                                                                                                   gAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTgg                                                    <  2:3803/247‑1 (MQ=255)                                                                                                                                                                                                     aaaTTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAggg                                              >  2:283511/1‑250 (MQ=255)                                                                                                                                                                                                      aaTTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTACGgtg                                             >  1:169388/1‑251 (MQ=255)                                                                                                                                                                                                                               aTCGGTACGGAAACGATCACGTTCCCCGTACTTGGCGGGGTGCTGACGCTGGCGAATCGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGATACGCCGGAACAGGCGCGATGGCTATTCGAACGCGGCGTTAATTTCAGTCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGGTTTTGTTCGCTTGAGAAGTAAACCGAATACGCCTCAGTGGTAAGGTGTGCTTACGTcgtttttaattcatag                    <  2:264471/251‑16 (MQ=255)                                                                                                                                                                                                                                             gATCACTTTCCCCGTACTTGACGCGGGGCTGCCGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGCCACGATGGGTATGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCTCTGGCTAAAGAAACATTATCCGCCGCATTGGTAAGTAGAGCTTACGTCACTTATTATTCATAGTGAAAGCAtttcgg      >  2:186381/1‑246 (MQ=255)                                                                                                                                                                                                                                                 aCTTTCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTg  >  1:197187/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       AAACTGCTCACTTCCCTGCCCGCACACCATTTCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAACAGGCGCGATGGCTAAGCGAACGCGGCGTTAATTTCATGCAAGGCTACTGGATTAGCCGCCCGTTACCGCTGGACGATTTTGTTCGCTGGCTAAAGAAACCGTATACGCCGCAGTGGTAAGGTGTGCTTACGTCCCTTATTATTCATAGTGAAAGCATGCCGGATTG  >  NC_000913/2271839‑2272329

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