SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to convert BED format to SAM/BAM? seq_newbie Bioinformatics 1 06-23-2011 09:11 AM
Looking process to convert gff3 format into ace format or sam format andylai Bioinformatics 1 05-17-2011 03:09 AM
how to convert from bam or sam to map rcorbett Bioinformatics 1 11-04-2010 11:49 AM
convert BAM -> SAM? thinkRNA Bioinformatics 4 06-30-2010 04:08 AM
Convert MAQ,SOAP alignment result to ACE format baohua100 Bioinformatics 0 03-02-2009 02:01 AM

Reply
 
Thread Tools
Old 01-24-2012, 02:05 AM   #21
Aman Mahajan
Member
 
Location: India

Join Date: Jan 2012
Posts: 22
Default

Yes. the command "head" works . shows me the initial contents of the file.
Aman Mahajan is offline   Reply With Quote
Old 01-24-2012, 02:13 AM   #22
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 801
Default

Quote:
shows me the initial contents of the file
And what are these contents? Could you please show the output?

This page suggests it should be a 9-column output, tab separated:

http://soap.genomics.org.cn/soap1/#Formatofoutput

To confirm tab separated, you can do this:

Code:
hexdump -C abc.soap | head
and look for '09' between fields in the left hand side of the output.
gringer is offline   Reply With Quote
Old 01-24-2012, 02:19 AM   #23
Aman Mahajan
Member
 
Location: India

Join Date: Jan 2012
Posts: 22
Default

>3 length 34 cvg_0.0_tip_0
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACTT
>5 length 34 cvg_0.0_tip_0
AGTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT

These are the first intial lines of the file abc.soap.


On running the command you mentioned, I get an output , not sure what you mentioned to look for means.
Aman Mahajan is offline   Reply With Quote
Old 01-24-2012, 02:24 AM   #24
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 801
Default

Code:
>3 length 34 cvg_0.0_tip_0
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACTT
>5 length 34 cvg_0.0_tip_0
AGTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT
That's a FASTA file, not a SOAP alignment file. Are you sure you ran SOAP (i.e. the SOAPaligner program) rather than SOAPdenovo to get this output?
gringer is offline   Reply With Quote
Old 01-24-2012, 02:30 AM   #25
Aman Mahajan
Member
 
Location: India

Join Date: Jan 2012
Posts: 22
Default

We ran soapdenovo on this...so we have to use soapaligner ??

I mean the output isn't same for soapdenovo and soapaligner??..

And lastly if I have to convert only this fasta file into SAM, what should I do?..is it possible?
Aman Mahajan is offline   Reply With Quote
Old 01-24-2012, 02:38 AM   #26
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 801
Default

Quote:
We ran soapdenovo on this...so we have to use soapaligner ??
I mean the output isn't same for soapdenovo and soapaligner??
Assembly and alignment are two different things. soap2sam.pl expects alignments in the format produced by SOAP (i.e. SOAPaligner). If you don't want to use SOAP for this, there are plenty of other programs for mapping that can produce SAM files directly as output.

Quote:
And lastly if I have to convert only this fasta file into SAM, what should I do?..is it possible?
You may be able to simulate reads using the coverage values in the sequence headers, but that would be silly because it wouldn't be a true representation of your data. Also, you would still need to work out a way to generate SAM files.
gringer is offline   Reply With Quote
Old 01-24-2012, 02:46 AM   #27
Aman Mahajan
Member
 
Location: India

Join Date: Jan 2012
Posts: 22
Default

Sorry to bother you, I am naive in this field.

Actually I did not use SOAPaligner as it requires a reference to run I think. So I used Soapdenovo. Now I want to run this output in Rseq for quantification of my data.

But Rseq requires SAM format as input so i have been trying to convert it to SAM.

Which softwares I must use , if you can suggest?
Aman Mahajan is offline   Reply With Quote
Old 01-24-2012, 02:52 AM   #28
Aman Mahajan
Member
 
Location: India

Join Date: Jan 2012
Posts: 22
Default

And I want to user mrfQuantifier in Rseq, so basically I need the SAM format as input for it or is there anyway to convert it to MRF directly without involving a SAM conversion step.
Aman Mahajan is offline   Reply With Quote
Old 01-24-2012, 03:29 AM   #29
arvid
Senior Member
 
Location: Berlin

Join Date: Jul 2011
Posts: 156
Default

Maybe you should read up on assembly, alignment and quantification before messing around too much with the different software...

The contigs generated by SOAPdenovo can be used as a reference for quantification. You can align the reads to those contigs with e.g. SOAPaligner to get a file with aligned reads (in the SOAP format) that you can convert into SAM for further analysis. Or you can use a software like RSEM for quantification.
arvid is offline   Reply With Quote
Old 01-24-2012, 03:38 AM   #30
arvid
Senior Member
 
Location: Berlin

Join Date: Jul 2011
Posts: 156
Default

Also you should possibly try a transcriptome-centric assembler such as Trinity, Trans-ABySS, Oases or SOAPdenovo-Trans instead of SOAPdenovo, which is optimized for genomic assembly.
arvid is offline   Reply With Quote
Old 01-10-2018, 08:59 PM   #31
hewm2008
Junior Member
 
Location: china,sz

Join Date: Jan 2018
Posts: 2
Default

Quote:
Originally Posted by KevinLam View Post
Dear All,
would you know if there's a convertor for SOAP alignment format to SAM?
see this software:
https://github.com/BGI-shenzhen/soap2bam
hewm2008 is offline   Reply With Quote
Old 01-10-2018, 09:05 PM   #32
hewm2008
Junior Member
 
Location: china,sz

Join Date: Jan 2018
Posts: 2
Default

This software can help all of you
https://github.com/BGI-shenzhen/soap2bam
hewm2008 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 02:36 PM.


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