Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,400,314	Δ1 bp	coding (1460/1527 nt)	abgT ←	p‑aminobenzoyl glutamate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,400,314	0	A	.	100.0%	116.8 / NA	25	coding (1460/1527 nt)	abgT	p‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/15);  total (10/15)

GGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCA  >  NC_000913/1400272‑1400350                                           |                                     ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                             <  1:883273/50‑1 (MQ=255) ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                             <  1:2668645/50‑1 (MQ=255)  gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                            <  1:525339/50‑1 (MQ=255)   acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                           <  1:292661‑M1/49‑2 (MQ=255)     cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                         <  1:965993‑M1/50‑4 (MQ=255)       aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                        >  1:2773924‑M1/1‑45 (MQ=255)        tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                       <  1:1170186‑M1/49‑6 (MQ=255)        tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                       <  1:2549686‑M1/49‑6 (MQ=255)             aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                 <  1:1646428‑M1/50‑12 (MQ=255)               aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa               >  1:1337393‑M1/1‑37 (MQ=255)                ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat              <  1:1040480‑M1/50‑15 (MQ=255)                ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat              <  1:54904‑M1/50‑15 (MQ=255)                 cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg             >  1:864076‑M1/1‑35 (MQ=255)                  gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc            >  1:2613536‑M1/1‑34 (MQ=255)                  gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc            <  1:893757‑M1/50‑17 (MQ=255)                  gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc            >  1:907883‑M1/1‑34 (MQ=255)                   aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:247791‑M1/1‑33 (MQ=255)                   aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:1848530‑M1/1‑33 (MQ=255)                   aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:846669‑M1/1‑33 (MQ=255)                    cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg          <  1:703964‑M1/50‑19 (MQ=255)                     aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt         <  1:924144‑M1/50‑20 (MQ=255)                       gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc        >  1:2358565‑M1/1‑29 (MQ=255)                           ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg   <  1:1225527‑M1/50‑26 (MQ=255)                           ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg   <  1:1943562‑M1/50‑26 (MQ=255)                            cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc  >  1:981657‑M1/1‑24 (MQ=255)                                           |                                     GGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCA  >  NC_000913/1400272‑1400350

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