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

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353                                              |                                        cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:5732600/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:6231918/1‑49 (MQ=38)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:3929342/1‑49 (MQ=38)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:1604581/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:6680530/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:3269027/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:6350263/50‑1 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:4510207/1‑50 (MQ=255)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:1984280/1‑49 (MQ=38)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:3343757/1‑49 (MQ=38)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:853413‑M1/50‑3 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtgt                            <  1:5788403‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:265441‑M1/1‑47 (MQ=255)         tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:267415‑M1/1‑46 (MQ=255)         tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:576360‑M1/1‑46 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:4869997‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:1305682‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:4510468‑M1/49‑6 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:6407075‑M1/48‑5 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:269934‑M1/1‑44 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:6523336‑M1/50‑11 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:4089220‑M1/1‑40 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:5474019‑M1/1‑39 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2622237‑M1/1‑39 (MQ=255)                  aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:3517751‑M1/1‑37 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:1798148‑M1/50‑15 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:6513939‑M1/1‑34 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:1534859‑M1/50‑17 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:662247‑M1/1‑32 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:2143671‑M1/1‑32 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:6699300‑M1/50‑20 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:5421765‑M1/50‑20 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:846219‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:2489624‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:5736016‑M1/50‑20 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:3159735‑M1/1‑30 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:5743698‑M1/1‑30 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:22478‑M1/50‑21 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4784193‑M1/50‑24 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        >  1:6636242‑M1/1‑27 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        >  1:4868051‑M1/1‑27 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:375055‑M1/50‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:2961969‑M1/50‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:1190658‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:5802704‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:3902026‑M1/49‑25 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:2931000‑M1/1‑24 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:1979716‑M1/1‑22 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:5689647‑M1/50‑29 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:326956‑M1/1‑22 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:3365594‑M1/50‑30 (MQ=255)                                              |                                        CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353

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