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
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.

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.