Error characters in fasta file
How do I remove invalid characters E, F, I, L, P, Q from a fasta file?
I get these invalid characters while running muscle command for multiple sequence alignment.
alignment
muscle
fasta
• 4,963 views
•
link
updated
by
zx8754
•
written
by
atitparajuli2018 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Multiple sequence aligments - parallelisation
written by Lada •I have a folder with approx 3.000 fasta files. Each fasta FILE corresponds to one gene (orthogroup) and it contains multiple sequences (orthologues from multiple …
-
matching PDB and PFAM sequences for contact mapping
written by Evan •I am trying to generate a contact prediction from PFAM MSAs but I need to reliably map a given protein family sequence (a specific sequence …
-
Multiple Sequence Alignment using muscle
written by atitparajuli2018 •I tried to align 112626 sequences of alfalfa using muscle command, but I get segmentation fault 11 error which is mainly due to low running …
-
Multiple Sequence Alignment using muscle
written by atitparajuli2018 •Hi, I keep getting an error segmentation fault 11 while running muscle command to align multiple sequences from a fasta file containing 112626 ESTs. Is …
-
Error in normalizing lines of sequence in a fasta file using NormalizeFasta (Picard)
written by RahelehHi, I want to run VarScan fpfilter for 3 .hc.vcf files (somatic, germline, LOH), however when I run this command bam-readcount -q 1 -b 20 …
-
samtools mpileup error
written by teresayahu •I am using samtools function mpileup to get reads in sites. The command I used is: samtools mpileup -l hg19.position -f hg19.fa -Q 20 -q …
-
AA prediction from CDR3 DNA sequence
written by bioinformatics.cancerHi, I was trying to figure out how CDR3 DNA sequence in TCR sequencing is converted to the corresponding AA sequence. From a data set …
-
Remove unwanted characters from FASTA file
written by setaHi everybody, I have some nucleotide sequence (fasta format) that contains many sequences with invalid characters, like 'e', 'q', 'i' and 'l', there is also …
-
HG19 Translation discrepancy or something else?
written by jacobsen.jeremy •Why is there Arginine where there should be a Tryptophan? Translation from XPasy of rna sequence for uc001abw.1 from knownGeneTxMrna sequence: ``` ... Q T …
-
Weird Problem Installing Biopython on OS X Mavericks
written by jxiang15 •Hello, I'm trying to install Biopython on my Mavericks laptop. When I run "python setup.py build" I get the following error: running build running build_py …
It's amino acid code. Does it assume the input is DNA for some reason?
Yes it does. Can I clean these characters from the file?
No, you should make sure your input is DNA
Please can you help me with the muscle command line for that?
Please read the manual, especially section 3. Asking us to do your work for you is not good etiquette.
I did use the -seqtype nucleo option in the muscle command still I was getting the error message so I was asking for the code that does not give an error message as: I* ERROR * Invalid parameter -SeqType nucleo
You already have an answer - you can only align multiple sequences of a single type (protein/DNA) using muscle. If you're sure your sequences are of the same type, we can help you with any error message you're seeing.
That doesn't change the fact that your sequences aren't nucleotides.
Dear, In a nucleotide fasta file, you won't be able to get this type of error. Make sure you have DNA file. And for muscle command please try this simple one and check you are getting the desired result
Asaf has already stated the point you're trying to make - what is the value you're adding to the discussion?