Pankaj Khurana <[email protected]> writes:
> I have a few 1000 fasta files. I would like to get the list showing the
> sequence name and their respective lengths.
If you don't mind using Haskell:
import Bio.Sequence
main = do
ss <- readFasta "file.fasta"
putStr (unlines [toStr (seqlabel s) ++ "\t" ++ show (seqlen s) | s <- ss])
You can probably do something similar using Bio{Perl,Python,Ruby,..}.
-k
--
If I haven't seen further, it is by standing in the footprints of giants
_______________________________________________
BBB mailing list
[email protected]
http://www.bioinformatics.org/mailman/listinfo/bbb
