Les germlines J+down ne sont pas gappées
Est-ce volontaire ou un bug qui pourrait avoir des conséquences sur les calculs de bio-productivity ?
En tout cas on pourrait peut-être le changer en inversant les deux blocs:
for key in key_downstream:
downstream_data[key][-1][1]['seq'] += l
(...)
if '>' not in l and current_files and feature == FEATURE_J_REGION:
l = gap_j(l)