Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,898,828	(T)5→4	coding (130/1878 nt)	bglF ←	PTS beta‑glucoside transporter subunit IIABC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,898,824	0	T	.	100.0%	86.7 / NA	21	coding (134/1878 nt)	bglF	PTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/9);  total (12/9)

AGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898743‑3898907                                                                                  |                                                                                    aGGAAGACATCGGCCACATGGTTACCTATGACCCCCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTa                                                                            <  1:329244/90‑1 (MQ=255)        cATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGc                                                                     >  2:92968/1‑90 (MQ=255)           cGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTg                                                                  <  2:281824/90‑1 (MQ=255)                acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCgctt                                                             <  1:379742/90‑1 (MQ=255)                             tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtt                                                >  1:14862/1‑90 (MQ=255)                                  accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt                                           >  2:253551/1‑90 (MQ=255)                                   ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          >  1:231920/1‑90 (MQ=255)                                             gCCACCGCTTTCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                <  1:486905/90‑1 (MQ=255)                                                      ttCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                       >  1:426643/1‑90 (MQ=255)                                                       tCCACCACCATAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      <  2:26858/90‑1 (MQ=255)                                                         caccaccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:49908/90‑1 (MQ=255)                                                             accaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCa                >  1:219231/1‑90 (MQ=255)                                                              ccaTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt               >  2:127060/1‑90 (MQ=255)                                                               caTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATc              <  1:253551/90‑1 (MQ=255)                                                               caTAATAATACCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATc              <  2:110335/90‑1 (MQ=255)                                                                      ataCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:14862/90‑1 (MQ=255)                                                                        aCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:229004/1‑90 (MQ=255)                                                                        aCCGGGGGTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:12365/1‑90 (MQ=255)                                                                           gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:113041/1‑90 (MQ=255)                                                                           gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:122550/1‑90 (MQ=255)                                                                           gggggTC‑TTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:195383/1‑90 (MQ=255)                                                                                  |                                                                                    AGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898743‑3898907

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