Questions tagged [reference-genome]

A reference genome (also known as a reference assembly) is a digital nucleic acid sequence database, assembled by experts as a representative example of a species' set of genes.

Filter by
Sorted by
Tagged with
2
votes
1answer
27 views

How do researchers account for heterozygosity in genome assembly?

this is a very general question about diploid genome assembly. I am wondering how people deal with heterozygous region in a genome when assembling ? Do they pick one of the haplotype to be in the ...
0
votes
0answers
27 views

Do I have to include the genome?

Let's say I'm using SBOL to describe an experiment where I've transformed a cell line by adding some plasmids. Do I have to include the whole genome of cell line in my representation? That's a lot ...
2
votes
1answer
327 views

How to find novel transcripts using GFFcompare?

I am trying to find novel transcripts from an RNA-seq database. Based on the advice I got, it seemed that using Stringtie for transcript assembly is a good way to go, and it supports novel transcript ...
0
votes
0answers
49 views

How should I deal with segmental duplications when aligning NGS reads to a reference genome?

This is a follow-up of my other question. I have been having trouble calling variants in the human SMN1 and SMN2 genes, because the human genome has a large segmental duplication there and these two ...
4
votes
0answers
88 views

Old versions of the reference genome and dbSNP

For benchmark purposes I’m trying to find an old version of the human reference genome (pre-GRCh37) and a matching version of dbSNP. Unfortunately it appears as though older versions of dbSNP aren’t ...
0
votes
2answers
82 views

Downloading Fasta for each chromosome of a complete genome

Whatever I am googling I am not able to find fasta file of each chromosome for hs37d5.fa. hs37d5.fa whole fasta file however is here ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/technical/reference/...
0
votes
0answers
32 views

What Ensembl genome version of wheat should I use for alignments? [duplicate]

There are many genome files available from Ensembl. Which one is the best to use/download for wheat RNA seq differential gene expression analysis? Also, tell me which tool or aligner or mapper to use? ...
4
votes
1answer
105 views

At what stage of a transcriptome assembly is it better to perform read contaminant filter?

I'm trying to assemble a bivalve transcriptome. Since bivalves are filter feeders, their transcriptomes tend to be highly contaminated by bacteria, algae and whatnot. Since I pooled several ...
2
votes
1answer
136 views

How to force pilon to correct a reference with low coverage mapped reads?

I have mapped some reads to a reference and want to change SNPs in the reference to SNPs found in the reads, even at low coverage. For example, here is a screenshot of one reference in IGV that I ...
0
votes
1answer
302 views

STAR Indexing Diference for Small and Large Genome File Output

I have a quick doubt on the output of the Genome Indexing, I have used the STAR program along with genome .fasta file and GFF file. Genome size is 3GB, here is the file output ...
1
vote
3answers
128 views

Genome scaffolding

I assembled a virus genome using Ray and got around 5000 contigs. Now I want to scaffold those contigs to get a full genome (I am aware of the fact that there might be some gaps). I found one tool ...
2
votes
1answer
28 views

GRCm37-designed exon target enrichment, which reference to use?

I have exomes from 24 individual mice. The exomes are the product of Roche's SeqCap EZ HyperCap target enrichment kit. I see that the mouse exome design comes from mm9/NCBI37, the previous major ...
2
votes
1answer
402 views

What process and input data is required for a cellranger reference transcriptome?

I'm analysing single-cell RNA-Seq data using the 10X Genomics cellranger platform. While they provide reference data for Mouse and Humans, other species require a ...
6
votes
2answers
135 views

Number of reference sequences in a SAM file

From a single record, I can get the reference sequence ID from the field rID, but how can I get the total number of different reference sequences stored in the whole SAM file? It can be simple as just ...
0
votes
1answer
79 views

Does the galGal5 chicken assembly have a chromosome 29?

The chromosome sizes at UCSC don't seem to contain chr29: ftp://hgdownload.soe.ucsc.edu/goldenPath/galGal5/bigZips/galGal5.chrom.sizes It has a chr28 and a chr30. Am I missing something or is there ...
4
votes
1answer
277 views

How to get results from Homo.sapiens package in bioconductor for a specific reference

