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

CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353                                             |                                        ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:1103586/1‑50 (MQ=255) ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:4161802/1‑50 (MQ=255) ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:5350720/1‑50 (MQ=255) ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:5513571/1‑50 (MQ=255)  cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:2995146/1‑49 (MQ=38)   ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:1643251/50‑1 (MQ=255)   ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:1874819/1‑50 (MQ=255)    gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:1523746/50‑1 (MQ=255)    gACCTATCGGCAGACCGACAAGATACCACGCAAGC‑ACATC‑GCAGCCATac                               <  1:1408163/50‑1 (MQ=255)     aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:506314/1‑49 (MQ=38)      ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             >  1:6657621‑M1/1‑48 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:815439‑M1/50‑4 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:187600‑M1/50‑4 (MQ=255)       cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:6740497‑M1/50‑4 (MQ=255)         aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:4536243‑M1/1‑45 (MQ=255)          tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:5082486‑M1/49‑6 (MQ=255)          tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:959589‑M1/1‑44 (MQ=255)           cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:5021718‑M1/1‑43 (MQ=255)           cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3542233‑M1/1‑43 (MQ=255)            ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:1501971‑M1/49‑8 (MQ=255)            ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:2245672‑M1/49‑8 (MQ=255)             gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:2804003‑M1/50‑10 (MQ=255)             gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:5794418‑M1/1‑41 (MQ=255)              cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:3425092‑M1/50‑11 (MQ=255)              cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:6580982‑M1/1‑40 (MQ=255)               aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:1289467‑M1/50‑12 (MQ=255)               aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1976125‑M1/1‑39 (MQ=255)                gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:250861‑M1/50‑13 (MQ=255)                gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:1997810‑M1/50‑13 (MQ=255)                gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:3011230‑M1/1‑38 (MQ=255)                 aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:754164‑M1/1‑37 (MQ=255)                 aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  <  1:1599663‑M1/50‑14 (MQ=255)                     aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:6117918‑M1/50‑18 (MQ=255)                     aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:2817727‑M1/50‑18 (MQ=255)                      cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:6695416‑M1/1‑32 (MQ=255)                      cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:5403562‑M1/1‑32 (MQ=255)                      cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:5365749‑M1/1‑32 (MQ=255)                       aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:702784‑M1/50‑20 (MQ=255)                       aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:2929885‑M1/1‑31 (MQ=255)                        aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:814068‑M1/50‑21 (MQ=255)                        aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:1719979‑M1/50‑21 (MQ=255)                         gATTCCACGCCAG‑AACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:2817821‑M1/50‑22 (MQ=255)                         gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:2715905‑M1/50‑22 (MQ=255)                         gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:4112865‑M1/1‑29 (MQ=255)                          aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         <  1:5159952‑M1/50‑23 (MQ=255)                           tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:986857‑M1/50‑24 (MQ=255)                           tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4297058‑M1/50‑24 (MQ=255)                            aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:671601‑M1/1‑26 (MQ=255)                             ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:5743848‑M1/50‑26 (MQ=255)                             ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:5384695‑M1/50‑26 (MQ=255)                             ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:4379097‑M1/50‑26 (MQ=255)                              cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:3616267‑M1/50‑27 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:6312843‑M1/50‑28 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:4867835‑M1/50‑28 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:3924433‑M1/50‑28 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:2337185‑M1/50‑28 (MQ=255)                               aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgaaggct    >  1:4588114‑M1/1‑23 (MQ=255)                                 gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:1013109‑M1/1‑21 (MQ=255)                                             |                                        CCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400270‑1400353

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