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

Permant error in the first step #134

Open
beginner984 opened this issue Aug 15, 2023 · 1 comment
Open

Permant error in the first step #134

beginner984 opened this issue Aug 15, 2023 · 1 comment

Comments

@beginner984
Copy link

Hi

What does this error means please?


fi1d18@RBGO-Server2:/data/Continuum/WES/vcf$ ./lofreq_star-2.1.0/bin/lofreq call -f resources_broad_hg38_v0_Homo_sapiens_assembly38.fasta -o vars.vcf NG-27280_CLTSS_LTS_001A_lib506241_7636_2_MarkedDup.bam
WARNING(lofreq_call.c|warn_old_fai): Index for fasta file (resources_broad_hg38_v0_Homo_sapiens_assembly38.fasta.fai) is older than fasta file! You should reindex (using faidx)!
[fai_fetch_seq] The sequence "NC_000001.11" not found
FATAL(plp.c|mplp_func:645): Can't compute BAQ or IDAQ without reference sequence
fi1d18@RBGO-Server2:/data/Continuum/WES/vcf$ 

I am just certain that the reference genome in this command and for fastq file alignment was the same.

Thanks for any help

@andreas-wilm
Copy link
Contributor

Well something is off with the reference. The index is out of date for one. Please reindex with samtools faidx. But it also seems like the reference used here and for mapping were different. You can compare what you get from samtools view -H and what you see in the .fai 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

2 participants