1518 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...1561 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...1908 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...2173 days ago
2996 days ago