Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,685,480	C→T	E47K (GAA→AAA)	pdeK ←	c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,685,480	0	C	T	100.0%	85.3 / NA	28	E47K (GAA→AAA)	pdeK	c‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/15);  total (13/15)

AACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATC  >  NC_000913/3685247‑3685721                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  aaCGACCAGCGAGATTTGAACCGGAAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACACAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGATTTTCAGTTGCGTAGCCgtgt                                                                                                                                                                                                                                  >  2:175479/1‑251 (MQ=255)           gagaTTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTCCCGGCACCGGGCGCGCTGGAATAAAACTTTTACGCAGCGACTGGAACTGGTTAGGCAGCACGACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGAtt                                                                                                                                                                                                                        <  2:49457/251‑1 (MQ=255)                  gAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCGTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTgg                                                                                                                                                                                                                 <  1:113776/251‑1 (MQ=255)                     ccGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGAcc                                                                                                                                                                                                              >  1:269254/1‑251 (MQ=255)                               aaaCAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACa                                                                                                                                                                                                    <  1:219331/251‑1 (MQ=255)                                                       cACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGATAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATaaaa                                                                                                                                                                             >  1:53031/1‑250 (MQ=255)                                                                ctctGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGCACAcc                                                                                                                                                                    >  2:81425/1‑250 (MQ=255)                                                                                                       aGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGtt                                                                                                                             >  2:229304/1‑250 (MQ=255)                                                                                                                        cACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCACTGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCACAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACgcgc                                                                                                           >  2:243254/1‑251 (MQ=255)                                                                                                                           ggCTGCCCACGCCTGCCGTTTTAATGCTGGCAAGAATAGCTTCCGCTGCGGGAATATCGCCTTGCAAAATAGCTGAAGATAAGGTTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTATCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGTTGCCCCCATTGCCATCTGCTTGATTGTGATCGAGCGACGTACTCGCaaa                                                                                                        <  1:252905/251‑1 (MQ=255)                                                                                                                                  cACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAACTAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGc                                                                                                                                      >  1:60753/1‑214 (MQ=255)                                                                                                                                  cACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGc                                                                                                                                      <  2:60753/214‑1 (MQ=255)                                                                                                                                         gCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCg                                                                                           <  1:49643/250‑1 (MQ=255)                                                                                                                                            ggTTTAATGCTGGCAAGAATAGCTTCCGCTGCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGGTGCGTAGCCGTGTTAGAGCGATTCGGCTGGACCAGATGGAACAGCAAAACGGTGCGAAAAATAAAAACTAACACCGGGACAACGGCTGCCACCATTGCCATCTGTTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGaaa                                                                                       <  1:160185/251‑1 (MQ=255)                                                                                                                                                tAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATccc                                                                                   >  2:83543/1‑251 (MQ=255)                                                                                                                                                              aTAGCTTCCGCTTCGGGCATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCg                                                                      >  1:67208/1‑250 (MQ=255)                                                                                                                                                                 gCTTCCGCTTAGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGCGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGCTTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGACACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGt                                                                   >  2:156121/1‑250 (MQ=255)                                                                                                                                                                    tCCGCTTCGGGAATATCGCCGTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGtata                                                                <  2:225966/250‑1 (MQ=255)                                                                                                                                                                           cGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATc                                                         <  1:208923/250‑1 (MQ=255)                                                                                                                                                                                atatCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCg                                                   <  2:212450/251‑1 (MQ=255)                                                                                                                                                                                     gCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACCGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGCGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAACAAATAAAAACGAACACCAGGACAAAGGCTGCCACCATTGCCATCTGCATGATTGTTAACAAGCGACTTACGCGCAAATGGTCTCCCCCCGGAACTCCCAGAtt                                                                              >  2:315984/1‑219 (MQ=255)                                                                                                                                                                                     gCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATGGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGAtt                                                                              <  1:315984/219‑1 (MQ=255)                                                                                                                                                                                         ttCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTaa                                            <  1:83543/249‑1 (MQ=255)                                                                                                                                                                                           cAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTGATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACCCCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTAagaga                                        <  1:253585/251‑1 (MQ=255)                                                                                                                                                                                                 aGCTGAAGATCAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAACATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTCTACCCGATCGCAGCGGCATTAAGAGAGGCGCt                                  >  1:54376/1‑251 (MQ=255)                                                                                                                                                                                                     gAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCTATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACTCCATGACAACGGCTGCCACCATTGCCATCTTCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATGGCCGCCCGAGTATACCCGATCGCAGCGTCATTAAGAGAGGCGCTATCt                              <  1:74573/251‑1 (MQ=255)                                                                                                                                                                                                            aGGGTTCGCGGACAGAGCGAGCAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACt                       <  2:114332/251‑1 (MQ=255)                                                                                                                                                                                                                                 cAATGCTTTTCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATc  >  1:91669/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  AACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACGAGCGACTTACGCGCAAATTGTCTCCCCCCGGAAATCCCAGATTGCCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATC  >  NC_000913/3685247‑3685721

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