I want to do a comparison of RNA-Seq datasets of same tissues in multiple species. I quantified the RNA-Seq using kallisto and now I want to read the kallisto files using tximport for DESeq2 analysis. My abundance files have identical transcript id for each species but the transcript lengths differ due to differences in gene lengths in diferent species. When I try to read the files using Tximport I get following error.
> sample <- read.table("samples.txt", header = T) > files <- file.path(dir, "kallisto_results",sample$samplenames, "abundance.tsv") > txt2gene <- read.table("Transcript-to-Gene-mapping.txt", header=T, sep = "\t") > txi.kallisto <- tximport(files, type="kallisto", tx2gene = txt2gene) reading in files 1 2 3 Error: all(txId == raw[[txIdCol]]) is not TRUE
The first two are for same species and the 3rd one with the error is for different species. Can someone please help me with this? I am really stuck as I cannot find out what this error mean. My guess is that this is due to difference in transcript length for the same Id in different species. But I am not sure!
Please let me know if you know what's going on here?
Thanks so much in advance!