User: Peter Chung

gravatar for Peter Chung
Peter Chung80
Reputation:
80
Status:
Trusted
Location:
Hong Kong
Last seen:
3 weeks, 3 days ago
Joined:
2 years, 6 months ago
Email:
q*************@gmail.com

Posts by Peter Chung

<prev • 15 results • page 1 of 2 • next >
0
votes
0
answers
108
views
0
answers
gatk4 GenomicsDBImport intervals with whole genome
... I am new in gatk tools, I would like to use GenomicsDBImport to merge GVCFs from multiple samples with whole genome. -L interval is a required option for GenomicsDBImport. For example, -L chr20 for contig chr20. -L chr20:1-100 for contig chr20, positions 1-100. However, I would like to import whol ...
gatk intervals genomicsdbimport written 4 weeks ago by Peter Chung80
0
votes
0
answers
178
views
0
answers
Comment: C: picard markduplicate output smaller file
... oh thanks. First I used bwa to align them and then use samtools sort to sort each bam files. Afterwards, I combined all the bam files into one bam file by samtools merge. After that, I used samtools addreplacerg to add readgroup. ###bwa and samtools sort REF="/data/data/reference/refs/ ...
written 8 weeks ago by Peter Chung80 • updated 8 weeks ago by finswimmer11k
0
votes
0
answers
178
views
0
answers
picard markduplicate output smaller file
... I am new in WGS analysis. First, I combine all my bam files into one and it's 157GB and then add read group on it to 159GB. Then I do the picard markduplicate step by using the following code: java -Xmx8g -Djava.io.tmpdir=${TMPFILE} -jar $PICARD MarkDuplicates \ INPUT=${FILE}.addRG.bam \ ...
gatk picard markduplicates bioinformatics written 8 weeks ago by Peter Chung80
1
vote
1
answer
183
views
1
answer
gatk HaplotyperCaller bam reference
... I am new in WGS sequencing. I am working on the gatk HaplotypeCaller step. The bam file is required for calling variants. I have generated a few bam files upstream starting from bwa. I have sample1.bam, sample1.addRG.bam for add read group and sample1.addRG.mkdup.bam for remove duplicates. Therefo ...
haplotypercaller gatk bam written 3 months ago by Peter Chung80
1
vote
0
answers
225
views
5 follow
0
answers
picard addorreplacereadgroups output issue
... I have a 189GB bam file and I need to add a readgroup in order to do downstream steps, so I used picard AddOrReplaceReadGroups to do it. Below is my command: java -Xmx8g -jar $PICARD AddOrReplaceReadGroups \ I=${name}.bam \ O=${name}.addRG.bam \ SORT_ORDER=coordinate \ RGID=${ ...
gatk readgroups picard java sequencing written 3 months ago by Peter Chung80
0
votes
0
answers
483
views
0
answers
use gatk phased vcf as input file of HapCut2
... I am new in sequencing analysis and would like to do haplotype blocks from the result of WGS. But I have some questions on the HAPCUT2. I generated a vcf from gatk and then I phased it using gatk ReadBackedPhasing into phased vcf. Commands below are from HAPCUT2: ./build/extractHAIRS [options ...
hapcut2 gatk vcf phasing written 10 months ago by Peter Chung80 • updated 10 months ago by RamRS20k
0
votes
2
answers
432
views
2
answers
bwa mem without -M option how to continue using picard
... I do the alignment without -M option and now I have all the bam file and would like to do the downstream steps using picard and GATK, what can I do to restore the -M option except realign all the fastq because I have more than hundreds fastq files. Thank you ...
assembly sequence sequencing written 11 months ago by Peter Chung80 • updated 11 months ago by dariober9.9k
1
vote
0
answers
974
views
0
answers
is it .impute2 format equal to .gen format
... I would like to apply qctool software from .impute2 format to .bgen format. However, qctool -g HELLO.impute2 -og HELLO.bgen Error comes out: !! Error (genfile::BadArgumentError): In argument(s) filetype_hint="guess", inferred type = "" to function genfile::SNPDataSource::create(): Unrecognised f ...
impute2 imputation qctool written 20 months ago by Peter Chung80 • updated 19 months ago by Biostar ♦♦ 20
0
votes
0
answers
1.3k
views
0
answers
Comment: C: Subsetting individuals from IMPUTE2 ref files: .hap and .legend
... I have that question too, https://www.biostars.org/p/239067/ do you know how to do it now? I would like to extract CHS CHB population out from .sample, .legend and .hap file of impute2 reference data phase3. ...
written 2.1 years ago by Peter Chung80
2
votes
0
answers
1.2k
views
6 follow
0
answers
imputation only choose one population as reference panel
... I have a gwas data to do the imputation by shapeit2 and impute2. I would like to use only east asians population in 1,000 Genomes haplotypes -- Phase 3 integrated variant set release in NCBI build 37 (hg19) coordinates from impute2 reference data: https://mathgen.stats.ox.ac.uk/impute/1000GP_Phase3 ...
shapeit impute2 population gwas written 2.1 years ago by Peter Chung80 • updated 22 months ago by ShirleyDai40

Latest awards to Peter Chung

Popular Question 3 months ago, created a question with more than 1,000 views. For TCGA UUID change back to TCGAbarcode
Popular Question 10 months ago, created a question with more than 1,000 views. For TCGA UUID change back to TCGAbarcode
Appreciated 10 months ago, created a post with more than 5 votes. For question about the impute2 output format
Popular Question 18 months ago, created a question with more than 1,000 views. For TCGA UUID change back to TCGAbarcode
Popular Question 18 months ago, created a question with more than 1,000 views. For looking for TCGA gene expression Normal Samples
Student 18 months ago, asked a question with at least 3 up-votes. For question about the impute2 output format
Supporter 2.1 years ago, voted at least 25 times.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1728 users visited in the last hour