Tweet
I want to to perform repeatmasking on ant genome. Ant is also included in the repeat library obtained from Repbase. The library RepeatMaskerLib.embl also contains my species. However, I also generated de novo repeat library by repeatmodeler (consensi.fa). How can I use both library for masking ant genome ? Should I concatenate RepeatMaskerLib.embl and consensi.fa or can I use both library with single command?
-
-
You can go both ways. Have not done this myself but the simplest way would be to pull put ant repeats from RepeatMaskerLib.embl and merge them with your custom repeat library (as fasta). You could then use the -lib option in RepeatMasker to supply the updated custom repeat library and mask the ant genome.
AFAIK, you cannot combine both libraries in a single command. -
Hey, I think this would be a great solution, but I have a problem...
I want to mask the repeats, but I want to classify them as well.
If I get the Repbase library in fasta format (to concatenate it with my own database and use it with -lib) then I lose the classificators (ex: #LINE/L1 --> ).
[I don't know if anyone came out with a script to put the "#" and "/" where they should be, but I don't know how, because usually there's only the Superfamily but not the Order (ex L1 but not LINE)]
I also tried to convert RepeatMaskerLib.embl to .fasta but I lose the classification info too.
I guess it won't work to put my repeats, in .embl, into RepeatMaskerLib.embl, but I haven't tried yet.
Did anyone find a way to solve this?
Thank you for any helpComment
-
or maybe...
Does anyone know how to use RepeatClassifier (part of RepeatModeler) only. Or make RepeatModeler classify a bunch of repeats already listed?
ThanksComment
-
RepeatClassifier can be run independently. If you hand it a fasta file it will generate a new file with a *.classified suffix. The output file will be a fasta file with all the id#class/subclass identifiers.Comment
-
Check out the util/buildRMLibFromEMBL.pl script included in the RepeatMasker distribution. This will convert the RepeatMaskerLib.embl file into a fasta file with the "#class/subclass" nomenclature.Comment
-
Thank you for your help.
I manage to get my library classified, but I still have another doubt.
When I use Repeatmasker with a custom library (-lib) the .tbl file seems to be made with the results of masking the genome with
Homo sapies repeats alone ("The query species was assumed to be homo"). My guess is that I can't obtain a .tbl output with a custom library (even if this library is made by RepeatModeler). Is this true?
Thank you again
NuriaComment
-
If the answer to my previous question is yes (so I can't get a real .tbl output when using RepeatMasker with the -lib option) is there any script that can produce something similar to the .tbl output with RepeatMasker .out or .gff files?
(My library is classified with #class/subclass identifiers).
Thank you for your time
NuriaComment
-
Yes. You can use the RepeatMasker/util/buildSummary.pl script to reprocess the *.out file and summarize the results. The buildSummary.pl creates a similar table to the *.tbl file but also includes tabulations of each individual repeat and each seq/chr of the input.Comment
-
Thank you, buildSummary seems to be what I was looking for, but I have another question.
I masked the same genome with different libraries and I used the buildSummary script to extract the repeat information from the *.out files.
But every .tbl file considers a different genome length. It is a number between the original total length and the excluding N/X-runs value. In my case 10 to 20 Mb smaller that the original length.
Consequently the % of repeats have important differences.
Can anyone help?
Thank you very much
NuriaComment
-
ok, I only had to read the buildSummary.pl script.
If someone else has the same problem, I solved using:
perl buildSummary.pl -useAbsoluteGenomeSize -genome tablewithscaffoldsnamesandlength.tsv -species yourspecies RMaskeroutput.out
NuriaComment
-
Hi Nuria,
I'm having the same problem. Just to be clear, you masked your genome using the RepeatModeler results and then used only the buildSummary.pl script to classify the repeats?
Thanks!Comment
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...Yesterday, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM -
Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.
Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...05-06-2026, 09:04 AM
Comment