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%	100.9 / NA	22	coding (1466/1527 nt)	abgT	p‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (12/10);  total (12/10)

GTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400262‑1400352                                               |                                             gTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                            >  1:982712/1‑47 (MQ=255)           ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                 >  1:2326216/1‑48 (MQ=37)             aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              <  1:317219/49‑1 (MQ=38)             aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              <  1:1876215/49‑1 (MQ=38)                tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           >  1:2212456‑M1/1‑46 (MQ=255)                  tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          <  1:249733‑M1/48‑5 (MQ=255)                  tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           <  1:1195296‑M1/47‑4 (MQ=255)                    ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                        <  1:260613‑M1/48‑7 (MQ=255)                    ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:1293715‑M1/1‑42 (MQ=255)                           cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                <  1:1969899‑M1/49‑15 (MQ=255)                            gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg               >  1:38669‑M1/1‑34 (MQ=255)                             aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc              >  1:621666‑M1/1‑33 (MQ=255)                              cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct             >  1:836836‑M1/1‑32 (MQ=255)                                aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           <  1:2669328‑M1/49‑20 (MQ=255)                                aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           >  1:1432615‑M1/1‑30 (MQ=255)                                aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           <  1:785059‑M1/49‑20 (MQ=255)                                 gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc          >  1:1569442‑M1/1‑29 (MQ=255)                                 gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc          <  1:28007‑M1/49‑21 (MQ=255)                                     ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      >  1:2567507‑M1/1‑25 (MQ=255)                                      cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg     >  1:544897‑M1/1‑24 (MQ=255)                                      cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg     >  1:732844‑M1/1‑24 (MQ=255)                                       aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc    <  1:2372135‑M1/49‑27 (MQ=255)                                         gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg  >  1:476419‑M1/1‑21 (MQ=255)                                               |                                             GTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400262‑1400352

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