remove last 76 letters from string [sorry, errata]

Jan Kaliszewski zuo at chopin.edu.pl
Thu Aug 6 02:24:22 CEST 2009


Dnia 06-08-2009 o 01:54:46 PeroMHC <macmanes at gmail.com> wrote:

> This snippet represents 3 individual DNA sequences. Each sequences is
> identified by the line starting with >
> The complete file has about 10 million individual sequences.
>
> A simple enough problem, I want to read in this data, and cut out the
> last 76 letters (nucleotides) from each individual sequence and send
> them to a new txt file with a similar format.

If I understand correctly you want sth like this:


with open(path_to_the_input_file) as fasta:
       with open(path_to_the_output_file) as nucleotides:
	for seq in fasta:
               print >>nucleotides, '> foo bar length=76'
               print >>nucleotides, seq[-76]


Cheers,
*j

-- 
Jan Kaliszewski (zuo) <zuo at chopin.edu.pl>



More information about the Python-list mailing list