Galaxy |

Protein Domains Filter (version 1.0.0)

Choose primary GFF3 file of all domains from Protein Domains Finder:

Minimum identity:

Protein sequence indentity threshold between input and mapped protein from db [0-1]

Minimum similarity:

Protein sequence similarity threshold between input and mapped protein from db [0-1]

Minimum alignment length:

Proportion of the hit length without gaps to the length of the database sequence [0-1]

Interruptions [frameshifts + stop codons]:

Tolerance threshold per every starting 100 amino acids of alignment sequence

Maximal length proportion:

Maximal proportion of alignment length to the original length of protein domain from database (including indels)

Select protein domain type:

Filter based on classification:

You can use preset options or enter an arbitrary string to filter a certain repetitive element type of any level. It must be a continuous substring in a proper format of Final_Classification attribute of GFF3 file. Classification levels are separated by | character. Filtering is case sensitive

WHAT IT DOES

This tool runs filtering on either primary GFF3 file of all domains, i.e. output of Protein Domains Finder tool or already filtered GFF3 file. Domains can be filtered based on:

Quality of alignment such as:

alignment sequence identity
alignment similarity
alignment proportion length
number of interruptions (frameshifts or stop codons) per 100 AA

Protein domain type

This filtration is based on "Name" attribute of GFF3 file

Repetitive element classification

In the text field you can specify a classification string you wish to filter. This filtration is based on "Final_Classification" attribute of GFF file, so it must be in the proper form (classification levels are separated by "|"). You can see which classifications occurs in your data taking a look into Classification summary table output. If you leave the field blank, domains of all classifications will be reported

All the records containing ambiguous domain type (e.g. RH/INT) are filtered out automatically. They do not take place in filtered gff file neither the protein sequence is derived from these potentially chimeric domains. Optimal results (for general usage) should be reached using the default quality filtering parameters which are appropriate to find all types of protein domains. Keep in mind that the results should be critically assessed based on your input data anyhow.

OUTPUTS PRODUCED:

Filtered GFF3 file
Translated protein sequences of the filtered domains regions of original DNA sequence in fasta format

Translated sequences are taken from the best alignment (Best_Hit attribute) within a domain region, however this alignment does not necessarily have to cover the whole region reported as a domain in gff file