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

ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTG  >  NZ_CP009273/4153744‑4153896                                                                         |                                                                                 aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                 <  1:172697/90‑1 (MQ=255)       gTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACa                                                           <  1:335848/90‑1 (MQ=255)            tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCGCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                      <  1:87550/90‑1 (MQ=255)             gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                     <  1:263141/90‑1 (MQ=255)                  gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                <  1:208938/90‑1 (MQ=255)                  gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                <  1:50606/90‑1 (MQ=255)                     gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGTCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCCGTCAtgt                                             >  1:492914/1‑90 (MQ=255)                        ccGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGtt                                          <  1:573707/90‑1 (MQ=255)                         cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                         <  1:349505/90‑1 (MQ=255)                           tCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgt                                       <  2:145252/90‑1 (MQ=255)                                    cGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAAttta                              >  2:332326/1‑88 (MQ=255)                                     gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                             <  2:482569/90‑1 (MQ=255)                                          tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                        <  1:332326/90‑1 (MQ=255)                                          tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                        <  1:444103/90‑1 (MQ=255)                                            tCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTATTTTTTTTCGCGGGACAAATGCCAGTCATGTGTTGGTGTTAATTTATTGCGTAc                      >  1:482569/1‑90 (MQ=255)                                                ccAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCt                  <  1:168443/90‑1 (MQ=255)                                                   aaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTTTCGGGACAAAAGCCCGGCATGTGTTGGGGTTAATTGATGGCGTACGCCTGaa               >  1:353367/1‑90 (MQ=255)                                                            ttCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggg      <  1:574243/90‑1 (MQ=255)                                                                ggTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggggggg  >  2:68407/1‑88 (MQ=255)                                                                         |                                                                                 ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTG  >  NZ_CP009273/4153744‑4153896

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