Questions tagged [sequence-analysis]

Umbrella term for understanding any given nucleic acid sequence code, either as a locus, loci or genome with itself or as a comparison to comparable nucleic acid sequence code either intraspecifically or interspecifically

Filter by
Sorted by
Tagged with
3
votes
1answer
66 views

Overlap in Paired-end Reads for Sequencing?

Regarding Sequencing: Do/can the paired-end reads have overlaps?
1
vote
1answer
52 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?
1
vote
1answer
42 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
52 views

SRA run data access

Trying to download run data of Run accession SRR2155174. Here's the link: https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR2155174 Please go to the link and go to the "Data Access" tab. Hereunder ...
0
votes
1answer
54 views

NCBI SRA database sample: control vs test

I was trying to download some data from NCBI SRA (SRA059451). There are 27 samples available for SRA059451. But i am unable to understand which samples are 'control' and 'test' samples. Please help me ...
2
votes
2answers
54 views

Tools or softwares for clustering evaluation

I am using a few clustering algorithms as well as my own tool for protein sequences. Also, I have benchmarked clustering results (i.e. ground truth) to compare with the results. Is there any tool/...
1
vote
2answers
45 views

Download sequences of isoforms

I want to collect all isoforms of all genes (Fasta/Fastq file of nucleotide and protein sequences) Arabidopsis Col-0. I am wondering if there is a straightforward way to download the file from any ...
2
votes
1answer
122 views

Any user friendly way to find rare mutations in whole genome raw?

Is there any user friendly way to find rare mutations in the individual human whole genome sequencing raw data? (from Dante, 30x coverage). To be more specific, I want to find mutations from this ...
1
vote
1answer
72 views

working with cutadapt

I'm working with ion torrent data where I apply the program cutadapt. I'm analyzing ITS seq data, as well as Matk. When I'm using cutadapt, I search the information for this genes, and built new ...
2
votes
1answer
492 views

How to convert VCF format file into Fasta(.fa) with the Reference data

I need a Software or Python Program for Converting VCF(.vcf) to FASTA(.fa) format with the help of reference (.fa)
0
votes
1answer
79 views

How to reproduce this specific sequence logo from 10 000 human mRNAs?

this is A sequence logo showing the most conserved bases around the initiation codon from 10000 human mRNAs. I would like to reproduce this sequence logo. How can I get the sequences? which database ...
0
votes
1answer
329 views

Is there a convention or standard of the highest bits in a sequence logo? [duplicate]

In bioinformatics, a sequence logo is a graphical representation of the sequence conservation of nucleotides (in a strand of DNA/RNA) or amino acids (in protein sequences). wiki says The total ...
1
vote
1answer
831 views

How to compute the Shannon entropy for a strand of DNA?

I'm confused by the computation of sequence logo. Wikipedia gives a process about this without a concrete example. Let's just consider DNA, so there are 4 bases (nucleic acids). The following data ...
2
votes
1answer
213 views

Why do ten rows (Figure_1) correspond to 2 bits (Figure_2) in a sequence logo?

Following this question, I'm confused with the computation of sequence logo Following data comes from the book "Machine Learning - A Probabilistic Perspective (Figure_1)" here is the corresponding ...
1
vote
1answer
376 views

What is "aligned sequences" and "consensus sequence" in the context of sequence logo? How to compute these?

In bioinformatics, a sequence logo is a graphical representation of the sequence conservation of nucleotides (in a strand of DNA/RNA) or amino acids (in protein sequences). so, What is "aligned ...
5
votes
2answers
435 views

How can I locate duplicated regions in a sequence?

I am facing an issue when trying to align short reads against a region in human chr5. The two Sensory Motor Neuron genes, (SMN1 and SMN2) are almost 100% identical and this causes the aligner to fail ...
0
votes
1answer
79 views

How to use extractfeat?

I am working on a script on bash linux to get the CDS (coding sequence) of a gene using extractfeat by EMBOSS. It gives me the error: Warning: No sequences written to output file I am unsure if my ...
1
vote
3answers
286 views

