**Hi everyone, **
This is my first time in Megan,
I and running shotgun data using paired-end reads
I used this
###1.Install Megan
MEGAN_Community_unix_6_18_0.sh
###2.This file in neccesary is we want meganizer
megan-map-Oct2019.db
###3.Download and install DIAMOND
###4.Download the NR database from here
ftp://ftp.ncbi.nih.gov/blast/db/FASTA/nr.gz
###5.Build a DIAMOND index
diamond makedb -d nr --in nr.gz
###6.This results in an index file
nr.dmnd
###7.Run DIAMOND (in query we use the name of our sequences) - ESTO SE DEMORÓ CORRIENDO 15 DÍAS
diamond blastx --db /vault2/bioinfDBs/diamondDB/nr.dmnd --out reads_shotgun.daa \ --outfmt 100 --query /vault2/homehpc/voteroj/SEQ_SHOTG/T4P4_T1_TCCGCGAA-GTACTGAC_L004_R2_001.fastq.gz
###8.This will produce a DAA file called
reads_shotgun.daa
###9. Before run MEGANIZER, we need the following database, we can download from Megan website:
###The following file maps NCBI-nr accessions to taxonomic and functional classes (eggNOG, InterPro2GO and SEED), unzip before use:
megan-map-Oct2019.db.zip
And then when I tried to meganizer using this script:
nohup /vault2/soft/bio/megan6/tools/daa-meganizer -i reads_shotgun.daa \ -mdb megan-map-Oct2019.db -v &
The result obtained; attached at the end of the message; I do not how to figure out, Do you have any suggestion, please.
**Thanks in advance, **
Best regards,
Vanessa O.J.
Version MEGAN Community Edition (version 6.18.9, built 4 Apr 2020)
Copyright © 2020 Daniel H. Huson. This program comes with ABSOLUTELY NO WARRANTY.
GLib-GIO-Message: Using the ‘memory’ GSettings backend. Your settings will not be saved or shared with other applications.
Java version: 11.0.1
Functional classifications to use: EGGNOG, INTERPRO2GO, SEED
Loading ncbi.map: 2,175,506
Loading ncbi.tre: 2,175,510
Loading eggnog.map: 30,875
Loading eggnog.tre: 30,986
Loading interpro2go.map: 13,501
Loading interpro2go.tre: 29,204
Loading seed.map: 13,662
Loading seed.tre: 21,085
Meganizing: reads_shotgun.daa
Meganizing init
Caught:
java.io.IOException: Invalid DAA file. DIAMOND run probably has not completed successfully.
at megan/megan.daa.io.DAAHeader.load(DAAHeader.java:163)
at megan/megan.daa.io.DAAModifier.removeAllMEGANData(DAAModifier.java:37)
at megan/megan.daa.DAAReferencesAnnotator.apply(DAAReferencesAnnotator.java:55)
at megan/megan.daa.Meganize.apply(Meganize.java:70)
at megan/megan.tools.DAAMeganizer.run(DAAMeganizer.java:246)
at megan/megan.tools.DAAMeganizer.main(DAAMeganizer.java:65)