Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	870,468	A→G	Q678Q (CAA→CAG)	pdeI →	EAL domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	870,468	0	A	G	100.0%	89.5 / NA	24	Q678Q (CAA→CAG)	pdeI	EAL domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/5);  total (19/5)

TGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAAGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGT  >  NZ_CP009273/870396‑870554                                                                         |                                                                                       tGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGAtt                                                                      >  1:375853/1‑91 (MQ=255)  gagaTCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATtg                                                                     >  1:178601/1‑91 (MQ=255)     aTCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATtgatg                                                                  <  1:209750/91‑1 (MQ=255)       cACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGAc                                                                >  2:632171/1‑91 (MQ=255)       cACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGAc                                                                >  1:785017/1‑91 (MQ=255)       cACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGAc                                                                >  1:59712/1‑91 (MQ=255)                       gAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGt                                                >  2:720565/1‑91 (MQ=255)                                        tttttATCTCAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGtttt                               >  1:444799/1‑91 (MQ=255)                                                tcAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAAt                       >  2:559775/1‑91 (MQ=255)                                                tcAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAAt                       >  1:604379/1‑91 (MQ=255)                                                tcAATATTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAAt                       >  1:334973/1‑91 (MQ=255)                                                   atatTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTaa                    >  1:802029/1‑91 (MQ=255)                                                     atTGAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATc                  <  2:49183/91‑1 (MQ=255)                                                        gAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTGTGAAATTAATCCAg               >  2:648075/1‑91 (MQ=255)                                                        gAACAACTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAg               >  1:304537/1‑91 (MQ=255)                                                          acaacTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGAc             <  1:720565/91‑1 (MQ=255)                                                          acaacTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGAc             <  2:785017/91‑1 (MQ=255)                                                          acaacTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGAc             <  2:334973/91‑1 (MQ=255)                                                             acTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGc          >  2:799366/1‑91 (MQ=255)                                                             acTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGc          >  1:653996/1‑91 (MQ=255)                                                             acTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGc          >  2:515717/1‑91 (MQ=255)                                                             acTGCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGc          >  1:217538/1‑91 (MQ=255)                                                                gCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATa       >  1:111235/1‑91 (MQ=255)                                                                gCGTAAACAGGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATa       >  1:147664/1‑91 (MQ=255)                                                                     aaCAGGGGGGACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGt  <  1:451940/91‑1 (MQ=255)                                                                         |                                                                                       TGAGATCACTGAAACCGTGATCGAGCGAATTAATGAACATTTTTATCTCAATATTGAACAACTGCGTAAACAAGGGGTACGGATATCGATTGATGACTTTGGCACCGGTTTGTCAAACCTGAAACGTTTTTATGAAATTAATCCAGACAGCATAAAGGT  >  NZ_CP009273/870396‑870554

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