Questions tagged [fasta]

should be used for questions specific to sequence file format `.fasta`. should **not** be used if you just used fasta file, or if the question is more general about sequence formats.

Filter by
Sorted by
Tagged with
-1
votes
0answers
17 views

How to identify where the mutation occurs frequently

I currently have 2 fasta files with 139 sequences. One that is already aligned and one that is not. Each file is Sars-cov2 sequence file from ncbi. I was wondering what I can do to find locations ...
0
votes
0answers
19 views

How to find the region with frequent mutation?

I am very new to bioinformatics and I was wondering if I can get some help regarding finding region that has high mutation rate. Currently, I have 2 fasta files, one that is aligned and one regular ...
0
votes
0answers
43 views

How can I understand FASTA sequence of Coronavirus? [closed]

In GenBank it is mentioned that the COVID-19 virus has an RNA genome, but here in the community someone wrote that it is cDNA and not RNA. For translation RNA must contain a start codon and poly A ...
1
vote
1answer
21 views

Cutting the sequence into several sequences with the information of a dataframe

I have a fasta file with several 120-concatenation protein sequences. I also have a data frame with 120 names of proteins in column 1 and their length in column 2. Using this dataframe I want to ...
1
vote
1answer
39 views

Identifying proteins in a concatenated protein sequences [closed]

I have a 120-concatenated protein sequence in a fasta file. I would like to split by proteins. How can I identify each protein?
1
vote
3answers
48 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
0answers
14 views

Retrieve FASTA sequences using sequence ids in cdbfasta

I'm using cdbfasta command to retrieve sequences using ids, all works well but I find some sequences repeated in my output and the command does not retrieve all the sequences I want.
1
vote
3answers
82 views

Edit all the fasta headers using awk

kindly your comment highly appreciated have more than 4k header sequences look like: ...
21
votes
5answers
7k views

Why does the FASTA sequence for coronavirus look like DNA, not RNA?

I'm looking at a genome sequence for 2019-nCoV on NCBI. The FASTA sequence looks like this: ...
0
votes
0answers
103 views

Using BaseRecalibrator in GATK4

GATK4's BaseRecalibrator uses a list of known variants to adjust the base quality scores in a BAM file. I would like to visualize the pre and post recalibration ...
149
votes
4answers
80k views

