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%	58.1 / NA	20	S84P (TCA→CCA)	btuB	TonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/9);  total (11/9)

CGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCG  >  NZ_CP009273/4153732‑4153890                                                                                     |                                                                           cGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCctct                                                                       >  2:510646/1‑90 (MQ=255)  gCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCctctt                                                                      >  2:418238/1‑89 (MQ=255)  gCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGTCAACTCCCCTCTa                                                                      >  1:404295/1‑90 (MQ=255)              ccTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgcg                                                          <  1:281366/90‑1 (MQ=255)                       atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATc                                                                        >  1:24312/1‑67 (MQ=255)                       atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATc                                                                        <  2:24312/67‑1 (MQ=255)                            gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                            >  2:550997/1‑90 (MQ=255)                                 gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGTCCGCTTCCCTCTTTTTTTATTCGCGGGACAAATGGCCGGCAtgt                                       >  1:143151/1‑90 (MQ=255)                                 gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTCGCGGTACAAAAGCCAGTCAtgg                                       >  1:521777/1‑89 (MQ=255)                                 gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAAAGCCAGTCAtgt                                       >  2:350193/1‑90 (MQ=255)                                   gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                     <  1:422146/90‑1 (MQ=255)                                    ccGTTTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGtt                                    <  1:36517/90‑1 (MQ=255)                                    ccGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGtt                                    <  2:207741/90‑1 (MQ=255)                                         ttCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtgt                               <  2:85955/90‑1 (MQ=255)                                         ttCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtgt                               <  1:98688/90‑1 (MQ=255)                                                 gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                       <  2:407027/90‑1 (MQ=255)                                                  tCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGGTg                      >  2:216851/1‑90 (MQ=255)                                                        tCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTTATGGCGTAc                >  1:39051/1‑90 (MQ=255)                                                               aaaCGGCGGTTCAGGTCAGGTCCCATCTATTTTTTTTTGCGGGACAAATGCCCGGCATGTGTTGGGGTTAATTGATTGGGGACGCCTGaa         >  1:38782/1‑90 (MQ=255)                                                                      ggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCg  <  1:86576/90‑1 (MQ=255)                                                                                     |                                                                           CGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCG  >  NZ_CP009273/4153732‑4153890

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