Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,286,815	T→C	H182H (CAT→CAC)	agaI →	galactosamine‑6‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,286,815	0	T	C	100.0%	61.8 / NA	17	H182H (CAT→CAC)	agaI	galactosamine‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/12);  total (5/12)

GATTGCGCGCAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACATGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGA  >  NC_000913/3286680‑3286933                                                                                                                                        |                                                                                                                       gATTGCGCGCAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACgag                                                                                                                     <  1:106545/139‑1 (MQ=255)        cgcAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCATGCCATATCAGTCAACTTGATGCCAGAACACATCAACCAGAGATGTTaa                                                                                                              >  1:538997/1‑139 (MQ=255)                         tttATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGcccc                                                                                             <  1:559030/139‑1 (MQ=255)                           tATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGt                                                                                           <  1:239221/139‑1 (MQ=255)                             tGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGa                                                                                         <  2:544790/139‑1 (MQ=255)                                gTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTc                                                                                      <  2:116182/139‑1 (MQ=255)                                gTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCAGATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTc                                                                                      <  2:381225/139‑1 (MQ=255)                                       gATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                               <  2:12356/139‑1 (MQ=255)                                                        cATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAgg                                                              <  1:130130/139‑1 (MQ=255)                                                         aTCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGa                                                                 >  1:551624/1‑135 (MQ=255)                                                         aTCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGa                                                                 <  2:551624/135‑1 (MQ=255)                                                                    aacgaacCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATg                                                  >  1:60904/1‑139 (MQ=255)                                                                                        gCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTg                              <  2:285599/139‑1 (MQ=255)                                                                                                                 cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTc                                                                                                  >  1:476382/1‑46 (MQ=255)                                                                                                                 cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTc                                                                                                  <  2:476382/46‑1 (MQ=255)                                                                                                                 cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAgg                                                              >  2:151142/1‑82 (MQ=255)                                                                                                                 cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAgg                                                              <  1:151142/82‑1 (MQ=255)                                                                                                                    gATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGa  >  2:43873/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                       GATTGCGCGCAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACATGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGA  >  NC_000913/3286680‑3286933

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