Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,392,826	T→C	Q18Q (CAA→CAG)	uspE ←	universal stress protein UspE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,392,826	0	T	C	100.0%	69.7 / NA	22	Q18Q (CAA→CAG)	uspE	universal stress protein UspE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/16);  total (6/16)

GAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAA  >  NZ_CP009273/1392746‑1392899                                                                                 |                                                                          gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  <  2:82354/90‑1 (MQ=255) gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  <  2:38827/90‑1 (MQ=255) gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  <  1:305032/90‑1 (MQ=255)     aGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGt                                                              >  1:531548/1‑90 (MQ=255)        cATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAg                                                           <  1:301675/90‑1 (MQ=255)                ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  <  1:333888/75‑1 (MQ=255)                ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  <  1:450187/75‑1 (MQ=255)                ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  >  2:333888/1‑75 (MQ=255)                ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                  >  2:450187/1‑75 (MQ=255)                     aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                              <  1:535950/90‑1 (MQ=255)                           tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                                        <  2:531548/90‑1 (MQ=255)                           tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                                        <  1:221027/90‑1 (MQ=255)                                ttttGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCt                                   <  1:403730/90‑1 (MQ=255)                                  ttGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGa                                 >  1:482360/1‑90 (MQ=255)                                     ccaccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATAc                              <  2:353599/90‑1 (MQ=255)                                       accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                            <  2:282840/90‑1 (MQ=255)                                       accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                            <  2:482360/90‑1 (MQ=255)                                                     gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc              <  1:76633/90‑1 (MQ=255)                                                     gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc              <  1:235420/90‑1 (MQ=255)                                                               aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCtt         <  2:326798/85‑1 (MQ=255)                                                               aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCtt         >  1:326798/1‑85 (MQ=255)                                                                 cTCGCCGCAATGCTGGCTGGTCGTCGTGGTTAGGATCGATAACAACGAGCATGTGCTGATACATAGCCATACAGGGGCTACTTActacaa  >  2:73214/1‑90 (MQ=255)                                                                                 |                                                                          GAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAA  >  NZ_CP009273/1392746‑1392899

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