SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Help with ABACAS kyoru Bioinformatics 1 07-12-2013 12:25 PM
Split a SAM file rahul Bioinformatics 6 12-20-2011 11:12 AM
split a fastq file lfaino Bioinformatics 4 04-14-2011 03:28 PM
FACS sorting chromosomes? k-gun12 General 0 03-19-2011 07:48 AM
All chromosomes and maq Layla Bioinformatics 1 02-13-2009 08:36 AM

Reply
 
Thread Tools
Old 01-30-2013, 11:18 AM   #1
yzzhang
Member
 
Location: florida

Join Date: Jan 2013
Posts: 67
Default split the abacas joined chromosomes

Hi, All. I just tried to align my contigs to a genome which has 9 chromosomes using ABACAS. I followed the PAGIT toolkit published on nature protocol, but after I joined the chromosomes together as required by abacas, I could not split the files, the script in the PAGIT "split2multifasta.pl did not work. Does some one have suggestions? Thank you.
yun
yzzhang is offline   Reply With Quote
Old 02-04-2013, 07:53 PM   #2
seb567
Senior Member
 
Location: Québec, Canada

Join Date: Jul 2008
Posts: 260
Default

Quote:
Originally Posted by yzzhang View Post
Hi, All. I just tried to align my contigs to a genome which has 9 chromosomes using ABACAS. I followed the PAGIT toolkit published on nature protocol, but after I joined the chromosomes together as required by abacas, I could not split the files, the script in the PAGIT "split2multifasta.pl did not work. Does some one have suggestions? Thank you.
yun
Can you provide details about the statement "PAGIT split2multifasta.pl did not work" ?
seb567 is offline   Reply With Quote
Old 02-05-2013, 04:20 AM   #3
yzzhang
Member
 
Location: florida

Join Date: Jan 2013
Posts: 67
Default

I looked through the split2multifasta.pl script, there were only few lines like that
use strict;


if (scalar(@ARGV)<5) {
die "usage: <Reference File> <ordered ABACAS query file> <resultName> <crunchfile> <Tag>\n\nThis script will split the query file accordinly to the reference, assuming the reference got processed with joinMultifasta.pl.\n";
}

my $spacer;
my $splitinfo = shift;
my $queryFile = shift;
my $crunchFile = shift;
my $resultName = shift;

print "Sorry, this script has not yet be included properly to the PAGIT distribution.\n";
.I email the author, and he kindly sent me the script and now I can use this pipeline now.
yzzhang is offline   Reply With Quote
Old 09-03-2013, 03:46 PM   #4
flipwell
Member
 
Location: Queensland

Join Date: Feb 2011
Posts: 14
Default

Has anyone had any trouble running this script? When I try to run it nothing happens. I don't get any errors it just seems to sit there and produce nothing. Does it take a long time to run?

Thanks
flipwell 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 08:36 AM.


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