How to fix FreeBayes error "unable to find FASTA index entry for bed files"
I'm trying to use freebayes to call variants through the following parameters using RNA-Seq data:
$freebayes --fasta-reference $ref \
--region $BED \
$bam \
> ${outdir}/${sample}.freebayes.vcf
and my bed files like this:
chr1:11868-12227
chr1:12612-12721
chr1:12974-13052
chr1:13220-14501
chr1:15004-15038
.....
but report an error:
unable to find FASTA index entry for '/path/to/GRCh38.exon.bed'
I'm pretty sure my fasta files have the index with the suffix .fai
Why did this happen?
• 1,389 views
•
link
updated
by
Ram
•
written
by
wlei091226 •
•
written
by
wlei091226 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
ATACseq Data Insert Size Plot
written by lingjun meng •
-
Help manually processing strand in paired-end reads when fishing for lariats in bulk RNA-seq
written by txema.heredia •
-
How to interpret the exon numbering for minus strand?
written by sebastian.swirski •
-
Unable to find FASTA index entry for 'unitig_2993_pilon error.
written by tmrhyd •
-
featureCount vs GENCODE gene lengths
written by user31888
-
Exon parsing from bed file
written by 1769mkc
-
Parsing GENCODE GTF to simpler BED files: am I reinventing the wheel ?
written by Charles Plessy
-
overlap gene and exon with bed file
written by bioguy24
-
Bed_intersect -s empty output
written by mileti.enrichetta •
I have the same issue. I check correspondence between fasta headers line and bed column. Both have chr1 as "prefix". May that fasta file has, other than Chr1 to Chr22, ChrM, Unknwon and HLA--### contigs while bed file not?