SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Plasma DNA Puzzle whw Sample Prep / Library Generation 9 08-03-2011 11:41 AM
BWA puzzle: Read fate is not deterministic? bhootnaath Bioinformatics 2 06-09-2011 12:50 PM
Looking process to convert gff3 format into ace format or sam format andylai Bioinformatics 1 05-17-2011 03:09 AM
puzzle in fq_all2std.pl biocc Bioinformatics 0 11-25-2010 04:33 AM
blast puzzle anyone1985 Bioinformatics 1 09-06-2009 02:45 AM

Reply
 
Thread Tools
Old 03-13-2009, 04:29 AM   #1
anyone1985
Member
 
Location: shanghai, chia

Join Date: Mar 2009
Posts: 67
Default format puzzle

Recently, I have finished Solexa assembly. I convert it to .ace by velvet's amos2ace. I want to transfer it to .phd file. But when I open it with Consed, it tell me that "BS line refers to a read that isn't in a AF line for this contig ". So I write a script like that :

#!/bin/perl -w
use Bio::SeqIO;
open TEMP, ">xp.phd" or die "can not create the file xp.phd";
close TEMP;
$in = Bio::SeqIO->new ( -file => "xp_4.ace", -format => 'ace');
$out = Bio::SeqIO->new (-file => ">xp.phd", -format => 'phd');
while( $seq = $in->next_seq())
{
print $seq->
$out->write_seq($seq);
}

But nothing happen. How should I do? or how .ace file transfer to .phd file without consed?
anyone1985 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:41 AM.


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