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

TGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGA  >  NC_000913/1400259‑1400353                                                  |                                              tGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:2714846/1‑50 (MQ=255)  gggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCa                                               >  1:1770444/1‑49 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             >  1:1171315/1‑48 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCaca                                             <  1:199514/50‑3 (MQ=255)   ggTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCac                                              <  1:3762025/49‑2 (MQ=255)            ccGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCa                                  <  1:2317862/50‑1 (MQ=255)             cGGACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCAt                                 <  1:4259749/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                <  1:1499551/50‑1 (MQ=255)              ggACCTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATa                                >  1:122654/1‑50 (MQ=255)                acctATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACt                              <  1:1400769‑M1/49‑2 (MQ=255)                 ccTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACta                             <  1:3878502‑M1/50‑3 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:1291021‑M1/1‑47 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            >  1:3629143‑M1/1‑47 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:1929464‑M1/50‑4 (MQ=255)                  cTATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtat                            <  1:756593‑M1/50‑4 (MQ=255)                   tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:3697446‑M1/1‑46 (MQ=255)                   tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:668464‑M1/1‑46 (MQ=255)                   tATCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2945629‑M1/1‑46 (MQ=255)                    aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2766752‑M1/1‑45 (MQ=255)                    aTCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           >  1:2212138‑M1/1‑45 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgca                         <  1:1734408‑M1/50‑7 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgc                          <  1:2371351‑M1/49‑6 (MQ=255)                     tCGGCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatg                           <  1:2381632‑M1/48‑5 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       <  1:4069488‑M1/50‑9 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:3818051‑M1/49‑8 (MQ=255)                       ggCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcag                        <  1:2850324‑M1/49‑8 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:2518028‑M1/1‑41 (MQ=255)                        gCAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcaga                       >  1:1554219‑M1/1‑41 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     <  1:2187922‑M1/50‑11 (MQ=255)                         cAGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:3197037‑M1/1‑40 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    <  1:1254151‑M1/50‑12 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:537240‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:417297‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:1356658‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2563936‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2264869‑M1/1‑39 (MQ=255)                          aGACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaa                     >  1:2239664‑M1/1‑39 (MQ=255)                           gACCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaac                    >  1:2672664‑M1/1‑38 (MQ=255)                            aCCGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaa                  >  1:3830491‑M1/1‑37 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                <  1:955897‑M1/50‑16 (MQ=255)                              cGACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatg                >  1:3119649‑M1/1‑35 (MQ=255)                               gACAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgc               <  1:949491‑M1/50‑17 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:498678‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              >  1:2311598‑M1/1‑33 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:344285‑M1/50‑18 (MQ=255)                                aCAAGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgct              <  1:3807052‑M1/50‑18 (MQ=255)                                  aaGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgt            <  1:3714408‑M1/50‑20 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:403171‑M1/1‑30 (MQ=255)                                   aGATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:2098084‑M1/1‑30 (MQ=255)                                    gATACTACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:1858935‑M1/1‑29 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcg          <  1:3564561‑M1/50‑22 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:2024064‑M1/1‑29 (MQ=255)                                    gATACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtc           >  1:906903‑M1/1‑29 (MQ=255)                                      tACCACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgat        <  1:1693972‑M1/50‑24 (MQ=255)                                         cACGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggc     >  1:1022214‑M1/1‑24 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:2837979‑M1/50‑28 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:577872‑M1/50‑28 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    >  1:635164‑M1/1‑23 (MQ=255)                                          aCGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggct    <  1:777547‑M1/50‑28 (MQ=255)                                           cGCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctg   >  1:124404‑M1/1‑22 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:913465‑M1/1‑21 (MQ=255)                                            gCCAGC‑ACATC‑GCAGCCATACtatgcagaaacaatgctgtcgatggctgc  >  1:1563149‑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 < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