1946 days ago
Comment on "Extract fasta sequence from a multifasta file with fasta header Ids"
perl -ne 'if(/^>(\S+)/){$c=grep{/^$1$/}qw(id1 id2)}print if $c' fasta.file If you have a large number of sequences that you want to extract, then you most likely have the sequence identifiers in a separate file. A...1953 days ago
Comment on "List of motif discovery tools !"
Motif algorithms This link on Nature explain: how does DNA sequence motif discovery work? Paper on BMC bioinformatics: survey of DNA motif finding algorithms This paper discusses some limitations and potentials of motif discovery algorithms (2005)2024 days ago
Comment on "Nanopolis: polish a genome assembly"
# Index the draft genome bwa index draft.fa # Align the basecalled reads to the draft sequence bwa mem -x ont2d -t 8 draft.fa reads.fa | samtools sort -o reads.sorted.bam -T reads...2142 days ago
Comment on "Ancestral sequence reconstruction steps !"
Parallelization of MAFFT for large-scale multiple sequence alignmentshttps://academic.oup.com/bioinformatics/article/34/14/2490/49160992164 days ago
2165 days ago
2180 days ago
Comment on "poRe: an R package for the visualization and analysis of nanopore sequencing data"
...des-3.6.2-Linux.tar.gz gunzip < SPAdes-3.6.2-Linux.tar.gz | tar xvf - Now, we are ready to go. First off, let’s extract the 2D sequence data as FASTQ from the MinION...2217 days ago
Comment on "Basic command-line to run BLAST"
...A/Human_vs_Avaga_protein/mouse.protein.faNew DB title: mouse.protein.faSequence type: ProteinKeep MBits: TMaximum file size: 1000000000BAdding sequences from FASTA; added 76216 sequ...2287 days ago
2295 days ago