Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,282,062	T→C	R29R (CGT→CGC)	agaS →	putative galactosamine‑6‑phosphate deaminase/isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,282,062	0	T	C	100.0%	70.7 / NA	21	R29R (CGT→CGC)	agaS	putative galactosamine‑6‑phosphate deaminase/isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/8);  total (13/8)

ACAACGCCCGACGTCAAGTTCATCAGACTAAGGATTGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCG  >  NC_000913/3281935‑3282177                                                                                                                                |                                                                                                                    aCAACGCCCGACGTCAAGTTCATCAGACTAAGGATTGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCtcact                                                                                                                >  1:517416/1‑133 (MQ=255)            cGTCAAGTTCATCAGACTAAGGATTGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCg                                                                                                     <  2:413052/133‑1 (MQ=255)                                ggATTGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTc                                                                                 >  2:49832/1‑133 (MQ=255)                                 gATTGAGTCATGCCACAATATTCCACCCCTTCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCa                                                                                <  1:185533/133‑1 (MQ=255)                                   ttGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAAt                                                                              <  2:98947/133‑1 (MQ=255)                                   ttGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAAt                                                                              <  1:394989/133‑1 (MQ=255)                                    tGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATa                                                                             >  2:209439/1‑133 (MQ=255)                                           tGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCt                                                                      >  2:6043/1‑133 (MQ=255)                                                  aaaTTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCg                                                               <  1:209439/133‑1 (MQ=255)                                                             ctgctgCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCTCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCaa                                                    >  2:317559/1‑133 (MQ=255)                                                             ctgctgCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCACAATAACATCCCTGCACCGTTCCTGCGCaa                                                    >  1:260452/1‑133 (MQ=255)                                                                      ccgcAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCt                                           <  1:49832/133‑1 (MQ=255)                                                                                tACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACCCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGatcatc                                 >  2:67247/1‑133 (MQ=255)                                                                                     ggACTGAAGACGAGATCCGCCATCTGCCTCGCGCATGGATCAGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGAc                            <  2:260452/133‑1 (MQ=255)                                                                                     ggACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCTCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGAc                            <  1:317559/133‑1 (MQ=255)                                                                                     ggACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGAc                            >  1:316188/1‑133 (MQ=255)                                                                                          gaagaaGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCTAAGAGAAACTGCGGGTCATCCTGACCGGAg                       >  1:227534/1‑133 (MQ=255)                                                                                          gaagaaGAGATCCGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAg                       >  1:451031/1‑133 (MQ=255)                                                                                                     ccGCCATCAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCg            >  2:571107/1‑133 (MQ=255)                                                                                                            cAGCCTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCGTCCTGACCGGATCCGGAACGTTGGCATTTa     >  1:196598/1‑133 (MQ=255)                                                                                                               ccTCGCGCATGGATCCGCTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCg  >  2:509720/1‑133 (MQ=255)                                                                                                                                |                                                                                                                    ACAACGCCCGACGTCAAGTTCATCAGACTAAGGATTGAGTTATGCCAGAAAATTACACCCCTGCTGCCGCCGCAACCGGTACATGGACTGAAGAAGAGATCCGCCATCAGCCTCGCGCATGGATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCG  >  NC_000913/3281935‑3282177

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