User: Peter Chung
Peter Chung • 80
- 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
... 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 ...
written 4 weeks ago by
Peter Chung • 80
0
votes
0
answers
178
views
0
answers
... 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 Chung • 80
• updated
8 weeks ago by
finswimmer ♦ 11k
0
votes
0
answers
178
views
0
answers
... 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 \
...
written 8 weeks ago by
Peter Chung • 80
1
vote
1
answer
183
views
1
answer
... 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 ...
written 3 months ago by
Peter Chung • 80
1
vote
0
answers
225
views
5 follow
0
answers
... 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=${ ...
written 3 months ago by
Peter Chung • 80
0
votes
0
answers
483
views
0
answers
... 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 ...
written 10 months ago by
Peter Chung • 80
• updated
10 months ago by
RamRS ♦ 20k
0
votes
2
answers
432
views
2
answers
... 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 ...
written 11 months ago by
Peter Chung • 80
• updated
11 months ago by
dariober ♦ 9.9k
1
vote
0
answers
974
views
0
answers
... 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 ...
written 20 months ago by
Peter Chung • 80
• updated
19 months ago by
Biostar ♦♦ 20
0
votes
0
answers
1.3k
views
0
answers
... 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 Chung • 80
2
votes
0
answers
1.2k
views
6 follow
0
answers
... 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 ...
written 2.1 years ago by
Peter Chung • 80
• updated
22 months ago by
ShirleyDai • 40
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.
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