Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

after running gimble parse, gimble.vcf.gz is null #29

Open
kuangzhuoran opened this issue Jun 7, 2024 · 7 comments
Open

after running gimble parse, gimble.vcf.gz is null #29

kuangzhuoran opened this issue Jun 7, 2024 · 7 comments

Comments

@kuangzhuoran
Copy link

my command is gimbleprep -f ./galili.fa -b ./galili.bam -v ./galili.golani.rename.vcf.gz -k -t 40

@GertjanBisschop GertjanBisschop transferred this issue from LohseLab/gIMble Jun 7, 2024
@GertjanBisschop
Copy link
Member

Thanks @kuangzhuoran for reporting this.
Could you provide a bit more information so we can try and reproduce the error?

Note though that the -b flag should be followed by a directory of bam files, not an actual bam file.

@kuangzhuoran
Copy link
Author

gimble.log.txt

@kuangzhuoran
Copy link
Author

Thanks for the reply, the -b parameter follows a folder。
Do I need to take a portion of the input file and send it to you.
I have now uploaded the log file

@kuangzhuoran
Copy link
Author

I used gatk to extract the SNPs, but it was running successfully before. bam file is using sambaba to remove PCR duplicates!

@GertjanBisschop
Copy link
Member

GertjanBisschop commented Jun 7, 2024

Can you be more specific about gimble.vcf.gz is null? What does the output look like? Obviously, gimbleprep is a filtering step. Is it possible that after applying filters there are no more SNPs left?

@kuangzhuoran
Copy link
Author

gimble vcf
Only the header

@kuangzhuoran
Copy link
Author

“Is it possible that after applying filters there are no more SNPs left”

It's possible.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants