SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Interpretation of bwa output - has it used all the reads gwilymh Bioinformatics 0 10-16-2014 08:38 AM
Quick Interpretation of FastQC output rnastar Bioinformatics 4 11-22-2013 11:43 AM
interpretation of output from HTSeq syintel87 Bioinformatics 1 01-06-2013 01:12 PM
Varscan-output interpretation bioman1 Bioinformatics 0 05-24-2012 12:50 PM
SAMtools flagstat output interpretation a2z@blr Bioinformatics 2 10-20-2011 01:23 PM

Reply
 
Thread Tools
Old 03-05-2015, 11:35 AM   #1
syintel87
Member
 
Location: Universe

Join Date: Dec 2012
Posts: 81
Default tblastx fmt1 Output Interpretation

Hello,
While locally running "blastx" with output option "-outfmt 1", I have a question.

The below is an example of my result. As seen, throughout the outputs, all pairs of anchored proteins are not matched. How can I interpret these non-identical alignment?

Quote:
Query_3 477 GSVEYVHMLNGTMCATTRTVCAL 409
medicagoRALF3 38 .GM.WI.QTKTAT.EGSIAD.M. 60

Query_4 34 PYHLLLQKFYKT 69
medicagoRALF11 104 ..NRGCS.Y.RC 115
Thank you in advance.
Attached Images
File Type: jpg blast.jpg (50.2 KB, 4 views)

Last edited by syintel87; 03-06-2015 at 07:56 AM.
syintel87 is offline   Reply With Quote
Old 03-05-2015, 11:58 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,049
Default

Are you searching with very short query sequences (like illumina reads)?
GenoMax is offline   Reply With Quote
Old 03-06-2015, 07:26 AM   #3
syintel87
Member
 
Location: Universe

Join Date: Dec 2012
Posts: 81
Default

My query sequences are contigs that are de novo assembled, some of which are short (e.g. 500) whereas others are very long (e.g. 100,000).

But database of blast is composed of several peptide proteins whose length is short (e.g. 44).
syintel87 is offline   Reply With Quote
Old 03-06-2015, 07:36 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,049
Default

Have you tried to do the search the other way around (using your peptides as query)?

Try using BLAT too. Especially if you know that you expect the peptides to be there in your data.
GenoMax is offline   Reply With Quote
Old 03-06-2015, 08:21 AM   #5
syintel87
Member
 
Location: Universe

Join Date: Dec 2012
Posts: 81
Default

In addition to "-outfmt 1", I tried other options as well.
There seem to be different ways of alignment.

0 = pairwise
1 = query-anchored showing identities
2 = query-anchored no identities
3 = flat query-anchored, show identities
4 = flat query-anchored, no identities

Quote:
[-outfmt 0]
Query_2 229 KMSFRYLFFAIKKYALSKF 173
thalianaRALF4 5 ...LTS...VSIVIV..L. 23

[-outfmt 1]
Query_2 229 KMSFRYLFFAIKKYALSKF 173
thalianaRALF4 5 ...LTS...VSIVIV..L. 23

[-outfmt 2]
Query_2 229 KMSFRYLFFAIKKYALSKF 173
thalianaRALF4 5 KMSLTSLFFVSIVIVLSLF 23

[-outfmt 3]
Query_2 229 KMSFRYLFFAIKKYALSKF 173
thalianaRALF4 5 ...LTS...VSIVIV..L. 23

[-outfmt 4]
Query_2 229 KMSFRYLFFAIKKYALSKF 173
thalianaRALF4 5 KMSLTSLFFVSIVIVLSLF 23
The results of [-outfmt 2] and [-outfmt 4] may be the results that I look forward to getting. However, I still cannot understand the principles and differences that distinguish output format 0 to 4.
syintel87 is offline   Reply With Quote
Old 03-06-2015, 08:33 AM   #6
syintel87
Member
 
Location: Universe

Join Date: Dec 2012
Posts: 81
Default

Ah, now I see! In formats "0, 1, and 3", dots stand for identities between query and target. And differences are shown with protein letters.

I should have posted after more consideration.
Thanks GenoMax! I am going to try "BLAT", too.
syintel87 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:47 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO