Wiebe Kooistra
- Total activity 11
- Last activity
- Member since
- Following 0 users
- Followed by 0 users
- Votes 0
- Subscriptions 2
Comments
Recent activity by Wiebe Kooistra Sort by recent activity-
Hey, Sure. Here is the VCF output line: SL4.0ch11 25548831 . T C,A 4.09785e+06 PASS AC=1,0;AF=0.523,0.006162;AN=2;BaseQRankSum=-0.295;DP=257674;ExcessHet=-0;FS=0;InbreedingCoeff=0.6408;MLEAC=819,8;...
-
Hi Genevieve Brandt (she/her) I did read about the local realignment. There are 173 ArtificialHaplotypeRG's marked sequences. That the numbers don't add up I understand, I was just wondering why my...
-
Hi Pamela Bretscher, I do not have a bamout image for this. The command that was run: gatk HaplotypeCaller -R ref.fasta -ERC GVCF --dont-use-soft-clipped-bases true --do-not-run-physical-phasing tr...
-
Hey Pamela Bretscher The full GVCF line is: Ls_98839 100 . A <NON_REF> . . END=100 GT:DP:GQ:MIN_DP:PL 0/0:2007:0:2007:0,0,0 The stack trace is: 15:11:08.430 INFO NativeLibraryLoader - Loading libgk...
-
Hey Pamela, For the GVCF file do you want to see multiple lines or just of the SNP location that I am expecting? For the full stack trace? Do you mean the output HaplotypeCaller is giving when runn...
-
Hey Pamela Bretscher, I am using HaplotypeCaller and not Mutect2, but I did change --pruning-lod-threshold to the same levels as the gentlemen did (1.3, 0.5 and 0.1). The genotype and the GQ score ...
-
Hi Pamela Bretscher, Thank you for you answer. I looked at all suggestions in both links. I found that if I turn off the downsampling, the GQ value becomes 0: Ls_98839 100 . A <NON_REF> . . END=100...