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

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353                                              |                                        cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:1278160/50‑1 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:4689124/1‑50 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:5302242/1‑50 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:4172801/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:718215/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:6590412/1‑50 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:5913133/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:5263406/50‑1 (MQ=255)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:3070160/1‑49 (MQ=38)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:1208081‑M1/50‑3 (MQ=255)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:1644247‑M1/50‑3 (MQ=255)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:2008427‑M1/50‑3 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:8231948‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:357690‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:2817904‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:1822935‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:7959563‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:7743560‑M1/48‑5 (MQ=255)            cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:4743986‑M1/1‑43 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:1675019‑M1/50‑9 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:1671943‑M1/50‑9 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:5122016‑M1/49‑8 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        >  1:507365‑M1/1‑42 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:4319073‑M1/49‑8 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:5043294‑M1/50‑10 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:3035390‑M1/50‑10 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:5751623‑M1/1‑40 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:8296040‑M1/1‑40 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1147662‑M1/1‑40 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2977602‑M1/1‑40 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3422366‑M1/1‑39 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:6782297‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:3742654‑M1/1‑38 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:6922686‑M1/1‑38 (MQ=255)                  aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:6592262‑M1/1‑37 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:1668550‑M1/50‑15 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:5379798‑M1/50‑15 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:1944826‑M1/50‑16 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:7582236‑M1/50‑16 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:7195941‑M1/50‑16 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:8016030‑M1/1‑34 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:730096‑M1/1‑34 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:1051119‑M1/50‑17 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:3970978‑M1/1‑34 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:2126170‑M1/50‑17 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:6397246‑M1/1‑34 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:5862817‑M1/1‑34 (MQ=255)                      aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:6480353‑M1/50‑18 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:2246164‑M1/50‑19 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:2743283‑M1/50‑19 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:7136819‑M1/50‑19 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:2502899‑M1/50‑19 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:4217778‑M1/50‑19 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:5463709‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:7791312‑M1/50‑20 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:5147774‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:3138568‑M1/50‑20 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:4025031‑M1/50‑21 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:7899386‑M1/1‑30 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:2202280‑M1/50‑21 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:302896‑M1/50‑22 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:7158014‑M1/50‑23 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:8224311‑M1/49‑25 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:5903850‑M1/1‑24 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:2568582‑M1/50‑27 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:5727216‑M1/1‑24 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:6168288‑M1/1‑24 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCTATACtatgcagaaacaatgctgtcgatggctg   <  1:3712947‑M1/50‑29 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:1970073‑M1/1‑22 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:3202104‑M1/50‑29 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:8106344‑M1/1‑22 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:5603231‑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: ‑