Mercurial > repos > jackcurragh > ribogalaxy_get_chrom_sizes

comparison get_chrom_sizes/calculating_chrom.sizes.py @ 17:f87a4b64e6ef draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author jackcurragh
date Tue, 17 May 2022 13:57:39 +0000
parents ab7b5928faea
children 06c2eaff31e5
comparison
equal deleted inserted replaced
16:ab7b5928faea 17:f87a4b64e6ef
2 # adapted from https://bioexpressblog.wordpress.com/2014/04/15/calculate-length-of-all-sequences-in-an-multi-fasta-file/ 2 # adapted from https://bioexpressblog.wordpress.com/2014/04/15/calculate-length-of-all-sequences-in-an-multi-fasta-file/
3 from sys import argv 3 from sys import argv
4 # python calculating_chrom.sizes.py genome_input.fa output.chrom.sizes 4 # python calculating_chrom.sizes.py genome_input.fa output.chrom.sizes
5 fasta_source = str(argv[1]) 5 fasta_source = str(argv[1])
6 prefix = str(argv[2]) 6 prefix = str(argv[2])
7 genome = str(argv[3]) 7 output = str(argv[3])
8 builtin = str(argv[4]) 8 genome = str(argv[4])
9 builtin = str(argv[5])
9 10
10 # genome = 'test-data/test.fasta' 11 # genome = 'test-data/test.fasta'
11 # output = "test-data/test_chrom.sizes" 12 # output = "test-data/test_chrom.sizes"
12 print(fasta_source, genome, builtin, prefix, output)
13 if fasta_source == 'builtin': 13 if fasta_source == 'builtin':
14 genome = builtin 14 genome = builtin
15 15
16 chromSizesoutput = open(output,"w") 16 chromSizesoutput = open(output,"w")
17 17