Question: Xenome Star Aligner
0
gravatar for lkmklsmn
4.9 years ago by
lkmklsmn870
United States
lkmklsmn870 wrote:

Hi,
after having classified reads from xenograft experiment using xenome. I try aligning human reads using STAR aligner.
Pseudo code looks like this: STAR --genomeDir genome --readFilesIn human_1.fastq human_2.fastq
Then I get the error:

EXITING because of FATAL ERROR in input reads: unknown file format: the read ID should start with @ or >

Apr 11 12:13:31 ...... FATAL ERROR, exiting [samopen] SAM header is present: 25 sequences. [sam_read1] reference 'SN:chrM LN:16571 66

What does this error mean? Are the human_*.fastq xenome outputs in a format unrecognizable by STAR?

Thanks

rnaseq • 2.6k views
ADD COMMENTlink modified 4.9 years ago by Irsan6.8k • written 4.9 years ago by lkmklsmn870

Can you post the first few (10 or so) lines of each of the fastq files?

ADD REPLYlink written 4.9 years ago by Devon Ryan88k
0
gravatar for Irsan
4.9 years ago by
Irsan6.8k
Amsterdam
Irsan6.8k wrote:

yup, you just need to add the @ character at the beginning of each seq id in your classified fastq.

cat human_1.fastq | paste - - - - | sed 's/^/@/g' | sed 's/\t/\n/g' > human_repaired1.fastq
ADD COMMENTlink modified 4.9 years ago • written 4.9 years ago by Irsan6.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 887 users visited in the last hour