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

IndexError: list index out of range #13

Open
SAADAT-Abu opened this issue Sep 17, 2023 · 3 comments
Open

IndexError: list index out of range #13

SAADAT-Abu opened this issue Sep 17, 2023 · 3 comments

Comments

@SAADAT-Abu
Copy link

I install fings and the test run went well but with actual data it is giving this error

File "/mnt/beegfs/asaadat/.local/bin/fings", line 8, in <module>
    sys.exit(main())
  File "/mnt/beegfs/asaadat/.local/lib/python3.9/site-packages/fings/FiNGS.py", line 223, in main
    applyfilters(tdata,ndata,sdata,pdict,resultsdir,vcfpath,referencegenome,args.PASSonlyin,args.PASSonlyout)
  File "/mnt/beegfs/asaadat/.local/lib/python3.9/site-packages/fings/filter_functions.py", line 83, in applyfilters
    ordered_fieldnames = flist[0].keys()
IndexError: list index out of range

Would be glad if someone can help

@cpwardell
Copy link
Owner

Can you upload the VCF file you're trying to filter and the log file from the failed run to some cloud storage for me to have a look at?

@SAADAT-Abu
Copy link
Author

Here is a link to the data. https://drive.google.com/drive/folders/1jUhxJeNU6gXq9016M29auCTLYI2tCrs5?usp=sharing
Thanks a lot.
A. Saadat

@cpwardell
Copy link
Owner

Your log and input file look good. Let's try some more things.

1.) Minor issue, shouldn't case any problems:
You've included a reference genome, but you're using the default filtering settings that don't need a reference genome.
Out of interest, can you tell me how many lines are in your reference genome?

wc -l Homo_sapiens_assembly38.fasta

2.) Are your bam files ok?
What sequencing technology was used to produce your data? Is it Illumina paired-end sequencing?
It looks like your VCF was produced using CRAM files; are you sure converting them to BAM format went ok? Can you share your bam files in the same google drive as your VCF? If you upload your CRAM files I might add CRAM support.
Alternatively, can you use samtools to give me some snippets from them? Please do this for both HK11.bam and WT11.bam.

How many reads are in the bam files?

samtools view HK11.bam | wc -l

What is the header like?

samtools view -H HK11.bam

I'm sure we can solve your problem. When that's done, I'm going to update FiNGS to clean it up, lots of the packages it depends on have changed and it's making installing/running it difficult for new users.

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