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%	95.8 / NA	29	E193G (GAA→GGA)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/16);  total (13/16)

TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCC  >  NZ_CP009273/1560734‑1560906                                                                                       |                                                                                       tAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTcccc                                                                                     >  1:126565/1‑90 (MQ=255) tAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTcccc                                                                                     <  2:559428/90‑1 (MQ=255) tAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTcccc                                                                                     >  2:1089076/1‑90 (MQ=255)   aaCCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCAt                                                                                   >  1:1087108/1‑90 (MQ=255)   aaCCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCAt                                                                                   <  2:65421/90‑1 (MQ=255)      cacaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaa                                                                                <  1:1089076/90‑1 (MQ=255)      cacaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaa                                                                                <  1:83033/90‑1 (MQ=255)       acaGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCAAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGaaa                                                                               <  2:649847/90‑1 (MQ=255)        caGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAg                                                                              >  2:992564/1‑90 (MQ=255)         aGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAga                                                                             >  2:993252/1‑89 (MQ=255)         aGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGt                                                                             >  2:70767/1‑90 (MQ=255)             cAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTg                                                                         >  2:345962/1‑90 (MQ=255)                    aGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCGTGAAAGTATTGAGGCTAt                                                                  >  2:905681/1‑90 (MQ=255)                    aGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTAt                                                                  >  1:71334/1‑90 (MQ=255)                    aGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTAt                                                                  >  1:183854/1‑90 (MQ=255)                     gCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTAtt                                                                 >  1:945833/1‑90 (MQ=255)                         ggCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGcc                                                             <  2:121339/90‑1 (MQ=255)                         ggCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGcc                                                             >  2:177698/1‑90 (MQ=255)                                     aaCTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTttctt                                                 >  2:695752/1‑90 (MQ=255)                                                   tCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTcc                                   <  1:933944/90‑1 (MQ=255)                                                   tCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTcc                                   <  1:289869/90‑1 (MQ=255)                                                     aGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAg                                 <  1:166840/90‑1 (MQ=255)                                                     aGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAg                                 <  1:345962/90‑1 (MQ=255)                                                     aGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAg                                 <  2:657671/90‑1 (MQ=255)                                                              cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa                        <  1:1222036/90‑1 (MQ=255)                                                              cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa                        <  1:745466/90‑1 (MQ=255)                                                              cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa                        <  1:212713/90‑1 (MQ=255)                                                              cagGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGagaa                        <  1:745467/90‑1 (MQ=255)                                                                         aGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTCCGGCGTTTTCCAGTAACGAGAAGATACGATAGt             <  1:880952/90‑1 (MQ=255)                                                                            taataTCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt          >  2:66823/1‑90 (MQ=255)                                                                                 tCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCt     <  1:843710/90‑1 (MQ=255)                                                                                    aTTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTcc  >  2:605846/1‑90 (MQ=255)                                                                                       |                                                                                       TAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATCTTCC  >  NZ_CP009273/1560734‑1560906

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