Question: short peptide mapping to Multiple seq alignment (MSA)
0
gravatar for dravinashvasait
6 months ago by
dravinashvasait0 wrote:

Hello all,

I have a set of short peptides (200) comprising 5-30 amino acids that l have to map (align) these peptides against MSA file order to detect polymorphisms present in short peptides. I tried fuzzpro (EMBOSS package) but it takes only one pattern (peptide seq) as input. Could you please recommend me some tools where I can give multiple peptides as input to map against multiple sequence alignment file?

ADD COMMENTlink modified 6 months ago by Elisabeth Gasteiger1.6k • written 6 months ago by dravinashvasait0

This is a Question not a Forum (Forum's are for extended discussion where there may not be a particular or singular right answer). I've changed it for you this time but please bear it in mind for future posts.

ADD REPLYlink written 6 months ago by jrj.healey13k

@jrj.healey. Thank you. I shall keep it in mind for next post. Thank you.

ADD REPLYlink written 6 months ago by dravinashvasait0
0
gravatar for jrj.healey
6 months ago by
jrj.healey13k
United Kingdom
jrj.healey13k wrote:

You can use Clustal for this. It is possible to align sequences against an existing HMM profile or MSA.

Take a look at the --help from clustalo:

Sequence Input:
  -i, --in, --infile={<file>,-} Multiple sequence input file (- for stdin)
  --hmm-in=<file>           HMM input files
  --dealign                 Dealign input sequences
  --profile1, --p1=<file>   Pre-aligned multiple sequence file (aligned columns will be kept fix)
  --profile2, --p2=<file>   Pre-aligned multiple sequence file (aligned columns will be kept fix)
  --is-profile              disable check if profile, force profile (default no)
  -t, --seqtype={Protein, RNA, DNA} Force a sequence type (default: auto)
  --infmt={a2m=fa[sta],clu[stal],msf,phy[lip],selex,st[ockholm],vie[nna]} Forced sequence input file format (default: auto)

Particularly, you'll be interested in the --profile options.

ADD COMMENTlink written 6 months ago by jrj.healey13k
0
gravatar for Elisabeth Gasteiger
6 months ago by
Geneva
Elisabeth Gasteiger1.6k wrote:

UniProt's PeptideSearch doesn't map peptides against a MSA, but can find exact peptide matches (multiple peptides can be entered simultaneously) in the UniProt Knowledgebase: https://www.uniprot.org/peptidesearch

ADD COMMENTlink written 6 months ago by Elisabeth Gasteiger1.6k
Please log in to add an answer.

Help
Access

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