Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,211,818	A→C	C73W (TGT→TGG)	yehY ←	glycine betaine ABC transporter permease YehY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,211,818	0	A	C	100.0%	78.4 / NA	21	C73W (TGT→TGG)	yehY	glycine betaine ABC transporter permease YehY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/9);  total (12/9)

CCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCACAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACC  >  NZ_CP009273/2211698‑2211931                                                                                                                          |                                                                                                                  ccagAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCtttttt                                                                                                  >  2:121707/1‑139 (MQ=255)   agAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTAcc                                                                                               <  1:408377/139‑1 (MQ=255)       aCCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGa                                                                                           <  1:476600/139‑1 (MQ=255)                                  cgcACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGgcgc                                                                  >  2:292250/1‑138 (MQ=255)                                         ccATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGAGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACg                                                          >  2:551189/1‑139 (MQ=255)                                                 cgcCAGTTGGGTCGCCGCCTTTCCAGCTCCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAg                                                  >  2:552006/1‑139 (MQ=255)                                                   ccAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccacc                                                    <  2:401222/134‑1 (MQ=255)                                                   ccAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccacc                                                    >  1:401222/1‑134 (MQ=255)                                                              gccgccTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccac                                                     <  1:321193/122‑1 (MQ=255)                                                              gccgccTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccac                                                     >  2:321193/1‑122 (MQ=255)                                                                        cAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCcaca                           <  1:323295/138‑1 (MQ=255)                                                                                      cAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGa             <  1:552006/139‑1 (MQ=255)                                                                                      cAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGa             >  2:430838/1‑139 (MQ=255)                                                                                        gCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccacca                                                   <  2:358948/99‑1 (MQ=255)                                                                                        gCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCccacca                                                   >  1:358948/1‑99 (MQ=255)                                                                                        gCAATACGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACg           >  1:124614/1‑139 (MQ=255)                                                                                              cGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAg                          <  1:532855/118‑1 (MQ=255)                                                                                              cGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAg                          >  2:532855/1‑118 (MQ=255)                                                                                              cGAAGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCt      >  1:389086/1‑138 (MQ=255)                                                                                                 aGACGAATTGCGCCAGAATGAGTGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCAcc  >  2:447049/1‑139 (MQ=255)                                                                                                 aGACGAATTGCGCCAGAATGAGAGCCCAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCAcc  <  1:551189/139‑1 (MQ=255)                                                                                                                          |                                                                                                                  CCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCT‑CCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCACAAATGCTGCCTTTTTTACCGGGAATAAAACAGGCCGTCAGCCAGGCGCAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACC  >  NZ_CP009273/2211698‑2211931

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