Skip to content

Commit

Permalink
update clean_fasta_cdna_cds #20
Browse files Browse the repository at this point in the history
  • Loading branch information
sinamajidian authored Apr 24, 2023
1 parent 57b10ea commit e10f629
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions archive/scripts/clean_fasta_cdna_cds.py
Original file line number Diff line number Diff line change
Expand Up @@ -67,7 +67,7 @@ def edit_record_write_faa(species_name_all_faa, faa_all, fiveLetter_species_dic,
# >lcl|AF092942.1_prot_AAC96311.1_11
prot_id_old_split= prot_id_old.split("_")
prot_id_old_split.remove("prot")
prot_id_edit = "_".join(prot_id_old_split)
prot_id_edit = ".".join(prot_id_old_split)

prot_id_new = fiveLetter_species+ prot_id_edit
prot.id = prot_id_new
Expand Down Expand Up @@ -97,7 +97,7 @@ def edit_record_write_fna(species_name_all_fna, fna_all, output_file_fna):
# >lcl|AF092942.1_prot_AAC96311.1_11
prot_id_old_split= prot_id_old.split("_")
prot_id_old_split.remove("cds")
prot_id_edit = "_".join(prot_id_old_split)
prot_id_edit = ".".join(prot_id_old_split)

prot_id_new = fiveLetter_species+ prot_id_edit
prot.id = prot_id_new
Expand All @@ -116,8 +116,8 @@ def edit_record_write_fna(species_name_all_fna, fna_all, output_file_fna):

if __name__ == '__main__':

input_folder_faa = sys.argv[1] # "data/"
output_folder_faa = sys.argv[2] # "DB/"
input_folder_faa = sys.argv[1]+"/" # "data/"
output_folder_faa = sys.argv[2]+"/" # "DB/"

output_file_fna = sys.argv[3] # "all_cdna.fa"

Expand Down

0 comments on commit e10f629

Please sign in to comment.