Planning to Write a script that reads a Fasta-formatted file and writes the number of sequences in it, followed by the accessions of the sequences. For example, if the input is
>Pig
ACGT
>Dog
AGTG
>Bunny
TATA
then the output should be
3
Pig
Dog
Bunny
any idea?
There is a full course on biopython from the Pasteur institute which includes reading Fasta format http://www.pasteur.fr/recherche/unites/sis/formation/python/ch11s03.html#exa_seq_string
I think you should read this course and use the Bio modules. Don't reinvent the wheel at first !
We're a friendly, industry-focused community of developers, IT pros, digital marketers, and technology enthusiasts meeting, networking, learning, and sharing knowledge.