Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,103,871	(A)7→6	intergenic (+71/+78)	BW25113_RS15410 → / ← yghD	tRNA‑Phe/GspM family type II secretion system protein YghD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,103,865	0	A	.	100.0%	107.9 / NA	26	intergenic (+65/+84)	BW25113_RS15410/yghD	tRNA‑Phe/GspM family type II secretion system protein YghD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (9/17);  total (9/17)

GATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103781‑3103949                                                                                     |                                                                                     gATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaa                                                                                 >  2:683389/1‑90 (MQ=255)            cGGGCACCACTAATTCTTAAGAACCCTCCCACAAGGCGGGTTTTTGCTTTTGGTTCTGACAATAACCTTCACG‑AAAAAATTAGCTTATaa                                                                     <  1:246692/90‑1 (MQ=255)             gggCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATaaa                                                                    >  2:407251/1‑90 (MQ=255)              ggCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAg                                                                   <  2:365977/90‑1 (MQ=255)                  ccacTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTg                                                               <  1:407251/90‑1 (MQ=255)                              aaGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTActc                                                   <  1:250971/90‑1 (MQ=255)                                  aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:428489/90‑1 (MQ=255)                                  aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:772635/90‑1 (MQ=255)                                   cccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGcc                                              >  2:65371/1‑90 (MQ=255)                                    ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:38909/90‑1 (MQ=255)                                            aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                                     <  1:683389/90‑1 (MQ=255)                                              ggCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACgag                                   <  2:193351/90‑1 (MQ=255)                                                 gggTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  1:99383/90‑1 (MQ=255)                                                 gggTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  1:95656/90‑1 (MQ=255)                                                 gggTTTTTGCTTTTGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCACGCCACGTTAACGAGAGt                                >  1:358847/1‑90 (MQ=255)                                                             ttGGATCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGat                    <  2:492940/90‑1 (MQ=255)                                                                 aTCTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTa                <  1:489655/90‑1 (MQ=255)                                                                   cTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  1:159924/90‑1 (MQ=255)                                                                   cTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:232011/90‑1 (MQ=255)                                                                   cTGACAATAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:79434/90‑1 (MQ=255)                                                                         aTAACCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTg        >  2:95921/1‑90 (MQ=255)                                                                           aaCCTTCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:81859/90‑1 (MQ=255)                                                                               ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:336064/1‑90 (MQ=255)                                                                               ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:346343/1‑90 (MQ=255)                                                                               ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  1:758679/1‑90 (MQ=255)                                                                               ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  1:509084/1‑90 (MQ=255)                                                                               ttCACG‑AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:600806/1‑90 (MQ=255)                                                                                     |                                                                                     GATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103781‑3103949

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