Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,626,297	A→G	pseudogene (330/1368 nt)	ydfJ ←	MHS family MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,626,297	0	A	G	100.0%	30.1 / NA	11	pseudogene (330/1368 nt)	ydfJ	MHS family MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/6);  total (5/6)

ACCCAACCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGATGGTAATAAAGAGCACTTTTTTACGCCCTATTCTGT  >  NZ_CP009273/1626234‑1626382                                                                |                                                                                      aCCCAACCCCTGAATAATACGCAACGTAACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAggt                                                             <  1:358904/90‑1 (MQ=255)           tGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATc                                                   <  2:96155/90‑1 (MQ=255)           tGAATAATACGCAACCTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATc                                                   <  2:107122/90‑1 (MQ=255)                 ataCGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCcatc                                             >  1:282655/1‑90 (MQ=255)                      cAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGt                                        >  2:117105/1‑90 (MQ=255)                       aaCGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTg                                       <  2:437731/90‑1 (MQ=255)                        aCGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATc                                                   <  1:470667/77‑1 (MQ=255)                        aCGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATc                                                   >  2:470667/1‑77 (MQ=255)                                        gggTGCAAAAACACCAATCTGTGCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGATGGTAATAAAGAGCAc                      >  1:416306/1‑90 (MQ=255)                                                  acacCAATCTGTGCGTAGGTCGCTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGATGGTAATAAAGAGCACTTTTTTACGc            <  2:439493/90‑1 (MQ=255)                                                            gtgCGTAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGATGGTAATAAAGAGCACTTTTTTACGCCCTATTCTGt  >  1:188902/1‑90 (MQ=255)                                                                |                                                                                      ACCCAACCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGATGGTAATAAAGAGCACTTTTTTACGCCCTATTCTGT  >  NZ_CP009273/1626234‑1626382

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