Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,566,914	A→G	V76A (GTG→GCG)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,566,914	0	A	G	100.0%	164.1 / NA	48	V76A (GTG→GCG)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/25);  total (23/25)

CGGATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCACGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATGCGCAACGGC  >  NC_000913/1566778‑1567050                                                                                                                                         |                                                                                                                                         cGGATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGt                                                                                                                                       >  2:263213/1‑140 (MQ=255)   gATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGCTTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTta                                                                                                                                     >  1:651109/1‑140 (MQ=255)   gATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTta                                                                                                                                     <  1:203708/140‑1 (MQ=255)   gATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTta                                                                                                                                     >  2:360741/1‑140 (MQ=255)      agagGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTtaata                                                                                                                                  >  2:147614/1‑140 (MQ=255)      agagGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTtaata                                                                                                                                  >  2:297519/1‑140 (MQ=255)        agGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTaataat                                                                                                                                >  1:841732/1‑140 (MQ=255)        agGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTaataat                                                                                                                                >  2:383256/1‑140 (MQ=255)        agGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTaataat                                                                                                                                >  1:307601/1‑140 (MQ=255)        agGATTTTTTTCAGAACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTaataat                                                                                                                                >  1:513826/1‑140 (MQ=255)         ggATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATc                                                                                                                               <  1:263213/140‑1 (MQ=255)         ggATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATc                                                                                                                               <  1:343375/140‑1 (MQ=255)              tttttCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCg                                                                                                                          <  1:864177/140‑1 (MQ=255)                         ccGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCAc                                                                                                               <  1:297519/140‑1 (MQ=255)                                 ccATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAAc                                                                                                       <  1:452806/140‑1 (MQ=255)                                  cATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACt                                                                                                      <  1:383256/140‑1 (MQ=255)                                       cGACAATTTCTACCGGAATTCCTATGTGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGc                                                                                                 >  2:1256588/1‑140 (MQ=255)                                        gACAATTTCTACCGGAATTCCTATGTGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCt                                                                                                <  1:1256588/140‑1 (MQ=255)                                        gACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCt                                                                                                <  2:307601/140‑1 (MQ=255)                                              ttCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACtt                                                                                          <  2:513826/140‑1 (MQ=255)                                               tCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTg                                                                                         >  2:161678/1‑140 (MQ=255)                                                         ttCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCa                                                                                     >  1:819650/1‑134 (MQ=255)                                                         ttCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCa                                                                                     <  2:819650/134‑1 (MQ=255)                                                                              ccGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCggg                                                          >  2:386544/1‑140 (MQ=255)                                                                                      gTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCg                                                  <  1:147614/140‑1 (MQ=255)                                                                                      gTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCACGCGGGTCGATGCg                                                  <  2:651109/140‑1 (MQ=255)                                                                                        aTGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGa                                                <  2:581533/140‑1 (MQ=255)                                                                                          gCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACa                                              >  1:196213/1‑140 (MQ=255)                                                                                           cTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACaa                                             <  1:103698/140‑1 (MQ=255)                                                                                            tGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAAt                                            <  1:161678/140‑1 (MQ=255)                                                                                                   gTATTAGCC‑TTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATaaa                                    >  2:613913/1‑140 (MQ=255)                                                                                                         gCCTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCg                               <  2:18395/140‑1 (MQ=255)                                                                                                          gccttCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGa                              <  1:613913/137‑1 (MQ=255)                                                                                                          ccTTTCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGa                              <  2:402756/140‑1 (MQ=255)                                                                                                            tttCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATa                            <  2:416814/140‑1 (MQ=255)                                                                                                              tCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACt                          >  2:75357/1‑140 (MQ=255)                                                                                                              tCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACt                          >  2:937268/1‑140 (MQ=255)                                                                                                              tCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACt                          >  2:265154/1‑140 (MQ=255)                                                                                                              tCGACATCGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACt                          <  1:265154/140‑1 (MQ=255)                                                                                                                     cGTCAACCTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATaa                   <  2:842926/140‑1 (MQ=255)                                                                                                                           ccTGGGCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCa             <  2:612391/140‑1 (MQ=255)                                                                                                                               ggCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCa                        >  1:768519/1‑125 (MQ=255)                                                                                                                               ggCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCa                        <  2:768519/125‑1 (MQ=255)                                                                                                                               ggCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATgcg         >  2:1021590/1‑140 (MQ=255)                                                                                                                               ggCAGAAAGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATgcg         >  1:347143/1‑140 (MQ=255)                                                                                                                                    aaaGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATGCGCAACg    >  1:1207876/1‑140 (MQ=255)                                                                                                                                     aaGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATGCGCAACgg   >  1:124037/1‑140 (MQ=255)                                                                                                                                      aGCGCGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATGCGCAACGGc  <  2:347143/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CGGATAGAGGATTTTTTTCAGCACCCGAAACCCCATCTCGACAATTTCTACCGGAATTCCTATGCGGGCATGCACTTCCGCGACGGTATGCTGGATTTGTATTAGCCTTTCGACATCGTCAACCTGGGCAGAAAGCACGTTAATAATCCAGCGTTCCATCGCACTCTTCAACTGCCGCTCAACTTGTTCATTACTCAAGAATTCTTCGGCATGCGGGTCGATGCGGACAATTCGATAAAACTCGATACTCAGATAATGAGCATGCGCAACGGC  >  NC_000913/1566778‑1567050

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: ‑