SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why isn't corset giving me output? aprice67 Bioinformatics 0 09-11-2015 12:00 PM
Why is this blast script giving me empty tabular output? hlyates Bioinformatics 10 06-16-2015 11:19 PM
coverageBed -split giving incorrect output? benwittner Bioinformatics 0 10-29-2012 10:34 AM
MIRA output for Illumina mapping giving 100% coverage! Kasycas Bioinformatics 1 09-10-2010 03:45 AM
CAP3 output *.ace in EagleView? mucku Bioinformatics 13 10-27-2009 05:24 AM

Reply
 
Thread Tools
Old 04-18-2016, 02:54 PM   #1
Alirez
Junior Member
 
Location: Durham, UK

Join Date: Sep 2014
Posts: 6
Default CAP3 not giving consensus and singlet files as output

Hi all,

I'm a bit confused as to why CAP3 is not giving .singlet and .contig output files as specified in the manual. The only output I am getting is the report of the assembly which I redirect from standard output.

There doesn't seem to be an option for toggling these outputs on/off. I have a script which can extract the consensus sequences from the report, but the .singlets file (with reads not incorporated into an assembly) is important.

Can anyone offer some help as to why I'm unable to generate these outputs?

Thanks.


ps. The command I am using is:

cap3 input.fa > output.report
Alirez is offline   Reply With Quote
Old 04-19-2016, 05:02 AM   #2
Alirez
Junior Member
 
Location: Durham, UK

Join Date: Sep 2014
Posts: 6
Default Problem solved...

In case anyone is interested, when cap3 is called from a python script (using os.system) all the desired output files are produced. I haven't tested this yet, but it may well work for bash and perl as well.

A
Alirez is offline   Reply With Quote
Reply

Tags
assebly, bioinfomatics, cap3

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:45 AM.


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