Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,491,710	A→G	intergenic (+304/+210)	BW25113_RS24825 → / ← insG	transposase domain‑containing protein/IS4‑like element IS4 family transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,491,710	0	A	G	100.0%	114.2 / NA	32	intergenic (+304/+210)	BW25113_RS24825/insG	transposase domain‑containing protein/IS4‑like element IS4 family transposase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/17);  total (15/17)

GAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACA  >  NZ_CP009273/4491626‑4491792                                                                                     |                                                                                   gagTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAg                                                                               >  2:249154/1‑90 (MQ=255)             tCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCaa                                                                   <  2:136419/90‑1 (MQ=255)              cATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAAc                                                                  <  2:326312/90‑1 (MQ=255)                      cTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGAt                                                          >  2:94054/1‑90 (MQ=255)                      cTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGAt                                                          >  2:254598/1‑90 (MQ=255)                        gTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGc                                                        >  2:78517/1‑90 (MQ=255)                        gTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGc                                                        >  2:426734/1‑90 (MQ=255)                        gTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGc                                                        >  2:201515/1‑90 (MQ=255)                          cGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCaa                                                      >  1:119557/1‑90 (MQ=255)                           ggTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAg                                                     <  1:180390/90‑1 (MQ=255)                           ggTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAg                                                     >  2:411440/1‑90 (MQ=255)                            gTGATGCGCTGAAAAATGGGTTAATGTTGAGTTCCTCAAACCTTTACCAATTAAGGGGCGTAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  2:46071/90‑1 (MQ=255)                            gTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  2:255780/90‑1 (MQ=255)                            gTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  1:110117/90‑1 (MQ=255)                            gTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  1:254598/90‑1 (MQ=255)                                  cgcTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTa                                              >  2:197965/1‑90 (MQ=255)                                   gcTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTaa                                             <  2:338532/90‑1 (MQ=255)                                     tAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATg                                           <  1:426734/90‑1 (MQ=255)                                                     gTTGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTgg                           <  1:26869/90‑1 (MQ=255)                                                       tGAGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGcc                         >  1:56373/1‑90 (MQ=255)                                                         aGTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTg                       <  2:5690/90‑1 (MQ=255)                                                          gTTCCTCGAACCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGc                      <  1:312569/90‑1 (MQ=255)                                                                  aaCCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  2:478922/52‑1 (MQ=255)                                                                  aaCCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    >  1:478922/1‑52 (MQ=255)                                                                   aCCCTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGAc             <  1:64197/90‑1 (MQ=255)                                                                     ccTTACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACaa           <  1:201515/90‑1 (MQ=255)                                                                       ttACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATg         >  2:38963/1‑90 (MQ=255)                                                                        tACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    <  2:466925/46‑1 (MQ=255)                                                                        tACCAATTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGc                                                    >  1:466925/1‑46 (MQ=255)                                                                             aTTAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAAc   <  1:249154/90‑1 (MQ=255)                                                                              ttAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACa  >  1:410745/1‑90 (MQ=255)                                                                              ttAAGGGGCATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACa  >  1:465018/1‑90 (MQ=255)                                                                                     |                                                                                   GAGTGCAATCGTTCATATACACTGTCGGTGATGCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACA  >  NZ_CP009273/4491626‑4491792

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: ‑