SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > 454 Pyrosequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
using sff_extract veronique Bioinformatics 24 04-28-2013 11:01 PM
sff_extract problems clostridium40 Bioinformatics 2 07-21-2011 05:08 AM
sff_extract question veena Bioinformatics 1 06-20-2010 10:31 PM
sff_extract Syntax Error Cresten Bioinformatics 3 12-21-2009 12:16 AM
sff_extract paired end extraction error coldturkey Bioinformatics 1 03-12-2009 04:48 AM

Reply
 
Thread Tools
Old 02-23-2010, 05:03 AM   #1
Mona
Member
 
Location: Uppsala

Join Date: Feb 2010
Posts: 27
Default Error using sff_extract.

Hi,
i m trying to extract fasta, qual and xml files from .sff file using sff_extract.

I have just all the provided commands for extraction, Following commands have been used.

sff_extract filename.sff

sff_extract -o filename.sff

$ sff_extract -s seq.fasta -q qual.fasta -x anci.xml my_454_file.sff

With all the above commands, i get the same error message: "struct.error: unpack requires a string argument of length 1600"
How i can remove this error?
Mona is offline   Reply With Quote
Old 02-23-2010, 11:38 PM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

This suggests either the SFF file is damaged (e.g. truncated due to a problem downloading - try redownloading it), or doesn't quite follow the standard sff_extract is expecting. Is this an original unmodified SFF file from your sequencing center?

Do you have the (Linux only) Roche tools to test with?
maubp is offline   Reply With Quote
Old 02-23-2010, 11:53 PM   #3
Mona
Member
 
Location: Uppsala

Join Date: Feb 2010
Posts: 27
Default

yes this is the original file from the sequencing centre.
and i dnt have the roche tools to test with.
Mona is offline   Reply With Quote
Old 02-24-2010, 12:03 AM   #4
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Getting the Roche "Newbler" tools is quite easy, assuming you have a suitable Linux machine to run them on. See:
http://seqanswers.com/forums/showthread.php?t=114

However, in the short term as mentioned above, I would strongly suggest you try re-downloading the SFF file from the sequencing center, or at least confirm you got the whole file by checking the filesize.

P.S. Seeing the whole error message might help.
maubp is offline   Reply With Quote
Old 02-24-2010, 12:23 AM   #5
Mona
Member
 
Location: Uppsala

Join Date: Feb 2010
Posts: 27
Default

Is it possible to get error in sff_extract when there is no vector used in the 454 sequencing.
As i know that the centre from which my data is coming they are not using any vector for sequencing. But A and B adapters are used. so is it still compulsory to screen the sequences and extract the fasta, qual and xml files from the sff file. and second thing do i have to perform the sequence masking for the tag sequence after the extraction of these three files? or only one of the above two steps is required.
Mona is offline   Reply With Quote
Old 03-18-2010, 05:13 AM   #6
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

Do you have an example SFF file with no vector used in the 454 sequencing that breaks sff_extract that you could share? I'd like to use it to check the SFF support in Biopython.

(Regarding the rest of your question, it probably depends on what you are going to do next - if you are using MIRA then it might be easier to ask on the mira_talk mailing list)
maubp is offline   Reply With Quote
Old 03-23-2010, 04:05 AM   #7
Mona
Member
 
Location: Uppsala

Join Date: Feb 2010
Posts: 27
Default

Yes i have the file that breaks sff_extract.
Can u send me ur email id?
Mona is offline   Reply With Quote
Old 03-23-2010, 06:37 AM   #8
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

After chatting to Mona privately, it does sound like the original problem was a partial SFF file. The full file works fine with sff_extract
maubp is offline   Reply With Quote
Old 03-23-2010, 06:59 AM   #9
Mona
Member
 
Location: Uppsala

Join Date: Feb 2010
Posts: 27
Default

Yes exactly.
Thnx maubp!!
Mona 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 01:22 AM.


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