Questions tagged [fasta]
To be used for questions specific tothe sequence file format `.fasta`. Please minimise usage if the question is more generally about sequence formats.
232
questions
191
votes
4
answers
123k
views
Why does the SARS-Cov2 coronavirus genome end in aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa (33 a's)?
The SARS-Cov2 coronavirus's genome was released, and is now available on Genbank. Looking at it...
...
- 1,725
42
votes
4
answers
58k
views
What is the difference between FASTA, FASTQ, and SAM file formats?
I'd like to learn the differences between 3 common formats such as FASTA, FASTQ and SAM. How they are different? Are there any benefits of using one over another?
Based on Wikipedia pages, I can't ...
- 1,313
34
votes
4
answers
9k
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:
...
- 443
31
votes
3
answers
23k
views
Uppercase vs lowercase letters in reference genome
I am using a reference genome for mm10 mouse downloaded from NCBI, and would like to understand in greater detail the difference between lowercase and uppercase letters, which make up roughly equal ...
- 2,133
28
votes
7
answers
9k
views
Read length distribution from FASTA file
I have a single ~10GB FASTA file generated from an Oxford Nanopore Technologies' MinION run, with >1M reads of mean length ~8Kb. How can I quickly and efficiently calculate the distribution of read ...
- 2,133
23
votes
4
answers
14k
views
What Ensembl genome version should I use for alignments? (e.g. toplevel.fa vs. primary_assembly.fa)
When you look at all the genome files available from Ensembl. You are presented with a bunch of options. Which one is the best to use/download?
You have a combination of choices.
First part options:
...
- 1,563
21
votes
10
answers
6k
views
What is the fastest way to calculate the number of unknown nucleotides in FASTA / FASTQ files?
I used to work with publicly available genomic references, where basic statistics are usually available and if they are not, you have to compute them only once so there is no reason to worry about ...
- 5,497
19
votes
7
answers
15k
views
16
votes
12
answers
9k
views
How to convert fasta file to tab delimited file
I have a fasta file like
>sample 1 gene 1
atgc
>sample 1 gene 2
atgc
>sample 2 gene 1
atgc
I want to get the following output, with one break between ...
- 945
15
votes
10
answers
3k
views
How to simulate NGS reads, controlling sequence coverage?
I have a FASTA file with 100+ sequences like this:
>Sequence1
GTGCCTATTGCTACTAAAA ...
>Sequence2
GCAATGCAAGGAAGTGATGGCGGAAATAGCGTTA
......
I also have a ...
- 2,739
13
votes
6
answers
392
views
Are there any databases of templates for common bioinformatic file formats?
I want some templates of different file formats that I can use to test my scripts and identify possible bugs in my code.
For example, consider nucleotide FASTA, a simple but often abused format, I ...
- 3,928
11
votes
1
answer
3k
views
Changing the record id in a FASTA file using BioPython
I have the following FASTA file, original.fasta:
>foo
GCTCACACATAGTTGATGCAGATGTTGAATTCACTATGAGGTGGGAGGATGTAGGGCCA
I need ...
- 496
11
votes
3
answers
6k
views
Converting a VCF into a FASTA given a reference with Python, R
I am interested in converting a VCF file into a FASTA file given a reference sequence with Python or R.
Samtools/BCFtools (Heng Li) provides a Perl script ...
- 1,681
10
votes
5
answers
20k
views
How to download FASTA sequences from NCBI using the terminal?
I have following accession numbers of the 10 chromosomes of Theobroma cacao genome.
...
- 331
10
votes
6
answers
7k
views
How do I find identical sequences in a FASTA file?
I want to create a database for a proteomics study. Therefore, the mapping from a given sequence to a protein ID has to be unique. I am wondering whether there is already a built-in function in ...
- 743
9
votes
5
answers
4k
views
How can longest isoforms (per gene) be extracted from a FASTA file?
Is there a convenient way to extract the longest isoforms from a transcriptome fasta file? I had found some scripts on biostars but none are functional and I'm having difficulty getting them to work.
...
- 301
9
votes
1
answer
6k
views
How to extract fasta from a blastdb
How to extract the sequence used to create a blast database. This is useful when you download a blastdb from somewhere else e.g. one of the databases provided by NCBI including the 16SMicrobial ...
- 322
8
votes
9
answers
22k
views
Remove/delete sequences by ID from multifasta
I have a fasta file like this:
>Id1
ATCCTT
>Id2
ATTTTCCC
>Id3
TTTCCCCAAAA
>Id4
CCCTTTAAA
I want to delete sequences that have the following IDs.
<...
- 375
8
votes
4
answers
1k
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 ...
- 1,681
8
votes
2
answers
3k
views
Subset FASTA file by species name
I have a problem: I've managed to download a massive fasta file of 1500 sequences, but now I want to split them into separate fasta files based on the genus.
EDIT
The fasta file looks like this:
...
- 303
8
votes
3
answers
147
views
How does one construct a cladogram of intraspecies relationships?
I have SNP data from several cultivars of rice which I have used to produce alignments, but I don't think that the usual models and algorithms used for generating phylogenetic trees are appropriate, ...
- 133
7
votes
5
answers
5k
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 ...
- 663
7
votes
4
answers
9k
views
How do I rename fasta headers?
I am trying to use the seqkit replace command to replace chromosome names in the format chr_I, chr_II, ... to I, II, .... I am ...
- 2,459
7
votes
1
answer
2k
views
Why is this makeblastdb command not working?
I am trying to make local databases for the ncbi-blast+ package (version 2.60). I am doing so for 4 T-cell receptor genes. 3 of the 4 (TRAV, TRAJ, TRBV) have worked fine, but I am having problems with ...
- 449
7
votes
1
answer
510
views
Split FASTQ and matching BAM into matching chunks
I am running a slow downstream analysis on a large set of nanopore reads (approx 3 million), and would like to split them into smaller chunks, run the analysis in massively parallel, and then ...
- 2,133
6
votes
3
answers
855
views
How to convert a binary matrix of gene presence or absence into a fasta sequence
I have a binary matrix of gene presence or absence which looks like: [roary output]
...
- 945
6
votes
2
answers
1k
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:
...
- 63
6
votes
4
answers
1k
views
Convert SRA to FastA
I'm trying to get the FastA files for some accessions (like NC_001416.1). I did not managed to find an FTP server or direct link to these files (I want to get it from command line with ...
- 221
6
votes
4
answers
612
views
multi-sequence alignment of samples with multiple contigs each
I want to perform a multi-sequence alignment on 12 samples that clustered based on cgMLST. Ultimate goal is to find out whether they differ by the presence of a few genes.
I have performed multi-...
- 283
6
votes
1
answer
1k
views
How to convert DNAbin to FASTA in R?
I am trying to convert DNAbin files to fasta format. Rationale: I want to use the fasta file to calculate non-synonymous/synonymous mutation rate.
my_dnabin1 is a ...
- 63
6
votes
1
answer
82
views
Detecting portions of human proteins with high degree of microbial similarity
I'm a newcomer to the world of bioinformatics, and in need of help solving a problem.
My goal is to take a list of human proteins, and identify segments (13-17aa in length) with a high degree of ...
- 421
6
votes
1
answer
338
views
Run gffread in multi-thread mode
Is there any possibility to run gffread in multi-thread mode? The answer seems to be 'no' from the manual (or gffread -h), as no multi-thread option is mentioned.
...
- 2,656
6
votes
1
answer
124
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 ...
- 859
5
votes
3
answers
2k
views
Removing duplicate FASTA sequences based on headers with Bash
I used the following command to remove duplicate FASTA sequences based on the header sequence:
...
- 51
5
votes
3
answers
2k
views
How to append numbers only on duplicates sequence names?
I have a reference database with contains 100s of sequences in fasta format. Some of these sequences have duplicate names like so:
...
- 945
5
votes
5
answers
9k
views
How to extract / convert gff3 CDS sequences to multifasta
I would like to extract all the CDS from a batch of genomes. I have found a perl script from BioStars but this does not seem to work for me. I would preferably like to have a script/ method which will ...
- 945
5
votes
2
answers
300
views
Rule Extraction from nnet results
I used a script in R language that uses nnet library to predict promoter bacteria and i would like to know how to extract rules from this neural network results.
As my input of the neural network i ...
- 161
5
votes
1
answer
2k
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
...
- 1,734
5
votes
3
answers
1k
views
finding unique headers in a fasta file using linux command line
I tried to use the following command
uniq -u reference.fasta >> reference_uniq.fasta
I'd like a count of the unique headers.
- 51
5
votes
1
answer
396
views
Can FASTA files have nucleotide and protein sequences within them; or must they only have 1 type?
Can FASTA files have nucleotide and protein sequences within them; or must they only have 1 type? For example, a FASTA file has 2 sequences. Can the first one encode amino acids while the second one ...
- 59
5
votes
4
answers
461
views
How to rename a file based on a directory name?
I have a lot of FASTA files that I recently scaffolded with RagTag. The output of RagTag includes a directory for each sample and the corresponding FASTA file written as "ragtag.scaffold.fasta&...
- 804
5
votes
4
answers
3k
views
How to map PDB chains to Uniprot IDs using API services
I have a lot of PDB IDs and I need to get uniprot fasta sequences of these PDB IDs special chains by API services. For example, imagine that I need to get fasta sequence of '1kf6' 'A' chain. The ...
- 767
5
votes
1
answer
723
views
Read and write FASTA files with more information than id and sequence
I am trying to read a fasta file, manipulate is in Python (using BioPython) and then write it back.
The format of my sequences is like:
...
5
votes
2
answers
1k
views
Fasta Sequence Identifier format?
I have been analyzing some virus DNA from the NCBI databases. The fasta sequences that I recieve from them have header lines that look like this:
...
- 399
5
votes
1
answer
892
views
Ancestry of the coronavirus 2019-nCov, WuHan city, China
In one of the answers to another question about the corona virus a link was given to this phylogenetic analysis of the virus.
Can somebody assist a non-bio type here? It seems to show that the ...
- 153
5
votes
1
answer
484
views
Can a data file in VCF format be converted into FASTA?
I'm considering purchasing the 'MyGenome' product by Veritas Genetics to analyze my genome for a project. I'd like the data to be in FASTA format, but Veritas only provides VCF data. Is it possible to ...
5
votes
3
answers
941
views
How can I only get the species name for fasta sequences from blast results?
I am trying to make a phylogenetic tree from sequences obtained with blast. I have files that contains hundreds of fasta protein sequences which are all named like ...
- 51
5
votes
2
answers
769
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 ...
- 9,442
5
votes
1
answer
176
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:
<...
5
votes
1
answer
165
views
Extracting sequences from FASTA beginning with common 5' end
I am trying to figure out the best way to extract sequences from a FASTA file which begin with a common 5' region of 43 nucleotides. Preferably, I would like to to allow for "fuzziness" in this region ...
- 53