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%	99.5 / NA	30	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/18);  total (12/18)

CTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGC  >  NC_000913/1894636‑1895073                                                                                                                                                                                                             |                                                                                                                                                                                                                                           cTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCAcc                                                                                                                                                                                             >  2:97211/1‑251 (MQ=255)      cGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCg                                                                                                                                                                                        <  1:207962/251‑1 (MQ=255)                           tCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTcc                                                                                                                                                                   >  1:131993/1‑251 (MQ=255)                                  cTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATg                                                                                                                                                            >  1:108280/1‑251 (MQ=255)                                         cGCTCTACAGCTTTTTGATGCTGTTCAAGGGTGCGTCAAGGAACACCTGCAAACATTGTTCACTCCTGCTAAATTGTTGACGCTAATTATTTCATGCTACCCGCCACATAGCCAGTAGAGTAAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTCCCCTGGCGTCAGGAAGCCGCTGCATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatc                                                                                                                                                      >  2:326566/1‑250 (MQ=255)                                             cGACAGCTTTGTGCTGCTATTCATGGGTGAGTGAAGGTAAACCTGCAAACATTTTTAACCCCGGCTAAATTGTTGGCGCTAATTATTTCAGGCTTCCCGGCACATAGCCAGTAGAGTCAGGACGGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGCACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCg                                                                                                                                                 <  1:161692/251‑1 (MQ=255)                                                cAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGtat                                                                                                                                              >  1:186155/1‑251 (MQ=255)                                                            gctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGat                                                                                                                                  <  2:295820/251‑1 (MQ=255)                                                               gTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTTCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCGTCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGatat                                                                                                                                <  2:329507/250‑1 (MQ=255)                                                                    tGGGTGAGTGAAGGTAAACCTGCAACCATTGTTAACTCCTGCTAACTTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTAACCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGTTCATCCGTATAGCCGCTTTGATATTGTGGt                                                                                                                          >  2:27488/1‑251 (MQ=255)                                                                      ggTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTGCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTTCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTAGTTCTCCCCCTTAAGCCACCGGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCg                                                                                                                        <  1:326566/251‑1 (MQ=255)                                                                               aGGTAAACCTGCAAACTTTGTTAAATCCTGCTAAATTGTCGGCGCTAATTATTTCATGCAACCCCGCACATAGCCAGTAGAGTAAGGACTGATGAAGAAGTTATCTCCCGCTGTG‑ATTACTTTAACCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGAGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTACAGCCGCTTTGATATTGTGGTCgccgaaccg                                                                                                                >  2:37509/1‑250 (MQ=255)                                                                                               aTTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCAc                                                                                               <  1:97211/251‑1 (MQ=255)                                                                                                    tAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACTACTTTc                                                                                           <  1:93968/250‑1 (MQ=255)                                                                                                            gCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACCTAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGa                                                                                   >  2:227663/1‑250 (MQ=255)                                                                                                                  ttgttgGCGCTAATTATTTCATGCTACCAGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgtt                                                                            >  1:96927/1‑251 (MQ=255)                                                                                                                    ggtgGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACAACTTTCGGTAAAGAAACCgttgt                                                                          >  2:204188/3‑251 (MQ=255)                                                                                                                         cgcTAATTATTTCATGCTACCCGGCACATTGCCATGAGAGTCAGGACTGATGAAGACGTTATCTCCCGC‑GTGGATTACTTTACCCGGGCGTCAGGACGCCTTTGAATTTTAGTTCTCCCGCTTAAGCCACCTGCCGGGGGCGATTCTTTTACACTCCGGCTATGCCTATCATCCTTATAGCCGCTTTGAGATTGTGGTCTCCGAGCCGATTTTCACTTTAACCACTTTCGGTAAAGAATCCGTTTTTAGTg                                                                     <  2:225449/251‑1 (MQ=255)                                                                                                                             aaTTATTTCATGCTTCCCGGCCCATAACCAGTAGAGTCAGGTCTTAAGAAAACTGTATCTCCCGCCGTG‑GTTTCTTTCCCCTGGCGTCAGGACGCCCCTTAATTTTATTTCTCCCTCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAg                                                                 <  1:199129/251‑1 (MQ=255)                                                                                                                                 attTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGaa                                                             <  1:223917/251‑1 (MQ=255)                                                                                                                                                    aTAGCCAGTAGTGTAAGGACTGATGAAGACGTTATCTCCCGCTTTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTGCGGTAAAGAAACCGTTGTTAGTTAAAGCGAAAAACGCACAACGACCACtg                                          <  2:137628/251‑1 (MQ=255)                                                                                                                                                    aTAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCTATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACg                                                 <  2:161763/244‑1 (MQ=255)                                                                                                                                                    aTAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACg                                                 >  1:161763/1‑244 (MQ=255)                                                                                                                                                             agagTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCGGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCg                                  <  2:138940/250‑1 (MQ=255)                                                                                                                                                                  cAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAg                            <  1:64376/251‑1 (MQ=255)                                                                                                                                                                   aGGACTTATGAAGACGTTATCTCCCGCTTTG‑ATTCCTTTACCCTGGCGGCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCTTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGTCCACTGATGACCCGCTACAgg                           <  2:238311/251‑1 (MQ=255)                                                                                                                                                                   aGGACTGATGAAGACGTTATCTCCCGTTGTT‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAgg                           <  2:186155/251‑1 (MQ=255)                                                                                                                                                                            gaAGACGTTATCTCCCGCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACGCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCcag                   <  1:33934/250‑1 (MQ=255)                                                                                                                                                                               aaCGTTATCTCCCGCTGTT‑ATTTCTTTTCCCTGGCGTCCGGACGCCGCTGAATTTTATTTCCCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAgg               <  2:244789/250‑1 (MQ=255)                                                                                                                                                                                            gCTGTG‑ATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGAGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCgcgc  >  2:163470/1‑251 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                           CTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTG‑ATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGC  >  NC_000913/1894636‑1895073

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: ‑