Examples: Monday, today, last week, Mar 26, 3/26/04
Avatar

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Badges

Third Anniversary250 LikesSecond Anniversary2500 Comments25 Agrees100 LikesFirst Anniversary1000 Comments500 Comments25 Likes5 Agrees5 Likes100 Comments10 CommentsPhotogenicName DropperFirst Comment

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community

Sheila admin

About

Username
Sheila
Location
Broad Institute
Joined
Visits
2,040
Last Active
Roles
Member, Broadie, Moderator
Points
2,753
Badges
17
Location
Broad Institute
Posts
4,872

Activity

  • bassu

    Hi Sheila,

    I have some question with regard on how to move forward to make my analysis more meaningful. As discussed in my previous question with VQSR. I have 10+ RNASeq data sets I have followed the GATK best practices using STAR2pass method. I have now 10+ vcf file using the following commands

    gatk -Xmx30048M -T HaplotypeCaller -R Homo_sapiens.GRCh38.dna.primary_assembly.fa -I recal.bam -dontUseSoftClippedBases -stand_call_conf 20.0 -stand_emit_conf 20.0 -o bqsr.vcf

    gatk -Xmx30048M -T VariantFiltration -R Homo_sapiens.GRCh38.dna.primary_assembly.fa -V bqsr.vcf -window 35 -cluster 3 -filterName FS -filter "FS'>'30.0" -filterName QD -filter "QD'<'2.0" -o filtered_bqsr.vcf

    What I currently do is use bcftools to merge the 10+ vcf file and do the downstream analysis. Should i do joint genotyping instead (my sample size will increase as I move)? I am quite confused which will be the best method.

    May 24