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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353                                                  |                                              tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:3032327/1‑50 (MQ=255)            ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:787755/1‑50 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:860185/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:408670/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:4591926/1‑49 (MQ=38)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:1180370/1‑50 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:4114982/1‑50 (MQ=255)                acctAACGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:3836625‑M1/49‑2 (MQ=255)                aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:3517174/1‑49 (MQ=38)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:8642‑M1/50‑3 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:390937‑M1/1‑47 (MQ=255)                    aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:1724726‑M1/1‑45 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:580995‑M1/48‑5 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:1658810‑M1/48‑5 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:3105264‑M1/48‑5 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:3261584‑M1/48‑5 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:205039‑M1/49‑8 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:4436471‑M1/49‑8 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        >  1:2479739‑M1/1‑42 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:3322528‑M1/50‑11 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:4414767‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1328672‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2685119‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:4342591‑M1/1‑38 (MQ=255)                             ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:4587609‑M1/1‑36 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:2303280‑M1/1‑35 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:393784‑M1/50‑16 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:3452716‑M1/50‑17 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:1937230‑M1/50‑18 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:337721‑M1/1‑32 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:1441283‑M1/1‑31 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:558241‑M1/1‑31 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:232397‑M1/1‑31 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:4692292‑M1/1‑31 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:22263‑M1/50‑20 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:1395664‑M1/1‑30 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:3385458‑M1/1‑30 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:3990000‑M1/50‑21 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:1400290‑M1/50‑22 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:45579‑M1/50‑24 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:2906267‑M1/50‑24 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:3936132‑M1/49‑25 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:4357961‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:4221142‑M1/1‑22 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:1666258‑M1/50‑29 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:2314857‑M1/50‑30 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:3840682‑M1/50‑30 (MQ=255)                                                  |                                              TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353

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