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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400259‑1400352                                                  |                                             tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                              <  1:948907/50‑1 (MQ=255)  gggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                              >  1:1257394/1‑49 (MQ=255)  gggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                              >  1:1984067/1‑49 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                 >  1:1311551/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                <  1:197622/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                               >  1:282043/1‑50 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                               >  1:2432396/1‑50 (MQ=255)                aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATact                             <  1:1769435/50‑2 (MQ=255)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                            >  1:634379‑M1/1‑48 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           <  1:2228341‑M1/50‑4 (MQ=255)                    aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:1199828‑M1/1‑45 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                         <  1:1413419‑M1/49‑6 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          <  1:2435608‑M1/48‑5 (MQ=255)                      cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:214524‑M1/1‑43 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                     <  1:1005334‑M1/50‑10 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                     <  1:596783‑M1/50‑10 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                   <  1:2023807‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                    >  1:233804‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                  <  1:614182‑M1/50‑13 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc          >  1:360737‑M1/1‑30 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc          >  1:1297395‑M1/1‑29 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga        >  1:1478419‑M1/1‑28 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga        >  1:1753389‑M1/1‑28 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat       <  1:771571‑M1/50‑24 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      <  1:2028854‑M1/49‑25 (MQ=255)                                         cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc    >  1:519637‑M1/1‑24 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct   >  1:2629168‑M1/1‑23 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg  <  1:851602‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg  <  1:1755275‑M1/50‑29 (MQ=255)                                                  |                                             TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400259‑1400352

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: ‑