Questions tagged [assembly]

Process of creating the original sequence from the read sequences that it generated during a sequencing experiment. Can refer to genome assembly, in which case the original sequence is a genome, or transcripts assembly, in which case the original sequences are RNA transcripts.

16 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
5
votes
1answer
65 views

What is the current state-of-the-art in assembling hybrid transcriptomes?

We are considering attempting de novo assembly of a species transcriptomes (i.e. without a reference genome) using the combined NGS outputs of Iso-seq and Illumina. One example I saw (Li et al 2017),...
4
votes
0answers
57 views

Can someone help me estimating the runtime of the pipeline applied by the vertebrate genome project?

The vertebrate genome project (VGP) has a lot of interesting publications such as this one. The rough pipeline is outlined below: Here the pipeline in more detail: While the paper describes all the ...
3
votes
0answers
207 views

How to submit a canu job on LSF high-performance computing cluster farm?

I am currently running canu on an LSF Linux server using the following script called assemble.sh: ...
2
votes
0answers
170 views

Platanus-allee phasing fail: Error(13): Error, SolveDBG exception!

I am using Platanus-allee 2.2.2 for heterozygous genome (~500mb) assembly with Illumina short reads and PacBio reads input data. I have the file contigs.fa from short reads but phasing step with ...
2
votes
1answer
44 views

polishing assembled genome to QV50 value

I am trying to assemble some genomes from pacbio reads using HGAP3 followed by quiver polishing. However, even after multiple rounds of quiver polishing, I cannot attain quality value 50 (QV50) for a ...
2
votes
0answers
41 views

Calculating alignment/mapping time

I am trying to assemble a plant genome using AWS resources using velvet. Plant genome is huge (> 10 times human genome) and coverage is around 30 x. We are planning for de novo assembly with Velvet (...
1
vote
1answer
28 views

Assembling all transcripts for an individual gene? (using single sequence to seed the assembly)

Let's say I have a candidate gene and I believe that in an individual sample, the genome sequence differs from the reference which then interferes with alignment. Is there a way for me to do a "...
1
vote
0answers
29 views

How to quantifiy of specific genes from shotgun metagenome?

I have googled a "lot", couldn't find any specific answer to the question. So, I am here seeking for your guidance. My question is similar to this. I have several metagenome (n=30). But for ...
1
vote
0answers
14 views

Repeat analysis on eukaryotic assemblies

I have a hundred insect genomes and I'm looking for repeated regions along these assemblies. First of all, I thought of using <...
1
vote
1answer
58 views

Interpreting contig alignments to a reference genome

I have applied two de novo genome assembly tools to data from the run SRR12707453, corresponding to a phage (I downloaded the data and have no relation with the authors of the study). Using rnaSPAdes ...
0
votes
0answers
50 views

How to filter a genome assembly consistsing of a large number of contigs?

I did some de novo genome assemblies with Illumina PE data using SPAdes, whereas most of them consisting of a large number of contigs(>1000). I have several questions below. Do we need to filter ...
0
votes
0answers
63 views

Resources to learn genome assembly workflow for small genomes (like viruses)

I have sequencing data of a few samples of a DNA genome virus. I'd like to learn de novo assembly of the ...
0
votes
0answers
66 views

Gffcompare issue: 0 reference transcripts loaded

I am trying to use gffcompare to compare my assembled transcriptome to a reference gtf that contains information about small open reading frames (sORFs). The reference gtf was obtained by processing ...
0
votes
0answers
46 views

How to obtain the extended unaligned bases of the read in reference based genome assembly?

I did illumina read mapping based on a particular gene sequence (reference sequence) using the following command, ...
-1
votes
2answers
48 views

Denovo, Stacks: Getting an “ambiguous redirect” error

I tried running the following command ...
-1
votes
1answer
72 views

Will using smaller kmers help get larger contigs? If not, then what?

I've been using MEGAHIT to assemble metagenomes, with particular focus on specific genomic areas. Sometimes all I get is gene fragments or pathway fragments (eg. if I know that genes A, B, C, D and E ...