SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
create scaffold file from contig file and agp file? Wallysb01 Bioinformatics 2 05-10-2016 02:53 AM
how to create output directory folder on velvet suhalsuhaimi Bioinformatics 5 12-28-2014 05:04 AM
Create .bed file from Genbank file ksw9 Bioinformatics 4 09-07-2014 03:52 PM
How to create count table for SOAP output? ahmadsam Bioinformatics 1 08-26-2014 09:51 AM
parsing blast output using filters using regular perl newbie2this Bioinformatics 0 11-10-2012 06:02 AM

Reply
 
Thread Tools
Old 06-30-2015, 01:34 AM   #1
vineetha
Junior Member
 
Location: India

Join Date: Jun 2015
Posts: 2
Default how to create an output file using perl

I have two files.One file contains blast result of contigs and the next file contains the contigs that are used for blast.I want an output text file of sequences in a separate folder which gives no hits in blast result from original fasta file with all contigs.I have no idea about this.Can you please help me to create perl script for this problem.

Help me pleaseee...
Attached Files
File Type: txt sequences.txt (8.5 KB, 6 views)
File Type: txt blastresult.txt (6.4 KB, 11 views)

Last edited by vineetha; 06-30-2015 at 10:22 PM.
vineetha is offline   Reply With Quote
Old 06-30-2015, 02:25 AM   #2
wyzandrea
Banned
 
Location: USA

Join Date: Apr 2015
Posts: 2
Default

I want to ask this question too.
wyzandrea is offline   Reply With Quote
Old 07-01-2015, 12:06 PM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

I don't think anyone wants to write a blast parser for you. If you are using blast from the command line, you can tell it to format the output a number of ways, one of those ways might be easier to parse for what you want.
swbarnes2 is offline   Reply With Quote
Old 07-01-2015, 03:51 PM   #4
SNPsaurus
Registered Vendor
 
Location: Eugene, OR

Join Date: May 2013
Posts: 521
Default

If you want an output file with the "no hits", you could just grep the blast file for "no hits", with -B to grab the query, then grep that to get all the querys only, then grep the blast file with the query file using -f and with grep -A to get the sequence itself. Need to convert the fasta file to single line fasta, though.
__________________
Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com
SNPsaurus 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 12:15 AM.


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