Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,332	T→A	E229V (GAA→GTA)	dgcQ ←	diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,332	0	T	A	100.0%	83.0 / NA	29	E229V (GAA→GTA)	dgcQ	diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (15/14);  total (15/14)

AGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGT  >  NC_000913/2027089‑2027556                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 aGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTg                                                                                                                                                                                                                           >  2:256507/1‑251 (MQ=255)      aaccAGGTCAAGAGGATCAACAATATTCCCTGTTGGCTGGTCAGGCTTGGCACAGGGCAACAATCTCAGCTTACTGT‑AATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTACTTAAAAATTGCTGCATAGTACGTACGTGCATACTCATCCCCAGCACACCATACCAGTAGTTATTACAGTCTACCGGAACAATGACGGTAACCTGCGGTTCAGTATTTCTGGCGTGTACCGGTTGcca                                                                                                                                                                                                                        >  2:54443/3‑247 (MQ=255)         cAGTTCCCGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAACAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGAGTGTACCGGTTGCGAAGTg                                                                                                                                                                                                                    >  2:279037/1‑250 (MQ=255)                  gggaacAAAAATATTCCCTGTTGGATGGTCGGTATTGTAAGAGGGCAAAAATCTCAGCTTACTTT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCTGTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCGACCCCAGCACGCCATTCCAGTAGTTATTACTGTCTACCGGAACGCTGTCGGGAACCTGCGGTTCAGTATTGCTGGCGTGTACCTGTGGCGAAGTTAACCAGCGTa                                                                                                                                                                                                          <  2:144271/246‑1 (MQ=255)                              atTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCGCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGAc                                                                                                                                                                                              <  2:319725/251‑1 (MQ=255)                                 cccTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAGCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGa                                                                                                                                                                                            <  1:222313/250‑1 (MQ=255)                                     gTTGGATGGTCAGTATTTGAAGAGGTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCCCCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTTACGAttt                                                                                                                                                                                         <  1:255176/249‑1 (MQ=255)                                                     tGGAAGAGGTCAAAAATCTCATCTGTC‑GTACATAGAGCGGATACTCGCCATCGAGGTTTTTATCGATGGCTTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACTGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTTCGAAGTGAACCAGGGTACCGCGCTGTTACGATTTTCTCGGTGTGAATGGcc                                                                                                                                                                        <  2:193463/250‑1 (MQ=255)                                                            ggTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAAcc                                                                                                                                                                >  1:148046/1‑251 (MQ=255)                                                               cAAAAATTTCAGCTTACTGT‑CATAGATCTGATACTCACCATCGAGGTTTTTACCGATGGCGTTTCTTAAAAAGTGCTGCATAGTACGCACGGGAATACTAATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCGGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAg                                                                                                                                                              <  2:75407/250‑1 (MQ=255)                                                                 aaaaTCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTAGCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATTCCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGt                                                                                                                                                           <  1:35936/251‑1 (MQ=255)                                                                              aCTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGc                                                                                                                                              <  1:121221/251‑1 (MQ=255)                                                                                       agagCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGtaata                                                                                                                                       >  2:138425/1‑250 (MQ=255)                                                                                            tGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGt                                                                                                                                 >  2:324699/1‑251 (MQ=255)                                                                                                               ttttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTaaa                                                                                                               >  1:49653/1‑250 (MQ=255)                                                                                                                tttATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGCCATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAAc                                                                                                              >  2:66975/1‑250 (MQ=255)                                                                                                                                                    aCGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGATATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCAcgcg                                                                         >  2:206591/1‑251 (MQ=255)                                                                                                                                                       cACGGGAATACGCTTCCCCTGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATGGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAAGAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCTGCAGGCgaga                                                                      <  2:142406/251‑1 (MQ=255)                                                                                                                                                           ggAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGacat                                                                   >  1:166794/1‑250 (MQ=255)                                                                                                                                                             aaTACTCATCCCCAGCAAGCCATACCAGTAGGTATTACTGTCTACCGAAACACTGACTGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGGACCGCGCGGTGACGATTTTATCCTTGTGAATGGCCAATAAAACCGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGgcatac                                                                 >  1:162807/1‑250 (MQ=255)                                                                                                                                                                        ccAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGTGATACATCGCCTGTTCa                                                     <  2:307963/251‑1 (MQ=255)                                                                                                                                                                                        gtagtTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCGGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGt                                     <  1:164670/251‑1 (MQ=255)                                                                                                                                                                                        gtagtTATTACTGGCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGt                                     <  2:265526/251‑1 (MQ=255)                                                                                                                                                                                                         ccGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGCGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGGCGGCTGCCTCGAAACGTAAAATCCGGCACGCGAGACATACACCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCgcagc                    >  2:249297/1‑251 (MQ=255)                                                                                                                                                                                                                  tgacgggcaCCTGCAGTTTAGTATTGCCGGCGTCCAACGGTTGCGAACGGAAACAGCGCAACGCGCGGAGACGATTTTCAAATTGGGAATGGCCACTATACCAGGGATGGGAGACATAGACGTAATAGCGCGTTGGAACATTGCGCGTACAAAAGGTAGGCTGCGTCGACACGTAAAATCCGGCACGCGAGACATACATCGCCAGTTCAACCATCGACGAGGAGTTGCGCGCTAATCGCAGCAAGAAAcc            >  2:20623/9‑250 (MQ=255)                                                                                                                                                                                                                    aCGGTAACCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAAc         >  1:238596/1‑251 (MQ=255)                                                                                                                                                                                                                          aCCTGCGGTTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAg   <  1:179949/251‑1 (MQ=255)                                                                                                                                                                                                                           ccTGCGGTTCAGTATTGCGGGAGTGTACGGGTTGAGAAGGGAACCAGCGTAGCGCGCGGTGACGAGTTGCTCGTTGTGTATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCGTGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGt  <  1:54443/251‑1 (MQ=255)                                                                                                                                                                                                                           ccTGCGGCTCAGTATTGCTGGCGTGTACCGGTTGCGAAGTGAACCAGAGTACCGCGCGGTGACGATTTTCTCGTTGTGCATGGCCAATAACCCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGt  >  2:32735/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 AGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT‑CATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCATACCAGTAGTTATTACTGTCTACCGGAACACTGACGGTAACCTGCGGTTCAGTATTGCTGGCGTGTTCCGGTTGCGAAGTGAACCAGCGTACCGCGCGGTGACGATTTTCTCGTTGTGAATGGCCAATAAACCAGGGTTGGGTGACATAGCCGTAATAACGCGTTGGAACATTGCGCGTAAACAAGGTCGGCTGCGTCGAAACGTAAAATCCGGCACGCGAGACATACATCGCCTGTTCAACCATCGACGAGGAGTTGTGCGCTAGTCGCAGCAAGTAACCAACTTCCAGT  >  NC_000913/2027089‑2027556

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