FastQTL fails to get genotype region

I am using fastQTL to do QTL mapping and got the following error:

ERROR: Failed to get region chr1:0-111675033

I think the tool is not very well documented and just do not know how to slove the problem.

fastqtl

• 2,434 views • link

written 7.6 years ago by Phoenix Mu •

What is the command-line you used to run FastQTL? Are you sure 111675033 bases falls within the length of chr1? Maybe region must be one-based indexed?

edit: it seems --region must be one-based:

The region of interest is specified with the standard chr:start-end format.

• 0 views • link

written 7.6 years ago by h.mon

Thank you for your reply. I already figured it out that the error was because I have different chromosome names in VCF files and phenotype files. One has "1" whereas the other has "chr1". There is no error after I used bcftools to rename the chromosomes in VCF.

• 0 views • link

written 7.6 years ago by Phoenix Mu •

0 answers

No answers yet.

More posts like this

Lifting over with rtracklayer vs UCSC
written 4 months ago by _quantum_girl_ •

When I lift this region: chr1 145686997 145808272 from hg38 to hg19 on UCSC browser I get: chr1 145626811 145748066. However if I do the …
R (tidyr::pivot_wider)
written 4.3 years ago by hosin •

Hi there, I just use tidyr::pivot_wider (https://tidyr.tidyverse.org/reference/pivot_wider.html) to get a matrix with following .csv: > Coordinates Numbers Values > chr1:110238914-110324454 NumB 1 > chr1:110238914-110324454 NumC …
Can't index .bed file with tabix
written 5.1 years ago by zjardyn •

Hello, I am trying to index .bed like files from Leafcutter output as described here: [Leafcutter sQTL][1] [1]: https://davidaknowles.github.io/leafcutter/articles/sQTL.html Leafcutter generates an .sh script to …
fastQTL phenotypes with no cis-variants
written 5.9 years ago by Cindy •

Hi, I've been using FastQTL to map cis-QTL for molecular phenotypes. Correct me if I'm wrong, but I just noticed that nominal pass is not …
Michigan imputation server quality control failed
written 7.4 years ago by Phoenix Mu •

I am using Michigan Imputation Server. I have performed a initial QC and removed monomorphic sites and sites with MAF > 0.5. However, the QC …
Few SNPs after merging my samples with 1000G phase 3
written 7.9 years ago by Phoenix Mu •

I used pruned SNPs from 1000G phase 3 to do a PCA, and got the following results: ![enter image description here][1] In the plot above, …
Please add an explicit type tag :NAME
written 8.3 years ago by mahdikhadem95 •

Dear all, I am totally new to exome-seq data analysis. Currently, I have got stuck in "base recalibration" step. When I use the following command: …
SHAPEIT error after run: x samples are included in several duos/trios
written 8.9 years ago by c_barbarosie •

Hi, I am new using SHAPEIT and I have a problem after running it. I am using the last version and I used the following: …
cwltool: Tool definition failed initialization
written 9.2 years ago by mepstein •

I had an old version of cwltool installed (i.e., the reference implementation, from https://github.com/common-workflow-language/cwltool), and when I tried to run it today I got: Tool …
How To Annotate A Human Dna Position
written 14.6 years ago by Free Man

<p>If I know a human DNA position and its code, for example: an "A" located in chr1:207510596-207510597(forward srand), how do I map this location to …