Question: Weird fastq format
gravatar for tw30282
15 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 • 540 views
ADD COMMENTlink modified 15 months ago by genomax51k • written 15 months ago by tw3028210
gravatar for mastal511
15 months ago by
mastal5111.9k 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 15 months ago • written 15 months ago by mastal5111.9k
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: 945 users visited in the last hour