Plot the density of genes in R
#column1 = chromosome name and column2 = start position of the gene # check if ggplot2 is installed, if so, load it, # if not, install and load it if("ggplot2" %i...2298 days ago
2089 days ago
Pack a perl program with their dependencies on Ubuntu !
...e medical records for easy access anytime, anywhere. To know more visit here: https://www.lazoi.com/Home Comparative Genomics !!This group is dedicated to Comparative genomics ... Dramatic results have em...1525 days ago
Create random 10000 SNPs in genome !
...] Generating output files .. Generating the correspondance map for genomic variants introduced during si...muSNP.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during si...1225 days ago
Create random 1000 INDEL in genome !
...] Generating output files .. Generating the correspondance map for genomic variants introduced during si...INDEL.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during si...1225 days ago
Create random 1000 CNVs in genome !
...] Generating output files .. Generating the correspondance map for genomic variants introduced during si...muCNV.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during si...1225 days ago
Create random 5 inversions in genome !
...] Generating output files .. Generating the correspondance map for genomic variants introduced during si...muINV.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during si...1225 days ago
Create random 2 translocations in genome !
...] Generating output files .. Generating the correspondance map for genomic variants introduced during si...TRANS.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during si...1225 days ago
Oneliner to convert lower-case to sequence masked with Ns
perl -pe '/^[^>]/ and $_=~ s/[a-z]/N/g' genomic.fna > genomic.N-masked.fna awk '{if(/^[^>]/)gsub(/[a-z]/,"N");print $0}' genomic.fna > genomic.N-masked.fna964 days ago
Extract the sequences with IDs !
#sed -i 's/\_/ /g' Delta_seqID_from_lineage_report.txt seqtk subseq genomic.fna Delta_seqID_from_lineage_report.txt > Delta.fasta #Split the fasta in 11 equal sequences subsets pyfasta split -n 11 Delta.fasta865 days ago