problem with sra toolkit fastq-dump sratoolkit.2.1.10-win64

[hui_shi]

Hello,

I was trying to use fastq-dump to convert a .sra file to .fastq format.

The data contains paired end reads, so when I type in commands,
fastq-dump --split-files .sra, it gives two files, *_1.fastq and *_2.fastq.

But the problem is the size of the two files are different, one is a lot bigger than the other. This doesn't seems to be right.

Does anyone know why this happens? and how to fix it? Shouldn't the two files contain the same number of reads, so that the size should be the same as well?

Thank you,

Hui

[vadim]

Could you post the run accession please?

[hui_shi]

Thanks for your reply. I didn't use --accession or -A, don't know if that's what you mean. Actually, I don't understand why using -A, except to modify the output name.

The full command I typed in was:

fastq-dump --split-files SRR443885.sra

it generates two files, SRR443885_1.fastq and SRR443885_2.fastq, but they are of different size.

Thanks,

Hui

[vadim]

The spot descriptor claims that the read length is 150bp, 75 forward and 75 reverse. But the actual read length is 100 bases only. I guess fastq-dump just obeys the rules by splitting the read at 75, so the reverse read becomes 3 times smaller than the forward.

The run metadata is here:

http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=run_browser&run=SRR443885

The experiment metadata is here:

http://www.ncbi.nlm.nih.gov/sra/SRX116341?&report=full

You can contact them at sra@ncbi about this.

[emp]

regarding small sra format file

I wanted to convert .sra file of size 289.98 Mb. when i am using sratool fastq-dump.. it is generating error message as:-



An error occurred during processing.
A report was generated into the file '.../ncbi_error_report.xml'.
If the problem persists, you may consider sending the file


kindly help me out for this ASAP

[dpryan]

Did you try looking at what was in the file? Given that it's an error report, perhaps it's informative...

[emp]

Respected dpryan,


The file is in xml format which on opening with google chrome is giving code. It does'nt have any error specified there in...

Here in is the file which i am trying to convert to fastq.

http://www.ncbi.nlm.nih.gov/sra/SRX370497

[dpryan]

Try upgrading your copy of the SRA toolkit, that usually fixes this sort of issue (I tried to extract the file you linked and couldn't until I upgraded my local copy of fastq-dump).

[emp]

kindly tell me how can i upgrade sra tool kit??


Is it through verbose??

[dpryan]

Try here, on NCBI

[tw7649116]

Dear hui_shi,

I have the same problem with convert the sra to pair-end fastq.
Do you remember how to solve this problem?

Thank you very much!!

[GenoMax]

Are you using the latest SRAtoolkit? Have you done the steps described in this post: http://seqanswers.com/forums/showpos...6&postcount=7?

[tw7649116]

Thank you GenoMax,

I used the latest version.
In that way we need to redownload the file?
Thank you! If I can't solve this, may be I need to download again!

Best wishes

[GenoMax]

If you have downloaded the .sra file have you tried the command this way?

c:\> fastq-dump.exe -F --split-files c:\path_to\SRAfile.sra

Sometimes people do submit unequal length data sets so if you have mismatched cycle numbers that is a possibility. What SRA# are you looking at?