Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,517,367	(A)5→4	coding (466/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,517,363	0	A	.	100.0%	93.9 / NA	21	coding (470/2190 nt)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/11);  total (10/11)

AAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGT  >  NC_000913/2517286‑2517439                                                                              |                                                                             aaaGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaa                                                                             >  1:449455/1‑79 (MQ=255)  aaGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaag                                                                          >  1:449462/1‑80 (MQ=255)   aGCATGTTGTAAATCAGCAAAGCGGTGAAAAGGCTTAGCAAATCAACGACCGAGAAAATGGCATCCGCATCGCCGaaaagg                                                                         <  1:449420/81‑3 (MQ=255)           gtAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATc                                                                >  1:449463/1‑81 (MQ=255)            tAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATc                                                                <  1:449421/80‑1 (MQ=255)            tAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCt                                                               >  2:449461/1‑81 (MQ=255)              aaTCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCttt                                                             <  1:449422/81‑1 (MQ=255)                  aGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATTGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGc                                                          <  1:449423/80‑1 (MQ=255)                  aGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCg                                                         <  2:449424/81‑1 (MQ=255)                          ggTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATc                                                  >  1:449471/1‑80 (MQ=255)                            tGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCa                                                 >  1:449449/1‑79 (MQ=255)                             gAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCa                                                 >  2:449469/1‑78 (MQ=255)                                  ggCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCt                                         >  2:449470/1‑81 (MQ=255)                                   gCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTc                                        <  1:449425/81‑1 (MQ=255)                                        gCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCca                                    <  2:449426/80‑1 (MQ=255)                                         cAAATCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCcaca                                  <  1:449427/81‑1 (MQ=255)                                            aTCAACGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGa                               <  1:449428/81‑1 (MQ=255)                                               aaCGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATa                             >  1:449464/1‑80 (MQ=255)                                               aaCGACCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATa                             >  2:449466/1‑80 (MQ=255)                                                   aCCGTGAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACAtgc                        <  1:449429/81‑1 (MQ=255)                                                        gAAAATGGCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCAt                   >  2:449474/1‑81 (MQ=255)                                                               gCATCCGCATCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATg             <  1:449430/80‑1 (MQ=255)                                                                       aTCGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGc    <  2:449431/81‑1 (MQ=255)                                                                         cGCCG‑AAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGt  >  1:449473/1‑81 (MQ=255)                                                                              |                                                                             AAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGT  >  NC_000913/2517286‑2517439

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