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.

Filter by
Sorted by
Tagged with
0
votes
1answer
20 views

How to de novo hybrid assemble with Pacbio CCS and Illumina PE reads

I would like to perform de novo genome assembly on a diploid microalgal strain. I have two datasets: PacBio CCS/HiFi reads, low coverage. Illumina PE 2x150 (standard shotgun) Does anybody have any ...
1
vote
1answer
35 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 "...
6
votes
1answer
70 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),...
0
votes
0answers
6 views

Mapping contigs back to samples from co-assembled metagenome

In an attempt to increase the quality of our metagenomic assembly and ensure we capture low abundance species we are co assembling ~20 samples. However after the assembly we need to perform some ...
-1
votes
2answers
25 views

How do you set the coverage in PacBio's Sequel II?

I am reading the Whole Genome Sequencing for de novo Assembly Best Practices Use the Sequel II or IIe System and SMRT® Cell 8M to sequence to desired coverage depth for complexity of genome 10- to 15-...
2
votes
3answers
174 views

What assembler is appropriate for High-Fidelity PacBio reads

What assembler is appropriate for High-Fidelity PacBio reads? For example, canu is good for high-error PacBio reads. But what algorithm to use for HiFi reads? Would it be OK to use canu without the ...
-1
votes
2answers
56 views
4
votes
0answers
58 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 ...
-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 ...
3
votes
2answers
146 views

Merge / Reconciliate several de novo transcriptome assemblies with different kmers

I am building a De Novo transcriptome reference assembly for an eukaryotic organism for which I have a genome. I've created several assemblies with rnaSpades using different kmer sizes (19 to 69 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 ...
0
votes
1answer
29 views

How to know if the DNA sequence has been assembled and why is it important to know how it was assembled?

I have downloaded my FASTA format files, that have the DNA sequences of the coding region of the genes and the DNA sequence of the complete genome, from NCBI. How can I recognize if these sequences ...
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 ...
1
vote
0answers
30 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
1answer
36 views

Trinity assembly from many samples

When you combine samples for de-novo transcriptome assembly with Trinity, do you suggest limiting the number of reads for each sample? I had read one of Matthew MacManes' papers awhile back suggesting ...
0
votes
2answers
66 views

Genome annotation

I have helped a lab sequence a mitochondrial genome. I used then MITOS to annotate the genome (warning them that I did it just for curiosity as I am not experienced). They submitted the sequence and ...
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: ...
3
votes
1answer
1k views

SPAdes error during assembly

I try to perform a hybrid assembly using Unicycler. Unicycler use SPAdes to assembly the illumina sequences, and the program crashes while in SPAdes. The spades.log file says something about an ...
3
votes
2answers
264 views

Is there any value in scaffolding the output contigs of MEGAHIT assembler given a metagenomic dataset?

As far as I understood, for most assembly programs, the scaffolding step takes into consideration paired-end information in order to get from contigs (contiguous sequences) to scaffolds (longer ...
1
vote
1answer
41 views

where to find sample contig data

Where or How can I find contigs? I am trying to learn Bioinformatics and I want to make a reference-based gene search. I also want to align contigs with a given reference sequence. From NCBI other ...
2
votes
1answer
198 views

Coverage required

I was came across a problem during an exercise in a book and I don't really know how to solve it. I feel like something's missing. "coverage, c = $NL/G$ (N=number of reads, L=read length, G=genome ...
0
votes
3answers
36 views

Genome QC + Assembly Pipeline semantics

I’m trying to create a pipeline for genome assembly. How best can I “redirect/pipe” from existing fasta files (or files in general) to other steps of the pipeline? I was thinking of going from the SRA ...
3
votes
1answer
21 views

Can SBOL encode DNA Assembly linkers?

Where in an SBOL file could you specify the linkers that would be placed as prefixes / suffixes between parts for assembly? Lets say you have a final construct of your design and want to assemble the ...
0
votes
2answers
66 views

creating a tab delimited file

I am working on a project using a fasta file. I am writing my command in nano within command-line and executing using python, also within my command-line. I would like my command to provide me with a ...
1
vote
2answers
99 views

Is it possible to filter contaminated reads for raw PacBio sequences (not HiFi reads) before assembly?

De novo genome assembly for non-model organisms face the issue of bacterial contamination. For assembled contigs with mostly bacterial-like sequences (based on BLAST search), the entire contig can be ...
7
votes
1answer
335 views

wtdbg2: practical implications of k-mer fsize and psize choice

I am using wtdbg2 2.3 to assemble a human genome (sequenced on PromethION from a cell line). I filtered out reads with low average quality, and now I am trying to determine the parameters that will ...
0
votes
1answer
32 views

Quast duplication ratio and mismatches percent

When analyzing Quast results it seems that it doesn't calculate mismatches and indels in a useful way if the "Duplication ratio" is over 1. For example, that's what I get for an assembly ...
2
votes
1answer
38 views

Gold standard benchmark

This page is claimed to contain a gold standard benchmark for viral genome assembly. https://github.com/cbg-ethz/5-virus-mix The claim is here: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5411778/ &...
18
votes
3answers
661 views

How to deal with heterozygosity during polishing of genome assembly based on long reads?

All the long-read sequencing platforms are based on single-molecule sequencing which causes higher per-base error rates. For this reason a polishing step was added to genome assembly pipelines - ...
1
vote
4answers
79 views

Tools for comparing/visualizing FASTAs?

I have two FASTAs/assemblies of the same species, but the bases are somewhat different. I would like to explore this. What tools/methods exist to compare two FASTAs and see the difference in ...
1
vote
2answers
75 views

How can I get or create a reference genome for Bacteria?

I am a computer engineer and nowadays trying to grasp some concepts of Bioinformatics particularly, reference genomes and genomic variants. My aim is to find the effect of sequence features on variant ...
1
vote
1answer
65 views

Which one is a more convenient assembly?

I have developed a software for de novo genome assembly. Its performance varies gradually according to how much data you employ. At initial stages it often produces contigs that look like that when ...
0
votes
0answers
51 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 ...
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 <...
4
votes
1answer
77 views

Viral genome assembly using broad viral ngs pipeline?

I am trying to assemble RNA virus genome using Broad Viral NGS pipeline BROAD VIRAL NGS PIPELINE. I am two questions : 1) As this pipeline requires unaligned bam format as input, how do I convert ...
2
votes
1answer
118 views

