Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	525,225	A→G	pseudogene (248/366 nt)	ylbG ←	helix‑turn‑helix domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	525,225	0	A	G	100.0%	42.8 / NA	13	pseudogene (248/366 nt)	ylbG	helix‑turn‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/6);  total (7/6)

TTCAGTATCGTCTTTTGAGATTCTCAAGGTTATGGAGTAACCAGACAGAACGGACACCACTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGAAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGT  >  NZ_CP009273/525100‑525348                                                                                                                              |                                                                                                                            ttCAGTATCGTCTTTTGAGATTCTCAAGGTTATGGAGTAACCAGACAGAACGGACACCACTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCg                                                                                                                  <  1:268342/137‑1 (MQ=255) ttCAGTATCGTCTTTTGAGATTCTCAAGGTTATGGAGTAACCAGACAGAACGGACACCACTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCg                                                                                                                  >  2:268342/1‑137 (MQ=255) ttCAGTATCGTCTTTTGAGATTCTCAAGGTTATGGAGTAACCAGACAGAACGGACACCACTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCGt                                                                                                                 <  2:183930/138‑1 (MQ=255)                                                           acTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCgcgg                                                      >  2:251281/1‑139 (MQ=255)                                                                               ccTGTTTACGCAGCTCGTTGCTGGCCCGGAGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTgcg                                   >  2:155425/1‑138 (MQ=255)                                                                                             tCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCtt                                                                      <  1:125919/89‑1 (MQ=255)                                                                                             tCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCtt                                                                      >  2:125919/1‑89 (MQ=255)                                                                                                           gTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGAc                                                  >  1:471198/1‑95 (MQ=255)                                                                                                           gTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGAc                                                  <  2:471198/95‑1 (MQ=255)                                                                                                           gTGCTGACCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTcg      >  2:502487/1‑139 (MQ=255)                                                                                                                 aCCATGTGTCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGt  <  1:155425/137‑1 (MQ=255)                                                                                                                      gtgtCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGAc                                                  <  1:204511/84‑1 (MQ=255)                                                                                                                      gtgtCGGGGAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGAc                                                  >  2:204511/1‑84 (MQ=255)                                                                                                                              |                                                                                                                            TTCAGTATCGTCTTTTGAGATTCTCAAGGTTATGGAGTAACCAGACAGAACGGACACCACTGTCGGAGATAAAAACCCCCTGTTTACGCAGCTCGTTGCTGGCCCGGTGCTGACCATGTGTCGGGAAAGCGACGGCGTAATCAACAACAGCCTGTTCAGTTGCCTCATCGGTACGGTTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGT  >  NZ_CP009273/525100‑525348

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