I'm not really sure i understand what you are asking for or why.
but here is my take;
cat contigs.fa | perl -p -e 's/^>.*/LINKER/g' | perl -p -e 's/\n//g' > Pseudo.fa
Announcement
Collapse
No announcement yet.
X
-
Well if you don't care about contig ordering then just put a bunch of 'N's (I suggest 100) between each contig. I'd probably use 'sed' for this. Something like:
sed -e 's/>.*/NNNNNNNNNNNNNNNNNN/' <input.fasta>output.fasta
Then reformat if desired (I used fasta_formatter) and put a header on top of the combined sequence.
Leave a comment:
-
Pseudo-molecule
Hi,
How may I concatenate contigs into a single pseudo-molecule using a linker sequence? Does anyone know of a script for this purpose?Tags: None
Latest Articles
Collapse
-
by seqadmin
The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
...-
Channel: Articles
11-27-2023, 01:15 PM -
-
by seqadmin
Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...-
Channel: Articles
11-09-2023, 07:02 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 02:24 PM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 02:24 PM
|
||
Started by seqadmin, Yesterday, 07:37 AM
|
0 responses
23 views
0 likes
|
Last Post
by seqadmin
Yesterday, 07:37 AM
|
||
Started by seqadmin, 12-04-2023, 08:23 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
12-04-2023, 08:23 AM
|
||
Started by seqadmin, 12-01-2023, 09:55 AM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
12-01-2023, 09:55 AM
|
Leave a comment: