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.
139
questions
3
votes
1
answer
53
views
compare fasta sequences in pairs and collect metrics
I have 96 fasta files (A1, A2, A3...) from one plasmid assembly pipeline, and I have another 96 fasta files (B1, B2, B3 ...) from another plasmid assembly pipeline.
I would like to compare pair ...
- 89
1
vote
0
answers
44
views
Why do XRAY (but not CryoEM) structures of ribosome in PDB have 2 assemblies?
When i started programming against PDB i had a mixture of confusion & frustration with the fact that certain cif files contain two actual structures aka ...
- 332
1
vote
0
answers
15
views
DNASTAR viral-host integration assembly keeps failing
I have two NGS files from an NGS company corresponding to the sequencing data from a tumor sample as follows:
TB_7710391_R1.FASTQ.gz
TB_7710391_R2.FASTQ.gz
I have downloaded the genome for MCPyV as ...
1
vote
0
answers
57
views
RagTag patch error--"Tuple index out of range"
This question was also asked on GitHub
I'm trying to correct a long-read assembly with a short-read scaffold; I'm hoping to fill in the short gaps in the scaffold with the matching long-read sections. ...
- 31
2
votes
1
answer
33
views
Velvet Optimizer automatically changes to hash-length 31
I'm trying to use Velvet Optimizer for a De Novo Assembly; I set my hash-lengths to be between 55 and 69
...
- 67
2
votes
1
answer
30
views
How to find specific types of assemblies for specific species using entrez tools?
How to find specific types of assemblies for specific species using entrez tools?
Task: Trying to specifically find transcriptomes and associated cDNA data for a list of speices.
I can use this ...
- 107
1
vote
1
answer
50
views
How to subset an SRA file for a single chromosome?
I used prefetch to get the Pacbio reads of chicken from the SRA database. I want to align these reads against a reference genome, but not all the reads. I am only interested in a particular region on ...
2
votes
2
answers
123
views
Improving prokaryotic assembly with other contig/scaffold-level data?
I have what at first sight appears to be a high-quality MAG (~10 pieces, high completion%) that I built from a hybrid assembly (Illumina + Nanopore data) from a cyanobacterium.
Workflow:
Quality ...
- 787
2
votes
1
answer
61
views
Sanger sequencing annotation error
I am a student in a Cancer lab. Working with sanger is new to me. While analyzing a report we found an insertion that has not been reported in any databases so far, we were working on checking if the ...
- 21
3
votes
1
answer
477
views
How to solve Nextflow error: "Trace file already exists"?
When trying to run epi2me-labs/wf-artic, I get the following error:
...
- 377
2
votes
2
answers
45
views
Calling isoforms from long read data generated from partially degraded RNA
What will be the best tool to call isoforms from long read data generated from partially degraded RNA. By mistake we processed some samples with poor quality RNA to generate long read. Now we are ...
- 193
2
votes
1
answer
33
views
Using very closely related strains to increase coverage for short read, de novo assembly
Do you think it's possible to combine short read Illumina libraries (WGS) from multiple closely related eukaryotic microbial strains (e.g. libraries from a re-sequencing study, >99% ITS1 sequence) ...
- 31
1
vote
1
answer
37
views
How is an X chromosome encoded into a fasta string?
A human "X" chromosome has a centromere and two "identical" chromatids. If the chromatids are not identical, this fact is not assembled, correct?
The fasta string for a chromomsome ...
- 239
4
votes
1
answer
51
views
How to promote assemblies into genomes in NCBI?
Note: I've never submitted an assembly/genome to NCBI, so excuse if my perspective is flawed.
I'm working with Drosophila subobscura. (spring fruit fly)
I see here https://www.ncbi.nlm.nih.gov/data-...
- 239
1
vote
3
answers
71
views
How can I assemble my genome from raw files?
I've had my whole genome sequenced (at 30x average coverage) by a lab, and they have provided the raw files to me (BAM, FASTQ, and VCF).
How can I assemble it?
And does assembly provide any further ...
- 111
2
votes
1
answer
71
views
need bam file for pilon
I just ran an assembly on yeast genomes using Flye and I want to polish those assemblies with Pilon but it requires a sorted BAM file.
How do I make a BAM file of the resulting assembled.fasta?
- 494
2
votes
1
answer
52
views
What is the best way to process yeast genomes?
I have obtained several hundred raw, unassembled yeast genomes from NCBI and I am looking for advice on how to process the genomes for downstream analysis.
I have a reference genome (S288C) to use for ...
- 494
3
votes
1
answer
65
views
How does one distinguish nuclear DNA from mitochondrial DNA when doing WGS?
I'm interested in doing de-novo sequencing but also phylogenetic analysis. In particular, after de-novo sequencing and annotating the genome, I need to align the CO1 gene and the nuclear 28S rRNA gene ...
- 257
1
vote
1
answer
36
views
MetaQuast for assembling samples from complex communities
I'm working with whole genome metagenomic samples from human skin, and I'm using MEGAHIT for assembly and MetaQuast for evaluation. However, MetaQuast requires a list of reference genomes for the ...
- 13
1
vote
1
answer
253
views
DNA genome string reconstruction from k-mer
I have the following quiz question, but the Pattern1 for both (ACC|ATA) and (CGA|ACT) are unique (just grep for ...
- 143
2
votes
1
answer
32
views
How good does the assembly of an NCBI prokaryotic genome have to be in order to argue gene loss?
NCBI has several labels for assembly completeness - Complete, Scaffold, Chromosome and Contig. Complete would be a circularized genome (or linear, rarely)
For a Complete genome it's fairly ...
- 787
0
votes
2
answers
31
views
Assemble Anaeroplasma species genome from metagenomic PacBio data
I have a fasta file containing reads generated by PacBio HiFi whole genome sequencing of a feces sample from mouse.
I would like to use this dataset to generate an assembled circularized genome for an ...
- 1
2
votes
2
answers
42
views
How can I improve or otherwise investigate an unreliable genome tree?
Summary My genome tree doesn't agree with my gene trees and I get the feeling that my genome tree might be wrong, possibly due to long branch attraction, but I don't know how to check/fix it.
...
- 787
0
votes
1
answer
137
views
Why must a maximal non-branching path be a contig?
The following is from Bioinformatics Algorithms:
Fortunately, we can derive contigs from the de Bruijn graph. A path in a graph is called non-branching if in(v) = out(v) = 1 for each intermediate ...
- 127
6
votes
2
answers
49
views
How to display novel genome assemblies or uncommon genome assemblies using the UCSC Genome Browser?
I want to display E.coli BW25113 (GenBank: CP009273.1) strain in UCSC browser. This strain is not listed in http://microbes.ucsc.edu/ browser. How can I display E.coli BW25113 assembly in the browser?
- 534
10
votes
2
answers
196
views
Extract sequence context of high-degree nodes in assembly graphs
I often use metaSPAdes to assemble short reads from human microbiomes. My simplified understanding of short-read de Bruijn graph assemblers is that they fail where ambiguous paths cannot be resolved. ...
- 563
0
votes
0
answers
90
views
How to manually curate a genome assembly for sequence variation or error?
I have a PacBio HiFi assembly of 1.1 Gb from a heterozygous species. I have aligned this assembly against a reference genome which is around 0.9 Gb. I can see that there are quite a few INDELs, ...
- 75
0
votes
1
answer
59
views
How to improve a genome assembly using Dovetail and PacBio assembly?
I have more of a conceptual question. I have two genome assemblies from the same plant, one from Dovetail technology (~998 Gb) and another is PacBio HiFi assembly (~1.1 Gb). The Dovetail assembly is ...
- 75
0
votes
1
answer
63
views
Length of Contigs in Transcriptome and Whole Genome Assembly
Why are there shorter contigs from transcriptome assembly than from a whole genome assembly? I know the difference between transcriptome and genome, but don't really understand what contigs are in the ...
- 13
1
vote
1
answer
103
views
Comparing homozygosity of k-mer plots
Attached are two kmer plots from two closely related species. Is that safe to say that the one on the left has higher homozygosity than the one in the right k-mer plot, due to a low to almost flat ...
0
votes
1
answer
195
views
Understand this Kmer plot from Merqury?
The attached figure is generated based on Illumina reads from multiple individuals compared to genome assembly. Looks like there are a lot of kmers are reads only (grey colored). Also, a blue peak (2x ...
0
votes
2
answers
90
views
What are some ways to check metagenomic bin quality?
I am new to metagenomic binning. I've used CheckM in order to estimate completeness and contamination values, and most of my bins of interest appear to have good values. My workflow was pretty ...
- 787
0
votes
1
answer
72
views
BUSCO results apparently inconsistent
I have two assemblies obtained from two different softwares for the same run data. In one of them, I get a BUSCO score of 69,4%. In the other one, I get 68,5%. The first assembly covers a 99,7% ...
- 349
0
votes
1
answer
361
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 ...
- 31
1
vote
1
answer
57
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,568
0
votes
2
answers
127
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-...
- 280
5
votes
1
answer
83
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 ...
- 280
0
votes
2
answers
161
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
1
answer
66
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 ...
- 111
1
vote
1
answer
318
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 ...
- 151
0
votes
2
answers
109
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 ...
- 349
1
vote
1
answer
285
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 ...
- 73
-1
votes
2
answers
129
views
0
votes
3
answers
68
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 ...
1
vote
1
answer
87
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 ...
- 349
3
votes
2
answers
82
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/
...
- 349
0
votes
0
answers
173
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 ...
- 81
1
vote
1
answer
81
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 ...
- 349
1
vote
0
answers
17
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 <...
- 123
1
vote
2
answers
486
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 ...