Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,716	G→T	G416C (GGT→TGT)	pdeC →	c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,276,716	0	G	T	100.0%	79.3 / NA	22	G416C (GGT→TGT)	pdeC	c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/9);  total (13/9)

ATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATT  >  NC_000913/4276477‑4276962                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aTTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTg                                                                                                                                                                                                                                             >  2:38383/1‑251 (MQ=255)        tATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAAtt                                                                                                                                                                                                                                       <  1:97458/250‑1 (MQ=255)                aTGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGAtttt                                                                                                                                                                                                                              <  1:229614/251‑1 (MQ=255)                      tCCGCGATCTGGGCGATTACCTGGCAAAACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTAtt                                                                                                                                                                                                                        >  1:120816/1‑251 (MQ=255)                                                                  gTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGatat                                                                                                                                                                            >  2:202568/1‑251 (MQ=255)                                                                     tCTATTATCCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAttt                                                                                                                                                                          <  2:199238/250‑1 (MQ=255)                                                                              cTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCg                                                                                                                                                                 >  1:319038/1‑250 (MQ=255)                                                                                          gATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTAtt                                                                                                                                                                                                                        <  1:118415/183‑1 (MQ=255)                                                                                          gATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTAtt                                                                                                                                                                                                                        >  2:118415/1‑183 (MQ=255)                                                                                              ttCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGa                                                                                                                                                 <  2:148970/250‑1 (MQ=255)                                                                                                                        atcaatcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTc                                                                                                                       >  2:196650/1‑250 (MQ=255)                                                                                                                         tcaatcaGCAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTATGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAt                                                                                                                     >  2:349919/1‑251 (MQ=255)                                                                                                                                                  ggTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGc                                                                                            >  1:11443/1‑251 (MQ=255)                                                                                                                                                       gTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGgcg                                                                                       >  2:64030/1‑251 (MQ=255)                                                                                                                                                                          tGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTaaaa                                                                     <  1:186225/250‑1 (MQ=255)                                                                                                                                                                             aGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGa                                                                 <  1:95420/251‑1 (MQ=255)                                                                                                                                                                                    gAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTg                                                           >  2:196158/1‑250 (MQ=255)                                                                                                                                                                                    gAGCATGCATGTCTTGATGTTGACAAAATGACGCCGATTAGTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGGCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGa                                                          <  1:29757/251‑1 (MQ=255)                                                                                                                                                                                                     tgatGACAAACTGACGCTGATTCTTCTGGCTTTCCGCCCGGAATGTTACGAAGTGGCAATTGCTGATTTTGGTCTTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTgagg                                         >  2:47978/4‑251 (MQ=255)                                                                                                                                                                                                           cAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAgg                                   <  2:316569/251‑1 (MQ=255)                                                                                                                                                                                                              aaTGAAGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGCTTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGtt                                 >  2:212386/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTgcg       >  1:54455/1‑251 (MQ=255)                                                                                                                                                                                                                                          caggcatgTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACGGAGGAGCAGGTTAACTTGCTGCGCAAACGCGGCGTGCGCTa    <  2:301496/251‑1 (MQ=255)                                                                                                                                                                                                                                            ggcatgTTACGAAGAGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCACAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAACTAGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGCTTGCGGAACACATCATCGCGCTGGAGCACAGCCTGGGGTTAACAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGATACCTGGCTGCGCAAACGCGTCGTGCGCTAtt  >  2:136328/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       ATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATT  >  NC_000913/4276477‑4276962

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