1524 days ago
Comment on "Perl script to count number of Ns in a multifasta file !"
If you want to count all 'N' in miltifasta file. (base) ? output_2_test git:(master) ? more scaffolds.fasta | grep -Ho N * | uniq -c 12 draf...1567 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 sequenc...1914 days ago
Comment on "Run miniasm assembler on nanopore reads !"
....gfa Convert GFA to FASTA: awk '/^S/{print ">"$2"\n"$3}' reads.gfa | fold > reads.fa And then count how many contigs: grep ">" reads.fa | wc -l...2178 days ago
3002 days ago