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

Unable to run Haptree #3

Open
vinaydeep26 opened this issue May 13, 2023 · 6 comments
Open

Unable to run Haptree #3

vinaydeep26 opened this issue May 13, 2023 · 6 comments

Comments

@vinaydeep26
Copy link

Hello I'm unable to Run Haptree for Some Reason, Please if you can reply soon that'll be great because we're working on a publication and will cite HaptreeX if it works.
I try to run the Command :
haptreex -v A_J_DBA_calls.vcf -d ../AJ_DBA_2J_sorted_using_Picard_Tools_dedup_withreadgroups_with_changed_sample_name.bam -o outputfile.txt &> log_test_hapree.txt &

but i always get an error like this:

Aborted                 (core dumped) haptreex -v A_J_DBA_calls.vcf -d ../AJ_DBA_2J_sorted_using_Picard_Tools_dedup_withreadgroups_with_changed_sample_name.bam -o outputfile.txt

ValueError: not BGZF compressed
raised from: sam_itr_next
/home/vagrant/.seq/lib/seq/stdlib/bio/c_htslib.seq:65:9

then i try to Bgzip the VCF file and tabix it using tabix -p vcf vcffile
then i get this error:

IndexError: list index out of range

raised from: _idx_check
/home/vagrant/.seq/lib/seq/stdlib/core/collections/list.seq:322:13
@vinaydeep26
Copy link
Author

vinaydeep26 commented May 13, 2023

Screenshot (295)
Screenshot (296)
Here are the screenshots. I'm not sure what am i doing wrong can you guys help me out please?
@inumanag @LillianZ @arshajii @HarisSmajlovic

@inumanag
Copy link
Collaborator

Can you please decompress VCF (via gunzip) and then see if it works with the uncompressed VCF?

@vinaydeep26
Copy link
Author

@inumanag In the first run i'm running with a decompressed file, then i compressed it and it started showing me the index error

@vinaydeep26
Copy link
Author

@inumanag @LillianZ

@pdimens
Copy link

pdimens commented Sep 26, 2023

@inumanag @LillianZ I see the same issue if using a BCF file or a bam file. It's not practical to be forced to decompress these large common datastores to perform this analysis

@pdimens
Copy link

pdimens commented Sep 26, 2023

@vinaydeep26 your issue is the sam file, I think, as the error called is raised from: sam_itr_next, suggesting it's the alignment file

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

3 participants