@@ -123,13 +123,15 @@ def write_fasta_from_json(self):
123123 % (j [i ]['model_id' ], j [i ]['model_name' ]))
124124 logger .info (
125125 "Please let the CARD Admins know! Email: card@mcmaster.ca" )
126+ snpList = [] # reset
126127
127128 try :
128129 for seq in j [i ]['model_sequences' ]['sequence' ]:
129- fout .write ('>%s_%s | model_type_id: 40293 | pass_bit_score: %s | SNP: %s | %s\n '
130- % (i , seq , pass_bit_score , ',' .join (snpList ), j [i ]['ARO_name' ]))
131- fout .write ('%s\n ' % (
132- j [i ]['model_sequences' ]['sequence' ][seq ]['protein_sequence' ]['sequence' ]))
130+ if len (snpList ) > 0 :
131+ fout .write ('>%s_%s | model_type_id: 40293 | pass_bit_score: %s | SNP: %s | %s\n '
132+ % (i , seq , pass_bit_score , ',' .join (snpList ), j [i ]['ARO_name' ]))
133+ fout .write ('%s\n ' % (
134+ j [i ]['model_sequences' ]['sequence' ][seq ]['protein_sequence' ]['sequence' ]))
133135 except Exception as e :
134136 logger .warning ("No model sequences for model (%s, %s). RGI will omit this model and keep running."
135137 % (j [i ]['model_id' ], j [i ]['model_name' ]))
@@ -154,13 +156,15 @@ def write_fasta_from_json(self):
154156 % (j [i ]['model_id' ], j [i ]['model_name' ]))
155157 logger .info (
156158 "Please let the CARD Admins know! Email: card@mcmaster.ca" )
159+ snpList = [] # reset
157160
158161 try :
159162 for seq in j [i ]['model_sequences' ]['sequence' ]:
160- fout .write ('>%s_%s | model_type_id: 41091 | pass_bit_score: %s | SNP: %s | %s\n '
161- % (i , seq , pass_bit_score , ',' .join (snpList ), j [i ]['ARO_name' ]))
162- fout .write ('%s\n ' % (
163- j [i ]['model_sequences' ]['sequence' ][seq ]['protein_sequence' ]['sequence' ]))
163+ if len (snpList ) > 0 :
164+ fout .write ('>%s_%s | model_type_id: 41091 | pass_bit_score: %s | SNP: %s | %s\n '
165+ % (i , seq , pass_bit_score , ',' .join (snpList ), j [i ]['ARO_name' ]))
166+ fout .write ('%s\n ' % (
167+ j [i ]['model_sequences' ]['sequence' ][seq ]['protein_sequence' ]['sequence' ]))
164168 except Exception as e :
165169 logger .warning ("No model sequences for model (%s, %s). RGI will omit this model and keep running."
166170 % (j [i ]['model_id' ], j [i ]['model_name' ]))
0 commit comments