Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,342,094	+G	coding (184/1356 nt)	melA →	alpha‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,342,090	1	.	G	100.0%	197.3 / NA	54	A60A (GCA→GCG)	melA	alpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (29/25);  total (29/25)

AAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCA‑GGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCA  >  NC_000913/4341967‑4342224                                                                                                                              |                                                                                                                                       aaaTATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCaaa                                                                                                                         >  2:230949/1‑140 (MQ=255)      tttttGGTGTTGTGTTCCATCGCGGGGGGCTAAAAAC‑GGCGAATTTTGCCCTGAGGGACATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAGGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCAc                                                                                                                    <  2:395335/137‑1 (MQ=255)      tttttGGTGATGTTTTCCACCGCGAGGCGCTAAAAAC‑GGCGCATATTCCCCTGATGGACATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCAc                                                                                                                    <  2:202014/137‑1 (MQ=255)      ttCTTGGTGATGTGTCCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCAc                                                                                                                    <  2:322317/140‑1 (MQ=255)              gATGTGTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCcaca                                                                                                            >  2:36000/1‑140 (MQ=255)                tgttTTCCATCG‑GGGGCCCTTAAAACGGGCGTATATTGCCCGGATGGACATTCACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTAAGCAGGGGGCCAGCGGCAAAATCACCTGCCACAcc                                                                                                          <  2:648159/140‑1 (MQ=255)                tgtgTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACAcc                                                                                                          >  1:613924/1‑140 (MQ=255)                tgtgTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACAcc                                                                                                          >  2:448761/1‑140 (MQ=255)                tgtgTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACAcc                                                                                                          >  2:737975/1‑140 (MQ=255)                   gtTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCaa                                                                                                       >  2:533440/1‑140 (MQ=255)                   gtTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCaa                                                                                                       >  1:679903/1‑140 (MQ=255)                       cATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTAGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACaga                                                                                                   >  1:202335/1‑140 (MQ=255)                         tCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaa                                                                                                 >  1:163746/1‑140 (MQ=255)                           gcgGGGCGCTAAAACC‑GGCGCATTTTGCCCTGATGGACATGGCCCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaaga                                                                                               <  1:642846/140‑1 (MQ=255)                            cgAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaagaa                                                                                              >  1:323409/1‑140 (MQ=255)                            cgAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaag                                                                                                >  2:808339/1‑138 (MQ=255)                            cgAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaag                                                                                                <  1:808339/138‑1 (MQ=255)                              aGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGCCCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaaga                                                                                               <  1:944078/137‑1 (MQ=255)                              aGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACagaaaga                                                                                               >  2:944078/1‑137 (MQ=255)                              aGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGc                                                                                            >  2:828180/1‑140 (MQ=255)                              aGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGc                                                                                            >  1:133943/1‑140 (MQ=255)                              aGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGc                                                                                            >  2:670583/1‑140 (MQ=255)                                 cgcTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCtt                                                                                         >  1:451728/1‑140 (MQ=255)                                   taaaaaaC‑GGCGATTTTTCCCTTGTTGAACATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTag                                                                                       <  2:451728/137‑1 (MQ=255)                                              cgcATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGa                                                                             <  2:323409/140‑1 (MQ=255)                                                            aTGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGtggtgg                                                               <  2:679903/140‑1 (MQ=255)                                                                 cATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCAttt                                                          >  1:163887/1‑140 (MQ=255)                                                                  aTTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTc                                                         >  1:114509/1‑140 (MQ=255)                                                                      cccccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGAt                                                     <  1:828180/139‑1 (MQ=255)                                                                      aCCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGAt                                                     >  1:260533/1‑140 (MQ=255)                                                                         ccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTgg                                                  <  1:856199/137‑1 (MQ=255)                                                                         ccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTgg                                                  <  1:448761/137‑1 (MQ=255)                                                                         ccACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGAAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTgg                                                  <  1:225584/140‑1 (MQ=255)                                                                          cACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                                                 <  2:260533/140‑1 (MQ=255)                                                                          cACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                                                 <  1:533440/140‑1 (MQ=255)                                                                           aCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGTAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCg                                                >  2:623384/1‑140 (MQ=255)                                                                           aCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCg                                                >  1:164129/1‑140 (MQ=255)                                                                            cccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCgg                                               >  1:270676/1‑140 (MQ=255)                                                                                ccTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTAt                                           >  1:355441/1‑140 (MQ=255)                                                                                          tCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCa                                 <  2:448642/140‑1 (MQ=255)                                                                                             cATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACgg                              <  2:163887/140‑1 (MQ=255)                                                                                              atatTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCACggg                             >  2:152071/1‑139 (MQ=255)                                                                                                atTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGa                           <  2:114509/140‑1 (MQ=255)                                                                                                 ttGTGGTGCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGAc                          <  1:737975/140‑1 (MQ=255)                                                                                                       tgCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACg                               >  2:457303/1‑129 (MQ=255)                                                                                                       tgCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACg                               <  1:457303/129‑1 (MQ=255)                                                                                                        gCGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTc                   <  2:163746/140‑1 (MQ=255)                                                                                                         cGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCg                  >  1:782493/1‑140 (MQ=255)                                                                                                         cGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCg                  <  1:670583/140‑1 (MQ=255)                                                                                                         cGTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGGGACTGATTTCg                  >  1:42203/1‑140 (MQ=255)                                                                                                          gTAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCGTTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGa                 >  1:300255/1‑140 (MQ=255)                                                                                                           tAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAg                <  2:164129/140‑1 (MQ=255)                                                                                                           tAAGCTGATGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTTCACCGGGACTGAtttttag                >  1:964276/1‑136 (MQ=255)                                                                                                                   tGGATTCAGCAGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTaa        <  2:782493/140‑1 (MQ=255)                                                                                                                         cagcagGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGTGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCa  <  1:36000/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                       AAATATTCTTGGTGATGTGTTCCATCGCGAGGCGCTGAAAAC‑GGCGCATATTGCCCTGATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCA‑GGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCA  >  NC_000913/4341967‑4342224

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