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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353                                                  |                                              tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:2372760/1‑50 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             <  1:1041666/50‑3 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             >  1:2678850/1‑48 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             <  1:4930142/50‑3 (MQ=255)            ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:1443684/50‑1 (MQ=255)            ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:6000012/1‑50 (MQ=255)            ccGGACCTATCGACAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:2498030/50‑1 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:3739459/50‑1 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:3393563/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:115443/1‑50 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:7059922/1‑50 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:4354076/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:2650179/1‑50 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:2227428/50‑1 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:4471581/50‑1 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:2909695/1‑50 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:2072724/1‑50 (MQ=255)                acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:1475640‑M1/49‑2 (MQ=255)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             >  1:6594507‑M1/1‑48 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:3248118‑M1/1‑47 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:5312403‑M1/50‑4 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:2742340‑M1/50‑7 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:6938186‑M1/50‑7 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:1803513‑M1/1‑44 (MQ=255)                      cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:5438842‑M1/1‑43 (MQ=255)                      cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:5654519‑M1/1‑43 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:2093472‑M1/50‑9 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:3490975‑M1/50‑9 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:507988‑M1/49‑8 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:6548747‑M1/50‑10 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:1421371‑M1/1‑41 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3990263‑M1/1‑40 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:586186‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:3583096‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:4132759‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:795703‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:4817719‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:6590300‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:6947452‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:4318063‑M1/1‑38 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:2749921‑M1/1‑37 (MQ=255)                             ccGACAAGATACCTCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:2161117‑M1/50‑15 (MQ=255)                             ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:5310739‑M1/1‑36 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:4758414‑M1/50‑16 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:1057455‑M1/50‑17 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:6356289‑M1/1‑34 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcggaaacaatgct              <  1:6073196‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:490171‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:2359251‑M1/50‑18 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:4315358‑M1/50‑19 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:4901602‑M1/1‑32 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:1367514‑M1/1‑31 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:4477953‑M1/1‑31 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:2721267‑M1/50‑21 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:1949612‑M1/1‑30 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:5247248‑M1/1‑30 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:4281314‑M1/1‑29 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:228403‑M1/1‑29 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:1608111‑M1/50‑23 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:1866248‑M1/50‑23 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:648666‑M1/50‑24 (MQ=255)                                       aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:494033‑M1/1‑26 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:3907564‑M1/50‑26 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:6191521‑M1/49‑25 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:2217076‑M1/50‑28 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:4633176‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:1230099‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:1281673‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:5561306‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:2734868‑M1/50‑29 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:3367385‑M1/1‑21 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:6814503‑M1/1‑21 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:3284841‑M1/50‑30 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:4100573‑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 < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