在分析转录组数据时,用bowtie2比对生成的bam文件,下一步call peak使用m6Aviewer,需要bam文件的index文件。所以我直接敲命令

samtools index xx.bam 
[E::hts_idx_push] NO_COOR reads not in a single block at the end 10 -1

报上面的错误。后来发现,samtools建索引时,bam或sam文件必须是排序好了的,而且必须使用samtools的默认排序方法,即坐标顺序排序。

time nohup samtools sort -@ 20 28.m6A.bam -o 28.m6A.sort.bam &

然后在进行index

相关文章:

  • 2021-07-03
  • 2022-12-23
  • 2021-05-04
  • 2021-08-25
  • 2021-11-05
  • 2021-07-15
  • 2022-12-23
  • 2021-05-26
猜你喜欢
  • 2019-09-04
  • 2022-12-23
  • 2021-05-26
  • 2021-09-08
  • 2022-12-23
  • 2022-12-23
相关资源
相似解决方案