Is there a way to do it? Sorry for the uninformative question, so I have downloaded an SRA file from NCBI and used included sratoolkit to split the file into two fastq sequences. I am trying to do a de novo assembly using these paired-end strand_specific reads. However, a required parameter is the average insert size. Does anyone know how to obtain this from an SRA file or fastq?
2 answers
Guessing an insert size length, assembling, mapping to the assembly, and then iterating with the improved insert size length (from the mappings) is a reasonable choice, and probably about the best you can do. You hopefully should have some rough idea from the library preparation method (size selection criteria or if it's jumping library or not).
In fact, Velvet does this automatically (from the 1.1 manual): "If the insert length of a library is unspecified, Velvet will attempt to measure it for you, based on the read-pairs which happen to map onto a common node." As they point out, it's critical to check the reported estimate to make sure it's sane.
Log in to answer this question.
-
Paired-end dataset from SRA database
written by melissachua90 •
-
How can i get a merge fastq file from sratoolkit like the file i directly download from ncbi?
written by Little.G •
-
Converting a bunch or SRA files using fastq-dump split-files
written by Mitra •
-
SRA toolkit (NCBI) - sra to fasta
written by t.kranenburg •
-
Problem in processing SRA file from NCBI
written by majeedaasim •
-
PolyA selection for transcriptome sequencing
written by Mehmet
-
Assembly pipeline De Novo | Reference
written by elbecerrasoto •
-
How can I change the default path to which SRA files are downloaded?
written by Cacau
-
fastq to NCBI sra file format
written by ewre
-
Determining Orientation Of Fastq Files
written by dan79 •
Please describe your question so people can help you. I think I understand what your asking, but without more information it is difficult to answer.
Edited, thanks.
You will need to align the reads (both pairs). Then you can find the insert lengths by parsing the SAM/BAM file.
Align the reads to a reference genome? This seems counterintuative considering the whole point of a de novo assembly is to not need a reference.
Good point. Sorry. I need to read more carefully. I don't know the answer. I look forward to seeing the best solution.