Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,894,969	C→T	C55C (TGC→TGT)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,894,969	0	C	T	100.0%	54.8 / NA	18	C55C (TGC→TGT)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/9);  total (9/9)

AACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGT  >  NC_000913/1894727‑1895211                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   aaCATTGTTAACTCCTGCTAAATGGTTGGCTCTAATTATTTCATGCTACCCGGCCCATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCAGCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAAc                                                                                                                                                                                                                                            <  1:96106/251‑1 (MQ=255)       gTTAACTCATGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACttt                                                                                                                                                                                                                                      >  1:54975/1‑251 (MQ=255)        ttAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACttt                                                                                                                                                                                                                                      <  1:101210/250‑1 (MQ=255)           aCTCTTGCTAAAGTTTTGGCGTTAATTATTTCATTCTACCCGTCACATGGCCAGTAGAGTCAGGACTGAGGAA‑TTGTTGACCCCCCGCAGTGATTCCTTTTCCCGGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCTCT‑AAGCCACCTTCAGGGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTCCTTTAACCACTTTCGGTa                                                                                                                                                                                                                                 <  2:224479/251‑1 (MQ=255)                              cgcTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTg                                                                                                                                                                                                               >  1:197200/1‑251 (MQ=255)                                     tattTCATGCTACCCGGCACATAGACAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTACAGAAACCGTTGTTAGTGAAAGCg                                                                                                                                                                                                         >  2:204725/1‑250 (MQ=255)                                                        caTAGCCAGTAGAGTCAGGACTGATTAAGACGTT‑ATCTCCCGCTGTGATTACTTTCCCCTGGCGTCAGGACGCCGCTGAATGTTATTTCTCCCGCTTAAGCCACCTGCCTTGGGCGATGCTGTTACACTCCGGCTATGCCGTTCATCCGTATAGCCGCTTTGATAGTGTGGTCGCCGATCCGATTGGTACTGTAGCCGCTTTCGGGAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACt                                                                                                                                                                                     <  1:13521/251‑1 (MQ=255)                                                         aTAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtg                                                                                                                                                                                    >  2:89506/1‑251 (MQ=255)                                                                              gatgaAGACGTT‑ATCTCCCGCGGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTcc                                                                                                                                                               <  1:228285/251‑1 (MQ=255)                                                                                  aaGACGTT‑ATCTCCCGCTGTGCTTCCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGAGAAAGCGACAAACGCACAACGACCACTGATGACCCGCTACAgtgtctccagc                                                                                                                                                            >  2:51153/1‑240 (MQ=255)                                                                                                                    tGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCTCTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGc                                                                                                                           >  1:115435/1‑250 (MQ=255)                                                                                                                                                         tAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgg                                                                                       <  1:126137/249‑1 (MQ=255)                                                                                                                                                                                   tgaaccaccggatatgacgtccatcTGTATAGCCGCATTGATATTTGGGGCGCCGATCCGATATGTACTTTATCCACTTCCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTTCAGCAGGTGCTGGAGCGCGCAGACAGTCGCCCAGCGCATAACGAAGAGTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTGTGCAACGATCTGGgccgcc                                                            <  2:2387/229‑1 (MQ=255)                                                                                                                                                                                               tATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                               >  1:9409/1‑251 (MQ=255)                                                                                                                                                                                                aTGCCGATCATCCGTATAGCCGATTTGATATTGTGGTCGCCGAGCCGACTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCGGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTg                                              >  1:168807/1‑251 (MQ=255)                                                                                                                                                                                                                                     cGCCGAGCCGATTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGata         <  1:89506/251‑1 (MQ=255)                                                                                                                                                                                                                                      gccgagccgaTTTGTACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGatat        >  1:114508/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccgaTTTGTACTTTAACCACTGTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGt  <  1:28840/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   AACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGT  >  NC_000913/1894727‑1895211

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