More posts like this
-
Adding read groups to mark duplicated bam file
written by rezahi I have a mark duplicated bam file and I want to use it to variant calling using GATK. GATK needs read groups and I …
-
For loop to parallel run commands
written by endretoth •After reading two articles/post (I put the links below) I decided to parallelize my commands to be more effective. I would like use a For …
-
Using GATK AddOrReplaceReadGroups
written by vctrm67 •Quick question regarding using AddOrReplaceReadGroups: if I have a BAM file that correspond to a single sample, but is comprised of multiple lanes, is it …
-
Mark duplicate "Picard"
written by Nelislam4 •Hi, I am wondering, How mark duplicate of Picard tools works? In other words. How Picard tools detect the duplicated read in BAM file? coz …
-
Extract consensus sequence reads (collapse PCR duplicates) from bam
written by Nicolas RosewickHi, I did an experiment where I designed a sequencing library focused on a specific genomic region of interest. To go very deep and to …
-
How do I find input arguments (RGPU, RGLB, etc) for AddOrReplaceReadGroups
written by MaxI created a .bam file without read groups in the header, and am trying to correct this by using picard tools' AddOrReplaceReadGroups. I've looked over …
-
Piping samtools output to picard
written by Yaseen Ladak •Hi All, I need some help I want to pipe the output from samtools to picard read group commands but I am failing: bowtie2 -x …
-
Adding Multiple Read Group Headers to a single BAM file
written by koen.vdberge •I have a BAM file which consists of sequencing reads derived from pooling 25 individuals and sequencing the pool of DNA. Variant calling, or calculating …
-
a VERY BASIC question about "add or replace groups" for BAM files
written by CrazyBQ: 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 …
-
GATK, SAM file doesn't have any read groups defined in the header
written by Chirag NepalHi all, I have been trying to use Mutect to compare results from Varscan and other tools. To run MuTect, pre-processing from GATK and Picard …
Replace the read group information or name or swap it for a completely different dataset? The former two can be done simply by modifying the BAM header (you can reheader the file with
samtools reheaderthen).