Bash script to get exon fragments from genome files !
#Exons are already defined in the GTF file, so we simply need to print lines that are marked exonic. gunzip -c genome_file.gtf.gz | awk 'BEGIN{OFS="\t";} $3=="exon" {print $1,$4-1,$5}' | bedtools sort | bedtools merge -i - | gzip > my_exon.bed.gz1365 days ago
1207 days ago
Set up WGD environment using conda !
#Wgd cannot be installed directly with bioconda at present, so it is slightly troublesome to install, because it #depends on a lot of software. wgd depends on the fol...1202 days ago
989 days ago
Commandline for paired end reads simulation with BBMap !
...genScaffoldInfo=true Writing chunk 1 Waiting for writing to finish. Finished. snpRate=0.0, max=0, unique=true insRate=0.0, max=0, len=(0-0) delRate=0.0, max=0, len=(0-0) sub...989 days ago
927 days ago
Bash script to simulate a genome !
...upper.diploid.fasta.gz \ ziplevel=6 \ ploidy=2 \ subrate=0.0192 \ indelrate=0.001 \ maxindel=20 \ no...+"` echo "best k=${k}" # make 30x haploid coverage for PacBio CLR reads # error rate from 13 -...917 days ago
Run Pango on your multifasta file !
...bined.fa (pangolin) [jnarayan@hn1 FASTA]$ pangolin --update pangolin already latest release (v3.1.16)...olearn updated to 2021-11-18 constellations updated to v0.0.24 scorpio already latest release (v0.3.14)...886 days ago
861 days ago
Installing SEVA environment in Conda !
...########################################################## | 100% r-modelr-0.1.8 | 218 KB | ###...########################################################### | 100% r-callr-3.7.0 | 440 KB | ##...837 days ago