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	85.7%	43.7 / 3.3	21	intergenic (‑361/+185)	tomB/acrB	Hha toxicity modulator TomB/efflux RND transporter permease AcrB
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (10/8);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTT  >  NZ_CP009273/476437‑476610                                                                                         |                                                                                      ccAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGat                                                                                      <  1:513283/90‑1 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAAt                                                                              >  1:225698/1‑90 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAAt                                                                              <  2:92255/90‑1 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAAt                                                                              >  2:553247/1‑90 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAAt                                                                              >  2:351177/1‑90 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAAt                                                                              >  2:444885/1‑90 (MQ=255)         gCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAAt                                                                              >  2:403907/1‑90 (MQ=255)              gAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTAtt                                                                         <  2:613829/90‑1 (MQ=255)                   tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTGTTAAATCTATTAgcgc                                                                    <  1:177197/90‑1 (MQ=255)                   tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                    <  1:646491/90‑1 (MQ=255)                   tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                    <  1:226849/90‑1 (MQ=255)                   tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                    <  2:139740/90‑1 (MQ=255)                   tGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAgcgc                                                                    <  2:210778/90‑1 (MQ=255)                        taataaCCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACtt                                                               >  1:455912/1‑90 (MQ=255)                                                        aaaTAATTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTa                               >  1:765465/1‑90 (MQ=255)                                                          ataatTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATg                             <  2:536995/90‑1 (MQ=255)                                                             atTAAGAATAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAAt                          >  1:767746/1‑90 (MQ=255)                                                                     tAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCg                  >  1:675543/1‑90 (MQ=255)                                                                     tAGCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCg                  <  1:329707/90‑1 (MQ=255)                                                                       gCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTg                >  2:142351/1‑90 (MQ=255)                                                                       gCCTCTAAATGATTATGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTg                >  1:44598/1‑90 (MQ=255)                                                                                     aTGGGTTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGttt  >  1:324594/1‑90 (MQ=255)                                                                                         |                                                                                      CCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTT  >  NZ_CP009273/476437‑476610

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