getting paired end datasets from SRA
I am searching SRA by keywords like "paired-end", but the sra-toolkit seems to only download one file (single-end) about 90% of the time. I just want to make sure my commands are all correct:
prefetch SRR13310323
fastq-dump -I --split-files --outdir fastq --gzip --skip-technical --readids --read-filter pass --dumpbase --split-3 --clip ~/ncbi/public/sra/SRR13310323.sra
The SRA entry says "Layout: PAIRED". It's totally possible that whoever uploaded entered incorrect metadata, but the SRA curation step could easily have validated whether it was paired-end or not, and generally prevented a lot of grief. so am i missing something in my command? can i filter my search in SRA differently to get the desired accession numbers?
• 1,962 views
•
link
updated
by
Kevin Blighe
•
written
by
from the mountains
•
written
by
from the mountains
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
scRNA-seq from SRA
written by ntsopoul •
-
Download and convert SRA file to fastq file
written by rescueson •
-
Paired-end dataset from SRA database
written by melissachua90 •
-
Downloading multiple studies SRA data through while loop
written by jobrien20 •
-
Problems with fastq-dump generating a single fastq file
written by pomodoro_sinensis •
-
faste-dump for single-end data
written by eva.lucarelli •
-
Batching fastq download w/sra-toolkit: Won't recognize IDs
written by harka14 •
-
From where to get Sub-telomeric Region Coordinates/ sequences for each chromosome
written by rohitsatyam102
-
SRA to FASTQ problem
written by m.t.lorenc •
-
Downloading single cell data from NCBI
written by V
'99% of the time', this works 'all the time':
It will download anything there on the record. I then prefer to inspect the files and retain whatever I need. It should split reads into R1 and R2