Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,257,606	C→T	Q254Q (CAG→CAA)	psuG ←	pseudouridine‑5'‑phosphate glycosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,257,606	0	C	T	100.0%	28.2 / NA	12	Q254Q (CAG→CAA)	psuG	pseudouridine‑5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/7);  total (5/7)

TGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGT‑TGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATCGGGTTCGCTACCACGAGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAA  >  NC_000913/2257392‑2257782                                                                                                                                                                                                                         |                                                                                                                                                                                 ngatgaGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGATCGCCAGAATGGCGG‑TGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGAACGCCGCATTGATAGTGTGTTCCGGCAt                                                                                                                                               >  1:149636/2‑250 (MQ=255)                       tCTGCCCGTAAAAAATTAACCCGCGAGACGCTTATATTCTTGGGCAATGTCGCTCGACAGAATGGAGT‑TGTTGAACACCAGCTGGATGTTGGATTGAAGGTTGG‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATCggn                                                                                                                         <  1:312193/250‑2 (MQ=255)                               tCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGCCTGGCGG‑AGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCACACTGTTCCGGAATCGGGTTCGCTAcc                                                                                                                >  1:290940/1‑251 (MQ=255)                                                          attgttgggCAGTTTCGCTCGCCAGATTGG‑GTGTGTTAAACGCCAGCTGGATGTTGGAGTTC‑GGGTGTTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTAGTTTCTTTACCAATGACGCCCTGAGCGTCAGCTTCATCTACCGCTTGATCGATCTCCGGATTGATAGTGGGTTCCGTCATAGCAAGCTGGTCCGGGATCGTGTTCGCTACCACGAGGCCACCGTTCAGACCGCTTTGc                                                                                     <  1:333920/244‑1 (MQ=255)                                                                                         gt‑tgttGAACACCAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCa                                                                                                                                                                   >  1:274556/1‑142 (MQ=255)                                                                                         gagacgtgtccttctccgatctggtGGATTTCAGGCGGT‑CCCCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATCGGGTTCGCTACCACGAGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTc                                                        <  1:264168/226‑1 (MQ=255)                                                                                                     ccAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGCa                                                                                                                                                <  1:426958/150‑1 (MQ=255)                                                                                                                                                  gCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCTTgatc                                                                                                                                                                             <  1:154893/77‑1 (MQ=255)                                                                                                                                                                                 cgttACCAATGACTCCCGGAGCTTCAGCTTCAGCTACCGCTTGATCGATCGCCGCTTTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATc                                                                                                                            <  1:506474/93‑1 (MQ=255)                                                                                                                                                                                             cTCCCTGAGCTTCAGCTTCAGCTACCGCTTGATCGATAGCCGCATTGATAGTGTGTTCCGGCATAGCAACTTGTTCCGGGATAGGGTTCGCTACCAAGCGGCAAACGTTAAGACCGCTTTGCCATTTCACCACCATTGAACGGGCACTTTCGCTGGCGCTGTCGCGACGAATGCTGACGTCACACGGGCTGGTGCGGCAGaaaaa  >  1:226436/1‑205 (MQ=255)                                                                                                                                                                                                            cttcagctACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGc                                                                                                                                                 <  1:8645/47‑1 (MQ=255)                                                                                                                                                                                                            cttcagctACCGCTTGATCGATCGCCGCATTGATAGTGTGTTCCGGc                                                                                                                                                 >  2:8645/1‑47 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                 TGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGT‑TGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGT‑CACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCTTCAGCTTCAGCTACCGCCTGATCGATCGCCGCATTGATAGTGTGTTCCGGCATAGCAAACTGTTCCGGGATCGGGTTCGCTACCACGAGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAA  >  NC_000913/2257392‑2257782

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