Why does the Wuhan coronavirus genome end in aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa (33 a's)?

The Wuhan coronavirus's genome was released, and is now available on Genbank. Looking at it... ...
0
votes
1answer
50 views

Interpreting MultiQC plot of GATK BaseRecalibration data?

I'm using GATK4 to build a somatic variant calling pipeline. The pipeline uses MultiQC to aggregate quality control data, and one of the QC measures reported is base quality score recalibration from ...
4
votes
4answers
103 views

Convert FASTA to FASTQ with dummy quality scores

I have a FASTA file which I would like to convert into FASTQ format as the tool I want to use my data in requires it in FASTQ format. So dummy quality scores are fine. Note: I am not using any ...
3
votes
1answer
39 views

How to generate two fasta files from samtools phase output?

samtools phase is an easy-to-use tool to phase variants from a bam file and a reference. For a sequence that has two haplotypes, like this: ...
0
votes
0answers
14 views

Obtaining HGDP project data in fasta format

I need to obtain sample data from modern humans in fasta format. I just need some megabytes of data from every individual. I actually use a script that obtains the cram file from here (ftp.1000genomes....
0
votes
1answer
19 views

SRA Toolkit and lebanese data

I am trying to extract data from this: https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR6245218 I installed SRA Toolkit, downloaded the SRR6245218 file and executed this: ...
2
votes
1answer
35 views

How to extract hypervariable region of mitochondrial dna from its fasta file?

Recently I've done some analysis on Human mitochondrial DNA and now I want to run the same analysis on only hypervariable regions of mitochondrial DNA to see whether the changes that I have observed ...
1
vote
1answer
29 views

Bedtools get fasta and ORF from a blastX run

Hello I made a blastX research (A query genome in nucleotide format translated into protein in 6 reading frames against a protein dabatase) And here is a head of ...
3
votes
5answers
865 views

How to download FASTA sequences from NCBI using the terminal?

I have following accession numbers of the 10 chromosomes of Theobroma cacao genome. ...
0
votes
0answers
31 views

Extract secondary structures and sequences from dssp files

I have a very messy problem and I really need your help. I have a lot of dssp files in a folder. Every dssp file contains informations about proteins. Also, I have a separate file containing two ...
1
vote
1answer
134 views

Filtering fasta files by ID

I have this code to filter some IDs from a fasta file: ...
0
votes
1answer
81 views

How to edit the headers of multiple fasta files from multiple folders? (recursively)

This is an expansion of my previous question, How to edit the headers of multiple fasta files from multiple folders? My directories are organized as follows: one main directory, in which I have ...
2
votes
1answer
64 views

How to edit the headers of multiple fasta files from multiple folders?

My directories are organized as follows: one main directory, in which I have multiple directories that end with a number ranging from 314 to 727, followed by .3 . For example, 'mgm4761314.3'. Within ...
2
votes
1answer
79 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)
1
vote
2answers
68 views

MiXCR: only create a single export file for all clonotypes

I am using following command from MiXCR to both align, assemble and export my input files: ...
3
votes
2answers
227 views

How to find common sequences among 6 multi-fasta files

I have 6 multi-fasta files, every of them contains ca 1500 sequences like that: ...
2
votes
1answer
123 views

Concatenate multiple fasta files into one file per unique entry

I have eight fasta files, one per gene, with 106 fasta entries per file, like this: gene1.fasta: ...
3
votes
1answer
55 views

For phylogenetic tree construction from core-genome which one is preferable: amino-acid based MSA or nucleotide based MSA?

The genomes are from same species. Is it true that, in phylogenetic tree constructed from amino-acid based MSA (multiple sequence alignment) some information are lost, so for phylogenetic ...
3
votes
1answer
168 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
0answers
22 views

return cds from a protein fasta based on gene ID in unix?

I'm trying to retrieve cds sequences for given genes. Which requires taking a list of genes gene27 gene28 Retrieving their Gene ID (LOC##) and using faa to ...
0
votes
1answer
63 views

Subset a protein fasta based on sequence length?

I have a strange protein fasta file. There are several entries for the same gene and I need to extract the entry with the longest width. ...
2
votes
1answer
45 views

how to merge two fasta based on full fasta header and remove the duplicate

I am wondering if anyone can help to merge two fasta file based on the full header as initially gene ID may be similar in some of the sequences and to remove the duplicate entry based on full fasta ...
0
votes
1answer
41 views

Function to find ATP binding sites

I'm currently doing a python project where I need to write code that extracts the ATP binding site from a FASTA file using a regular expression. I've tried the following code but it does not work. It ...
2
votes
2answers
175 views

Issue with visualising cladogram/phylogenetic tree with multiple sequence alignment data in R?

I would like to visualize tree with multiple sequence alignment. My din.newick file is shown below, ...
2
votes
3answers
392 views

How to reverse complement the DNA sequences for given inverse/reverse coordinates?

I have the series of coordinates in id.txt file, whose coordinates sequences are in genome.fasta file. The coordinates of id.txt ...
0
votes
0answers
63 views

Reconstructing reference FASTA from CRAM file

I would like to reconstruct (or download) the original FASTA reference file that was used to create a CRAM file in my possession. I realize that I can download the contigs specified in the CRAM header ...
2
votes
1answer
55 views

Alignment for predicting DNA hybridization?

I am currently working on a Computer Science project where we are trying to build a large set of orthogonal single-stranded DNA sequences. The goal would be to ensure that when put in solution, the ...
5
votes
1answer
81 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 ...
0
votes
2answers
107 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/...
4
votes
1answer
945 views

Splitting fasta file into smaller files based on header pattern

I have to split this fasta files into smaller files and write them into individual files my files ...
3
votes
3answers
187 views

GTF from consensus sequence

I am new to using bioinformatic tools and I was hoping this community could help clear something up for me. I need to generate a gtf file. My data are a set of complete HA genes for influenza B ...
1
vote
2answers
40 views

Get & Annotation of isomiR (isoform of miRNA) profile from RNA-seq

I should get and analyse isomiR (isoform of miRNA) profile from RNA-seq with FASTA format. The best tool for this is "DeAnnIso". But here is a problem with this tool, it only allows you to upload ...
2
votes
4answers
1k views

Can a gff file be converted to a fasta file?

I downloaded an annotated genome file in gff format here. I would like to use it for proteomics. Though I need it in fasta format. Is there any tool that converts gff to annotated fasta? I see the ...
6
votes
2answers
418 views

How do I re-name the headers of my Fasta file?

I appologize for asking this, but I really am really bad with regex... Can someone help me transform the headers of my fasta files from this: ...
3
votes
1answer
394 views

scRNA-seq, 10x cellranger pipelines

I'm starting to do sc-rnaseq using 10x cellranger pipelines, and i add TdTomato sequence to mouse reference genome and add an entry in the gtf. My code is But when I makref, it remindered that: ...
5
votes
1answer
132 views

How to find all variable-length seqs with an exact 5' and 3' match in a FASTA file

Context I am interested in finding all of the promotors specific to a particular sigma factor. I have identified the -35 and -10 sites from the literature, bold denotes -10, -35, binding sites: <...
2
votes
1answer
120 views

Question about research in bioinformatics for Computer Science student [closed]

I am a computer science student, who was interested in taking research in the field of bioinformatics. The first research idea I had was to make a model that could predict the possibility of type 2 ...
2
votes
5answers
2k views

How to create a .bed file from .fasta? [duplicate]

I have some problems with creating a .bed file for hg19 which I will be able to load on IGV. The fasta file contains rows of this form: ...
4
votes
2answers
139 views

How to align output of grep --color=always? (To QC fasta/fastq files)

Grepping out short sequences from a fasta or fastq file is a really useful way to look at sequencing data. Using the option --color=always makes this even more ...
1
vote
2answers
101 views

convert tRNA coordinates Result to Fasta format

I believe many of the experts here have gone through this issue, I have run tRNAscan-SE on my genome to predict all the tRNAs, and have successfully gotten all the tRNA coordinates like start, end ...