More posts like this
-
Does read group id need to be unique between multiple bam files?
written by antmantras 8Hi. I want to do a variant calling analysis with some WGS data. I did not provide read groups in the alignment step so I …
-
Add SM group STAR
written by Lucy 17Hi, I have samples sequenced over three lanes that I am mapping together (paired end reads, three forward and three reverse) using STAR as follows: …
-
Extract VCF of tumor variants from a list of SNPs
written by vctrm67 9I have a list of germline SNPs. I want to get the read level information, etc. of these SNPs in the tumor BAM file and …
-
How to split multi run bam files and add read groups
written by MAPK 214I have a few hundred WGS samples that came without read groups. I need to run GATK on these samples, but before running GATK, I …
-
Set bam read group from existing tag
written by jbalberge 17In bam file from single-cell RNA-seq experiments processed with 10X's cellranger 3.0 pipeline, the @RG:SM tag is set by sample while the @CB is computed …
-
GATK BaseRecalibrator and platform?
written by vctrm67 9Posting here because the GATK forums are slow to respond. To my understanding, BaseRecalibrator requires certain fields, of which includes the platform information. I received …
-
Using GATK AddOrReplaceReadGroups
written by vctrm67 9Quick question regarding using AddOrReplaceReadGroups: if I have a BAM file that correspond to a single sample, but is comprised of multiple lanes, is it …
-
How treat a Bam file with mutiple read groups (RG)
written by zamani_2012 0Hi everyone, I am pretty new to the NGS data analysis. I have downloaded a WES dataset at BAM file format from SRA database. I …
-
a VERY BASIC question about "add or replace groups" for BAM files
written by CrazyB 28Q: What should I put in the read groups in my BAM files? Yes, I've read about the "read groups" thing on biostars (e.g. Picard …
-
Picard Addreplacereadgroup Error: Rg Id Not Found In Header
written by William 540<p>I have a bam file which is missing read groups in the header. </p> <p>Picard ValidateSamFile.jar gives : ERROR: Record 1, Read name M0VN6:402:2486, RG …
I would also like to know how to do this!