comment on

Well I have unfortunately encountered another error. I read the qual file in like a fastA file and put it into a hash and I can't write it like a qual file for some reason.. is there some sort of parsing that I have to do that the BioPerl modules won't do for me?

--------------------- WARNING ---------------------
MSG: seq doesn't validate, mismatch is 4040404040404040404040404040404
+040404040404040404040404040404040404040393737373737373737373737373737
+373737373535353535353030302828283030323235373737373737373737373737373
+737373737373737373737373737373535353737343232323232323230161515152427
+2626131313232313131317212121212119211618171918
---------------------------------------------------

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: Attempting to set the sequence to [404040404040404040404040404040
+404040404040404040404040404040404040404039373737373737373737373737373
+737373737353535353535303030282828303032323537373737373737373737373737
+373737373737373737373737373737353535373734323232323232323016151515242
+72626131313232313131317212121212119211618171918] which does not look 
+healthy
STACK: Error::throw
STACK: Bio::Root::Root::throw /usr/lib/perl5/site_perl/5.8.8/Bio/Root/
+Root.pm:357
STACK: Bio::PrimarySeq::seq /usr/lib/perl5/site_perl/5.8.8/Bio/Primary
+Seq.pm:270
STACK: Bio::PrimarySeq::new /usr/lib/perl5/site_perl/5.8.8/Bio/Primary
+Seq.pm:221
STACK: Bio::Seq::new /usr/lib/perl5/site_perl/5.8.8/Bio/Seq.pm:484
STACK: tagfinder.pl:227
[download]

In reply to Re^2: Reading in Qual File Using Bio::SeqIO by twaddlac
in thread Reading in Qual File Using Bio::SeqIO by twaddlac

Are you posting in the right place? Check out Where do I post X? to know for sure.
Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
<code> <a> <b> <big> <blockquote> <br /> <dd> <dl> <dt> <em> <font> <h1> <h2> <h3> <h4> <h5> <h6> <hr /> <i> <li> <nbsp> <ol> <p> <small> <strike> <strong> <sub> <sup> <table> <td> <th> <tr> <tt> <u> <ul>
Snippets of code should be wrapped in <code> tags not <pre> tags. In fact, <pre> tags should generally be avoided. If they must be used, extreme care should be taken to ensure that their contents do not have long lines (<70 chars), in order to prevent horizontal scrolling (and possible janitor intervention).
Want more info? How to link or How to display code and escape characters are good places to start.


Keep It Simple, Stupid
	PerlMonks