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%	159.2 / NA	33	coding (1460/1527 nt)	abgT	p‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (15/18);  total (15/18)

CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353                                             |                                        ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:1969525/1‑50 (MQ=255)  cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:1201375/1‑49 (MQ=38)  cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:408110/1‑49 (MQ=38)   ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:443757/1‑50 (MQ=255)   ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:2563646/1‑50 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:787935‑M1/50‑4 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:3722896‑M1/1‑47 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:3416329‑M1/50‑4 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:2878592‑M1/50‑4 (MQ=255)        tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2328874‑M1/1‑46 (MQ=255)          tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:256766‑M1/50‑7 (MQ=255)          tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:3259557‑M1/48‑5 (MQ=255)           cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2087253‑M1/1‑43 (MQ=255)             gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:2389698‑M1/50‑10 (MQ=255)             gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:1382131‑M1/50‑10 (MQ=255)             gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:372896‑M1/1‑41 (MQ=255)              cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:2289826‑M1/50‑11 (MQ=255)               aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:3647589‑M1/50‑12 (MQ=255)               aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:3291792‑M1/50‑12 (MQ=255)                gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:1370689‑M1/1‑38 (MQ=255)                 aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  <  1:3878157‑M1/50‑14 (MQ=255)                  ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:2270981‑M1/50‑15 (MQ=255)                   cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:1146291‑M1/1‑35 (MQ=255)                    gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:1891263‑M1/1‑34 (MQ=255)                    gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:3946089‑M1/50‑17 (MQ=255)                       aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:2274230‑M1/50‑20 (MQ=255)                        aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:3581522‑M1/50‑21 (MQ=255)                         gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:2142520‑M1/50‑22 (MQ=255)                            aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:2767018‑M1/1‑26 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:120742‑M1/1‑23 (MQ=255)                                cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:414077‑M1/50‑29 (MQ=255)                                 gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:2893739‑M1/1‑21 (MQ=255)                                 gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:1976710‑M1/50‑30 (MQ=255)                                             |                                        CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353

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