Question: How to extract the mobile genetic elements
gravatar for anabio86
3 months ago by
anabio860 wrote:

Hi folks, I need to do alinhamento using just the chromosomal DNA, without the mobile genetic elements, I would like to know if there is a bioinformatic tools that can do it from fasta files? thx

sequencing alignment assembly • 226 views
ADD COMMENTlink modified 10 weeks ago • written 3 months ago by anabio860

I guess you could subset your reference FASTA file to include only the sequences you wish to include. What have you tried?

ADD REPLYlink written 3 months ago by RamRS20k

Are you referring to masking parts of the genome that are repetitive before aligning? repeatmasker is the program used to do this. If you are working with human genome then repeat-masked versions are available.

ADD REPLYlink written 3 months ago by genomax60k

Actually, I am working with bacterial DNA, I need to perform a phylogenetic analysis using references from database, but before to compare with these references I need to remove the sequences that belonged to mobile genetic elements such as plasmids, phages this things

ADD REPLYlink written 3 months ago by anabio860

mobile genetic elements

I think this can lead to some confusion. As a human geneticist I think of transposons when I read mobile genetic elements, and not phages/plasmids/....

ADD REPLYlink written 3 months ago by WouterDeCoster35k

We need lots more information, such as what bacteria, whether it has completed reference genomes etc.

Detecting mobile elements is non-trivial. If your genome is not well curated already the task will be almost impossible. There is also not one single tool which will do all of this. Plasmid finding has its own set of tools (e.g. PlasFlow, PlasmidTron and others), prophages can be detected with PHAST/PHASTER (e.g.). I'm no expert on detection of transposons etc, but there are reviews out there which compare and benchmark many tools (e.g ).

You will largely have to try things and see what the results look like.

ADD REPLYlink written 3 months ago by jrj.healey9.7k

I just gave up, seems impossible to do that

ADD REPLYlink written 10 weeks ago by anabio860

Nothing is impossible. Why Portuguese ('alinhamento')?

ADD REPLYlink written 10 weeks ago by Kevin Blighe35k

If I understand correctly, you are looking for FASTA files of chromosomes but not plasmids, right? If so, you can download them from NCBI using edirect as follows:

esearch -db assembly -q 'GCF_000008865.2' | elink -db assembly -target nuccore -name assembly_nuccore_refseq | efilter -q 'NOT plasmid[filter]' -source refseq | efetch -format fasta
ADD REPLYlink written 10 weeks ago by vkkodali860
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1657 users visited in the last hour