I want to use the Homo.sapiens package in Bioconductor to retrieve the chromosome location start and end for each gene symbol in a specific reference (e.g. hg19 or hg38). Right now, I am using the ...
3
votes
2answers
143 views

Two variants associated with same chromosomal position in genotype data for one individual

I'm trying to play around with my family's raw genetic data from Ancestry.com so that I can do some genetic analyses, and I noticed in the AncestryDNA raw data files that some regions of the genome ...
1
vote
1answer
256 views

How can I subset a reference based on only the first chromosome?

I have the GATK FASTA reference for hg19. I would like to subset this based on the first chromosome. Normally what I do is open the FASTA as a text file, and then delete all other chromosomes except ...
1
vote
1answer
94 views

population admixture

I am new to Admixture and I would like to use fastNGSadmix to look for admixture in one individual (modern). Admixture is supposed to look for proportions of ancestral populations in a sample (either ...
2
votes
1answer
186 views

Java error when launching Pilon

I am trying to use Pilon to improve a reference based on some Illumina data I have got. So, I aligned the Illumina reads to the reference using bwa. Then, I want to use Pilon to improve the reference ...
5
votes
2answers
110 views

Creating a new reference genome B by changing genome A with mismatches from BAM file without long reads

This question is slightly related to this one: Improve a reference genome with sequencing data Only in my case, I have a starting genome reference A, of a short phage genome, and reads from a related ...
6
votes
1answer
89 views

Bacterial genome annotation of a clinical isolate strain?

So I'm basically so new to this that I'm just trying to find out what tools, methods, and keywords I should go look up by myself. I have a unique strain of a bacteria. I was given RNAseq data for ...
4
votes
2answers
359 views

Is it possible to do alignment within Python? Check variants against reference?

I'm currently adding a few SNPs randomly into a FASTA within python using BioPython. In the following example from BioPython, I add an SNV at location "5" http://biopython.org/DIST/docs/api/Bio.Seq....
4
votes
1answer
205 views

Given a Genomic Ranges of SNPs, how to inject these SNPs in genome via BSGenome?

Let's say I have the genome hg19 loaded into R via BSGenome ...
8
votes
4answers
793 views

How to manipulate a reference FASTA or bam to include variants from a VCF?

I have some software which takes fastas as the input. I need to include SNVs and InDels from a VCF into the reference hg38 and then use this. The problem is, I don't know of an algorithmically sound ...
2
votes
1answer
101 views

Is there a way to retrieve ENSEMBL IDs from a search query?

I would like to be able to batch download FASTA files from ENSEMBL. I normally would use the API to download them from the ENSEMBL gene IDs. However, the IDs aren't ones I know, only ones I would get ...
6
votes
5answers
632 views

What is the standard way to work with a diploid reference genome? Complementary strands?

At the moment, the standard reference genomes (e.g. hg19, hg38) are haploid genomes. We know that the human genome is diploid. Naturally, the latter would be the respectively correct representation of ...
6
votes
1answer
174 views

Can a customized GRCh38 .gtf file be used with any of the GRCh38 released patches?

I got a customized GRCh38.79 .gtf file (modified to have no MT genes) and I need to create a reference genome out of it (for 10xGenomics CellRanger pipeline). I suspect that the .79 part is the ...
8
votes
1answer
673 views

What are all the reference files produced by bwa index, and are these dependent upon whether the reference is zipped?

I have indexed a gzipped reference with bwa: bwa index reference.fa.gz, which produces a series of other files ...
3
votes
3answers
75 views

Which reference to use for read mapping for popular model organisms

What is the "best" assembly for the popular model organisms: human (GRCh37 and GRCh38 are obvious, I'd pick whatever bwakit uses) mouse (GRCm37/GRCm38, OK) but what about non-human/mouse ones? ...
11
votes
2answers
4k views

Downloading a reference Genome for Bowtie2

How do I download a reference genome that I can use with bowtie2? Specifically HG19. On UCSC there are a lot of file options.
16
votes
3answers
2k views

Convert a BAM file from one reference to another?

I have a set of BAM files that are aligned using the NCBI GRCh37 human genome reference (with the chromosome names as NC_000001.10) but I want to analyze it using a BED file that has the UCSC hg19 ...