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%	100.2 / NA	22	coding (1460/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)

GGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400272‑1400352                                           |                                       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)                                           |                                       GGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400272‑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: ‑