Question: Weird fastq format
gravatar for tw30282
22 months 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 • 686 views
ADD COMMENTlink modified 22 months ago by genomax62k • written 22 months ago by tw3028210
gravatar for mastal511
22 months 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 22 months ago • written 22 months 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: 744 users visited in the last hour