Mobile element discovery with MELT brings different results when chromosomes are analysed separate
0
0
Entering edit mode
6 months ago

I get around 5x less mobile elements when calling each chromosome separate compared to when monolith genome .bam is supplied. What is the reason for that? Here is my example command:

java -jar ${meltDir}/MELT.jar Single -h ${refFasta} -bamfile ${id}.${chr}.bam -n ${refDir}/hg38.genes.bed -t mei_list.txt -w . -d 10000
melt • 351 views
ADD COMMENT
0
Entering edit mode

A quick look at the manual indicates that MELT is designed to be run on genome(s).

MELT is optimized for performing discovery in a large number of individual sequencing samples using the Sun Grid Engine (SGE). MELT also has two additional workflows: analysis without SGE (for adaptability to other parallel computing platforms) and single genome analysis.

So running it separately on each chromosome probably treats them as different genomes. So you are not using the tool as designed/intended.

ADD REPLY
0
Entering edit mode

It's slower when computing whole genome. I also haven't tried running multiple genomes at the same time, that could bring a different results.

ADD REPLY

Login before adding your answer.

Traffic: 2458 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6