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

CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400269‑1400352                                              |                                       cccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGcc                                  <  1:2551700/50‑1 (MQ=255)  ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                 >  1:3219366/1‑50 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                 >  1:3557427/1‑49 (MQ=38)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                <  1:2031803/50‑1 (MQ=255)   cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                <  1:1721210/50‑1 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                               >  1:7575975/1‑50 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                               >  1:7277319/1‑50 (MQ=255)    ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                               >  1:7211781/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              >  1:1813702/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              >  1:2374701/1‑50 (MQ=255)     gACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              <  1:595714/50‑1 (MQ=255)      acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                             <  1:2976521‑M1/49‑2 (MQ=255)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              >  1:4873264/1‑49 (MQ=38)      aCCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATac                              >  1:1329482/1‑49 (MQ=38)       ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                            <  1:6029910‑M1/50‑3 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           <  1:3253955‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           <  1:4306849‑M1/50‑4 (MQ=255)        cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                           >  1:1484682‑M1/1‑47 (MQ=255)          aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:3302125‑M1/1‑45 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                        <  1:6946190‑M1/50‑7 (MQ=255)           tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:5360711‑M1/1‑44 (MQ=255)            cGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                          >  1:1354926‑M1/1‑43 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                       >  1:1197209‑M1/1‑42 (MQ=255)             ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                       >  1:4169316‑M1/1‑42 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                     <  1:6215131‑M1/50‑10 (MQ=255)              gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaa                     <  1:1619167‑M1/50‑10 (MQ=255)               cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                    <  1:4403701‑M1/50‑11 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                   <  1:6093026‑M1/50‑12 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                    >  1:2954868‑M1/1‑39 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                    >  1:1301639‑M1/1‑39 (MQ=255)                aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                    >  1:823355‑M1/1‑39 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaaca                  <  1:388236‑M1/50‑13 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                   >  1:6628537‑M1/1‑38 (MQ=255)                 gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                   >  1:4915892‑M1/1‑38 (MQ=255)                  aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                 >  1:1998824‑M1/1‑37 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                <  1:3515176‑M1/50‑15 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                >  1:5949502‑M1/1‑36 (MQ=255)                   ccGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaat                >  1:4609278‑M1/1‑36 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg               >  1:1030455‑M1/1‑35 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg               <  1:4907451‑M1/50‑16 (MQ=255)                    cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg               >  1:5078302‑M1/1‑35 (MQ=255)                     gACAAGATGCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc              <  1:240453‑M1/50‑17 (MQ=255)                     gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc              >  1:7116924‑M1/1‑34 (MQ=255)                      aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct             >  1:1702945‑M1/1‑33 (MQ=255)                      aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct             <  1:1917952‑M1/50‑18 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg            <  1:3815992‑M1/50‑19 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg            >  1:7087484‑M1/1‑32 (MQ=255)                       cAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctg            >  1:7404365‑M1/1‑32 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           >  1:5579987‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           <  1:4369375‑M1/50‑20 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           >  1:2301769‑M1/1‑31 (MQ=255)                        aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt           >  1:4055755‑M1/1‑31 (MQ=255)                         aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc          >  1:6757277‑M1/1‑30 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga        >  1:1659662‑M1/1‑28 (MQ=255)                           aTACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcga        >  1:7376512‑M1/1‑28 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat       <  1:5528765‑M1/50‑24 (MQ=255)                            tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat       <  1:7031265‑M1/50‑24 (MQ=255)                             aCCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      >  1:147504‑M1/1‑26 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatgg     >  1:508120‑M1/1‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      <  1:1089792‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      <  1:2841582‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      <  1:427438‑M1/49‑25 (MQ=255)                              ccACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatg      <  1:3832599‑M1/49‑25 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc    >  1:2080933‑M1/1‑24 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc    <  1:5642299‑M1/50‑27 (MQ=255)                               cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc    <  1:2560938‑M1/50‑27 (MQ=255)                                aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct   <  1:1283829‑M1/50‑28 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg  <  1:5156555‑M1/50‑29 (MQ=255)                                 cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg  >  1:3552762‑M1/1‑22 (MQ=255)                                              |                                       CCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAG  >  NC_000913/1400269‑1400352

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