Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,515,903	T→C	T644A (ACC→GCC)	pdeA ←	putative c‑di‑GMP phosphodiesterase PdeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,515,903	0	T	C	100.0%	117.8 / NA	35	T644A (ACC→GCC)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/19);  total (16/19)

CGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGTGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTG  >  NC_000913/2515766‑2516032                                                                                                                                          |                                                                                                                                  cGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGgcgc                                                                                                                                 >  2:477590/1‑140 (MQ=255) cGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGgcgc                                                                                                                                 >  1:606717/1‑140 (MQ=255)      tCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAg                                                                                                                            <  2:433356/140‑1 (MQ=255)       cGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGt                                                                                                                           >  2:516510/1‑140 (MQ=255)             aTCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCaa                                                                                                                     >  1:176711/1‑140 (MQ=255)                  tAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCa                                                                                                                <  1:614928/140‑1 (MQ=255)                  tAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCa                                                                                                                <  1:477590/140‑1 (MQ=255)                       gATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCg                                                                                                           >  2:463192/1‑140 (MQ=255)                          ctcACAATCATCGCATCGAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGaa                                                                                                        >  1:433898/1‑140 (MQ=255)                            cacaATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAAc                                                                                                      <  1:463192/140‑1 (MQ=255)                            cacaATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAAc                                                                                                      <  2:606717/140‑1 (MQ=255)                                     cGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAAc                                                                                                      >  2:705670/1‑131 (MQ=255)                                     cGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAAc                                                                                                      <  1:705670/131‑1 (MQ=255)                                           cAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAgct                                                                                       >  2:603572/1‑140 (MQ=255)                                             gCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAg                                                                                         <  2:529328/136‑1 (MQ=255)                                             gCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAg                                                                                         >  1:529328/1‑136 (MQ=255)                                              cGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAgctgct                                                                                    <  2:428826/140‑1 (MQ=255)                                               gtgtTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTc                                                                                   <  1:53354/140‑1 (MQ=255)                                                  ttCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGAt                                                                                >  1:588564/1‑140 (MQ=255)                                                  ttCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGAt                                                                                >  1:292642/1‑140 (MQ=255)                                                       gACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGAtgttgt                                                                           >  1:334032/1‑140 (MQ=255)                                                                      aaaGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTg                                                            >  1:305018/1‑140 (MQ=255)                                                                        aGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCGCAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCg                                                          >  1:27441/1‑140 (MQ=255)                                                                        aGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCg                                                          <  1:675581/140‑1 (MQ=255)                                                                        aGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCg                                                          >  2:660527/1‑140 (MQ=255)                                                                              cATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGaa                                                    <  2:292642/140‑1 (MQ=255)                                                                              cATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGaa                                                    <  2:334032/140‑1 (MQ=255)                                                                                   aTTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGc                                               <  1:172174/140‑1 (MQ=255)                                                                                     tttGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCt                                             <  2:433898/140‑1 (MQ=255)                                                                                               tCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGt                                   <  1:516510/140‑1 (MQ=255)                                                                                                  gCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAGACGCCTGCTCCTCGGTGAt                                <  1:748261/140‑1 (MQ=255)                                                                                                                 aaCCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGc                 <  2:588564/140‑1 (MQ=255)                                                                                                                  aCCGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGcc                >  2:223820/1‑140 (MQ=255)                                                                                                                   ccGTTCGTAGTTGGCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCt               <  1:603572/140‑1 (MQ=255)                                                                                                                         gTAGTTGGCATATCCGGCGCCGAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCgggg         <  2:27441/140‑1 (MQ=255)                                                                                                                                gCATATCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGAtgtg  >  2:70577/1‑140 (MQ=255)                                                                                                                                     tCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGAt     <  1:704862/132‑1 (MQ=255)                                                                                                                                     tCCGGCGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGAt     >  2:704862/1‑132 (MQ=255)                                                                                                                                          |                                                                                                                                  CGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCCTGCAAACGCTTTAACCGTTCGTAGTTGGCATATCCGGTGCCAAAGTCATCAATCGCAATCCGGAAACCAAACTTATGCAGCTGCTCGATGTTGTACATGCTGCTTTCTGCGTTAGAAAACGCCTGCTCCTCGGTGATCTCAAGAATGACCGCCTGCGGGGAGATGTG  >  NC_000913/2515766‑2516032

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