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

CCGGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGA  >  NC_000913/1894592‑1895076                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              ccGGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCAGCTGTTCATGGGTGAGGGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTg                                                                                                                                                                                                                                            <  1:443167/251‑1 (MQ=255)              cACGCCGCTCTTCTGCCACCAGCGCGATTGCCTAGCCGCTGCTCATCACTTGTGCGATCAGTTCCTGGATCCGCTCGACCGCTTTTTGCTGCTGTTCATGGGTGAGAGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTACATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTCTCTCCCGCTGTGATTACTTTACCCTGGCGTCAGgccgcc                                                                                                                                                                                                                               >  1:325657/1‑251 (MQ=255)                    nctctTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGa                                                                                                                                                                                                                          >  1:133368/2‑250 (MQ=255)                                   gcgcGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCCGTTCCTGGATCGGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCGCACATTGTTAACTCCTGCTAATTTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTCGAGTCAGGACTGATGCAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGACTTTTATTTCTCCCGc                                                                                                                                                                                                          >  1:204316/1‑251 (MQ=255)                                   gcgcGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACAGAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCTCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGc                                                                                                                                                                                                          <  1:218107/251‑1 (MQ=255)                                        aTTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAAGTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAg                                                                                                                                                                                                     <  1:373595/251‑1 (MQ=255)                                            ccTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCAc                                                                                                                                                                                                 >  1:398453/1‑251 (MQ=255)                                                                                              gCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCAAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGACAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGtat                                                                                                                                                 >  1:57524/1‑249 (MQ=255)                                                                                                      ntgctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTCCCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatc                                                                                                                                                         >  1:53718/2‑233 (MQ=255)                                                                                                        gctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTg                                                                                                                                       <  1:366317/249‑1 (MQ=255)                                                                                                                  ngtgagtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCAGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCACGCTTAAGCCACCTGCCGTGGGCGCTGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCg                                                                                                                           >  1:38252/2‑251 (MQ=255)                                                                                                                       gtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTACATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGCTGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccg                                                                                                                        >  1:393097/1‑249 (MQ=255)                                                                                                                                     gCAAACATTGTTAACTCCTGCTAAAGTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGTCTGATGAAGACTTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCGCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATTCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACtt                                                                                                         <  1:225430/250‑1 (MQ=255)                                                                                                                                                              gctggcgaaatAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCCGGACTGGTGAAGACGTAATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGCCGCCGCTGAATTTTATTTATCCCGCTTAAACAACCTGCCGTGGGCGCTGCTTTTACACTCCGGCTATGCCGGTCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAAc                                                                                                     >  1:211903/11‑229 (MQ=255)                                                                                                                                                                                tCATGCTACCCGGCACATAGCCAGTTGAGTCAGGACTGATGCAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGACGTGGGCGATGATTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGc                                                                   >  1:391153/1‑245 (MQ=255)                                                                                                                                                                                   ngCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCAc                                                                                                  >  1:266431/2‑211 (MQ=255)                                                                                                                                                                                               caTAGCCAGTAGAGTCAGGACTGATGATGACGTTATCTCACGATGTGATTACTTTACCCTGGCGAAAGGACGCCGCTGAATTTTATTTATCCCGCTTAAGCAACATGCCGTGGACGATGCTTTTACACTCCGGCTATGCCGATCCTCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGc                                                                   >  1:241531/1‑230 (MQ=255)                                                                                                                                                                                                     cAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTAACCTGGAATCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACATGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCg                                                                                                                                                            >  1:64108/1‑135 (MQ=255)                                                                                                                                                                                                       gTAGAGTCAGGACTGATGAAGACGTTNTCTCCCNNTGNGATTACTTTACCCTGGCGTNAGNNNGNNNCNGAATTNTANNTCTCCCGCTNAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGc                                                                                                                                                                   >  2:97808/1‑126 (MQ=255)                                                                                                                                                                                                       gTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGc                                                                                                                                                                   <  1:97808/126‑1 (MQ=255)                                                                                                                                                                                                              cAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCTTGGGCGATGCTGTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTg                                                                                                              <  1:392736/172‑1 (MQ=255)                                                                                                                                                                                                               aGGACGGAGGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAAGTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATAGTGTGGTCGCCGAGCCGATTTGCACTTTAACCACGTTCgg                                                                                            <  1:445126/189‑1 (MQ=255)                                                                                                                                                                                                                gtgcTGATGAAGTCGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTan                                                                                         <  1:268511/188‑2 (MQ=255)                                                                                                                                                                                                                ggACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTNTAGNNNNnttnn                                                                                                                                      >  2:268511/1‑144 (MQ=255)                                                                                                                                                                                                                   cTGATGAAGACGTTATCTCCCGCTGTGATTTCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgttn                                                                              <  1:270521/199‑2 (MQ=255)                                                                                                                                                                                                                   cTGATGAAGACGTTATCTCACGCTGTGATTTCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATANCCGNNNNNatnnnnnnnnnn                                                                                                                           >  2:270521/1‑144 (MQ=255)                                                                                                                                                                                                                       tgaAGACGTTATCTCCCGCTGTGATTACTGTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATgn                                                                                                                                                              <  1:47557/115‑2 (MQ=255)                                                                                                                                                                                                                           gACGTTATCGCCCGCTGTGATTTCTTTCCCCTGGCGTCACGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTATTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCcagcagn                  <  1:48899/251‑2 (MQ=255)                                                                                                                                                                                                                                tATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCTCTTAAGCCACCTGCCGTTGGCGATGCTTTTACACTCCGGCTATGCCTATCATCCGTAGAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAg                                                                                       <  1:426730/177‑1 (MQ=255)                                                                                                                                                                                                                                       cGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTAt                                                                                                                                                                                                                   >  1:29544/1‑46 (MQ=255)                                                                                                                                                                                                                                       cGCTGTGATTACTTTACCCGGGCGTCAGGACGCCGGTGAATTTTAt                                                                                                                                                                                                                   <  2:29544/46‑1 (MQ=255)                                                                                                                                                                                                                                           ntgATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCTACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCAGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAACGCGACAAACGCACCACGACCACTGATGCCCCGCTACAGGTGCTCCAGCAGGTgtcaccttgcagcat  >  1:54462/2‑237 (MQ=255)                                                                                                                                                                                                                                             gATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACTCTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgcn                                                                                                                         <  1:82049/130‑2 (MQ=255)                                                                                                                                                                                                                                                nnnnnnacnNNNNCGTNAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCAc                                                                                                  <  2:266431/144‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              CCGGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGA  >  NC_000913/1894592‑1895076

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