Hi,
I was wondering whether there's a way to perform protein BLAST searches allowing for particular positions to be degenerate. Nucleotide BLAST allows the usage of the IUPAC degenerate namings (N = A/C/G/T; K = G/T; and so on). For instance, say I have the following canonical sequence:
KRSCATTMG
And I want to perform a BLAST search with the first position being either K or R, and positions 6 and 7 being either T or S.
I'm aware BLAST allows for certain variability within the sequence, what I need to do is biasing such allowance in particular positions towards particular amino acids, while not allowing variability in other positions (in my example, for instance, I need position 4 to be a C).
Maybe a faster approach would be to just download BLAST databases and conditional filter through a python script?
1 answer
Even if BLAST can do what you want, it is most likely not the best choice for it, especially if you have short search strings.
You may want to check out regular expressions like the ones used in PROSITE. For what you want, the pattern would be something like [KR]-R-S-C-A-[ST](2)-MG. There are programs that can search databases with these patterns:
Log in to answer this question.
-
On retrieving SNPs from dbSNP
written by Agenor Neto •
-
What's best practice when determining a degenerate base when creating a primer?
written by O.rka
-
Subsetting BAM by positions
written by בת אל •
-
bcftools merge heterozygous calls
written by evelyn
-
Annotate genomic positions with dbSNP rsIds
written by Jimbou
-
How to identify contigs containing the degenerate primer sequences using python?
written by xupbuy •
-
HOW EDIT VCF FILE COLUMNS
written by misbahabas •
-
generate consensus sequence from unphased heterozygous calls
written by atul_shah1 •
-
Non-standard nucleotides in my fast converted from a bam
written by biogirl
-
Mapping Genomic Position To Protein Position
written by John Knowlg •
