Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,895,043	T→C	L80P (CTA→CCA)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,895,043	0	T	C	100.0%	118.0 / NA	33	L80P (CTA→CCA)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/16);  total (17/16)

CGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTG  >  NC_000913/1894912‑1895177                                                                                                                                    |                                                                                                                                       cGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGc                                                                                                                                >  1:172499/1‑140 (MQ=255)       tGCCGATCATCCGTATAGCCGTTTTGATTTTGTGGTCGCCGAGCCGTTTTGCACTTTAACCACTTTCGGTAAAGAACCCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          <  1:177502/140‑1 (MQ=255)       tGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          <  2:20193/140‑1 (MQ=255)       tGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          <  1:370864/140‑1 (MQ=255)       tGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          >  1:60517/1‑140 (MQ=255)       tGCCGATCACCCGTATAGCCGTTTTGATTTTGTGGTCGCCGAGCCGATTTGCACTTTACCCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          <  2:76696/140‑1 (MQ=255)       tGCCGATCACCCGTATAGCCGTTTTGATATTGTGGCCGCGGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAACCGCACAACGACCACTGATGACCCGCCACAGGTGCTCcagc                                                                                                                          <  2:60517/140‑1 (MQ=255)            atcatcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTg                                                                                                                     >  1:215559/1‑140 (MQ=255)                       aGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCgcgc                                                                                                          <  2:172499/140‑1 (MQ=255)                        cccGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCa                                                                                                         <  1:90805/139‑1 (MQ=255)                         ccGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAg                                                                                                        >  2:3693/1‑140 (MQ=255)                            cTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACa                                                                                                     >  1:197795/1‑140 (MQ=255)                             tttGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCa                                                                                                         <  1:81012/135‑1 (MQ=255)                             tttGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCa                                                                                                         >  2:81012/1‑135 (MQ=255)                                       tgGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACg                                                                                          >  1:372604/1‑140 (MQ=255)                                         gTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCa                                                                                        >  2:173788/1‑140 (MQ=255)                                            gccgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATaa                                                                                     >  2:161156/1‑140 (MQ=255)                                                 gccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAg                                                                                >  1:295212/1‑140 (MQ=255)                                                                    aCTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGggcgg                                                             <  1:161156/140‑1 (MQ=255)                                                                       ttCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGgcgc                                                          >  1:374255/1‑140 (MQ=255)                                                                                   aCCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGgttgtt                                              <  1:252987/140‑1 (MQ=255)                                                                                        tgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                         >  1:369617/1‑140 (MQ=255)                                                                                        tgttAGTAAAAGCAAAAAACGCACAACGCCCCCTGAGGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                         <  2:374255/140‑1 (MQ=255)                                                                                         gttAGTGAAGGCGAAAACCCCACAACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTa                                        <  2:256091/140‑1 (MQ=255)                                                                                                      aaaaaCCCACAACGCCCCCTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccgc                           <  2:215559/140‑1 (MQ=255)                                                                                                             cacaACGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                    >  2:81977/1‑140 (MQ=255)                                                                                                                aaCGACCACTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                 >  1:260804/1‑140 (MQ=255)                                                                                                                       aCTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGcc          <  1:3693/140‑1 (MQ=255)                                                                                                                        cTGATGACCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCa         >  1:246333/1‑140 (MQ=255)                                                                                                                             gaCCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCa               <  1:388061/129‑1 (MQ=255)                                                                                                                             gaCCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCa               >  2:388061/1‑129 (MQ=255)                                                                                                                              aCCCGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAtt   >  2:2818/1‑140 (MQ=255)                                                                                                                               cccGCCACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTg  <  2:295212/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       CGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTG  >  NC_000913/1894912‑1895177

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