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%	146.6 / NA	39	H182H (CAT→CAC)	agaI	galactosamine‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/24);  total (15/24)

CGCAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACATGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGACAGATCGTT  >  NC_000913/3286687‑3286942                                                                                                                                 |                                                                                                                                cgcAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTaa                                                                                                                       >  2:14267/1‑139 (MQ=255)              tGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCg                                                                                                          >  1:59256/1‑139 (MQ=255)                         gTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTc                                                                                               <  2:1006099/139‑1 (MQ=255)                         gTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTc                                                                                               <  1:413819/139‑1 (MQ=255)                         gTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTc                                                                                               <  1:817982/139‑1 (MQ=255)                                gATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  2:13707/139‑1 (MQ=255)                                gATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  1:146777/139‑1 (MQ=255)                                gATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  1:985065/139‑1 (MQ=255)                                gATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  1:382036/139‑1 (MQ=255)                                     gggAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCt                                                                                   >  1:230928/1‑139 (MQ=255)                                     gggAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCt                                                                                   <  1:1070657/139‑1 (MQ=255)                                       gAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTa                                                                                 >  2:868966/1‑139 (MQ=255)                                       gAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTa                                                                                 <  1:310131/139‑1 (MQ=255)                                            aCGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCt                                                                            <  1:868966/139‑1 (MQ=255)                                                tCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGa                                                                                                  >  1:1089000/1‑113 (MQ=255)                                                tCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGa                                                                                                  <  2:1089000/113‑1 (MQ=255)                                                             aacgaacCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACgag                                                                                                                              >  1:819068/1‑72 (MQ=255)                                                             aacgaacCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACgag                                                                                                                              <  2:819068/72‑1 (MQ=255)                                                                    cGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGa                                                    <  2:528851/139‑1 (MQ=255)                                                                    cGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGa                                                    <  2:230928/139‑1 (MQ=255)                                                                           aaGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTg                                             <  1:14267/139‑1 (MQ=255)                                                                            aGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        >  1:179644/1‑95 (MQ=255)                                                                            aGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  2:179644/95‑1 (MQ=255)                                                                            aGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGt                                            >  2:850410/1‑139 (MQ=255)                                                                                      cGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGt                                                                                                            >  1:274863/1‑65 (MQ=255)                                                                                      cGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGt                                                                                                            <  2:274863/65‑1 (MQ=255)                                                                                              cATATCAGTCAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAg                          <  2:59256/139‑1 (MQ=255)                                                                                                       cAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGa                                                                          <  1:888284/82‑1 (MQ=255)                                                                                                       cAACTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGACCTGa                                                                          >  2:888284/1‑82 (MQ=255)                                                                                                          cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        >  2:313104/1‑65 (MQ=255)                                                                                                          cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGAt                                                                                        <  1:313104/65‑1 (MQ=255)                                                                                                          cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGa                                                    <  2:122564/101‑1 (MQ=255)                                                                                                          cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGa                                                    >  1:122564/1‑101 (MQ=255)                                                                                                          cTTGATGCCAGAACACAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATg              >  2:148453/1‑139 (MQ=255)                                                                                                                      acacAGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGACAGATCGtt  <  1:173706/139‑1 (MQ=255)                                                                                                                          aGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGa                                                                          >  2:532279/1‑63 (MQ=255)                                                                                                                          aGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGa                                                                          <  1:532279/63‑1 (MQ=255)                                                                                                                          aGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGaa                       <  1:696606/114‑1 (MQ=255)                                                                                                                          aGCAACACGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGaa                       >  2:696606/1‑114 (MQ=255)                                                                                                                                 |                                                                                                                                CGCAAAGGCGGTCTGGATTTATGCGTTCTCGGATTGGGGAAAAACGGTCATCTTGGGCTGAACGAACCGGGAGAAAGCCTGCAACCGGCCTGCCATATCAGTCAACTTGATGCCAGAACACAGCAACATGAGATGTTAAAAACCGCGGGTCGCCCCGTGACTCGTGGGATCACCTTAGGCCTGAAGGATATTCTCAATGCCCGCGAAGTTTTGTTACTGGTGACTGGCGAAGGAAAGCAGGATGCGACAGATCGTT  >  NC_000913/3286687‑3286942

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