Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,400,309	Δ1 bp	coding (1465/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,308	0	A	.	100.0%	117.7 / NA	25	coding (1466/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)

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCA  >  NC_000913/1400259‑1400350                                                  |                                           tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                            >  1:1632182/1‑50 (MQ=255) tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                            >  1:1735071/1‑50 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                          <  1:983967/50‑3 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCac                                           <  1:741333/49‑2 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCac                                           <  1:584432/49‑2 (MQ=255)              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:907883‑M1/1‑34 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc            <  1:893757‑M1/50‑17 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:247791‑M1/1‑33 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:846669‑M1/1‑33 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct           >  1:1848530‑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)                                                  |                                           TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCA  >  NC_000913/1400259‑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: ‑