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($_ =~ /^>...2161 days ago
2154 days ago
2151 days ago
2103 days ago
Generates a genome coverage plot with R
library(CoverageView) ##draw a coverage plot for a test case BAM file #get a BAM test file treatBAMfile2083 days ago
Generate simulated polyploid genome !
#Generate 3% divergence msbar -point 4 -count 16558 toy.fasta > toyheterozygous3percent.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheteroz...2081 days ago
Perl script to create a consensus of nucleotide sequences !
use strict; use warnings; my @instances = qw ( AAAAA ATCGA ATAAA ); my @instances2 = qw ( AAAAA AACGA ATAAA AGA...{$b} $prefOrder{$a} } qw/A T G C/ ]->[0]; } return @cons; } #refere...2060 days ago
Setting up falconUnzip conda environments for genome assembly !
➜ Analysis_Results conda create -n denovo_asm Solving environment: done ## Package Plan ## environment location: /home/urbe/anaconda3/envs/denovo_asm Pro...2018 days ago
1753 days ago
Perl script to run in parellel !
#!/usr/bin/perl use strict; use warning...my ($sequence_data_ref) = parse_genome_files($ARGV[0]); my %genome=%{$sequence_data_ref}; my...; foreach my $chr_set (keys %genome) { $count++...e FILE or die $!; } sub parse_genome_files { my $file=shift;...1714 days ago