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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353                                                  |                                              tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:4903598/1‑50 (MQ=255) tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:4625866/1‑50 (MQ=255) tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:4133202/1‑50 (MQ=255)  gggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:342227/1‑49 (MQ=255)   ggTAAATACCCGGTCCTATCGGCAGACCGACAAGATACCACGCCAGCac                                              <  1:3806415/49‑2 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             >  1:1217593/1‑48 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             >  1:1765436/1‑48 (MQ=255)           cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:1623417/50‑1 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:3639758/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:4691359/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:1639589/1‑49 (MQ=38)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:5609826/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:5351463/50‑1 (MQ=255)               gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:3173060/1‑50 (MQ=255)                acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:4068422‑M1/49‑2 (MQ=255)                aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:5273063/1‑49 (MQ=38)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             >  1:4701009‑M1/1‑48 (MQ=255)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:4067862‑M1/50‑3 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:5796359‑M1/50‑4 (MQ=255)                   tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3741209‑M1/1‑46 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:4927422‑M1/49‑6 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        >  1:2655237‑M1/1‑42 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:4539486‑M1/49‑8 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:3543694‑M1/1‑41 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:4835815‑M1/1‑41 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3867876‑M1/1‑40 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3314839‑M1/1‑40 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:3332341‑M1/50‑11 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:1719052‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:5246533‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:1645476‑M1/1‑38 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:207584‑M1/1‑38 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:5115465‑M1/1‑38 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  <  1:487088‑M1/50‑14 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:5139728‑M1/1‑37 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:2805993‑M1/1‑37 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:3047129‑M1/1‑37 (MQ=255)                             ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:72647‑M1/50‑15 (MQ=255)                             ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 <  1:1079895‑M1/50‑15 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:4155577‑M1/1‑35 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:4511570‑M1/1‑35 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:6427293‑M1/50‑16 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:5243555‑M1/50‑17 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:4319607‑M1/50‑17 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:5159962‑M1/50‑17 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:4042668‑M1/1‑34 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:3313351‑M1/50‑17 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:3695981‑M1/50‑17 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:2513670‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:3206377‑M1/1‑33 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:4885528‑M1/50‑18 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagagacaatgctg             <  1:1713462‑M1/50‑19 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             <  1:4546393‑M1/50‑19 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:3120599‑M1/1‑32 (MQ=255)                                 cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:2992398‑M1/1‑32 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:677754‑M1/50‑20 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:5714242‑M1/50‑20 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:4866024‑M1/50‑22 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:686715‑M1/1‑28 (MQ=255)                                     aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:5669971‑M1/1‑28 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4770583‑M1/50‑24 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4543871‑M1/50‑24 (MQ=255)                                       aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:5011233‑M1/1‑26 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:4353825‑M1/1‑25 (MQ=255)                                        ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:1493001‑M1/49‑25 (MQ=255)                                         cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:282778‑M1/1‑24 (MQ=255)                                         cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:2090408‑M1/50‑27 (MQ=255)                                         cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:2130460‑M1/50‑27 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:2003232‑M1/1‑23 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:3942887‑M1/50‑29 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:4075487‑M1/1‑22 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtggatggctgc  <  1:6343727‑M1/50‑30 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:6288201‑M1/1‑21 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:6311597‑M1/50‑30 (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 < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