SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Sam 2 ace jdiezperezj Bioinformatics 2 05-24-2012 01:08 PM
454 .ace to .bam conversion issue pmiguel Bioinformatics 14 06-24-2011 12:30 AM
NEw to Chip-seq and have .bam/.sam/.bam.bai files... then what? NGS newbie Bioinformatics 11 05-25-2011 08:48 AM
How do I convert 454 ace to a regular ace? lskatz Bioinformatics 5 11-22-2010 08:31 AM
Convert ACE to SAM (and SAM to ACE)? dan Bioinformatics 0 05-04-2009 03:55 AM

Reply
 
Thread Tools
Old 05-14-2010, 12:30 PM   #1
JueFish
Member
 
Location: Connecticut

Join Date: May 2010
Posts: 42
Default ACE to SAM/BAM pile-up

Hi all,

I'm trying to generate predicted SNP sites using a new program called ACCUSA, but need to get my data into SAM/BAM pile-up format first. I'm going from a 454-data generate ACE file to this format, but first I need to convert that to SAM/BAM format (before I can use SAMTOOLS to do the pileup) and have gotten myself a bit confused about what program to use to do this conversion. Any thoughts out there about what would be a good utility to do this conversion from ACE to SAM/BAM?

Thanks
JueFish is offline   Reply With Quote
Old 06-16-2010, 04:23 AM   #2
ulz_peter
Senior Member
 
Location: Graz, Austria

Join Date: Feb 2010
Posts: 219
Default

I would be interested in this conversion too, ACE is by far too big to visualize easily.
I know there is a way of parsing the ACE file in Biopython, but the AlignIO package yet doesn't support SAM-File Output as far as I know.
Really looking forward to standardizing file formats in NGS...
ulz_peter is offline   Reply With Quote
Old 06-16-2010, 05:40 AM   #3
imilne
Member
 
Location: JHI, Dundee, UK

Join Date: Jan 2010
Posts: 68
Default

If you use the AMOS tools, I believe you can go from ACE to AFG/BANK, and from there to SAM/BAM.

Iain
__________________
Our software: Tablet | Flapjack | Strudel | CurlyWhirly | TOPALi
imilne is offline   Reply With Quote
Old 06-16-2010, 09:39 AM   #4
jkbonfield
Senior Member
 
Location: Cambridge, UK

Join Date: Jul 2008
Posts: 146
Default

Gap5 can import ACE too and can output SAM. Although the SAM flags for 1st/2nd pair are a bit broken in the current release; fixed in the SVN tree. I've no idea how it compares to AMOS for validity or efficiency.

Both of these are around-the-houses solutions though so I'm sure someone must have written a dedicated perl hack for it by now. Anyone?

James
jkbonfield 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:05 PM.


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