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%	67.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 (7/14);  total (7/14)

ACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153730‑4153894                                                                                       |                                                                               aCCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCAt                                                                             <  1:584517/90‑1 (MQ=255)       ggCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTAttt                                                                       >  1:534374/1‑90 (MQ=255)           gTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTGAGGGCAAGGCCCAGCTAttttttt                                                                   >  1:698463/1‑88 (MQ=255)                ccTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgcg                                                              <  2:264542/90‑1 (MQ=255)                     gTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCCCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACa                                                         <  1:195031/90‑1 (MQ=255)                      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                        <  2:74906/90‑1 (MQ=255)                         atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGGAGGCCCCATCTATTTTTTTTTGCGGGAaaaaag                                                     >  2:137799/1‑88 (MQ=255)                         atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGTCAGCTCCCCCCTTTTTTTTTTTGCGGGACaaaac                                                     >  1:277650/1‑88 (MQ=255)                          tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                    <  2:198545/90‑1 (MQ=255)                           gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                   <  1:322018/90‑1 (MQ=255)                              gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGGACAAATGCCAGt                                                >  2:422427/1‑90 (MQ=255)                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGATCAGGTCAGCTCCCCTCTATTTTTTTTCCCCGGACAAAAGCCAGGCAtgt                                           >  2:519890/1‑90 (MQ=255)                                     gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                         <  1:510879/90‑1 (MQ=255)                                       cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                       <  1:711168/90‑1 (MQ=255)                                          cTTCCGGGCGTCGATATCCCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtg                                    <  2:666125/90‑1 (MQ=255)                                           ttCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGgtgt                                   <  1:130196/90‑1 (MQ=255)                                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                           <  2:247283/90‑1 (MQ=255)                                                          tCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTAc                    <  2:534374/90‑1 (MQ=255)                                                          tCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATCGATGGCGTAc                    <  2:402799/90‑1 (MQ=255)                                                          tCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAAAGCCAGTCATGTGTTGGGGTTAATTGATGGCGTAc                    >  2:762033/1‑90 (MQ=255)                                                                            cAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggggg  <  2:356995/90‑1 (MQ=255)                                                                                       |                                                                               ACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153730‑4153894

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