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

GGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400260‑1400353                                                 |                                              gggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:1955614/1‑49 (MQ=255)  ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             >  1:689048/1‑48 (MQ=255)          cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:4699242/50‑1 (MQ=255)           ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:5504137/50‑1 (MQ=255)           ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:3126937/50‑1 (MQ=255)            cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:3015850/1‑49 (MQ=38)            cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:1090090/50‑1 (MQ=255)            cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:562229/50‑1 (MQ=255)             ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:41326/1‑50 (MQ=255)             ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:3755480/1‑50 (MQ=255)             ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:988071/1‑50 (MQ=255)              gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:4454720/50‑1 (MQ=255)                   aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3748153‑M1/1‑45 (MQ=255)                    tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:55864‑M1/50‑7 (MQ=255)                    tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:1851836‑M1/49‑6 (MQ=255)                    tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:3341970‑M1/49‑6 (MQ=255)                    tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:852999‑M1/48‑5 (MQ=255)                      ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:3222809‑M1/50‑9 (MQ=255)                      ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:1268232‑M1/50‑9 (MQ=255)                      ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:426863‑M1/49‑8 (MQ=255)                      ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        >  1:1898464‑M1/1‑42 (MQ=255)                      ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:586394‑M1/49‑8 (MQ=255)                       gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:5357230‑M1/50‑10 (MQ=255)                        cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2943771‑M1/1‑40 (MQ=255)                        cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:2658761‑M1/50‑11 (MQ=255)                         aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:1291206‑M1/50‑12 (MQ=255)                         aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2391586‑M1/1‑39 (MQ=255)                         aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1906320‑M1/1‑39 (MQ=255)                          gACCGACAAGATACCACGCCAGC‑ATATC‑GCAGCCATACtatgcagaaaca                   <  1:3538842‑M1/50‑13 (MQ=255)                           aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:778963‑M1/1‑37 (MQ=255)                            ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:3205876‑M1/50‑15 (MQ=255)                            ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:181378‑M1/1‑36 (MQ=255)                             cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:3113719‑M1/50‑16 (MQ=255)                             cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:5101949‑M1/1‑35 (MQ=255)                             cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:3394928‑M1/50‑16 (MQ=255)                              gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:549443‑M1/50‑17 (MQ=255)                              gACAAGATACCACGCCAGC‑ACAAC‑GCAGCCATACtatgcagaaacaatgc               <  1:1116931‑M1/50‑17 (MQ=255)                                cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:3913075‑M1/50‑19 (MQ=255)                                 aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:1960409‑M1/50‑20 (MQ=255)                                   gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:4185901‑M1/50‑22 (MQ=255)                                   gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:374285‑M1/50‑22 (MQ=255)                                    aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:5367153‑M1/50‑23 (MQ=255)                                    aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:1258771‑M1/50‑23 (MQ=255)                                    aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:922931‑M1/50‑23 (MQ=255)                                     tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:3368008‑M1/50‑24 (MQ=255)                                       ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:4164860‑M1/1‑25 (MQ=255)                                       ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:1419617‑M1/1‑25 (MQ=255)                                       ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:646777‑M1/49‑25 (MQ=255)                                       ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:816171‑M1/49‑25 (MQ=255)                                        cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:2227478‑M1/1‑24 (MQ=255)                                        cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:3492430‑M1/50‑27 (MQ=255)                                         aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:3140771‑M1/50‑28 (MQ=255)                                           gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:1804440‑M1/50‑30 (MQ=255)                                           gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:1617438‑M1/1‑21 (MQ=255)                                                 |                                              GGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400260‑1400353

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