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($_ =~ /^>...2134 days ago
2128 days ago
2127 days ago
Long reads mapper bash script !
#!/bin/bash #only for LONG READS mapping #USAGE: runMapper.sh minimap2...echo "Name of the file used for mapping $fileName, present at $refFas...lName == "bwa" ]; then echo "Mapping with $toolName" $bwaMemLo...ame == "ngmlr" ]; then echo "Mapping with $toolName" if [ $reads...2123 days ago
2124 days ago
2076 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 treatBAMfile2056 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...2054 days ago
Installing pb-assembly on Linux !
[jnarayan@hmem00 ~]$ module avail -------------------------------------------------------- /usr/share/Modules/modulefiles ------------------------------------------...2035 days ago
2027 days ago