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%	60.4 / NA	19	Q18Q (CAA→CAG)	uspE	universal stress protein UspE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (3/16);  total (3/16)

CGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTT  >  NZ_CP009273/1392741‑1392892                                                                                      |                                                                   cGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggctggt                                                                <  2:104066/90‑1 (MQ=255)      gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                           <  1:394708/90‑1 (MQ=255)      gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                           <  2:279578/90‑1 (MQ=255)      gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                           <  2:28087/90‑1 (MQ=255)           gTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGtt                                                      >  2:390526/1‑90 (MQ=255)                        aaaaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaaca                                         <  1:451619/90‑1 (MQ=255)                          aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                       <  2:182060/90‑1 (MQ=255)                          aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                       <  1:108121/90‑1 (MQ=255)                                tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                           <  2:354769/64‑1 (MQ=255)                                tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                           >  1:354769/1‑64 (MQ=255)                                  tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                               <  1:224852/90‑1 (MQ=255)                                  tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                               <  2:327391/90‑1 (MQ=255)                                            accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                     <  1:390526/90‑1 (MQ=255)                                                   cGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATa              <  2:484959/90‑1 (MQ=255)                                                   cGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATa              <  2:218306/90‑1 (MQ=255)                                                          gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc       <  2:282443/90‑1 (MQ=255)                                                             aataaaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTcc    <  1:239900/90‑1 (MQ=255)                                                               taaaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCtt  <  2:87438/90‑1 (MQ=255)                                                                    aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                  <  2:182175/53‑1 (MQ=255)                                                                    aGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                  >  1:182175/1‑53 (MQ=255)                                                                                      |                                                                   CGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTT  >  NZ_CP009273/1392741‑1392892

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