Skip to content
Snippets Groups Projects
Commit 50a21532 authored by Samuel Mondal's avatar Samuel Mondal
Browse files

reducing the scope of cli.py to work only after bedtools

parent 22146b82
No related branches found
No related tags found
1 merge request!57reducing the scope of cli.py to work only after bedtools
...@@ -9,20 +9,18 @@ parser = argparse.ArgumentParser( ...@@ -9,20 +9,18 @@ parser = argparse.ArgumentParser(
prog = 'transcript_sequence_extractor', prog = 'transcript_sequence_extractor',
description = 'extracts transcript sequences from genome sequence and ouputs transcripts with PolyA tail added to them') description = 'extracts transcript sequences from genome sequence and ouputs transcripts with PolyA tail added to them')
parser.add_argument('--input_fasta_file', parser.add_argument('--input_fasta_file',
help='genome fasta file') help='fasta file obtained from bedtools')
parser.add_argument('--input_gtf',
help='gtf file')
parser.add_argument('--output_file_name', parser.add_argument('--output_file_name',
help='output fasta file') help='Name of the output fasta file')
args = parser.parse_args() args = parser.parse_args()
def main(): def main():
LOG.info("sequence_extractor begins") LOG.info("sequence_extractor begins")
exon_extraction_from_gtf() fasta_list = exon_concatenation(args.input_fasta_file)
fasta_list = exon_concatenation(args.)
final_list = poly_a_addition_to_fasta_list(fasta_list) final_list = poly_a_addition_to_fasta_list(fasta_list)
list_to_file(final_list,args.output_file_name) with open(args.output_file_name, 'w') as fasta_out:
fasta_out.write('\n'.join('%s\n%s' % x for x in final_list))
LOG.info("sequence_extractor ends") LOG.info("sequence_extractor ends")
if ___name__ == 'main': if ___name__ == 'main':
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment