Mapping with BWA-mem or BWA-sampe in one go with python script !
BAM files and mapping BESST requires sorted and indexed BAM files as input. Any read aligner + samtools can be used to obtain such files. Read pairs needs to be aligned...2177 days ago
Perl subroutine for reading multifasta file !
sub readfasta { (my $file)=@_; my %sequence; my $header; my $temp_seq; #suppose fasta files contains multiple sequences; open (IN, "2197 days ago
Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here...2193 days ago
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($_ =~ /^>/...2186 days ago
2180 days ago
2179 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...1590 days ago
2176 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)...2159 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;...2149 days ago