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

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353                                              |                                        cccGGACCTAGCGGCAGACCGACAGGATACCACGCCAGC‑ACATC‑GCAGcc                                   <  1:4476640/50‑1 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCCCGCCAGC‑ACATC‑GCAGCCa                                  >  1:2575587/1‑50 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:4825841/1‑50 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:3518901/50‑1 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  >  1:6483883/1‑50 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:1678928/50‑1 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:5389889/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               <  1:4630208/50‑1 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:6041241/1‑50 (MQ=255)      acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:2491297‑M1/49‑2 (MQ=255)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                               >  1:788324/1‑49 (MQ=38)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             >  1:782755‑M1/1‑48 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:1493176‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:7276932‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:2391534‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:4706018‑M1/1‑47 (MQ=255)         tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:5961426‑M1/1‑46 (MQ=255)         tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2337505‑M1/1‑46 (MQ=255)         tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:1891103‑M1/1‑46 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3245784‑M1/1‑45 (MQ=255)           tCGGCAGACCGACAAGATACCGCGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:2928307‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:6998529‑M1/48‑5 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:5005951‑M1/48‑5 (MQ=255)            cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2672907‑M1/1‑43 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                      <  1:2181986‑M1/50‑10 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:4905985‑M1/1‑41 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:4430467‑M1/1‑40 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:6898853‑M1/1‑39 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3793953‑M1/1‑39 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:6278816‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:3240119‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                   <  1:1976760‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:3093849‑M1/1‑38 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:3354006‑M1/1‑36 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:1696832‑M1/1‑36 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                 >  1:6177840‑M1/1‑36 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:3297480‑M1/50‑16 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:1945909‑M1/1‑35 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:118094‑M1/50‑17 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               >  1:5454326‑M1/1‑34 (MQ=255)                      aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:6019971‑M1/1‑33 (MQ=255)                      aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:7022148‑M1/1‑33 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg             >  1:4928197‑M1/1‑32 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:4626578‑M1/50‑20 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            >  1:927936‑M1/1‑31 (MQ=255)                        aaGATACCACGCAAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:3402954‑M1/50‑20 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:906434‑M1/50‑21 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           <  1:6574027‑M1/50‑21 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:5729617‑M1/50‑22 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:2876722‑M1/50‑22 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:5768826‑M1/1‑29 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:942843‑M1/1‑29 (MQ=255)                          gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:4722456‑M1/1‑29 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:4699057‑M1/1‑28 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga         >  1:4034702‑M1/1‑28 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:4839426‑M1/50‑24 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        >  1:23669‑M1/1‑27 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:3710099‑M1/50‑24 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        >  1:6775154‑M1/1‑27 (MQ=255)                             aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:568087‑M1/1‑26 (MQ=255)                             aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:5806526‑M1/1‑26 (MQ=255)                             aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       >  1:6024381‑M1/1‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:10421‑M1/50‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:1934771‑M1/1‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:1960830‑M1/50‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      <  1:2703610‑M1/50‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg      >  1:3577261‑M1/1‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:3232945‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:6386411‑M1/49‑25 (MQ=255)                              acaCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg       <  1:3015077‑M1/48‑25 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:1152204‑M1/1‑24 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:6153460‑M1/50‑27 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     <  1:2530260‑M1/50‑27 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:5132885‑M1/1‑23 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:1157696‑M1/1‑23 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:5921518‑M1/50‑28 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:2118915‑M1/1‑23 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:4135875‑M1/50‑28 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:6869924‑M1/50‑28 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:6501793‑M1/50‑28 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:5955505‑M1/50‑29 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   <  1:5851169‑M1/50‑29 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:6740831‑M1/50‑30 (MQ=255)                                  gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  <  1:3206336‑M1/50‑30 (MQ=255)                                              |                                        CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400269‑1400353

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