Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

picard pacbio

Hello,everyone,
I have done the Pacio reads alignment with BWA-SW,and generate a .sam file,how can I tranfer the .sam to .bam with header?
I used samtools and picard ,but I could not gain a .bam with header.What the Addorreplacereadgroups parameter RGID,RGLB,RGPL,RGPU,RGSM should be with Pacbio data??

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    With Picard tools, you just need to specify the output file name with ".bam" as suffix, and the tool will automatically write a bam file instead of a Sam file.

    Regarding the read group information, please look up the sam/bam file specification to learn what these fields describe.

Sign In or Register to comment.