Open antoineszatkownik opened 2 years ago
:exclamation: Make to check out our User Guide.
consensus sequence should look like this : LHCCGVTNYrDWxDTNYYSEKGFPKSCCKLEdxxCxPQRDxxxxxxxADKVNNE
consensus sequence made by hhmake look like this : LxCCGvxnyxdWxxTxYyxexGxPkSCCxxxxxxxxxdxxkVxxe
the hhmake command line is the following : hhmake ./msa_s_exon_2_0.a2m -o ./msa_s_exon_2_0.hhm -v 0 -name ENSG00000000003#2_0 -id 100.0 -M a2m -qsc 0
I've tried this hhmake command with several parameters values and in .fasta format (with M=50 no qsc) with .a2m and no qsc or different qsc values.
In .fasta format it just choose the first sequence in the query MSA as consensus
Include as many relevant details about the environment you experienced the issue in.
:exclamation: Make to check out our User Guide.
Expected Behavior
consensus sequence should look like this : LHCCGVTNYrDWxDTNYYSEKGFPKSCCKLEdxxCxPQRDxxxxxxxADKVNNE
Current Behavior
consensus sequence made by hhmake look like this : LxCCGvxnyxdWxxTxYyxexGxPkSCCxxxxxxxxxdxxkVxxe
Steps to Reproduce (for bugs)
the hhmake command line is the following : hhmake ./msa_s_exon_2_0.a2m -o ./msa_s_exon_2_0.hhm -v 0 -name ENSG00000000003#2_0 -id 100.0 -M a2m -qsc 0
Context
I've tried this hhmake command with several parameters values and in .fasta format (with M=50 no qsc) with .a2m and no qsc or different qsc values.
In .fasta format it just choose the first sequence in the query MSA as consensus
Your Environment
Include as many relevant details about the environment you experienced the issue in.