1573 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...1616 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...1963 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...2228 days ago
3052 days ago