Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,560,820	T→C	E193G (GAA→GGA)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,560,820	0	T	C	100.0%	50.2 / NA	16	E193G (GAA→GGA)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/7);  total (9/7)

TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATC  >  NZ_CP009273/1560734‑1560902                                                                                       |                                                                                   tAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTcccc                                                                                 <  2:342242/90‑1 (MQ=255)      cacaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATTAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaa                                                                            <  1:676166/90‑1 (MQ=255)      cacaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaa                                                                            <  1:509867/90‑1 (MQ=255)       acaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaaa                                                                           >  2:192688/1‑90 (MQ=255)        caGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAg                                                                          >  2:449080/1‑90 (MQ=255)             cAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTg                                                                     >  1:175481/1‑90 (MQ=255)                                   cAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTttc                                               >  2:584062/1‑90 (MQ=255)                                    aaaCTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTttct                                              >  1:648901/1‑90 (MQ=255)                                           tGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCgg                                       <  1:390436/90‑1 (MQ=255)                                              cTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGt                                    <  1:158284/90‑1 (MQ=255)                                                 aaTCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGtttt                                 >  1:274945/1‑90 (MQ=255)                                                   tCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTcc                               <  2:162578/90‑1 (MQ=255)                                                      gAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGttt                                  <  1:446612/84‑1 (MQ=255)                                                      gAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGttt                                  >  2:446612/1‑84 (MQ=255)                                                      gAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGt                            >  1:18602/1‑90 (MQ=255)                                                                            taataTCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt      >  2:79227/1‑90 (MQ=255)                                                                                atCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATc  <  2:648901/90‑1 (MQ=255)                                                                                       |                                                                                   TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATC  >  NZ_CP009273/1560734‑1560902

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