get read length from sra
0
0
Entering edit mode
3.0 years ago
Nobody ▴ 30

Hi everyone!

I need a method "api" to retrieve reads length from SRA.

Thanks in advance

SRA rna-seq • 1.6k views
ADD COMMENT
1
Entering edit mode

You will only get average length of reads (which may or may not be very useful) using Entrezdirect.

$ esearch -db sra -query "SRR13310323" | efetch -format runinfo | cut -d ',' -f 1,7
Run,avgLength
SRR13310323,109
ADD REPLY
1
Entering edit mode

The suggested methods works. Here is another method using the SRA toolkit

fastq-dump -X 1 -Z SRR5790106 | sed -n '2p' | awk '{ print length }'

Read 1 spots for SRR5790106

Written 1 spots for SRR5790106

100

this assumes that read length is the same for all reads

ADD REPLY

Login before adding your answer.

Traffic: 877 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6