Questions tagged [picard]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1
vote
1answer
50 views

How to demultiplex a mix of single-indexed and dual-indexed samples

The problem If I have a sample sheet that contains both single-indexed and dual-indexed samples, I can split it up into two sample sheets and then run bcl2fastq on each one. However, when doing this, ...
3
votes
2answers
126 views

Is there a safe catch-all adapter sequence for trimming?

I would like to trim/mark adapters using trimmomatic or picard MarkIlluminaAdapters from a series of Illumina Paired-End read fastqs. The fastq files may have been done using different kits or ...
1
vote
1answer
28 views

Including Picard tool in galaxy

I am running the Picard tool CollectSequencingArtifactMetrics in galaxy (locally in my pc). The output files are saved into disk and not passed to stdout and (I guess) they can't be passed as output ...
0
votes
0answers
84 views

Permanent error with picard

I am sure I set the path right but whatever I am trying the command not working Any help please? ...
1
vote
0answers
54 views

CollectHsMetrics base coverage output has overlapping targets

I am trying to get base coverage information, and am using the --PER_BASE_COVERAGE output from Picard tools to get a text file of base level coverage. However, some ...
0
votes
0answers
63 views

picard CollectSequencingArtifactMetrics

I am trying to run the Picard tool CollectSequencingArtifactMetrics without success. I am working with DATA and REFERENCE GENOME. I aligned the data to the reference genome with bowtie2. Next, I ...
3
votes
1answer
30 views

Determining Read Groups

Which Read Groups are correct: ...
5
votes
1answer
148 views

How to validate that BAMs have been downloaded correctly?

I currently have several hundred BAM files which were downloaded by someone else. These have remained untouched---before working with them, I would like to double-check that these BAMs have been fully ...
4
votes
2answers
306 views

GATK CombineVariants complains the contig order in the VCF files

I have called variants on two strains of C. elegans separately. I now want to merge the VCF files into one using the following code: Create a sequence dictionary of the reference sequence Sort the ...
1
vote
1answer
191 views

How can I subset a reference based on only the first chromosome?

I have the GATK FASTA reference for hg19. I would like to subset this based on the first chromosome. Normally what I do is open the FASTA as a text file, and then delete all other chromosomes except ...
2
votes
1answer
68 views

Picard CollectGcBiasMetrics ignoring certain chromosomes/sequences

What's the easiest way to run Picard GCBias ignoring certain chromosomes/sequences in the reference? Looking at the CollectGcBiasMetrics, it seems there isn't a bed file option that can be passed. <...
5
votes
1answer
593 views

How to remove all BAM read groups from all reads (not just the header)?

I have problem with one my BAMs---it appears to have invalid read groups. Normally when I have such a problem, I remove all the read groups from the BAM header as follows: ...
6
votes
1answer
311 views

low-memory high-speed replacement for Picard MarkDuplicates

I am running Picard MarkDuplicates with the following parameters below. On the file described, it takes about 41.6Gb of RAM memory and about 20-25 minutes to compute (only uses 1 core AFAICS). ...