Perl subroutine to read genome/reads fasta file !
sub readSeqFromFasta{ my $file = $_[0]; my (%map,$id,$seq); open(IN,$file) or die("Cannot open file for reading $file:$!\n"); while(){ chomp; if($_ =~ /^>/...2127 days ago
2121 days ago
2120 days ago
Install Parrot Virtual Machine !
#Parrot is a virtual machine designed to efficiently compile and execute bytecode for dynamic languages jit@jit-HP-Pro-3335-MT:~/Downloads/parrot-8.1.0$ wget ftp://ft...1531 days ago
Long reads mapper bash script !
#!/bin/bash #only for LONG READS mapping #USAGE: runMapper.sh minimap2 ref.fa reads.fa/fq 48 ont echo -e "This is a bash scrip to map your long reads and make it v...2116 days ago
2117 days ago
Installing Porechop on Ubuntu !
➜ Tools git:(master) ✗ git clone https://github.com/rrwick/Porechop.git Cloning into 'Porechop'... remote: Counting objects: 1579, done. remote: Total 1579 (delta 0)...2099 days ago
2092 days ago
Perl script to extract sequence by Ids from multifasta file !
#!/usr/bin/perl -w use strict; my $idsfile = "$ARGV[0]"; my $seqfile = "$ARGV[1]"; my %ids = (); open FILE, $idsfile; while() { chomp; $ids{$_} += 1;...2089 days ago
2069 days ago