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...1577 days ago
Comment on "Extract fasta sequence from a multifasta file with fasta header Ids"
...){$c=grep{/^$1$/}qw(id1 id2)}print if $c' fasta.file If you have a large number o...ely have the sequence identifiers in a separate file. Assuming that you have one s...$c=$i{$1}}$c?print:chomp;$i{$_}=1 if @ARGV' ids.file fasta.file1924 days ago
2113 days ago
Comment on "LAMSA: fast split read alignment with long approximate matches"
...-soft-clip Use soft clipping for supplementary alignment. [false] -C --comment Append FASTQ comment to SAM output. [false] -o --output [STR] Output file (SAM format). [stdout] -h --...2137 days ago
2137 days ago
2151 days ago
Comment on "Run miniasm assembler on nanopore reads !"
...Bio data. Minimap only accepts two FASTQ files and you need to map your FASTQ file against itself. So, if you ha...iasm: This uses your concatenated FASTQ file and the PAF file output from...FASTQ reads back to the contigs (the PAF file) created in the fist step. Th...2181 days ago
2189 days ago
Comment on "Perl one-liner for bioinformatician !!!"
Convert file to lowercase dd $ dd if=input.txt of=output.txt conv=lcase tr $ tr '[:upper:]' '[:lower:]' < input.txt > output.txt awk $ awk '{ print tolower($0...2228 days ago
Comment on "Basic command-line to run BLAST"
...1:52:20New DB name: /home/urbe/DATA/Human_vs_Avaga_protein/mouse.protein.faNew DB title: mouse.protein.faSequence type: ProteinKeep MBits: TMaximum file size: 1000000000BAdding seque...2258 days ago