Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,811	A→G	E147E (GAA→GAG)	pdeL →	DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,811	0	A	G	100.0%	191.1 / NA	54	E147E (GAA→GAG)	pdeL	DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/27);  total (27/27)

AATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGG  >  NC_000913/332674‑332947                                                                                                                                          |                                                                                                                                         aaTCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCa                                                                                                                                        >  2:22355/1‑140 (MQ=255)  aTCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCAt                                                                                                                                       >  1:319196/1‑140 (MQ=255)  aTCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCAt                                                                                                                                       >  1:218304/1‑140 (MQ=255)           aCCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACg                                                                                                                              >  1:501036/1‑140 (MQ=255)             cATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACggg                                                                                                                            >  2:213821/1‑140 (MQ=255)               tatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGaa                                                                                                                          >  1:258119/1‑140 (MQ=255)                 tCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAAtt                                                                                                                        >  2:63890/1‑140 (MQ=255)                     cACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATcc                                                                                                                    >  1:513451/1‑140 (MQ=255)                      aCGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATccc                                                                                                                   >  1:133620/1‑140 (MQ=255)                       cGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCGGCTGGGAGCATCCACAAACGGGAATTATCCCa                                                                                                                  >  1:230795/1‑140 (MQ=255)                       cGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCa                                                                                                                  >  2:210248/1‑140 (MQ=255)                       cGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCa                                                                                                                  <  2:454056/140‑1 (MQ=255)                       cGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCa                                                                                                                  >  2:500977/1‑140 (MQ=255)                       cGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCa                                                                                                                  >  1:184831/1‑140 (MQ=255)                        gCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCAc                                                                                                                 <  2:100337/140‑1 (MQ=255)                             aaGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGAt                                                                                                            <  2:268758/140‑1 (MQ=255)                             aaGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGAt                                                                                                            <  2:513451/140‑1 (MQ=255)                                  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGtt                                                                                                       >  1:267135/1‑140 (MQ=255)                                  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGtt                                                                                                       >  2:295829/1‑140 (MQ=255)                                  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGtt                                                                                                       >  2:409448/1‑140 (MQ=255)                                   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGttt                                                                                                      <  1:22355/140‑1 (MQ=255)                                    cAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCGGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTa                                                                                                     <  2:230795/140‑1 (MQ=255)                                               tGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCg                                                                                          <  1:63890/140‑1 (MQ=255)                                                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                                                                                         >  1:178803/1‑140 (MQ=255)                                                 gAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                                                                                        <  2:258119/140‑1 (MQ=255)                                                    aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGTGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                                                                                     <  2:353055/140‑1 (MQ=255)                                                       cACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATc                                                                                  <  1:12387/140‑1 (MQ=255)                                                              tCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCtt                                                                           <  1:183651/140‑1 (MQ=255)                                                               cAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTa                                                                          <  1:324972/140‑1 (MQ=255)                                                               cAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTa                                                                          <  1:295829/140‑1 (MQ=255)                                                                      tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCAt                                                                   >  1:359892/1‑140 (MQ=255)                                                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAAt                                                                <  1:359680/140‑1 (MQ=255)                                                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAAt                                                                <  1:339926/140‑1 (MQ=255)                                                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGAc                                                             <  1:200770/140‑1 (MQ=255)                                                                              aCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGAccc                                                           >  1:365481/1‑140 (MQ=255)                                                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCa                                                       <  1:535249/140‑1 (MQ=255)                                                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCa                                                       <  1:210248/140‑1 (MQ=255)                                                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCa                                                       <  1:409448/140‑1 (MQ=255)                                                                                       cTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACtga                                                  <  2:133620/140‑1 (MQ=255)                                                                                            cgcAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGaaa                                             <  1:485423/140‑1 (MQ=255)                                                                                                 aCTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGAc                                        <  2:184831/140‑1 (MQ=255)                                                                                                   tGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTg                                      <  2:501036/140‑1 (MQ=255)                                                                                                      cGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCgg                                   <  2:359892/140‑1 (MQ=255)                                                                                                         aCTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGata                                >  2:65368/1‑140 (MQ=255)                                                                                                           tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATAtt                              >  2:466887/1‑140 (MQ=255)                                                                                                             aCGGGCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCt                            >  1:108929/1‑140 (MQ=255)                                                                                                               gggCTGTGAGGTGCTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGa                          <  2:493468/140‑1 (MQ=255)                                                                                                                            cTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATAtt                              <  2:135731/123‑1 (MQ=255)                                                                                                                            cTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATAtt                              >  1:135731/1‑123 (MQ=255)                                                                                                                            cTTGTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACa             >  1:273838/1‑140 (MQ=255)                                                                                                                               gTCCGCTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACAttt          >  2:124890/1‑140 (MQ=255)                                                                                                                                    cTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgctgc     >  2:232611/1‑140 (MQ=255)                                                                                                                                    cTGGGAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgctgc     <  2:218304/140‑1 (MQ=255)                                                                                                                                       ggAGCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCgg  >  1:231531/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                         AATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGG  >  NC_000913/332674‑332947

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