Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,647,969	G→A	Y265Y (TAC→TAT)	oprB‑II ←	carbohydrate‑selective porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,647,969	0	G	A	100.0%	65.7 / NA	20	Y265Y (TAC→TAT)	oprB‑II	carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/8);  total (12/8)

CTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTG  >  NC_002947/1647846‑1648106                                                                                                                            |                                                                                                                                           ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCttt                                                                                                                  >  1:404811/1‑150 (MQ=255)                 aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑t                                                                                                 <  1:426488/150‑1 (MQ=255)                 aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑t                                                                                                 <  2:702726/150‑1 (MQ=255)                 aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGACTTGGATACG‑t                                                                                                 <  1:751833/150‑1 (MQ=255)                             tgcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc                                                                                                                                     >  1:183314/1‑103 (MQ=255)                             tgcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc                                                                                                                                     <  2:183314/103‑1 (MQ=255)                               cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAg                                                                                   <  2:404811/150‑1 (MQ=255)                                     cGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGaa                                                                             <  1:456179/150‑1 (MQ=255)                                                          gccgcTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCAGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACGACTTCGAACAGCATTAGGGCACACTtggg                                                        >  1:612055/1‑149 (MQ=255)                                                                           gcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCt                                       >  2:793998/1‑150 (MQ=255)                                                                                           cTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                       >  2:177561/1‑150 (MQ=255)                                                                                            ttGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAt                      >  2:180255/1‑150 (MQ=255)                                                                                             tGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATCAATCTTTAACATTGTGTGCACGTTGTGGTCC‑TACCAGTTCGACCGGAATTCCCGCAACCTCCTTGCTATTGCCTCTGAGCTTGAAACCATTGCCGGAtt                     <  1:889258/150‑1 (MQ=255)                                                                                                gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACTGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTTGa                  >  2:819633/1‑150 (MQ=255)                                                                                                gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa                  >  2:723161/1‑150 (MQ=255)                                                                                                gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa                  >  1:706759/1‑150 (MQ=255)                                                                                                           gCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGAGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAggg       <  1:793998/150‑1 (MQ=255)                                                                                                            ccTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGt      >  2:101098/1‑150 (MQ=255)                                                                                                             cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt     >  2:20286/1‑150 (MQ=255)                                                                                                                ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTg  >  2:169354/1‑150 (MQ=255)                                                                                                                            |                                                                                                                                           CTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACG‑TGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTG  >  NC_002947/1647846‑1648106

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