Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,027,131	T→C	H296R (CAT→CGT)	dgcQ ←	putative diguanylate cyclase DgcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,027,131	0	T	C	100.0%	133.7 / NA	38	H296R (CAT→CGT)	dgcQ	putative diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/20);  total (18/20)

AGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCAT  >  NC_000913/2026996‑2027266                                                                                                                                        |                                                                                                                                        aGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTc                                                                                                                                     <  2:451848/140‑1 (MQ=255) aGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTc                                                                                                                                     <  1:479171/140‑1 (MQ=255)  gCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCa                                                                                                                                    >  2:405191/1‑140 (MQ=255)             aaaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCACTTCACGAGGATCAAAACTATTCCCTGTTGGTCGGTCAGGATTGGAaga                                                                                                                         >  2:263580/1‑140 (MQ=255)                       gACGTTCCCAGCTAACATAGCGACTGTCCATACGAA‑GCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGGTTGGAAGAGGTCAAAAAtc                                                                                                              <  2:635065/140‑1 (MQ=255)                        aCGTTCCCAGCTAACTTAGCGACTGTCCAAACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAAtc                                                                                                              <  1:306855/140‑1 (MQ=255)                        aCGTTCCCAGCTAACCTAGCGATTGACCAAACGAATGGCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAAtc                                                                                                              <  1:263580/140‑1 (MQ=255)                        aCGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAAtc                                                                                                              <  1:638817/140‑1 (MQ=255)                              ccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCt                                                                                                        <  1:667303/140‑1 (MQ=255)                                 gCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTAc                                                                                                     <  1:100526/140‑1 (MQ=255)                                        tAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATa                                                                                              <  1:650492/140‑1 (MQ=255)                                         aGCGACTGTCCATACGATTGCCTCCCCGCGTATCAGGTTACATTGCCTTCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATag                                                                                             <  2:245841/140‑1 (MQ=255)                                           cGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCGAAAATCTCAGCTTACTGTCATagag                                                                                           <  2:565451/140‑1 (MQ=255)                                              cTGTCCATACGAATGCCGCCCCCCGTGTCATTTTCCATTCCCTGCGCCCGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTg                                                                                        <  2:212113/140‑1 (MQ=255)                                               tGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGa                                                                                       >  1:89319/1‑140 (MQ=255)                                                           tGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                                           >  1:124048/1‑140 (MQ=255)                                                           tGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                                           >  1:373869/1‑140 (MQ=255)                                                           tGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCACTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCg                                                                           >  1:354273/1‑140 (MQ=255)                                                             ccgccCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAg                                                                         >  1:45591/1‑140 (MQ=255)                                                                                           gcCAGCAAGGCCAGTTCCCGAGTATCAAAAAGATTCCCTGTTGGACGGGCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtg                                           <  2:238826/140‑1 (MQ=255)                                                                                           gcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtg                                           >  1:483320/1‑140 (MQ=255)                                                                                           gcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtg                                           >  1:35550/1‑140 (MQ=255)                                                                                           gcCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtg                                           <  2:513868/140‑1 (MQ=255)                                                                                             cAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATtgct                                         <  2:141895/140‑1 (MQ=255)                                                                                                cAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGAGACTCACCAGCGAGGTTATTATCGATAGCGTTTCTTAAAAATTGCTGCa                                      >  2:249915/1‑140 (MQ=255)                                                                                                      cAGTTCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTAc                                >  2:100843/1‑140 (MQ=255)                                                                                                         ttCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCa                             >  2:592639/1‑140 (MQ=255)                                                                                                         ttCACGAGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCa                             <  2:402740/140‑1 (MQ=255)                                                                                                               aGGATCAAAAATATTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGaa                       >  1:544832/1‑140 (MQ=255)                                                                                                                          tatTCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATcccc            >  2:632056/1‑140 (MQ=255)                                                                                                                            ttCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAg          >  1:86062/1‑140 (MQ=255)                                                                                                                            ttCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAg          >  2:592761/1‑140 (MQ=255)                                                                                                                            ttCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAg          <  2:159163/140‑1 (MQ=255)                                                                                                                            ttCCCTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGATTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAg          >  2:57855/1‑140 (MQ=255)                                                                                                                              cccTGTTGGACGGTCAGGATTGTAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCa        <  1:249915/140‑1 (MQ=255)                                                                                                                               ccTGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTATCTTAAAAATTGCTGCTTAGTACGCCCGGGTATACTCATCCCCAGCAc       >  2:192536/1‑140 (MQ=255)                                                                                                                                 tGTTGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGc     <  1:48236/140‑1 (MQ=255)                                                                                                                                    tGGACGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCAt  <  1:632056/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        AGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGTCATAGAGCTGATACTCACCATCGAGGTTTTTATCGATGGCGTTTCTTAAAAATTGCTGCATAGTACGCACGGGAATACTCATCCCCAGCACGCCAT  >  NC_000913/2026996‑2027266

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