Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,896,063	T→C	F420S (TTC→TCC)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,896,063	0	T	C	100.0%	103.1 / NA	31	F420S (TTC→TCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/13);  total (18/13)

TATCGACGAACTGGAACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTAAGACGTGGA  >  NC_000913/1895938‑1896174                                                                                                                               |                                                                                                                tATCGACGAACTGGAACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTcctgctct                                                                                                          >  2:482363/1‑133 (MQ=255)    cgacgaACTGGAACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACAGCCATTATCGGTC‑AGATTTCCGGCTCTgcg                                                                                                       >  1:140670/1‑133 (MQ=255)          aCTGGAACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGt                                                                                                 >  2:25580/1‑133 (MQ=255)                    cAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCg                                                                                       >  2:291303/1‑133 (MQ=255)                    cAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGATCTGCGGGCGGTGAAATTGTCg                                                                                       >  1:229130/1‑133 (MQ=255)                       cgtcGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCg                                                                                    <  2:56255/130‑1 (MQ=255)                         aCGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGAt                                                                                  <  1:291303/133‑1 (MQ=255)                                    ggTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGgaaga                                                                       <  1:469761/133‑1 (MQ=255)                                                     tATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGaaa                                                      >  2:407840/1‑133 (MQ=255)                                                             cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                                              <  1:361266/133‑1 (MQ=255)                                                             cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                                              <  1:407840/133‑1 (MQ=255)                                                             cTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGAt                                              <  1:482363/133‑1 (MQ=255)                                                                 tGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAg                                          >  1:305493/1‑133 (MQ=255)                                                                  gCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGt                                         <  2:305493/133‑1 (MQ=255)                                                                     gCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGAGGAATATAAGGAAACTTTTGAAAAAGTTaa                                      >  1:155447/1‑133 (MQ=255)                                                                       aaCATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGTGGAATATAAGGAAACTTTTGAGAAAGTTAATc                                    >  2:214388/1‑133 (MQ=255)                                                                                aCCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGa                           >  1:20999/1‑133 (MQ=255)                                                                                     tATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCaa                      >  1:485549/1‑133 (MQ=255)                                                                                             tCCGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGaa              >  2:116025/1‑133 (MQ=255)                                                                                              ccGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAATGTTAATCCTATACTGAAGCAACTGGAGAAg             >  1:206254/1‑133 (MQ=255)                                                                                               cGCACGCTGACTGCCATTAACGGACAAAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAACGTTAATCATATCCTGAAGCAAATGGAGAAg             >  1:176660/1‑133 (MQ=255)                                                                                               cGCACGCTGACTGCCATTAACGGACAAAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAg             >  1:128087/1‑133 (MQ=255)                                                                                               cGCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGt            >  2:88335/1‑133 (MQ=255)                                                                                                gCACGCTGACTGCCATTAACGGACGAAACTTCCTGCTCTGCCGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGt            <  2:176660/133‑1 (MQ=255)                                                                                                gCACGCTGACTGCCATTAACGGACAAAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGt            <  2:128087/133‑1 (MQ=255)                                                                                                gCACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTa           <  2:20999/133‑1 (MQ=255)                                                                                                 cACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTAAATCGTATCCTAAAACAACTGTAGAAGTaa          >  1:197975/1‑133 (MQ=255)                                                                                                 cACGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTATTGATAAAATTAGTCGTATCCTGAACCAACTCGAGAAGTaa          >  2:226593/1‑133 (MQ=255)                                                                                                  aCGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGga                 <  1:305164/125‑1 (MQ=255)                                                                                                  aCGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGga                 >  2:305164/1‑125 (MQ=255)                                                                                                  aCGCTGACTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTAAg         <  1:507744/133‑1 (MQ=255)                                                                                                     gtGACTGCCAATAACCGAC‑AAATTTCGTGCTCTGCGGGCGGAGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTAAGACg      <  2:229130/132‑1 (MQ=255)                                                                                                         cTGCCATTAACGGAC‑AAATTTCCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTAAGACGTGGa  <  1:25580/133‑1 (MQ=255)                                                                                                                               |                                                                                                                TATCGACGAACTGGAACCGCAGCGACGCAATGCCTGGTGCGGCAGCATTGGCTATTTGAGCTTTTGCGGCAACATGGATACCAGTATTACTATCCGCACGCTGACTGCCATTAACGGAC‑AAATTTTCTGCTCTGCGGGCGGTGGAATTGTCGCCGATAGCCAGGAAGAAGCGGAATATCAGGAAACTTTTGATAAAGTTAATCGTATCCTGAAGCAACTGGAGAAGTAAGACGTGGA  >  NC_000913/1895938‑1896174

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