Questions tagged [genome]

for questions related to whole genome analyses. For questions related to individual genes, use tag gene instead. Questions specific to format (like fasta) should be tagged with a format-specific tag.

Filter by
Sorted by
Tagged with
0
votes
0answers
37 views

STAR aligner multiple fastq files

I’m using STAR to align fastq files from SMART-seq2. I have raw data folder containing sub-folders with samples names the sub-folders each contain fastq file. How can I make a bash command in order ...
2
votes
0answers
22 views

Do RepeatModeler results contain functional domains?

The repeat families predicted by RepeatModeler contain known TEs and unknown ones. How do we know whether some of these may actually be a functional domain of a gene? It's predicted as repeats may ...
0
votes
0answers
8 views

Can we use reference strain as an out-group strain to root SNP based species tree?

I have generated a species tree based on whole-genome SNP using CSI_SNP_Phylogeny, wherein I used Reference strain as an outgroup to root my species tree (SNP fasta were downloaded and RaxML used to ...
0
votes
3answers
26 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 ...
2
votes
1answer
64 views

What tools can I use to look up my alleles, genotypes, or phenotypes from my sequenced DNA (WGS)?

(I am not a Bioinformatics expert, please forgive and educate me if I've used any wrong terms or assumptions here) I bought a "Whole Genome Sequence" kit, which gave me the following VCF ...
0
votes
0answers
10 views

Question: Combine BedTools Closest and Window Command

I am new to using Bedtools so please excuse my possible naive question. Let's say my file A is Snps and file B is genes. The Closest command with -k flat let one set how many genes that are closest to ...
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/ &...
0
votes
0answers
37 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
40 views

How to extract certain genes including non-coding regions?

I would like to look at the non-coding regions of some genes. I have around 1000 full genome assemblies and I was able to extract the nucleotide sequences for certain proteins with ...
1
vote
1answer
24 views

Mapping protein refseq to Gene ID

I have a protein refseq (eg, NP_000029). How can I get the corresponding gene ID (ag, APC) from NCBI using an R package?
1
vote
1answer
41 views

how to train a gene dataset with a nearest shrunken centroid classifier?

I have a data file named "geneexp.csv". the data contains information about gene expression of three different cell types (CD4 and CD8, CD19) I want to classify cells by performing the ...
1
vote
2answers
44 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 ...
2
votes
1answer
53 views

Merging several fq.gz files or R1 and R2 classes into a single one

I have the following files: ...
2
votes
2answers
84 views

How can I subset WGS data to the level of WES variants?

I would like to compare mutational signatures1 in patients from different studies, however some studies are based on exome seq (i.e. ~20,000 coding variants) and some are from whole genome seq (i.e. ~...
0
votes
1answer
19 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 ...
0
votes
1answer
90 views

kraken2 to OTU table creation

I did run kraken2 and I get two kinds of output How do I generate OTU table from these outputs for phyloseq usage?
0
votes
0answers
137 views

How to fix “expired or not present file ~/.gm_key!” error in braker?

I'm trying to run a braker with genome file only to see if it's working or not before running with all the datasets. But I'm getting this error again and again. ...
1
vote
2answers
43 views

Why does Coronavirus genome contain poly A region at the 3' end in GenBank? [duplicate]

If you look at the Coronavirus genome below, it contains poly-A region at the end. I know that the viral RNA contains a poly-A tail. But, the genome in GenBank is stored in cDNA format. That means the ...
1
vote
1answer
42 views

Get raw genome sequence

I'm currently working on testing some matching algorithms for strings. I would like to do some tests on raw genomic data as I expect different results from random strings given the lower entropy. How ...
1
vote
1answer
37 views

How do I download the mitochondrial haplogroup datasets for human genetics online?

I seem to have landed on the mitomap.org site, but I don't know what to make of it or what do with it / how to get the genomes onto my computer. It sounds like the genomes are stored in GenBank, but ...
1
vote
1answer
168 views

How to fix RepeatMasker.lib.nsq missing error when running RepeatModelor?

I have run RepeatModelor and I don't see the classified output file. I got this error. ...
0
votes
0answers
33 views

Why I get Singularity error when I try to run MAKER Genome Annotation on HPC?

I need to run MAKER genome annotation for my assembly file on HPC and I'm new to this. When I try to do that I get this error. ...
2
votes
1answer
50 views

Truble to run a multiprocessing kmer count script

Hi there I have this code: ...
1
vote
2answers
67 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 ...
0
votes
0answers
36 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
56 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
0answers
8 views

How do I include repeat purity, default slippage, default stutter, and minimum flanking (left and right) in Tandem Repeat Finder's output?

I am attempting to create a markerInfoFile for use in a program called popSTR (GitHub Documentation: https://github.com/DecodeGenetics/popSTR). The marker info file should contain information about ...
0
votes
1answer
1k views

Samtools Index: Chromosome Blocks not Continuous

I am working with short-read whole-genome sequences from the NCBI's SRA. I have aligned and sorted all of my short-read sequences and am attempting to index each sequence into .bai format using ...
0
votes
1answer
60 views

Seurat DE t.test

I am new to the Seurat package and was looking at the tutorials on the Seurat website and am just curious to know if there is any way to perform a t.test for the expression of a specific gene between ...
0
votes
0answers
15 views

Select for a specific proportion of a genomic feature

Currently have a bed file with exon coordinates for a list of about 20 genes, and I want to select the exons that fall within the first 3/4 or 2/3 of the coding region for each gene. Does anyone know ...
0
votes
0answers
37 views

Coding vs non-coding DNA length

I am trying to calculate the total exonic length (in bp) in order to see where the coding-noncoding ratio of roughly 1% to 99% comes from. We know all chromosomes total about 3 billion base pairs. So ...
0
votes
0answers
11 views

What is the origin of the 'source Gene ID' references given in the 'gene_presence_absence.csv' output of Roary?

I am learning to use Roary for preparing a pan genome for some lactobacillus strains. In the 'gene_presence_absence.csv' output of Roary (which i view in excel), a 'source Gene ID' is given for each ...
1
vote
1answer
42 views

BWA: Detecting Variation between Reference Genome and Short-Read Sequences

I need to identify all loci in the short-read sequence at which the number of microsatellite repeats (i.e. number of copies of "AA," "GTC," etc.) differ from the reference genome, ...
1
vote
1answer
31 views

How to use the DI-tector software for defective genomes analysis?

I am trying to use a free .py script from this site http://www.di-tector.cyame.eu/ while unfortunately there's neither manual nor any documentation presented. ...
1
vote
1answer
40 views

How can I convert codon coordinates to genomic position?

I am looking for a given mutation in IGV, which accepts coordinates in the form of chr<X>:<Y> Where X is the chromosome number and Y is the base ...
1
vote
1answer
22 views

Get a single number representing the contact probability between a pair of genomic loci using HiC data

Hello to the experts in analyzing chromosome structure data in HiC format, I have a very basic question. I have a HiC file (specifics are mentioned below), using which I would like to obtain contact ...
3
votes
0answers
51 views

What does it mean if aligned genomes are not of the same length?

I am aligning 4 genomes of Psuedomonas chlororaphis using progressive mauve software. Command used: ...
0
votes
0answers
27 views

Finding regions between ORFs

I'm doing work on mRNA and need to get predicted mRNA sequences from a genome. I'm working with this genome https://www.ncbi.nlm.nih.gov/nuccore/U00096.3 (E. coli k-12 mg1655). So have a genome ...
1
vote
2answers
89 views

how to use list of gene id to get cds sequence(cds fasta file have many annotation, only gene id: is same to query id)

i have a question when i want to extract cds sequence using gene id. but cds file is not just start with >gene is, it has many other annotation. the only same is star with gene: cds fasta: ...
6
votes
2answers
654 views

BERT Language Model and Gene Sequences - How Do I Relate Clusters of Sequences?

I hope you'll indulge a question from a computer scientist with limited bioinformatics knowledge. I've been working with the Google tool for language modeling called BERT. It's generally regarded ...
1
vote
1answer
24 views

How can I find the position of every mutation where the Allele Frequency is greater than X, in regards to a reference such as Hg19?

I have a Human genome sequenced in a BAM file (along with other files with the indels, snps, cnvs). I want to find every mutation with regards to the reference Human genome. However the majority of "...
1
vote
1answer
118 views

Getting sequences of one transcript per gene out of annotation and genome

I got a genome (data/genome/genome.fasta) and braker-based genome annotation (data/genome/annotation.gff3), now I would like to ...
3
votes
2answers
1k views

How can I classify the 3 clades(S, G, V) of the coronavirus without using protein data?

On GISAID they classified the coronavirus using 4 clades(S, G, V, Other). I downloaded around 1,000 complete genomes of the coronavirus from GISAID and I would like to classify each one as belonging ...
1
vote
0answers
35 views

Is it possible to determine the genomic context of a gene in a Whole Genome Shotgun project?

I performed several BLAST searches and identified several genes of interest in a Whole Genome Shotgun (wgs) project. I know that gene X is located on scaffold 1234, from nucleotide 1 - 2250, while ...
3
votes
1answer
57 views

Overlap in Paired-end Reads for Sequencing?

Regarding Sequencing: Do/can the paired-end reads have overlaps?
1
vote
1answer
50 views

Location of Reads in Sequencing

I have questions regarding sequencing: Are paired-end reads from different stands? What about single reads, are they coming from both strands?
20
votes
8answers
5k views

Since every human has a different DNA (different combinations of C, G, A, T) what does it mean to have the genome done? [closed]

I'm confused about the difference between genome and DNA. Is it correct to say that the same type of bacteria has the same DNA? But my understanding is that it is not correct to say that the same type ...
1
vote
1answer
38 views

Warning in fastqc

I am checking the quality control of my sequences using the Fastqc tool. For some steps, like "Per base GC content", etc., I received a warning. So, I was wondering whether I Should also take care of ...
1
vote
1answer
106 views

Per Base Sequence Content in fastqc

I have a question regarding "Per Base Sequence Content" plot for "fastqc": In the fastqc documentation, it is written: "In a random library you would expect that there would be little to no ...
2
votes
0answers
27 views

What Cellbase resources are 0-indexed, and which are 1-indexed?

The UCSC Genome Browser Team is clear about that the genome browser GUI indexing which is 1-based, closed interval. Genomic sequence retrieval via the UCSC REST API is 1-based, closed interval. All ...