Deep learning RNA sequences

Currently I'm working on a project, which combines deep learning with RNA sequences. I'll try to predict pseudotorsion angles [1] from raw rna sequence. The ideas is to train a neural network with raw ...
4
votes
1answer
90 views

Generating DNA sequences with constraints

I would like some advice on potential strategies to address the following problem. I want to write a program that will generate DNA sequences that are optimized on two constraints based on an input ...
3
votes
1answer
211 views

Protein SCA and DCA

Could you please suggest suitable R packages for statistical coupling and direct coupling analysis of protein sequences? I know there are Python packages for those: SCA and DCA. I was looking for ...
2
votes
3answers
1k views

CDS length for each human gene

Does anyone know where and how could I download a list of all human genes and the length of the coding sequence for each gene? Is it possible to do this on the NCBI site, ensembl?
1
vote
0answers
37 views

Building cancer gene database

I am trying to build a database of different cancer genes i.e lung cancer, throat cancer etc. Searched through NCBI website by searching for nucleotide sequences, but I am unable to find appropriate ...
8
votes
3answers
3k views

What is a quick way to find the reverse complement in bash

I have a DNA sequence for which I would like to quickly find the reverse complement. Is there a quick way of doing this on the bash command line using only GNU tools?
1
vote
1answer
82 views

Genome annotation of a Bacillus strain

I finished the sequencing of the genome of a biocontrol activity strain Bacillus amyloliquefaciens by Illumina hiseq. I have already gone through some papers but there were no details of processing in ...
1
vote
1answer
56 views

How to allow "half windows" in biopython?

I am using Biopython to examine protein sequences including the Kyte and Doolittle hydrophobicity. However, for short sequences, this returns an empty value; the value is not allowed if the window ...
0
votes
1answer
60 views

How to find the position-specific weight matrix (PSWM) for a list of genes having different sizes?

Wiki has the algorithm for finding position-specific weight matrix (PSWM) for a list of nucleotide sequences having an equal number of nucleotide. How do we find the PSSM for a list of nucleotide ...
1
vote
1answer
44 views

Download program for ΔG for TM insertion

I am trying to calculate the ΔG of TM insertion for many protein sequences. Has anyone come across the source code or a programmatic version for the von Heijne ΔG for TM insertion? The web server is ...
1
vote
2answers
60 views

Retrieve RNA sequencing data for human p53 colon cancer cell lines

I want to download the next generation RNA sequencing data in different cell lines. I want the data (alongwith normal) for p53 wild type, p53 knock-out and p53 null cell lines (human colon cancer). ...
3
votes
2answers
4k views

Hardware Requirements (specs) for Bioinformatics-dedicated desktop [closed]

I know this is a somewhat general or vague question, but I’m interested in your opinions. I must build a desktop for general bioinformatics activities with human genomes. I will work with Python and R ...
3
votes
1answer
69 views

De novo motif discovery in protein sequences

I am trying to build a features matrix to be used for Random Forest based classification. I'd like to add, as features, short motifs which are common to all the protein sequences belonging to a ...
4
votes
1answer
65 views

Codon usage analysis for whole genomes

I am new to bioinformatics. So if these questions seem you to a bit childish please forgive me. I have two queries. I am intending to perform a codon usage analysis followed by correspondence ...
4
votes
1answer
140 views

Tool to show DNA sequence and allowing upload of own graph data

Background We want to be able to load (or request) data for a genome including the sequence and gene annotation (bacteria). Then, we want to load our own annotation which should be displayed as a ...
0
votes
1answer
274 views

Fastqc- Per Base Sequence Quality

I am trying to figure out how to interpret the "Per Base Sequence Quality"? What does Position in reads (bp) mean? Also in order to draw this box-plot graphics, more than one quality scores are needed ...
0
votes
2answers
949 views

Python API for working with ENSEMBL genomes

