Mercurial > repos > cpt > cpt_wig_bigwig

annotate cpt_wig_bigwig/fasta_len.py @ 0:4939ae73e79b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author cpt
date Fri, 13 May 2022 05:42:27 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
4939ae73e79b Uploaded
cpt
parents:
diff changeset
1 #!/usr/bin/env python
4939ae73e79b Uploaded
cpt
parents:
diff changeset
2 import argparse
4939ae73e79b Uploaded
cpt
parents:
diff changeset
3 from Bio import SeqIO
4939ae73e79b Uploaded
cpt
parents:
diff changeset
4 import sys
4939ae73e79b Uploaded
cpt
parents:
diff changeset
5
4939ae73e79b Uploaded
cpt
parents:
diff changeset
6
4939ae73e79b Uploaded
cpt
parents:
diff changeset
7 if __name__ == "__main__":
4939ae73e79b Uploaded
cpt
parents:
diff changeset
8 parser = argparse.ArgumentParser()
4939ae73e79b Uploaded
cpt
parents:
diff changeset
9 parser.add_argument("fasta", type=argparse.FileType("r"))
4939ae73e79b Uploaded
cpt
parents:
diff changeset
10 args = parser.parse_args()
4939ae73e79b Uploaded
cpt
parents:
diff changeset
11
4939ae73e79b Uploaded
cpt
parents:
diff changeset
12 for record in SeqIO.parse(args.fasta, "fasta"):
4939ae73e79b Uploaded
cpt
parents:
diff changeset
13 sys.stdout.write("%s\t%s\n" % (record.id, len(record)))