Question: Weird fastq format
gravatar for tw30282
2.6 years ago by
tw3028210 wrote:

Hi all,

I downloaded a RNAseq raw data and found it in weird format. The data is generated by ABI_SOLID platform. Data ID is SRR179593.fastq and here are a few lines inside it:

@SRR179593.1 mendel_20110320_FRAG_BC_Ryan_RNA_Seq_2_58_240/1
@SRR179593.2 mendel_20110320_FRAG_BC_Ryan_RNA_Seq_2_58_300/1
@SRR179593.3 mendel_20110320_FRAG_BC_Ryan_RNA_Seq_2_58_638/1

I tried to map it to the genome using STAR:

cd .
/usr/local/star/2.4.1c/bin/STAR --runThreadN 4 \
--genomeDir /escratch4/tw30282/tw30282_Mar_17/STAR/ENSEMBL.homo_sapiens.release-75 \
--readFilesIn SRR179593.fastq.gz \
--outFileNamePrefix Wei_NPC_rep1_STAR_SE. \
--outSAMtype BAM SortedByCoordinate \
--readFilesCommand zcat

But an error returned and said input file in wrong format.

Is there anyway to transform or fix the format?

Thanks, Tianming

rna-seq format fastq sra • 878 views
ADD COMMENTlink modified 2.6 years ago by genomax74k • written 2.6 years ago by tw3028210
gravatar for mastal511
2.6 years ago by
mastal5112.0k wrote:

ABI SOLiD data is in colorspace rather than basespace, which is why you have 0123 instead of ACGT for the sequence. You need to use an aligner that works with SOLiD data. Normally, if you weren't getting the data from SRA the SOLiD data would have 2 separate files, one with the sequence reads and one with the read qualities.

Tophat2 works with colorspace data, you have to specify the options that your data and genome index are in colorspace, but I don't see anything in the Hisat manual so I don't think it does.

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by mastal5112.0k
Please log in to add an answer.


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