Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,894,839	T→C	L12P (CTC→CCC)	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,839	0	T	C	100.0%	49.7 / NA	16	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/7);  total (9/7)

TTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGC  >  NC_000913/1894632‑1895057                                                                                                                                                                                                                 |                                                                                                                                                                                                                           ttGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGc                                                                                                                                                                                 >  2:65132/1‑251 (MQ=255)                             cATCAGTTCATGGATCCGCTCGACAGCTATTTGCTGCTGTTCATGGGTGGGTGAAGGAAAACCTGCAAACATTGTTCACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGACTTTTATTTCTACCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACt                                                                                                                                                     >  2:251697/1‑251 (MQ=255)                                    tCCTGGCTCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGGAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGGGATTACTTTACCCTGGCGTCAGGAcg                                                                                                                                                                                                              >  1:72140/1‑187 (MQ=255)                                    tCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATGATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGAcg                                                                                                                                                                                                              <  2:72140/187‑1 (MQ=255)                                                  gACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCGGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGt                                                                                                                                <  1:251697/251‑1 (MQ=255)                                                                                          ccTGCAAACATTGTTAACTCATGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCCAGACGATTTgtac                                                                                        >  1:51171/1‑248 (MQ=255)                                                                                                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)                                                                                                                                                                                                                 |                                                                                                                                                                                                                           TTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTT‑ATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGC  >  NC_000913/1894632‑1895057

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