Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,342,680	(A)5→4	coding (770/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,676	0	A	.	100.0%	137.4 / NA	33	coding (766/1356 nt)	melA	alpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (17/16);  total (17/16)

CTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAG  >  NC_000913/4342596‑4342748                                                                                 |                                                                         cTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCa                                                                          <  1:893019/81‑1 (MQ=255)    aTGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaa                                                                        <  2:893020/80‑1 (MQ=255)     tGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaa                                                                       <  1:893021/80‑1 (MQ=255)     tGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaa                                                                       <  1:893022/80‑1 (MQ=255)      gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaagc                                                                     >  2:893061/1‑79 (MQ=255)      gAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaagc                                                                     >  1:893073/1‑79 (MQ=255)        aGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCaaaagctg                                                                   <  2:893023/81‑5 (MQ=255)           agggcaggCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGc                                                               >  1:893063/1‑81 (MQ=255)              gcaggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTAt                                                            <  1:893024/81‑1 (MQ=255)              gcaggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTAt                                                            <  2:893025/81‑1 (MQ=255)              gcaggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTAt                                                            <  1:893026/81‑1 (MQ=255)                 ggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTc                                                         <  2:893027/81‑1 (MQ=255)                 ggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTc                                                         >  1:893086/1‑81 (MQ=255)                    aCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGt                                                       >  1:893070/1‑80 (MQ=255)                    aCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTc                                                      <  2:893028/81‑1 (MQ=255)                     ccGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCa                                                     >  2:893082/1‑81 (MQ=255)                        aaGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACgg                                                  <  1:893029/81‑1 (MQ=255)                        aaGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACg                                                   >  2:893071/1‑80 (MQ=255)                          gCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGaa                                                >  1:893075/1‑81 (MQ=255)                            cGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTAGTCACGGAATc                                              >  2:893090/1‑81 (MQ=255)                                     aTGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACa                                     >  2:893058/1‑81 (MQ=255)                                        gCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAttt                                  >  1:893096/1‑81 (MQ=255)                                        gCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACAttt                                  <  2:893030/81‑1 (MQ=255)                                          aaTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTg                                >  1:893072/1‑81 (MQ=255)                                          aaTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTg                                <  2:893031/81‑1 (MQ=255)                                          aaTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTg                                <  2:893032/81‑1 (MQ=255)                                                cGCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAg                           >  1:893089/1‑80 (MQ=255)                                                 gCTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGt                          <  2:893033/80‑1 (MQ=255)                                                  cTGCCAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTac                        <  2:893034/81‑1 (MQ=255)                                                     ccAGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTacac                      <  1:893035/80‑1 (MQ=255)                                                       aGAATATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACAcc                   <  2:893036/81‑1 (MQ=255)                                                         aaTATTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGt                 <  1:893037/81‑1 (MQ=255)                                                           tatTGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTgg               >  1:893074/1‑81 (MQ=255)                                                             ttGTGCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGt              <  2:893038/80‑1 (MQ=255)                                                              tgtgCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGtt             >  1:893080/1‑80 (MQ=255)                                                               gtgCGCTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCTGAACATTTTGCTGAGTACACACCGTGGTtta           >  2:893095/1‑81 (MQ=255)                                                                 gcgcTACGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTttat          >  2:893093/1‑80 (MQ=255)                                                                       cGAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGc     >  2:893088/1‑79 (MQ=255)                                                                        gAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAg  >  2:893101/1‑81 (MQ=255)                                                                        gAAATGTTC‑AAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAg  <  1:893039/81‑1 (MQ=255)                                                                                 |                                                                         CTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAG  >  NC_000913/4342596‑4342748

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