Dear bioinformaticians, I could find two python APIs to mine ENSEMBL genome databases. There is a pyensembl which is more popular and then there is ...
3
votes
3answers
162 views

Effect of mutation in DNA sequence on transcription factor binding sites

How much does a single mutation/alteration of a nucleotide affect the presence of a transcription factor binding site (TFBS)? I am from computer science background(Obviously). I want to make a ...
12
votes
1answer
4k views

Better aligner than bowtie2?

Bowtie2 is probably the most widely used aligner because of it's speed. Burrow-wheeler (BW) algorithms (including bwa) tend to ...
3
votes
1answer
48 views

Is there any tools that can generating FBA model from genome sequences or proteome sequences?

Suppose I have all sequences of genes or proteins of certain species, are there any free tools I can use to generate a flux balance analysis model out of them? What I was thinking is more like to ...
4
votes
2answers
2k views

What is cellranger doing in comparison to other methods?

I've recently started working with the 10X-Genomics platform with Illumina (MiSeq and HiSeq) for single-cell RNA-Seq. I've been recommended the "cellranger" (version 2.1.0) which I understand handles ...
1
vote
1answer
139 views

How to predict the distance between two residues in a protein sequence

I have some protein sequences and I need to predict the physical distance between each two residues in the protein when folded correctly in 3 dimensions; I need to predict this distances just by ...
2
votes
1answer
134 views

SNP data Download Help

I am new to the field of Bioinformatics. Single Nucleotide Polymorphism or SNP data are available on the the Internet to download.SNP related data are valuable for variety of reasons ranging from ...
0
votes
1answer
47 views

Error while Importing fastq sequences to jMHC software

I am trying to analyze my .fastq files (MiSeq exon reads) with the jMHC software, but output files I have obtained have no results. I think I might be doing something wrong during the first step: "...
8
votes
3answers
345 views

Books on bioinformatics algorithms

I'm looking for a book about bioinformatics algorithms, such as alignment, BLAST search, and variant calling. I'm hoping reading about this subject will give me a deeper understanding of the ...
2
votes
1answer
558 views

Combine Fastq by writeFastq is not working properly

Currently I have 4 fastq file for each patient as WES data. For WES pipeline, I want to combine those fastq files. I am using writeFastq from ShortRead package in Bioconductor.But I think this ...
1
vote
0answers
25 views

How to use DCA to deduce the contacts between 2 proteins ensembles?

I would like to use DCA (direct coupling analysis) to analyze the interface between 2 protein ensembles (how the 2 ensembles interacts). Do you have access to a tutorial on this matter? A source code? ...
0
votes
1answer
177 views

What is the --rev_comp_mapping_barcode parameter in the QIIME1 script extract_barcodes.py?

I'm using the QIIME1 scripts extract_barcodes.py to extract barcodes from a dual-indexed MiSeq library, and then demultiplex my reads using split_libraries_fastq.py. I have been having trouble ...
5
votes
7answers
2k views

calculating nucleotide frequency per column

I have some sequences shown below ...
1
vote
0answers
50 views

How to rebuild the workflow for mining in oomynete genomes in UGENE

I want to recreate the UGENE workflow proposed in Convergence of retrotransposons in oomycetes and plants. After being able to understand the data and hmmsearch Here and Here and getting most of the ...
2
votes
1answer
243 views

Mini-bulk RNA-seq and scRNA-seq analysis differences?

We have RNA-seq data from libraries prepared using the Smartseq2 single-cell protocol on 500 cells (mini-bulk) / library. The complexity is much better than with single-cells (~14k genes for 1.5M ...
3
votes
2answers
67 views

Identifying relevant SNPs from a list

I have a list of all SNPs (a list of RsIDs) that fall in the coding sequence of several thousand human genes, and I'm looking for a convenient method to come up with a list of known SNPs that alter an ...
8
votes
3answers
949 views

Least present k-mers in the human genome

What are the least present k-mers in the human genome at different sizes? Starting with k=4 and going up in size until k=10, what are the k-mers least seen (or not at all) in the human genome? I am ...