Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,517,170	T→C	A221A (GCA→GCG)	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,170	0	T	C	100.0%	199.5 / NA	55	A221A (GCA→GCG)	pdeA	putative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/23);  total (32/23)

ATTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCTGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCA  >  NC_000913/2517041‑2517306                                                                                                                                  |                                                                                                                                         aTTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCcagccagc                                                                                                                                >  2:129792/1‑140 (MQ=255)  tttAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCt                                                                                                                               <  2:315263/140‑1 (MQ=255)                      cTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTtct                                                                                                           >  2:351167/1‑140 (MQ=255)                      cTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTtct                                                                                                           >  2:217437/1‑140 (MQ=255)                       ttACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTtctc                                                                                                          >  1:483826/1‑140 (MQ=255)                         aCCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCtt                                                                                                        <  1:640596/140‑1 (MQ=255)                         aCCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCtt                                                                                                        >  1:567595/1‑140 (MQ=255)                               cccGAGGGTAAAGATGATGAAGAAAACGGGTACAAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCa                                                                                                  >  1:809244/1‑140 (MQ=255)                                ccGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCaa                                                                                                 <  1:363377/140‑1 (MQ=255)                                   aGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAacga                                                                                              <  2:938819/140‑1 (MQ=255)                                    gggTAAAGATGATGAAGAAAACGGGTACAAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAacgac                                                                                             <  2:809244/140‑1 (MQ=255)                                      gTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACgg                                                                                           <  1:217437/140‑1 (MQ=255)                                               atgaAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATc                                                                                  >  1:730185/1‑140 (MQ=255)                                                tgaAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATcc                                                                                 <  1:454788/140‑1 (MQ=255)                                                tgaAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATcc                                                                                 >  2:789715/1‑140 (MQ=255)                                                   agaaAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTg                                                                              <  2:483826/140‑1 (MQ=255)                                                         cGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCcac                                                                        >  1:1053856/1‑140 (MQ=255)                                                            gTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAat                                                                     >  1:533865/1‑140 (MQ=255)                                                                  ggTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTgc                                                               >  2:580235/1‑140 (MQ=255)                                                                          gCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTg                                                       >  1:790086/1‑140 (MQ=255)                                                                            aaTAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTggg                                                     <  1:789715/140‑1 (MQ=255)                                                                            aaTAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTggg                                                     >  1:570449/1‑140 (MQ=255)                                                                            aaTAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTggg                                                     <  1:580235/140‑1 (MQ=255)                                                                                     gTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTgcgc                                                             >  1:213069/1‑123 (MQ=255)                                                                                     gTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTgcgc                                                             <  2:213069/123‑1 (MQ=255)                                                                                     gTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACaa                                            >  2:164897/1‑140 (MQ=255)                                                                                         tCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCAt                                        <  1:346240/140‑1 (MQ=255)                                                                                         tCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCAt                                        <  1:129792/140‑1 (MQ=255)                                                                                         tCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCAt                                        <  1:127615/140‑1 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGTACTTACAATCATGCGAGTgaga                              >  1:767184/1‑140 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTgaga                              >  2:535776/1‑140 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTgaga                              <  2:533865/140‑1 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTgaga                              >  2:457660/1‑140 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTga                                <  2:544682/138‑1 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTga                                <  2:212926/138‑1 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTga                                >  1:212926/1‑138 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTga                                >  1:544682/1‑138 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACACTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAg                                   >  2:36722/1‑135 (MQ=255)                                                                                                   tgCACAACAGAAGTAGCAGCACACTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAg                                   <  1:36722/135‑1 (MQ=255)                                                                                                     cacaACAGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGAta                            >  1:967868/1‑140 (MQ=255)                                                                                                          caGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGa                       >  2:10567/1‑140 (MQ=255)                                                                                                           aGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGaa                      <  1:351167/140‑1 (MQ=255)                                                                                                           aGAAGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGaa                      >  2:325290/1‑140 (MQ=255)                                                                                                             aaGTAGCAGCACGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGaaaa                    >  2:675277/1‑140 (MQ=255)                                                                                                                  gcagcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtg               >  2:207135/1‑140 (MQ=255)                                                                                                                    agcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgtt             >  1:809873/1‑140 (MQ=255)                                                                                                                    agcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgtt             >  2:606916/1‑140 (MQ=255)                                                                                                                     gcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgttg            <  2:567595/140‑1 (MQ=255)                                                                                                                     gcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgttg            <  1:10567/140‑1 (MQ=255)                                                                                                                     gcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgttg            >  1:309372/1‑140 (MQ=255)                                                                                                                     gcaCGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAtgttg            <  1:325290/140‑1 (MQ=255)                                                                                                                        cGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTaa         >  2:1031290/1‑140 (MQ=255)                                                                                                                        cGCTTAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTaa         >  2:825747/1‑140 (MQ=255)                                                                                                                           ttAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAg                   <  2:272715/127‑1 (MQ=255)                                                                                                                           ttAACGCCGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAg                   >  1:272715/1‑127 (MQ=255)                                                                                                                               cgccgcCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGca  <  2:730185/140‑1 (MQ=255)                                                                                                                               cgccgcCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGca  <  2:309372/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         ATTTAAAAACGGATAGCGAAGCTTACCGACCCCGAGGGTAAAGATGATGAAGAAAACGGGTACCAGGTAACCGGCAATAAAGTCGTTTTCATAAGGTGTGCACAACAGAAGTAGCAGCACGCTTAACGCTGCCAGCCAGCTTAAGGTAAATGCGCGTTTCTCTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCA  >  NC_000913/2517041‑2517306

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