Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,566,497	G→A	A215V (GCT→GTT)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,566,497	0	G	A	100.0%	87.9 / NA	26	A215V (GCT→GTT)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

CTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATCGCGATATCAATACTGTTAA  >  NC_000913/1566258‑1566719                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               cTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTaa                                                                                                                                                                                                                      >  1:292275/1‑250 (MQ=255)             aTACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTg                                                                                                                                                                                                         >  2:203310/1‑251 (MQ=255)                        aaCAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGa                                                                                                                                                                                              <  1:224121/251‑1 (MQ=255)                                            aataatTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTtat                                                                                                                                                                             >  1:9975/1‑248 (MQ=255)                                               aatTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCCCCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGataa                                                                                                                                                                       <  1:282009/251‑1 (MQ=255)                                                                    tATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGt                                                                                                                                                   <  1:167837/250‑1 (MQ=255)                                                                               aaaaTTTAACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGACTACTTCGGCAATAACACTAAAATAATGTCGAACTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGAAGGCTTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTAtt                                                                                                                                       >  1:249955/1‑251 (MQ=255)                                                                                            cAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCttttt                                                                                                                          <  1:255543/251‑1 (MQ=255)                                                                                             aGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGAttt                                                                                                                                                                                >  2:640/1‑196 (MQ=255)                                                                                             aGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGAttt                                                                                                                                                                                <  1:640/196‑1 (MQ=255)                                                                                                                 tCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGt                                                                                                      >  1:283601/1‑250 (MQ=255)                                                                                                                      tgtTACGAATGGTTTGATTGAAAATACCGTCGAAATACTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACg                                                                                                  >  2:49817/1‑249 (MQ=255)                                                                                                                                     gattgaAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGttt                                                                                 <  1:385126/251‑1 (MQ=255)                                                                                                                                      attgaAAAAACCATCGAAAACCTGAAGCAGACGGGAGATATGGCCTACTTCGGCAATACCTCTAAAATAATGTGGACCTTTAAGGGTAAGCCACAGGGCAAAATCAACCTGGCTTAAAGGCTAACTACTGCCTAAATCAGAAGCCAGCAGGATTATATAGATAATATCTATTTCCCACGAAAGCAGTGAGGCTATTTGCCGTTATTTTTCTACTTCGGCGTTTCCCAGTAACGAGAAGATACGGTAGtttt                                                                                <  2:127809/251‑1 (MQ=255)                                                                                                                                                 cggcGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCttt                                                                     >  1:113810/4‑251 (MQ=255)                                                                                                                                                 cGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAAGGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCttt                                                                     <  2:230969/251‑1 (MQ=255)                                                                                                                                                  gTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGAAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCNTTg                                                                    >  2:365236/1‑251 (MQ=255)                                                                                                                                                              aaTCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGt                                                        <  1:49817/251‑1 (MQ=255)                                                                                                                                                                         agaTATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAgga                                              >  1:163893/1‑249 (MQ=255)                                                                                                                                                                                      ttCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATCCGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGtc                                 >  1:270311/1‑250 (MQ=255)                                                                                                                                                                                             aTACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTc                         <  1:228910/251‑1 (MQ=255)                                                                                                                                                                                                 cACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGACAGCATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATc                     >  1:361296/1‑251 (MQ=255)                                                                                                                                                                                                           atGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATCGCGATATCaa           <  1:391783/251‑1 (MQ=255)                                                                                                                                                                                                                  ccTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTGCTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACGGAAGGTAAACGCGCGAGTCATCCCTTCCATCGCGAGGTCAATACTGtt    <  2:37389/251‑1 (MQ=255)                                                                                                                                                                                                                   cTTTATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATCGCGATATCAATACTGtt    >  1:292027/1‑250 (MQ=255)                                                                                                                                                                                                                     ttATGGTTAAACCACAGGCCAAAATCAACCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATCGCGATATCAATACTGTTaa  <  1:365236/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               CTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATTTGCGATACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCATGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCCTTTGAGGCACTACTGTCACTAAAGGTAAACGCGCGAGTCATCACTTCCATCGCGATATCAATACTGTTAA  >  NC_000913/1566258‑1566719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