a protein failed to be aligned to genome by tblastn
Hi all:
I have a protein sequence which I know the genomic region where it came from. also I have the CDS.
by blasting the CDS to the genome, I can find the exact position of that transcript. however, the first 200 AA of the protein could not be aligned to the corresponding region by tblastn.
Could any one help me to understand what happened ? Thanks in advance! Best,
2 answers
I solved this problem by turing off two options of tblastn: -seg and -comp_based_stats.
the former was for filtering low-complexity sequence.
Try tblastx. Usually it is much more successful than tblastn.
• 1 views
•
link
written
by
natasha.sernova
Log in to answer this question.
More posts like this
-
Visualize and annotate protein domains within a genomic context, given complement
written by Madde •
-
how to get genomic position from protein pdb ID+AA residue position within the protein
written by lait
-
Extracting a Genomic Region From Nucleotide Accession
written by 13dsc •
-
tophat: the Searching for junctions via segment mapping was skipped.
written by Pei
-
How to convert exon genomic coordinates to protein coordinates
written by carlota.rubio •
-
Representative AA change from point mutation given by genomic lcus
written by sak042 •
-
tophat, the parameter "Expected inner distance between mate pairs"?
written by Pei
-
How to annotate a sequence with its genomic coordinates ?
written by Pei
-
how to annotate a genomic internval with its gene ID, sequence, seq-type
written by Pei
-
Obtaining Aa Position From Snpeff Output:
written by Angel
I think their might be a large intronic region after 200 AA as BLAST looks for local alignment and the region for the 200 AA is little bit far upstream BLAST might have missed it. try tblastn using only the 200 AA and it may give alignment upstrem of your present alignment.
Hope this works.
Thanks! I have tried this but found that the first 200 AA could not be aligned to genome by tblastn when using this 200AA as query alone. as a control, another 200AA in the middle of the protein could be aligned to the correct position with a %identity < 60%.