hh-suite icon indicating copy to clipboard operation
hh-suite copied to clipboard

reformat.pl puts #=GC at top, should be at bottom

Open dstern opened this issue 2 years ago • 0 comments

:exclamation: Make to check out our User Guide.

Expected Behavior

#=GC line is recommended at bottom of alignment (https://en.wikipedia.org/wiki/Stockholm_format##=GC) Does not seem like big deal, but causes alphafold2 to crash, because af2 expects #=GC at bottom (in their module remove_empty_columns_from_stockholm_msa - https://github.com/deepmind/alphafold/blob/0be2b30b98f0da7aecb973bde04758fae67eb913/alphafold/data/parsers.py) This is a problem when using custom multiple sequence alignments as input.

Current Behavior

in conversion to Stockholm format, #=GC RF placed at top of alignment

Steps to Reproduce (for bugs)

Please make sure to execute the reproduction steps.

reformat.pl a3m sto file.a3m file.sto

HH-suite Output (for bugs)

Please make sure to post the complete output of the tool you called. Please use gist.github.com.

Context

Providing context helps us come up with a solution and improve our documentation for the future.

Your Environment

Include as many relevant details about the environment you experienced the issue in.

  • Version/Git commit used:
  • Server specifications (especially CPU support for AVX2/SSE and amount of system memory):
  • Operating system and version:

dstern avatar Nov 19 '21 23:11 dstern