What is the meaning of these misaligned reads in a sequencing run?

I am analyzing some SARS-CoV-2 sequencing runs abd often find read alignments like the one in the image. ...
1
vote
2answers
89 views

Bacterial DNA at the tail of transcriptome reads. What does that mean?

I am assembling a transcriptome obtained from the Internet. The transcriptome was extracted from a human cancer tissue that had been previously grafted into a mouse. I have detected that many ...
3
votes
1answer
86 views

Viral genome finishing

I have assembled poxvirus genome using Ray. The assembly is good. Out of several thousand contigs I got, I was able to get one scaffold using Contiguator tool, which is about 90% of my genome. I have ...
1
vote
1answer
44 views

Is there a simple command for outputting a tab delimited columns?

I am working on a fasta file and am writing my command in nano within command-line and executing using python, also within a command line. My objective is to get my command to provide me with a tab ...
3
votes
2answers
115 views

Contamination on genome assembly

I had a question for the community. I have a genome of a new species that has been sequenced via 150pb Illumina paired-end. To verify the quality of the assembly I used the ...
0
votes
1answer
23 views

Is it possible to convert BAM file from one genome assembly to the other?

I Have multiple BAM files that are referenced to UCSC genome assembly GRCh37/hg19 that are read in different time frames. Now, I am planning a different studies that require assembling all the data ...
3
votes
1answer
60 views

How to find all WGS assemblies accessions of a species

Some background Similar to the OP of https://www.biostars.org/p/377840/, I would like to programmatically BLAST a sequence to a local database of all WGS assemblies. Since this isn't feasible for the ...
2
votes
2answers
166 views

Genome assembly of SRR12196449 with SPAdes

I am trying to assemble the run SRR12196449 with SPAdes. The description of their project is: This project expected to standardize a method for amplification and ...
1
vote
1answer
34 views

Finding simple sequence from reads with significant overlap

I wrote a "script" to pull out reads from a huge fastq file in an iterative manner, by finding homology to the previous sequence. It should be relatively easy to overlap them and assemble ...
2
votes
1answer
82 views

PacBio long-reads impact in transcriptome de novo assembly?

We are strongly interested in assembly a good transcriptome of reference for a non-model organism and build a local database. We have sequenced the same individual with Illumina (150 millions of pair-...
1
vote
1answer
200 views

Difference between paired-end, mate-pair and long read

I writing here because I have some questions for you. I wondered what the essential differences were between paired-end, ...
0
votes
1answer
40 views

After KEGG and GO analysis, how to make tables+phylogenetic trees

hope everyone is ok. I used Trinity to do a de novo transcriptome assembly, then blastp/blastx and then used Blast2GO software to do KEGG and GO analysis. So i got some txt files with header : for GO ...
1
vote
1answer
101 views

RCSB API: limited returned result

Is there a limit of query results on RCSB API? Interestingly, no matter how to adjust my query criteria, only 10 assemblies are returned. Here is an example query json: (truncated from a longer query ...
1
vote
1answer
63 views

Gap-fill assembly with PacBio reads

I would like to gap-fill (or correct) a mammalian-sized assembly for which we have BACs that have been sequenced with PacBio CLR data. I have seen there are pbm22/gcpp tools available for correcting ...
0
votes
0answers
66 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 ...