Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,188	(T)5→4	coding (718/1587 nt)	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,184	0	T	.	100.0%	164.0 / NA	39	coding (714/1587 nt)	pdeC	c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/19);  total (20/19)

GCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTT  >  NC_000913/4276047‑4276323                                                                                                                                          |                                                                                                                                            gcagcaACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAttt                                                                                                                                           >  2:451886/1‑140 (MQ=255)     caACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTCCAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAttttata                                                                                                                                       <  2:775113/140‑4 (MQ=255)     caACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAttttata                                                                                                                                       >  1:233618/1‑137 (MQ=255)     caACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCAttttata                                                                                                                                       >  1:752975/1‑137 (MQ=255)       aCGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAat                                                                                                                                    <  2:26034/140‑1 (MQ=255)               tcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTAt                                                                                                                            >  2:727828/1‑140 (MQ=255)                      aTTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCg                                                                                                                     >  2:17179/1‑140 (MQ=255)                         cATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTg                                                                                                                  <  1:9253/140‑1 (MQ=255)                         cATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTg                                                                                                                  <  1:303798/140‑1 (MQ=255)                           tttGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGcc                                                                                                                >  1:336334/1‑140 (MQ=255)                                aGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCccgccg                                                                                                           >  1:97466/1‑140 (MQ=255)                                      tAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAttt                                                                                                     >  1:257756/1‑140 (MQ=255)                                          ggTACAAAGAAATGGCTTTTTATATGCGACAGTTTATTCAACAAAACCCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATACCCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg                                                                                                 <  2:581757/139‑1 (MQ=255)                                          cGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg                                                                                                 <  1:4590/140‑1 (MQ=255)                                             aCAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAg                                                                                              <  2:336334/140‑1 (MQ=255)                                              cAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg                                                                                                 >  2:629014/1‑136 (MQ=255)                                              cAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg                                                                                                 <  1:629014/136‑1 (MQ=255)                                               aaaGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                                                            <  1:725252/140‑1 (MQ=255)                                                     aTGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCGACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTGGCCCGCCGGTATTTTGGGGAGTCTTGTTc                                                                                      <  1:785424/140‑1 (MQ=255)                                                     aTGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTc                                                                                      <  2:233618/140‑1 (MQ=255)                                                      tGGCTATTTATATGCGACTGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCt                                                                                     >  1:168215/1‑140 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTg                                                                                                   <  2:585850/126‑1 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTg                                                                                                   >  1:585850/1‑126 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                                                            >  2:312076/1‑133 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                                                            <  1:312076/133‑1 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTcc                                                                                     >  1:581757/1‑139 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCt                                                                                     >  1:614232/1‑140 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCt                                                                                     >  1:62723/1‑140 (MQ=255)                                                      tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCt                                                                                     >  1:622673/1‑140 (MQ=255)                                                           aTTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt                                                                                >  2:286776/1‑140 (MQ=255)                                                           aTTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt                                                                                <  2:752975/140‑1 (MQ=255)                                                              tataTGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTAc                                                                             <  1:17179/140‑1 (MQ=255)                                                               atatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTAct                                                                            <  1:286776/140‑1 (MQ=255)                                                                      aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTa                                                                     >  2:734511/1‑140 (MQ=255)                                                                      aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTa                                                                     >  2:707797/1‑140 (MQ=255)                                                                      aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTa                                                                     >  1:320841/1‑140 (MQ=255)                                                                      aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTa                                                                     >  2:260714/1‑140 (MQ=255)                                                                          tttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTa                                                                 >  1:517806/1‑140 (MQ=255)                                                                               tCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGa                                                            <  2:168215/140‑1 (MQ=255)                                                                                cGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGAc                                                           <  2:622673/140‑1 (MQ=255)                                                                                  aCAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaa                                                         <  1:451886/140‑1 (MQ=255)                                                                                                     ccATTGTTGCGACTTCATATCAACGTCTTATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAAc                                      <  1:727828/140‑1 (MQ=255)                                                                                                                                ttATAACCCA‑TTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAAc           >  1:105728/1‑140 (MQ=255)                                                                                                                                          tttttATAATCATCTTTTTTTTGCTTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACttt  <  2:517806/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                            GCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTT  >  NC_000913/4276047‑4276323

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