Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,516,893	G→C	R314G (CGC→GGC)	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,516,893	0	G	C	100.0%	112.2 / NA	31	R314G (CGC→GGC)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/17);  total (14/17)

TTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTC  >  NC_000913/2516659‑2517130                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              tttgggtgtCCTGCATCAACGGCAGCAGGGAACGGGAAAGGGAGCGGATAAAGTTAACGCTCATCATTAAGCCGGAATGACGCCCCATAAACTCAAGATTGTCAGTGCGCAGGCAGCAAAAACTCTTCCCCGCCTCTTGCTCCTGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACTGATCTGTCATCGCCTTCAAATGCCAGCGGCGATTACGCCATTCACTGGCAGGATATATGCGCa                                                                                                                                                                                                                                 <  2:30671/241‑1 (MQ=255)   tCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCAACat                                                                                                                                                                                                                             >  2:199977/1‑251 (MQ=255)         tCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAAAGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAATCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATag                                                                                                                                                                                                                         >  2:176669/1‑249 (MQ=255)               aTCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGAGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGc                                                                                                                                                                                                                  >  1:220485/1‑250 (MQ=255)                       caggggtgTACGGAAAATTGAGCGGATAAAGTGAACGCGCATCATTAAGCCGTAATGACGACGCAAAAACTCAAGATTAGAAAGGCACAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGTGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCTATTCACTTCCTTGATAAATGCTCACCATATAGAGTAGGCAAACGCTaa                                                                                                                                                                                                         <  2:348704/244‑1 (MQ=255)                                          gAGCGGAGACAGTGATCGAGCAACATAGATCCGGAATGACGTCTACTAAACGCAAGATTATCAATTCGCGGGCAGCAAAAACTCTGTCTCGCGTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCTTGCAAATGCTACCGGCTATTAAGCCATTCTCTGCCATGATAAATGCGCACCTTATAGAGCAGGCAAACGCTAAAGTAAATCTGTACCGCaa                                                                                                                                                                                       <  2:24561/250‑1 (MQ=255)                                                       gatcgcgCATCATTAAGCCTTAATTACGACTCATAAACCCAAGATTATCAACGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCGTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGGCAGCGCCTGAAAATGCCAGCGGCGATTAAGCCATTTACTGCCATGATAAAAGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAg                                                                                                                                                                          <  2:435019/247‑1 (MQ=255)                                                            cgcATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAata                                                                                                                                                                    <  1:286311/251‑1 (MQ=255)                                                                      aGCCGTAATGACGACTCATAAACTCAAGATTAACAATGCTCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTcc                                                                                                                                                          <  2:187911/251‑1 (MQ=255)                                                                          gTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCAccc                                                                                                                                                      >  1:459227/1‑251 (MQ=255)                                                                               gacgacTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACGCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATAGAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTCCCACCCCTTg                                                                                                                                                  <  1:98452/250‑1 (MQ=255)                                                                                 cgacTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGAAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCAACCCCTTGCa                                                                                                                                                >  1:294943/1‑250 (MQ=255)                                                                                  gacTCCTAACCTCAAGAATATCAATGCGCAGGCAGCACAACCTCTTGCCAGCCTCTTGCTCAGGCGCTTATCCCACCGCACGAGGGTTGGGTAGGAGCGTTAACGGATGTGTAAGAGCCTGCAAATGCCAGCGGCGCTTAAGCCATTCCCTGACATGATAAATACGCACCATATAGCGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAATTGCCAGCGAATATTCGGTTTACACCCCATGCaa                                                                                                                                               >  2:512730/1‑250 (MQ=255)                                                                                         aaaCTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCt                                                                                                                                       <  1:513730/251‑1 (MQ=255)                                                                                               aaGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCAGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGACTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACAATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTACACCCCTTGCAAAAAGTTATGGTTGt                                                                                                                                 >  1:473929/1‑251 (MQ=255)                                                                                                aGATTATCAATGCGAAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGATAAAGGAAATCAGCACCCCAAGAATAAATGCCAGCGAATATTCCGTTTCCACCCCTTGCAAAAAGTTCTGGTTGt                                                                                                                                 >  2:53226/1‑250 (MQ=255)                                                                                                 gATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCGAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTa                                                                                                                                >  1:397980/1‑250 (MQ=255)                                                                                                      tCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAACACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCAACCCTTGCAAAAAGTTCTGGTTGTAATTc                                                                                                                            >  2:89159/1‑249 (MQ=255)                                                                                                         aTGCGCAGTCAGCAAAAACTCTTGCCCGCCTCTTGCACCGGCGCTTGTTCCAACGCACGAAAGTTGGTTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCGACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAg                                                                                                                       <  2:437985/251‑1 (MQ=255)                                                                                                           gcgcAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAAGTCAGAAGGc                                                                                                                     <  1:420759/251‑1 (MQ=255)                                                                                                               ggccagcaAAACGTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTGGGGTAGGAGCGATAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATACGGTTTCCTCCCCTTGCAAAAAGTTCTGGTGGTAATTCAGAAGGCAAAg                                                                                                                 <  2:85875/248‑1 (MQ=255)                                                                                                                    gcaAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTg                                                                                                            <  1:441729/251‑1 (MQ=255)                                                                                                                     caAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGa                                                                                                           <  1:113709/251‑1 (MQ=255)                                                                                                                            tctTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGcc                                                                                                    >  1:293549/1‑251 (MQ=255)                                                                                                                                         ttGCTCCGCCGCTTGATCCAACGCACGGAAGTTGGGTAGGAGCGTTAACGGATCTGTCCACGTCTGCACATACCAGCGGCGATTAAAACATTCACTGCCATGGTAAATGCGCACCATATAGAGCAGGCAAACGCTCAAGGAAATAAGCACAGCAAGAATAAATGCCAGCGCATATTCGGCTTCCACCACTTGCAAAAAGTTCTGGTTGGAATTCAGAAGGCAAAGCGTTGAAAAAGCCGAGGCGAGAttt                                                                                        >  2:434771/1‑250 (MQ=255)                                                                                                                                           gCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCTTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTaa                                                                                      <  1:8620/250‑1 (MQ=255)                                                                                                                                                  cgcTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGAt                                                                              <  1:271205/251‑1 (MQ=255)                                                                                                                                                       gTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGATAGCg                                                                          >  1:316150/1‑250 (MQ=255)                                                                                                                                                             aCGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCCATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGATAGCGAAGCTTa                                                                   <  1:135663/251‑1 (MQ=255)                                                                                                                                                                                                     cgaatGCAAATGCAAGCAGAGATTCAGCGATTCACTGCCATGGTAAATGCGCAACCTATAGAGCAGGCAAACGCTAAAGGAAATCAGCAACGCCAGTATAAATGCCAGCAAATATTCGGTCTCCACCCCATGCCACAAGTGCCGGCTGTAATTCAGAAAGCAAGGCGTTGTAACCGGCGAGGAAAGATTTAAAAACGGATAGCGCAACTTACCGACCCCCACGGACAAGGTGCTGAAGAAAACggacacc                            >  2:360536/5‑245 (MQ=255)                                                                                                                                                                                                                              aGCCATTCACTGCCATGATAAATGCTCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTTGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTTCCAGTTAACCTGCAATAAAGTCGTTTTc  <  2:274138/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              TTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTTGCCCGCCTCTTGCTCCGGCGCTTGTTCCAACGCACGAAAGTTGGGTAGGAGCGTTAACGGATCTGTCAGCGCCTGCAAATGCCAGCGGCGATTAAGCCATTCACTGCGATGATAAATGCGCACCATATAGAGCAGGCAAACGCTAAAGGAAATCAGCACCGCAAGAATAAATGCCAGCGAATATTCGGTTTCCACCCCTTGCAAAAAGTTCTGGTTGTAATTCAGAAGGCAAAGCGTTGAAACAGCCCAGGTGAGATTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTC  >  NC_000913/2516659‑2517130

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