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
1
vote
1answer
29 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, ...
1
vote
1answer
48 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-...
2
votes
2answers
137 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 ...
0
votes
1answer
45 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
86 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 ...
0
votes
0answers
34 views

Which of the Transcriptome assembly method is best for identifying novel lncRNAs?

I'm working with human samples and I'm trying to identify novel lncRNAs from tumor samples of Prostate cancer. I'm using reference based transcriptome assembly with ...
0
votes
0answers
11 views

Filtering out nonarthropod sequences using Blobtools2

I am using Blobtools2 and am working my way through the sample dataset that they provide before using my own. The filtering tutorial I am following is --> https://blobtoolkit.genomehubs.org/...
0
votes
0answers
46 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 ...
1
vote
1answer
25 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 ...
3
votes
1answer
94 views

Difference between genome assembly and genome sequence alignment to a reference to find structural variants

I'm trying to determine what the difference and benefits of genome assembly and genome sequence alignments are when trying to identify structural variants or transposons in populations. I've been ...
2
votes
2answers
66 views

How to fragment genomes into non-overlapping sequences of differing sizes?

I have bacterial chromosomes/plasmids that I have downloaded off refseq in a multi-fasta that I would like to fragment into non-overlapping fragments of between 1 kb and 15 kb. I have tried googling ...
1
vote
2answers
78 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 ...
1
vote
1answer
25 views

Comparing gene abundances between metagenomes

My workflow until now: Find fragments of a marker gene in unassembled metagenomes > download and assemble metagenomes > recover the gene neighborhood / gene set of interest Right now I have a ...
0
votes
0answers
2 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 ...
1
vote
1answer
71 views

Decontaminating RNA seqs for de novo transcriptome assembly and annotation of novel eukaryotes

I have raw paired-end RNA-seq reads for two novel eukaryotic species. Some background: the reads represent a copepod (arthropod) species each. The mRNA for each read set was obtained by extracting ...
0
votes
0answers
35 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
1answer
64 views

“megahit: command not found” after correct installation with miniconda

I want to use MEGAHIT (git: https://github.com/voutcn/megahit) to assemble metagenome shotgun sequencing reads (illumina). I install the package with the command <...
0
votes
1answer
41 views

Methods to predict one connected structure from two known separate structures

Currently, I am facing such a question: how to predict one connected protein structures from two separate known structures? It should be different from the protein-protein docking method, which is ...
2
votes
1answer
32 views

Assemble reads for a specific gene

I have a lot of unassembled sequencing data and I want to build phylogeny for some genes from these data. I can retrieved the sequene by mapping reads to a reference sequence but if the two species ...
0
votes
1answer
34 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
64 views

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

I guess that's a basic question but I can't find an answer. I am analyzing some SARS-CoV-2 sequencing runs. I often find read alignments like the one in the image. ...
2
votes
0answers
88 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
2answers
74 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 ...
0
votes
0answers
36 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 ...
1
vote
1answer
26 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
33 views

Optimum parameters for genome assembly via canu

I have been trying to assemble my genome via canu. Canu gave me about 200 contigs using default parameters. My genomesize is ~9Mbp. Can you recommend some parameters I can tweak to decrease the ...
1
vote
3answers
51 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 ...
0
votes
1answer
37 views

Can I get longer contigs by changing MEGAHIT settings?

My current settings are megahit -r input.fq --num-cpu-threads 32 --min-contig-len 300 --presets meta-large -o output. I picked the 'presets meta-large' because I am ...
1
vote
2answers
126 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 ...
0
votes
1answer
30 views

mapping and de novo assembly of plant without reference genome

I'm studying on a plant that has no reference genome and only has one scaffold assembly and one gff3 annotation file. Can I create an index with the same assembly and gff3 in STAR and do the mapping? ...
0
votes
1answer
97 views

nanopore - where to retrieve information from the basecaller used

Where can I find information about the basecaller used in a specific nanopore run (we have a Gridion machine)? I know already that we have Guppy v. 3.0.3, which I ...
1
vote
1answer
30 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
2answers
33 views

fixing chromosome with new contig

I want to fix a chromosome of around 50-100Mb (chr.fasta) for a region that is missing a few kilobases of data. I have the missing sequence with a window upstream ...
0
votes
0answers
44 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, ...
0
votes
2answers
159 views

Determining the quality of a publicly available genome assembly

I am looking for homologous genes in various protozoan species using BLAST. The genome sequences of these species are deposited in NCBI's WGS database. NCBI’s webpage includes the global statistics ...
2
votes
1answer
167 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 ...
4
votes
1answer
60 views

Authoritative source on human cytogenetic regions?

I am looking for a database that would keep track of human cytogenetic regions and genomic coordinates per genome assembly. I had expected the Genome Reference Consortium to have it, or Ensembl, or ...
2
votes
1answer
248 views

Consensus sequence from SAM or BAM file?

I am trying to perform reference-based assembly. Most of the tutorials teach how to create a bam file and view alignemnts in IGV or Tablet. But, I want a assembled genome sequence in fasta format. How ...
5
votes
1answer
84 views

Is there a way to assemble contigs starting from a specific sequence?

My work involves searching for marker genes/fragments in metagenomic databases (like the Sequence Read Archive). Once I find these sequences, I would like to know more about the neighboring genomic ...
3
votes
1answer
60 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
48 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 ...
1
vote
1answer
30 views

Is there a resource where I can donwload E. coli assemblies grouped by their pathotypes?

I will preface this by saying I'm not a microbiologist, so I apologise for any haziness and I'll be happy to expand. I would like to download the assemblies of different E. coli genomes and group ...
0
votes
1answer
264 views

Assembly by stringtie

I run this cmd ./stringtie G1_sorted.bam -B -o G1.gtf -G Triticum_aestivum.IWGSC.42.gtf -p 4 -C G1.refs.gtf -A G1.abund.tab Error is: ...
2
votes
1answer
334 views

Is nanopolish worth it since faster polishing software is available?

For Oxford Nanopore contigs produced by any long-read assembler has anyone performed any benchmarks to compare the polishing tools racon, ...
3
votes
2answers
127 views

Separation of mixed plasmid DNA sequences post whole-plasmid sequencing

Imagine a DNA sample containing a mixture of different intact plasmids. These samples are sequenced using either MiSeq or HiSeq sequencing. Would it possible to assemble these plasmids post-sequencing ...
5
votes
3answers
2k views

How to get fasta alignment file from SAM/BAM file?

I am not talking about consensus sequence, I know how to get consensus sequence using mpileup in samtools/bcftools. As I understand , SAM/BAM files are basically sequence alignment format so it's ...
7
votes
1answer
285 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 ...
2
votes
1answer
108 views

How is consensus alignment for OLC assembly usually implemented?

Background In C++ I want to implement a simple class that performs overlap layout consensus assembly but I can not figure out the most logical data structure to use for the consensus alignment step. ...
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 (...
2
votes
1answer
145 views

E.coli Sequencing & Analysis

I have been given the task of assembling a 'new' Ecoli genome and analysing the genes present etc. The Ecoli is a new strain, and has been taken and run on a Nextseq 500 in high-output mode with ...