Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	510,152	Δ1 bp	coding (164/459 nt)	ybbJ ←	NfeD family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	510,151	0	C	.	100.0%	85.5 / NA	21	coding (165/459 nt)	ybbJ	NfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/10);  total (11/10)

TCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTT  >  NZ_CP009273/510074‑510232                                                                              |                                                                                  tCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCa                                                                      >  2:976146/1‑90 (MQ=255) tCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCa                                                                      >  2:126856/1‑90 (MQ=255)          cgccgcGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAcca                                                             <  1:126856/90‑1 (MQ=255)                gACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGa                                                       >  2:1087490/1‑90 (MQ=255)                    aCCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGaccac                                                   >  2:54101/1‑90 (MQ=255)                         ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                              <  1:1087490/90‑1 (MQ=255)                         ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                              <  1:54101/90‑1 (MQ=255)                         ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                              <  1:904408/90‑1 (MQ=255)                                    ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                   <  1:129321/90‑1 (MQ=255)                                    ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                   <  1:328110/90‑1 (MQ=255)                                      aGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGc                                 >  2:686856/1‑90 (MQ=255)                                               gcagcGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                        <  2:736283/90‑1 (MQ=255)                                                  gcGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCt                     >  2:41145/1‑90 (MQ=255)                                                    gTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTcc                   >  2:999845/1‑90 (MQ=255)                                                        ggATGGCAAACATTCCCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaa               <  1:249917/90‑1 (MQ=255)                                                            ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa           <  2:379594/90‑1 (MQ=255)                                                            ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa           <  2:633110/90‑1 (MQ=255)                                                            ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa           >  1:776811/1‑90 (MQ=255)                                                            ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa           >  1:1150182/1‑90 (MQ=255)                                                             gCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt          >  2:30893/1‑90 (MQ=255)                                                                     tACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAttt  >  1:126242/1‑90 (MQ=255)                                                                              |                                                                                  TCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTT  >  NZ_CP009273/510074‑510232

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