Bash script to handle Multifasta files
...FASTA sequences according to their length $ awk '/^>/ {printf("%s%s\t",(N>0?"\n":""), $0);N++;next;} {printf("%s",$0);} END {printf("\n");...nning of headers in a FASTA file $ awk '{if ($0 ~/_/) {printf ">";} print $0; }' input.fast...1367 days ago
1226 days ago
1225 days ago
1213 days ago
Install and set up i-adhore for synteny and wgd analysis ! -- step by step --
...usr/include/x86_64-linux-gnu/bits/stdio2.h:34:43: note: ‘__builtin___sprintf_chk’ output between 8 and 26...usr/include/x86_64-linux-gnu/bits/stdio2.h:34:43: note: ‘__builtin___sprintf_chk’ output between 8 and 26...1213 days ago
Parse the NCBI taxonomy database with Perl !
...se Bio::DB::Taxonomy; use warnings; my $sps="human"; # Get one from a NCBI taxonomy database my $dbh = Bio::DB::Taxonomy->new(-source => 'flatfile',...1182 days ago
Perl script to get all the descendent in a tree !
...e warnings; use Bio::DB::Taxonomy; #Rotifera 10190 # Get one from a NCBI taxonomy database my $dbh = Bio::DB::Taxonomy->new(-source => 'flatfile',...1174 days ago
Commandline for paired end reads simulation with BBMap !
...axinsert=600, gaussian] Writing reference. Executing dna.FastaToChromArrays2 [mixed.fa, 1, writeinthread=false, genscaffoldinfo=true, retain, waitforwriting=false, gz=true, maxl...1000 days ago
Python script to read FASTA and FASTQ file !
# !/usr/bin/env python3 # -*- coding: utf-8 -*- from pysam import FastxFile def read_fasta_q_file(fasta_q_file): """Parse FASTA/Q file using `pysam.FastxFile`....986 days ago
972 days ago