I’ve got blast output after performing blastx of my sequences against Cazy database
This is how the output looks:
I would like to import this blast.tab output into Megan and visualize the functional entities on EC. But I am not sure which mapping file I should use as CAZY has got its own mapping file
So im confused on which mapping file to use, should I use CAZY mapping file or megan-map-Jul2020-2-X.db.zip? Will I be able to get visualize both the enzyme family and class on EC?
Did you try megan-map-Feb2022.db.zip? That contains mappings of accessions to EC numbers.
That mapping won’t work straight out of the box because MEGAN uses single integers as IDs.
For example, the CAZY mapping file shows:
AA2 NP_001237785.1 220.127.116.11
Opening a MEGAN mapping file with sqlite3, you can see how this is represented there:
SQLite version 3.26.0 2018-12-01 12:34:55
Enter “.help” for usage hints.
sqlite> select EC from mappings where Accession=“NP_001237785”;
If numbers is an array containing the EC numbers, I use the following calculation:
id = 100 * 100 * 10000 * numbers + 100 * 10000 * numbers + 10000 * numbers + numbers
So, if you want to create your own mapping of accessions to EC ids used by MEGAN, apply this formula.