samtools 生成bcf文件报错:[bcf_sync] incorrect number of fields (0 != 5) at 0:0

起源

用samtools 生成bcf文件,命令如下

samtools mpileup -ABuf ref.fa L.bam >L.raw.bcf

错误

输出文件一直增大。但是bcf文件内容只有表头,bcf文件有一行以下信息

[bcf_sync] incorrect number of fields (0 != 5) at 0:0

提示如下:

/opt/gridengine/default/spool/node179/job_scripts/373553: line 4: syntax error: unexpected end of file

解决

尝试好久,换了一个版本的samtools,既然成功了,
之前版本 samtools-1.3.1 换成 0.1.18
原因暂时不明,怀疑如下:
1:版本
2:安装编译问题

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这个错误通常是由于在生成BAM或VCF文件的索引文件时,索引文件的排序顺序与原始文件不匹配导致的。为了解决这个问题,你可以按照以下步骤进行操作: 1. 确认你的原始文件和索引文件是否是匹配的。你可以使用以下命令确认: ``` samtools view -H your_file.bam | grep '^@SQ' | cut -f 2 | sed 's/SN://' | sort | uniq ``` 这将列出原始文件中存在的所有染色体名称。然后,使用以下命令检查索引文件中的染色体名称: ``` bcftools index -s your_file.vcf.gz ``` 确认两个列表是否完全匹配。 2. 如果你的原始文件和索引文件不匹配,你需要重新生成索引文件。你可以使用以下命令重新生成: ``` samtools index your_file.bam ``` 或者 ``` bcftools index your_file.vcf.gz ``` 然后再次确认两个文件的染色体名称是否匹配。 3. 如果你的文件已经正确排序并且你已经重新生成了索引文件,但仍然出现此错误,请检查你的文件是否存在连续的染色体块。你可以使用以下命令检查: ``` samtools view -H your_file.bam | grep '^@SQ' | cut -f 3 | sed 's/LN://' | awk '{sum+=$1} END {print sum}' ``` 这将列出原始文件中所有染色体的长度总和。然后,使用以下命令检查索引文件中的染色体块长度: ``` bcftools index -s your_file.vcf.gz ``` 如果你的文件存在间隔或断裂的染色体块,请尝试重新生成文件或使用其他工具来修复文件
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值