SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
ask for help on soap xujie De novo discovery 0 12-13-2010 07:30 AM
soap semna Bioinformatics 0 12-02-2010 12:53 AM
Need help on soap denovo sundar De novo discovery 4 11-29-2010 04:03 AM
SAM to SOAP bioenvisage Bioinformatics 2 09-24-2010 12:48 AM
SOAP and quality Insane Bioinformatics 3 03-26-2010 02:20 AM

Reply
 
Thread Tools
Old 08-09-2008, 06:49 PM   #1
lemon027
Junior Member
 
Location: HongKong

Join Date: Aug 2008
Posts: 9
Default how to use soap?

I can't understand how to use it!

I don't kown
$soap-a is a database fa. or your query fa?

who can help me

I'M only a undergraduate student!
lemon027 is offline   Reply With Quote
Old 08-09-2008, 06:56 PM   #2
owomo
Junior Member
 
Location: china

Join Date: May 2008
Posts: 8
Default

-a <str> your query file,
-d <str> reference sequences file, is a database file
owomo is offline   Reply With Quote
Old 08-09-2008, 08:04 PM   #3
lemon027
Junior Member
 
Location: HongKong

Join Date: Aug 2008
Posts: 9
Default

but I fetch a sequence from my database to be a query file,then I use
soap -a query file -d mydatabaes file -o ......
no result!!!!!

so I converse the database to query file,just like
soap -a mydatabases file -d query file -o......
6 sequences in the result!

Why?
I can't understand!
lemon027 is offline   Reply With Quote
Old 08-11-2008, 12:08 AM   #4
zee
NGS specialist
 
Location: Malaysia

Join Date: Apr 2008
Posts: 249
Default

Your query sequence should be one or more short reads. The database should contain one or more very large sequences e.g. contigs, chromosomes, etc.

So :

soap -a <short-reads> -d <contigs/chromosome>

For -a sequences in FASTA or FASTQ formats are accepted. Database should be in FASTA only.
zee is offline   Reply With Quote
Old 08-29-2008, 06:06 AM   #5
amstisla
Junior Member
 
Location: Berlin

Join Date: Aug 2008
Posts: 4
Default

Hello.
Commonly we are using ELAND for Illumina data alignment. But for the some projects it's too short frame length available by ELAND because of its restriction (32bp as max). I need to align 40bp fragments and tried to use SOAP. The problem is that I couldn't force SOAP to taking in account quality's data (*prb files). It's working only with default '40'. I didn't find any information in official SOAP documentation or here, except notification about fastq. So the questions:
1. Is it possible to directly attach set of *.prb files to SOAP alignment process?
2. If not, how to convert *.prb and *.seq files to fastq? Are some tools available?

Thanks in advance!
Slava, MPIMG Berlin.

Last edited by amstisla; 08-29-2008 at 06:16 AM.
amstisla is offline   Reply With Quote
Old 08-29-2008, 09:49 AM   #6
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
2. If not, how to convert *.prb and *.seq files to fastq? Are some tools available?
Yup, that's what you need to do. Maq may come with a script that'll do the job. Otherwise, it's pretty straightforward to write one yourself.

I don't think that the alignment part actually pays attention to the quality scores, though. I think they only come into play in the columns of the output where the SNPs are listed, the alternate letter is listed with its quality score.

So if you aren't looking for SNPs, it might not matter much.

Last edited by swbarnes2; 08-29-2008 at 09:52 AM.
swbarnes2 is offline   Reply With Quote
Old 08-31-2008, 10:30 AM   #7
spirit
Member
 
Location: Canada

Join Date: Feb 2008
Posts: 11
Default

Maybe you can try ZOOM. It handles *_seq.txt and *_prb.txt automatically.


Quote:
Originally Posted by amstisla View Post
Hello.
Commonly we are using ELAND for Illumina data alignment. But for the some projects it's too short frame length available by ELAND because of its restriction (32bp as max). I need to align 40bp fragments and tried to use SOAP. The problem is that I couldn't force SOAP to taking in account quality's data (*prb files). It's working only with default '40'. I didn't find any information in official SOAP documentation or here, except notification about fastq. So the questions:
1. Is it possible to directly attach set of *.prb files to SOAP alignment process?
2. If not, how to convert *.prb and *.seq files to fastq? Are some tools available?

Thanks in advance!
Slava, MPIMG Berlin.
spirit is offline   Reply With Quote
Old 10-29-2008, 07:35 AM   #8
nservant
Junior Member
 
Location: france

Join Date: Oct 2008
Posts: 3
Default

Does anybody ever use SOAP in color space ?

- First Alignment : a bank against itself
/opt/soap_1.11/soap -a test.fa -d test.reference.fa -s 8 -o soap.out
-> OK. All sequences matched

- First Alignment : a bank against itself
/opt/soap_1.11/soap -a test.csfasta -d test.reference.csfasta -s 8 -o soap.out
->0 alignments !

Any ideas ?
thanks.
nservant is offline   Reply With Quote
Old 10-29-2008, 08:27 AM   #9
amstisla
Junior Member
 
Location: Berlin

Join Date: Aug 2008
Posts: 4
Default

Quote:
Originally Posted by nservant View Post
Does anybody ever use SOAP in color space ?
How did you convert your data to color space? Usually I use solid2fastq.pl script provided with MAQ to generate pseudo base reads. Beside you need special reference. It's mean less to try align on both strands of typical reference by the way of SOAP (because of color space!). Just obtain reference sequence in reverse order and after concatenate "+" and "-" to one fasta file. When starting SOAP point -n 1 option to preserve attempts of program to look matches on "complimentary" strand.
Cheers, Slava
amstisla is offline   Reply With Quote
Old 10-29-2008, 08:47 AM   #10
nservant
Junior Member
 
Location: france

Join Date: Oct 2008
Posts: 3
Default

I made a mistake ! I don't remember why, but i was convinced SOAP could work in color space.
And it should not be able to do it.

Just for answer your question Slava, I use the convert_to_dibase_encoding.pl program from ABI to convert my sequences from base to color space.
nservant 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 05:03 AM.


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