Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,516,215	Δ1 bp	coding (1618/2190 nt)	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,212	0	G	.	100.0%	90.2 / NA	22	coding (1621/2190 nt)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (11/11);  total (11/11)

TGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACC  >  NC_000913/2516138‑2516287                                                                           |                                                                            tgtTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAgggcaga                                                                       <  2:398307/81‑5 (MQ=255)     gccagccaCTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACt                                                                  >  2:398340/1‑81 (MQ=255)      ccagccaCTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACt                                                                  >  2:398339/1‑80 (MQ=255)       cagccaCTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTa                                                                <  2:398308/81‑1 (MQ=255)                 aacaaGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTc                                                       <  2:398309/80‑1 (MQ=255)                  acaaGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCAt                                                     <  2:398310/81‑1 (MQ=255)                        aTTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATgcc                                               >  2:398344/1‑81 (MQ=255)                              gCACATGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGACAGAAACTTATCCGGGGTCATAATGCCGCCGTc                                         <  2:398311/81‑1 (MQ=255)                                aCTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCat                                       <  1:398312/81‑1 (MQ=255)                                   tGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATAtt                                    >  1:398336/1‑81 (MQ=255)                                      aaaTCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATAttt                                   >  2:398356/1‑79 (MQ=255)                                      aaaTCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATAttt                                   >  2:398348/1‑79 (MQ=255)                                      aaaTCAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCa                                 <  2:398313/81‑1 (MQ=255)                                          cAAAACGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTcg                             >  1:398338/1‑81 (MQ=255)                                            aaaCGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTcgcg                           <  2:398314/81‑1 (MQ=255)                                              aCGCGCGCTAAGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTcgcgcg                         <  2:398315/81‑1 (MQ=255)                                                        aGGTTAAACTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCtcat               <  1:398316/81‑1 (MQ=255)                                                             aaaCTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATaa          >  1:398357/1‑81 (MQ=255)                                                              aaCTGAGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAAc         >  2:398345/1‑81 (MQ=255)                                                                   aGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTc     <  2:398317/80‑1 (MQ=255)                                                                   aGCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCa    >  2:398341/1‑81 (MQ=255)                                                                    gCAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCAc   <  1:398318/81‑1 (MQ=255)                                                                     cAATAA‑GGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCAcc  >  1:398358/1‑81 (MQ=255)                                                                           |                                                                            TGTTGCCAGCCACTTCAACAAGGATTCCAGCACTTGCAAATCAAAACGCGCGCTAAGGTTAAACTGAGCAATAAGGGGCAGAAACTTATCCGGGGTCATAATGCCGCCGTCATATTTCAGTCGCGCGAGGATCTCATCATAACCTTCACC  >  NC_000913/2516138‑2516287

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