Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,286,411	A→C	N48H (AAT→CAT)	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,411	0	A	C	100.0%	167.7 / NA	47	N48H (AAT→CAT)	agaI	galactosamine‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/25);  total (22/25)

AATGGAACGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGCCAG  >  NC_000913/3286269‑3286558                                                                                                                                               |                                                                                                                                                    aaTGGAACGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTg                                                                                                                                             <  2:303338/151‑1 (MQ=255)     gAACGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGAttt                                                                                                                                         <  1:415961/151‑1 (MQ=255)       aCGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGc                                                                                                                                       >  2:85007/1‑151 (MQ=255)        cGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGcc                                                                                                                                      <  2:348821/151‑1 (MQ=255)         gAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCt                                                                                                                                     <  2:296513/151‑1 (MQ=255)                 gCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGAc                                                                                                                               <  2:654809/149‑1 (MQ=255)                    tCTGGGGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGGGATTTGCCTGGCGAc                                                                                                                               >  1:350149/1‑146 (MQ=255)                       ggtggtGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCa                                                                                                                       <  2:915234/151‑1 (MQ=255)                             gCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa                                                                                                                  <  2:610861/150‑1 (MQ=255)                                  tttACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACAGGAGCCACGCCATTACTg                                                                                                            >  1:528688/1‑151 (MQ=255)                                    tACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGAc                                                                                                          >  2:876607/1‑151 (MQ=255)                                     aCTTAAGGAATTTGATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACg                                                                                                         >  1:226875/1‑151 (MQ=255)                                     aCTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACg                                                                                                         >  2:90005/1‑151 (MQ=255)                                         aaGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATc                                                                                                     >  1:266962/1‑151 (MQ=255)                                             aaTTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTa                                                                                                 <  1:310549/151‑1 (MQ=255)                                                   aTCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGt                                                                                           <  2:757243/151‑1 (MQ=255)                                                     ccTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTa                                                                                          <  1:489456/150‑1 (MQ=255)                                                     ccTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAg                                                                                         <  1:900862/151‑1 (MQ=255)                                                      cTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGGGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGt                                                                                           >  2:306177/1‑148 (MQ=255)                                                          gCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGaaaaa                                                                                    <  2:796906/151‑1 (MQ=255)                                                          gCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGaaaaa                                                                                    <  2:929526/151‑1 (MQ=255)                                                            aaaCCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAAt                                                                                  <  2:304244/151‑1 (MQ=255)                                                                 cTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATccacc                                                                             >  1:506712/1‑151 (MQ=255)                                                                    cAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACcagc                                                                          <  2:757083/151‑1 (MQ=255)                                                                       cAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACcagcagc                                                                       >  1:569382/1‑151 (MQ=255)                                                                         cgTTGACAACTATACGGGGTTAAGGGAACGTGCCAGCCAATATTTATAGCACGTAATCCATAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCAACGCCATTACTGACGCATCATTATCTGGTAGAAAAAATCCACcagcagc                                                                       >  1:177253/2‑149 (MQ=255)                                                                             gAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGCGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGtt                                                                  <  1:798615/150‑1 (MQ=255)                                                                               aaaCTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGAt                                                               >  1:901433/1‑151 (MQ=255)                                                                                               aGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCt                                               >  2:592624/1‑151 (MQ=255)                                                                                                    aCGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTg                                          >  2:489804/1‑151 (MQ=255)                                                                                                        gCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAg                                       <  2:514383/150‑1 (MQ=255)                                                                                                          cAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTc                                    <  1:799016/151‑1 (MQ=255)                                                                                                              gAATATTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTcgacg                                >  1:216017/1‑151 (MQ=255)                                                                                                                atatTTATTGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGaa                              <  1:300449/151‑1 (MQ=255)                                                                                                                 tatttattGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAAt                             >  1:51361/1‑151 (MQ=255)                                                                                                                 tatttattGGCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAAt                             >  2:761719/1‑151 (MQ=255)                                                                                                                  atttattGGCCGTGGTCCGTAGCAAACCGCTTGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCATCTCACCTTCGTGAAGCTCGACGAAt                             <  2:498153/150‑1 (MQ=255)                                                                                                                          gCCGTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATc                    >  1:274369/1‑151 (MQ=255)                                                                                                                             gTGATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGc                 >  2:405115/1‑151 (MQ=255)                                                                                                                               gATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGcc                <  1:298579/150‑1 (MQ=255)                                                                                                                               gATCCGTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGcc                <  2:693533/150‑1 (MQ=255)                                                                                                                                    gTAGCAAACCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAAc          <  1:81519/151‑1 (MQ=255)                                                                                                                                         aaaCCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGc     >  2:362743/1‑151 (MQ=255)                                                                                                                                           aCCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGcc    >  2:697813/1‑150 (MQ=255)                                                                                                                                           aCCGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGCCa   <  1:90003/151‑1 (MQ=255)                                                                                                                                            ccGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGcc    >  2:146364/1‑149 (MQ=255)                                                                                                                                            ccGCATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGCCAg  <  2:21817/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    AATGGAACGAGGCACTGCGTCTGGTGGTGCCTCTTTACTTAAGGAATTTCATCCTGTGCAAACCCTTCAGCAAGTTGAAAACTATACGGCGTTAAGTGAACGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAGCAGGTTGATGTCAGCCAGCTCACCTTCGTGAAGCTCGACGAATGGGTGGATCTGCCATTAACGATGCCAG  >  NC_000913/3286269‑3286558

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