Is anyone aware of a conversion process to take gsMapper output (newbler) and turn it into the more standard SAM/BAM format?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
The development version of bank2contig can take the -S option to output an AMOS bank as SAM. To get from gsMapper to an AMOS bank you'll have to use the outputted .ace file that gsMapper gives you. I haven't gotten it to to work properly though as the resulting SAM file doesn't seem to contain the base qualities.
I'm seriously considering scripting it myself if Roche doesn't deliver soon
Comment
-
I've scripted a GS Mapper to SAM format converter in Ruby.
You can find it here:
At the moment it only handles single end reads.
The code might be a bit ugly as I used this as an educational project to get into Ruby. Instructions can be found in the readme.
Caveat: The sam file compiles to bam, but I haven't thoroughly checked it for errors, so there might still be some bugs. Please feel free to report any if you see them
Comment
-
Asifrim,
The bank2contig wtih the -S option doesn't give a SAM format. According to the documentation online, conversion to SAM should be using the -s option but as of version 2.0.8 that is not available. The -S option in version 2.0.8 give a Simple Layout Style.
I also tried your Ruby script GstoSam but it truncated prematurely.
MQ.
Comment
-
Hi MQ,
I seem to have mistyped it, it should have have been downcase: -s for SAM
(http://sourceforge.net/apps/mediawik...le=Bank2contig)
Did the ruby script return any kind of error? What version of Ruby are you using? Are you outputting the standard GS mapper output or the tabbed version of the output? If I remember correctly the script works on the default output type (which is the multiline output per read).
Comment
-
Originally posted by asifrim View PostI've scripted a GS Mapper to SAM format converter in Ruby.
You can find it here:
At the moment it only handles single end reads.
The code might be a bit ugly as I used this as an educational project to get into Ruby. Instructions can be found in the readme.
Caveat: The sam file compiles to bam, but I haven't thoroughly checked it for errors, so there might still be some bugs. Please feel free to report any if you see them
So I am using your script on gsMapper output. It seems to work to a first approximation. (Which is more than I can say for any of the other methods I have tried...) But it does not seem like I am quite there. Your script creates a .sam file -- but there is no header created. The sam specification claims http://samtools.sourceforge.net/SAM1.pdf
The SAM format consists of one header section and one alignment section.The whole header section can be absent, but keeping the header is recommended.
Code:samtools view -b Synechocystis-1.sam [bam_header_read] EOF marker is absent. [main_samview] fail to read the header.
Version: 0.1.7 (r510)
Is there a reason why for my purposes a header section is required? Perhaps because I was running gsMapper, rather than gsAssembler?
--
Phillip
PS Yes, I agree Roche programs should offer SAM/BAM output as well as .ace
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 08:47 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Today, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
57 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Comment