Hello All - I am looking for help on how to solve a re-occuring problem. I have a file with certain sequences in it that need to be removed. The sequences are always different but the fix is always the same remove those sequences and leave the rest. Another team ID's the bad sequences and then I... (3 Replies)
Hey,
I've been trying to break a massive fasta formatted file into files containing each gene separately. Could anyone help me? I've tried to use the following code but i've recieved errors every time:
for i in *.rtf.out
do
awk '/^>/{f=++d".fasta"} {print > $i.out}' $i
done (1 Reply)
Hi,
I am having a file of dna sequences in fasta format which look like this:
>admin_1_45
atatagcaga
>admin_1_46
atatagcagaatatatat
with many such thousands of sequences in a single file. I want to the replace the accession Id "admin_1_45" similarly in following sequences to... (5 Replies)
I have two files. File1 is shown below.
>153L:B|PDBID|CHAIN|SEQUENCE
RTDCYGNVNRIDTTGASCKTAKPEGLSYCGVSASKKIAERDLQAMDRYKTIIKKVGEKLCVEPAVIAGIISRESHAGKVL
KNGWGDRGNGFGLMQVDKRSHKPQGTWNGEVHITQGTTILINFIKTIQKKFPSWTKDQQLKGGISAYNAGAGNVRSYARM
DIGTTHDDYANDVVARAQYYKQHGY
>16VP:A|PDBID|CHAIN|SEQUENCE... (7 Replies)
Hello,
I have 10 fasta files with sequenced reads information with read sizes from 15 - 35 . I have combined the reads and collapsed in to unique reads and filtered for sizes 18 - 26 bp long unique reads. Now i wanted to count each unique read appearance in all the fasta files and make a table... (5 Replies)
I have a fasta file as follows
>sp|O15090|FABP4_HUMAN Fatty acid-binding protein, adipocyte OS=Homo sapiens GN=FABP4 PE=1 SV=3
MCDAFVGTWKLVSSENFDDYMKEVGVGFATRKVAGMAKPNMIISVNGDVITIKSESTFKN
TEISFILGQEFDEVTADDRKVKSTITLDGGVLVHVQKWDGKSTTIKRKREDDKLVVECVM
KGVTSTRVYERA
>sp|L18484|AP2A2_RAT AP-2... (3 Replies)
Hi
How can I extract sequences from a fasta file with respect a certain criteria? The beginning of my file (containing in total more than 1000 sequences) looks like this:
>H8V34IS02I59VP
SDACNDLTIALLQIAREVRVCNPTFSFRWHPQVKDEVMRECFDCIRQGLG
YPSMRNDPILIANCMNWHGHPLEEARQWVHQACMSPCPSTKHGFQPFRMA... (6 Replies)
I could calculate the length of entire fasta sequences by following command,
awk '/^>/{if (l!="") print l; print; l=0; next}{l+=length($0)}END{print l}' unique.fasta
But, I need to calculate the length of a particular fasta sequence specified/listed in another txt file. The results to to be... (14 Replies)
I have a fasta file as follows
>sp|Q8WWQ8|STAB2_HUMAN Stabilin-2 OS=Homo sapiens OX=9606 GN=STAB2 PE=1 SV=3
MMLQHLVIFCLGLVVQNFCSPAETTGQARRCDRKSLLTIRTECRSCALNLGVKCPDGYTM
ITSGSVGVRDCRYTFEVRTYSLSLPGCRHICRKDYLQPRCCPGRWGPDCIECPGGAGSPC
NGRGSCAEGMEGNGTCSCQEGFGGTACETCADDNLFGPSCSSVCNCVHGVCNSGLDGDGT... (3 Replies)
Hi,
I have to add 7 bases of specific nucleotide at the beginning and ending of all the fasta sequences of a file. For example, I have a multi fasta file namely test.fasta as given below
test.fasta
>TalAA18_Xoo_CIAT_NZ_CP033194.1:_2936369-2939570:+1... (1 Reply)
Discussion started by: dineshkumarsrk
1 Replies
LEARN ABOUT DEBIAN
ntdpal
NTDPAL(1) Primer3 User Manuals NTDPAL(1)NAME
ntdpal - Provides Primer3's alignment functionality
SYNOPSIS
ntdpal [-g gval] [-l lval] [-m mval] [-f1, f2, f3] [-p] [-s] [-e] {seq1} {seq2} {mode}
DESCRIPTION
Ntdpal (NucleoTide Dynamic Programming ALignment) is a stand-alone program that provides Primer3's alignment functionality (local, a.k.a.
Smith-Waterman, global, a.k.a. Needleman-Wunsch, plus "half global").
OPTIONS -g gval
gval is a (positive) float (.01 precision) specifying penaltiy for creating a gap respectively (the penalties are subtracted from the
output score)
-l val
lval is a (positive) float (.01 precision) specifying penaltiy for lengthening a gap respectively (the penalties are subtracted from
the output score)
-a
Causes the scoring matrix to be modified by dpal_set_ambiguity_codes.
-e
Causes the end postion of the alignment in both sequences to be printed. Do not confuse with the 'e' mode.
-f1, -f2, -f3
Force specific implementations. -f2 forces use an implementation that might provide more informative error messages, possibly at the
expense of some speed.
-h
Use a different scoring matrix: G and C matches = 3, A and T = 2, and mismatches = -0.5. (The default scoring matrix assigns 1 to a
match, and -1 to a mismatch.)
-p
Causes the alignment to be displayed on stderr.
-s
causes only the score to printed.
-m mval
is the maximum allowed gap (default is 3).
seq1 and seq2
are the sequences to be aligned.
mode
is one of g, G, l, or L specifying a global, global end-anchored, local, or local end-achored alignment respectively. For backward
compatibility e is equivalent to G.
REFERENCE
Please cite Rozen, S., Skaletsky, H. "Primer3 on the WWW for general users and for biologist programmers." In S. Krawetz and S. Misener,
eds. Bioinformatics Methods and Protocols in the series Methods in Molecular Biology. Humana Press, Totowa, NJ, 2000, pages 365-386.
SEE ALSO primer3_core(1)oligotm(1)COPYRIGHT
Copyright (C) 1996,1997,1998,1999,2000,2001,2004,2006,2007,2008 Whitehead Institute for Biomedical Research, Steve Rozen
(http://jura.wi.mit.edu/rozen), Helen Skaletsky
All rights reserved. On Debian-based systems, please consult /usr/share/doc/primer3/copyright to read the licence of ntdpal.
This manual page was written by Charles Plessy <plessy@debian.org> for the Debian(TM) system (but may be used by others). Permission is
granted to copy, distribute and/or modify this document under the same terms as oligotm itself.
ntdpal 1.1.4 05/09/2008 NTDPAL(1)