Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,153,816	T→C	S84P (TCA→CCA)	btuB →	TonB‑dependent vitamin B12 receptor BtuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,153,816	0	T	C	100.0%	68.1 / NA	21	S84P (TCA→CCA)	btuB	TonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/16);  total (5/16)

CTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGT  >  NZ_CP009273/4153746‑4153902                                                                       |                                                                                       cTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCgg                                                                     <  2:543429/90‑1 (MQ=255)      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                                <  1:120600/90‑1 (MQ=255)      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                                <  2:177090/90‑1 (MQ=255)         atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTTAGGTCAGCTCCCCTCTATTTTTTTTTGCGGGAaaaaag                                                             >  1:23999/1‑88 (MQ=255)          tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                            <  1:682014/90‑1 (MQ=255)                gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                      <  1:568768/90‑1 (MQ=255)                     gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                                 <  1:494799/90‑1 (MQ=255)                       cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                               <  2:163879/90‑1 (MQ=255)                            tCCGGGCGTCGATATCCCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGtt                                          <  1:461088/90‑1 (MQ=255)                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                   <  2:104693/90‑1 (MQ=255)                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                   <  2:47861/90‑1 (MQ=255)                                        tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                              <  2:665187/90‑1 (MQ=255)                                        tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                              <  2:647994/90‑1 (MQ=255)                                        tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                              <  2:343083/90‑1 (MQ=255)                                         atCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTa                             <  2:706033/90‑1 (MQ=255)                                               cAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTCGCGGGACAAATGCCAGTCATGTGTTGGGGTTAATTGATGGCGGACGCCTg                       >  1:26645/1‑90 (MQ=255)                                                aaaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtg                                            >  2:145226/1‑68 (MQ=255)                                                aaaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtg                                            <  1:145226/68‑1 (MQ=255)                                                       cggTTCAGGTCAGCTCCCTTCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGc               <  2:437002/90‑1 (MQ=255)                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg   >  2:154648/1‑90 (MQ=255)                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg   >  2:622263/1‑90 (MQ=255)                                                                    ctccCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGt  <  2:199997/90‑1 (MQ=255)                                                                       |                                                                                       CTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGT  >  NZ_CP009273/4153746‑4153902

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