gossamer icon indicating copy to clipboard operation
gossamer copied to clipboard

Published 20 hours ago • verified publisher icon data61

Xenome classify output file

Open ssinghk opened this issue 3 years ago • 2 comments

Hello, I am using xenome classify to process xenograft RNA-seq data. The output species specific .fastq files are missing '@' or '>' and the STAR aligner is unable to map.

ssinghk avatar Aug 08 '21 20:08 ssinghk

Which version are you using? I see this bug in xenome version 1.0.1 downloaded from https://github.com/data61/gossamer/issues/9#issuecomment-402958742

I use sed '1~4s/^/@/' each_output.fastq to get around it.

RamRS avatar Mar 03 '22 16:03 RamRS

Hi,

I would also suggest to introduce a + in the 3rd line of each read with sed -i 's/^$/+/' each_output.fastq. I found some problems in downstream analyses using the "FASTQ" format from Xenome.

dgomezpere avatar Jun 19 '22 20:06 dgomezpere