Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,892,551	C→T	T113I (ACC→ATC)	queD →	6‑carboxy‑5,6,7,8‑tetrahydropterin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,892,551	0	C	T	100.0%	36.9 / NA	15	T113I (ACC→ATC)	queD	6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/8);  total (7/8)

GTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCG  >  NC_000913/2892346‑2892795                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     gTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATTGATTTGGGATCAGGTTAAACCCGGTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTcggcgg                                                                                                                                                                                                         <  1:192199/251‑1 (MQ=255)                                 cGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAggg                                                                                                                                                                          <  1:7574/250‑1 (MQ=255)                                                   gTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCAAGGTCTGGAAAACCCAACCAGCGCGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCt                                                                                                                                                       >  1:241359/1‑251 (MQ=255)                                                                                                    ccAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGt                                                                                                      <  1:7938/251‑1 (MQ=255)                                                                                                        gTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGGTAAACCTGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAAt                                                                                                                                                                                                    <  1:126235/153‑1 (MQ=255)                                                                                                                  cccAACCAGCGAGGTTTTAGCCAAATGGATTTGGGATCCGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGCGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGGCAATTTCCCTTTAGTAACGCTATGTTAGGGGGTTGTGTTCTGGGTATCTGGGGCATGACATGGACGACGACTGCAACATTATTATTATTGGTGCCGGTATTGCCGGCACCGc                                                                                        >  1:103186/1‑251 (MQ=255)                                                                                                                        cAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTc                                                                                                                                                                                                 <  1:241092/140‑1 (MQ=255)                                                                                                                                     gCAAAATGGATTTGGGATCAGCTTAAATCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGACGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGGGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGGCATTATTATTATTGGTGCCGGTATTGCAGGCAACGCTTGCGCGTTACGCTgcgcg                                                                     >  1:180198/1‑251 (MQ=255)                                                                                                                                                              aaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTCAAAGACACCTGCATCGCAGCTTGTATCTATCGCGGCGAATGATAAGCGTGAGTCGGCGGACAAATTCCCTTTAGTAACGCTCTGTTCGGGTGTTGTGTTATGGCTCTCTGGGGAATGACATGGAACACGCCTGCGACATTATTATTCTTGGTGACGGTATTGCAGGAACCGCTTGAGCATTACCCTGCGCGCGAGCGGGTTTATCAGTTTTGTTAc                                            >  1:122205/1‑251 (MQ=255)                                                                                                                                                                     tgtgCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTTGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTTGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTAGTATTGGTGCCGGCATTGCATGCACCGCTTGCGCGTGACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTgn                                         <  1:47994/247‑2 (MQ=255)                                                                                                                                                                      gtgCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGAGAATGATAAGAGTGTGTCGGCGGCAAATTTCCCTTAAGTAACGCTATGTTAGGGGGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAAcgc                                    >  1:108659/1‑251 (MQ=255)                                                                                                                                                                                      cGGTGATGGTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAAcn                                                                                                                                                                                     <  1:50906/90‑2 (MQ=255)                                                                                                                                                                                             ggTAAAAGAAACCTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGGCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCAtn                                                                                                                                               <  1:24669/121‑2 (MQ=255)                                                                                                                                                                                                        ncTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTACGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGCCATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTgg                                         >  1:25789/2‑212 (MQ=255)                                                                                                                                                                                                        ccTGCATCGCAGGTTGTATCTATCGCGGCGAATGATAAGCGTGTGTCGGCGGTCAATTTCCCTTCAGTAACGATATGCTAGGGTGTTGTGTTCTGGATATATGGGGCATGACATGGAAGACGACTGCGCCATTATTATTATTGGTGCCGGTATTGCAGGCACCGATTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGCACGCGCTGAAATCCCCGGCAGCAAATATTTTTCcggcg  >  1:171160/1‑251 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     GTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCG  >  NC_000913/2892346‑2892795

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