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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353                                                  |                                              tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               <  1:1629223/50‑1 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             <  1:2093354/50‑3 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCac                                              <  1:2625468/49‑2 (MQ=255)           cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:420865/50‑1 (MQ=255)            ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:236385/1‑50 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:2144420/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:3083874/1‑50 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:1678222/1‑50 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:1957588/50‑1 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:4209130/1‑50 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:2703060/1‑50 (MQ=255)                acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:3804469‑M1/49‑2 (MQ=255)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:1605816‑M1/50‑3 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:2365767‑M1/50‑4 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:4328982‑M1/1‑47 (MQ=255)                   tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:413608‑M1/1‑46 (MQ=255)                    aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:4390052‑M1/1‑45 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:2952587‑M1/50‑7 (MQ=255)                      cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:4151454‑M1/1‑43 (MQ=255)                      cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3276459‑M1/1‑43 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:2127102‑M1/50‑10 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:2596176‑M1/50‑11 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:2101965‑M1/50‑11 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:2793754‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2351742‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:3644358‑M1/50‑13 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:2652776‑M1/1‑35 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:2114585‑M1/50‑16 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:1437963‑M1/50‑16 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:3274924‑M1/1‑34 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:3606890‑M1/1‑33 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:1739643‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:379532‑M1/1‑33 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:93668‑M1/50‑19 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:895469‑M1/1‑32 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:4335467‑M1/1‑32 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:3340617‑M1/50‑20 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:3657579‑M1/1‑31 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:2782569‑M1/50‑21 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:774813‑M1/50‑21 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:679156‑M1/50‑21 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:4172218‑M1/50‑22 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4174789‑M1/50‑24 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:1999653‑M1/50‑24 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        >  1:2752868‑M1/1‑27 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:3897882‑M1/50‑24 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:3793765‑M1/50‑26 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:820080‑M1/50‑26 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:2636007‑M1/49‑25 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:3683022‑M1/49‑25 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:1142365‑M1/50‑29 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:2427783‑M1/50‑30 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:4385313‑M1/50‑30 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:653350‑M1/1‑21 (MQ=255)                                                  |                                              TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑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: ‑