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
1
vote
1
answer
64
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 ...
- 15
0
votes
2
answers
1k
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
0
votes
1
answer
149
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
2
answers
169
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 ...
- 199
3
votes
1
answer
141
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 ...
- 261
1
vote
1
answer
145
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 ...
- 349
1
vote
1
answer
750
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, ...
- 199
2
votes
1
answer
132
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
2
answers
295
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
...
- 349
1
vote
2
answers
159
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 ...
- 13
1
vote
1
answer
121
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 ...
- 185
0
votes
1
answer
127
views
Resources to learn genome assembly workflow for small genomes (like viruses)
I have sequencing data of a few samples of a DNA virus.
I'd like to learn de novo assembly of 'short read' data, to construct a scaffold and then count the abundance of each strain in the data.
I ...
- 169
1
vote
1
answer
95
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 ...
- 2,264
4
votes
1
answer
611
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 ...
- 41
2
votes
2
answers
120
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 ...
- 53
1
vote
2
answers
104
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 ...
- 349
1
vote
1
answer
127
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 ...
- 777
3
votes
1
answer
29
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
1
answer
164
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 ...
- 111
1
vote
0
answers
128
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 ...
- 31
0
votes
1
answer
316
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 <...
- 123
0
votes
1
answer
65
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
1
answer
37
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 ...
- 127
0
votes
1
answer
55
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
...
- 11
2
votes
1
answer
191
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.
...
- 349
2
votes
0
answers
305
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 ...
- 21
3
votes
2
answers
580
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 ...
- 777
0
votes
1
answer
143
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 ...
- 777
2
votes
1
answer
63
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 ...
- 21
0
votes
1
answer
74
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 ...
- 21
1
vote
4
answers
317
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,373
0
votes
1
answer
186
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 ...
- 777
2
votes
3
answers
229
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 ...
- 2,449
0
votes
1
answer
74
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? ...
- 1
0
votes
1
answer
563
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 ...
- 273
1
vote
1
answer
36
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 ...
- 113
2
votes
2
answers
37
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 ...
- 2,264
0
votes
0
answers
49
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,
...
- 109
0
votes
2
answers
859
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 ...
- 113
2
votes
1
answer
216
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 ...
- 21
4
votes
1
answer
69
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 ...
- 165
2
votes
1
answer
444
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 ...
- 709
6
votes
1
answer
118
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 ...
- 777
3
votes
1
answer
93
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 ...
- 709
4
votes
1
answer
88
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 ...
- 709
1
vote
1
answer
32
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 ...
- 439
0
votes
1
answer
780
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:
...
- 31
2
votes
1
answer
422
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,111
3
votes
2
answers
245
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 ...
- 145
5
votes
3
answers
4k
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 ...
- 357