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

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353                                              |                                        cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:4295242/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:4251463/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:3383619/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:5634929/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:1236497/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:424196/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:5170206/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:4430207/50‑1 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:1580295/50‑1 (MQ=255)      acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:3573507‑M1/49‑2 (MQ=255)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:179478/1‑49 (MQ=38)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:5590591/1‑49 (MQ=38)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:1290289‑M1/50‑3 (MQ=255)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             >  1:3669376‑M1/1‑48 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:4690193‑M1/1‑45 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3808139‑M1/1‑45 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:1590344‑M1/1‑45 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:607402‑M1/1‑45 (MQ=255)            cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:1643894‑M1/1‑43 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        >  1:306183‑M1/1‑42 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:5085353‑M1/50‑10 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:1573621‑M1/50‑10 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:2438525‑M1/50‑12 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1646487‑M1/1‑39 (MQ=255)                aGACCGACAAGAAACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:4781932‑M1/50‑12 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:1841901‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:4055439‑M1/1‑38 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:4781894‑M1/1‑38 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:1139564‑M1/1‑38 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:1776933‑M1/50‑15 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:737829‑M1/50‑15 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:71306‑M1/1‑35 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:905881‑M1/1‑35 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:811993‑M1/50‑16 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:279162‑M1/1‑35 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:1757018‑M1/50‑16 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:1822402‑M1/1‑34 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:2469610‑M1/1‑32 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:1453509‑M1/50‑19 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:2022091‑M1/50‑20 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:5055879‑M1/1‑30 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:516230‑M1/50‑21 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:3201499‑M1/50‑22 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:3197120‑M1/1‑28 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:5260989‑M1/1‑28 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:1606093‑M1/50‑24 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:5785273‑M1/1‑23 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:5883874‑M1/1‑23 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:4213024‑M1/50‑30 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:367340‑M1/1‑21 (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 < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