http://qs1969.pair.com?node_id=1076180


in reply to Re: A program to extract the reads and modify the seq ID
in thread A program to extract the reads and modify the seq ID

Your script only prints the sequences and not the header with the OP's wanted appended info.

  • Comment on Re^2: A program to extract the reads and modify the seq ID

Replies are listed 'Best First'.
Re^3: A program to extract the reads and modify the seq ID
by robby_dobby (Hermit) on Feb 26, 2014 at 04:02 UTC

    Yes, didn't I say that I fixed the program "mostly for clarity/readability"? I also mentioned this in my last point. :-)

      The OP says that the program ...neither prints out the output nor gives me no errors.... In fact, it doesn't print anything. I'm not sure that you've clarified it or made it more readable--nor did the OP request such--but your refactoring now makes it print only the sequences.

      Please don't get me wrong: your bullet points for the OP were very well done. I suppose I just assumed too much thinking that the code you posted was a corrected version of the OP's code. Instead, the OP must now fix your code--or his/her own--to get it to work correctly. The latter, however, is certainly not a bad thing, and your bullet points should help...

        Ah, I should have read the description before diving into code. :-)

        But yes, my approach when giving code is not to give everything, but get OP to fix (some remaining) issues - I've found that learning is better reinforced that way.