Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	476,525	A→G	intergenic (‑361/+185)	tomB ← / ← acrB	Hha toxicity modulator TomB/efflux RND transporter permease AcrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	476,525	0	A	G	100.0%	110.3 / NA	30	intergenic (‑361/+185)	tomB/acrB	Hha toxicity modulator TomB/efflux RND transporter permease AcrB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/17);  total (13/17)

AGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGT  >  NZ_CP009273/476439‑476608                                                                                       |                                                                                    aGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgtta                                                                                  >  2:2191239/1‑90 (MQ=255) aGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGgtta                                                                                  >  1:2306831/1‑90 (MQ=255)       gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAAt                                                                            >  1:2108058/1‑90 (MQ=255)               tatGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgc                                                                    <  2:1100884/90‑1 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  <  2:506379/90‑1 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  <  2:1550328/90‑1 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  <  1:2544782/90‑1 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  <  2:1115168/90‑1 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  >  1:295344/1‑90 (MQ=255)                 tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  >  2:110590/1‑90 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  >  1:1198225/1‑85 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                  <  2:1198225/85‑1 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACtt                                                             >  2:83918/1‑90 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACtt                                                             >  2:778633/1‑90 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACtt                                                             >  2:417797/1‑90 (MQ=255)                      taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACtt                                                             >  2:2105609/1‑90 (MQ=255)                         taaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAat                                                          <  1:1380170/90‑1 (MQ=255)                               aGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAAt                                                    <  2:2108058/90‑1 (MQ=255)                                 gaatgaatAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATta                                                  <  1:486254/90‑1 (MQ=255)                                          aaGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAAt                                         <  2:2945041/90‑1 (MQ=255)                                            gAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCa                                       <  2:295344/90‑1 (MQ=255)                                               ttACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTa                                    >  1:2370396/1‑90 (MQ=255)                                                tACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGttatta                                                                               <  1:2753889/46‑1 (MQ=255)                                                tACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGttatta                                                                               >  2:2753889/1‑46 (MQ=255)                                                       aataatTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTAt                            <  1:1086532/90‑1 (MQ=255)                                                                   tAGCCTCTAAATGATTATGGGTTATTGAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCg                <  2:133037/90‑1 (MQ=255)                                                                   tAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCg                <  2:722680/90‑1 (MQ=255)                                                                     gCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTg              >  2:2468211/1‑90 (MQ=255)                                                                      ccTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTgg             <  1:417797/90‑1 (MQ=255)                                                                                aTTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTgg   >  1:2141623/1‑90 (MQ=255)                                                                                aTTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTgg   <  1:410412/90‑1 (MQ=255)                                                                                 ttATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGt  >  1:3002337/1‑90 (MQ=255)                                                                                       |                                                                                    AGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGT  >  NZ_CP009273/476439‑476608

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