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

CCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153731‑4153901                                                                                      |                                                                                      ccGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATc                                                                                   <  2:242559/90‑1 (MQ=255)            cGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTc                                                                        <  1:375264/90‑1 (MQ=255)              aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  1:400628/90‑1 (MQ=255)               ccTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgcg                                                                     <  1:90985/90‑1 (MQ=255)                    gTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCCGCCCCCATCTATTTTTTTTCGCGGGACa                                                                >  2:342791/1‑90 (MQ=255)                        atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATg                                                            >  1:157668/1‑90 (MQ=255)                          gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:182386/90‑1 (MQ=255)                          gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:74317/90‑1 (MQ=255)                             gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                                       <  1:358964/90‑1 (MQ=255)                             gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGTACAAAAGCCAgg                                                       >  2:48718/1‑89 (MQ=255)                               gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                     <  1:158291/90‑1 (MQ=255)                                  gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTTCGGGACAAAAGCCAGTCAtgt                                                  >  2:182386/1‑90 (MQ=255)                                    gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                                <  1:131176/90‑1 (MQ=255)                                      cGTCTTCCGGGCGTCGATATCCCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  2:263973/90‑1 (MQ=255)                                      cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  2:396701/90‑1 (MQ=255)                                      cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  1:229753/90‑1 (MQ=255)                                      cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  1:84350/90‑1 (MQ=255)                                      cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  1:153970/90‑1 (MQ=255)                                                  gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                  <  1:53129/90‑1 (MQ=255)                                                  gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                  <  2:383617/90‑1 (MQ=255)                                                       tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                             <  1:48718/90‑1 (MQ=255)                                                       tatCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGt                             <  1:353396/90‑1 (MQ=255)                                                           aCCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGACAAAAGCCAGTCATGTGTTGGGGTTAATTTATGGCGTACGc                         >  2:79112/1‑90 (MQ=255)                                                                     gcggTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGACAAAAGCCAGTCATGTGTTGGGGTTAATTGATGGCGTACGCCTGAATCTgg               >  2:35117/1‑90 (MQ=255)                                                                     gcggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTgg               >  1:158241/1‑90 (MQ=255)                                                                                  gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  2:176228/1‑90 (MQ=255)                                                                                  gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  1:320740/1‑90 (MQ=255)                                                                                      |                                                                                      CCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153731‑4153901

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