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

Error reading file issues while indexing #10

Open
srbehera opened this issue Sep 9, 2022 · 5 comments
Open

Error reading file issues while indexing #10

srbehera opened this issue Sep 9, 2022 · 5 comments

Comments

@srbehera
Copy link

srbehera commented Sep 9, 2022

I used excord to create the bed.gz file and then run giggle index. But I am getting the following error
giggle: Error reading file "A94_Exp5_6/cache.0.idx": End of file

The command that is used:
giggle index -i A94_Exp5_6.bed.gz -o A94_Exp5_6 -s -f

@ryanlayer
Copy link
Owner

What was your excord command?

@ryanlayer
Copy link
Owner

Can you try putting A94_Exp5_6.bed.gz into a folder by itself then running
giggle index -i "folder/*gz" -o A94_Exp5_6 -s -f

@srbehera
Copy link
Author

I tried that initially, but it gave me the same error. So I was checking only one file at a time. Here is my excord script

#!/bin/bash
input="/DIR_PATH/P70_06.txt"
excord="/DIR_PATH/excord"
ref="/DIR_PATH/genome_mainchrs.fa"
out="/DIR_PATH/output/"
while IFS= read -r line
do
  f="$(basename -- $line .bam)"
  samtools view -b $line | ${excord} --discordantdistance 500 --fasta ${ref} /dev/stdin | LC_ALL=C sort --buffer-size 2G -k1,1 -k2,2n -k3,3n | bgzip -c > ${out}${f}.bed.gz
done < "$input"

@ryanlayer
Copy link
Owner

ryanlayer commented Sep 14, 2022 via email

@srbehera
Copy link
Author

Sure. I am sending you the link to the box folder.

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