Perl script to break the contigs by 'N'
#!/usr/bin/perl -w use Bio::SeqIO; use strict; my $fasta = Bio::SeqIO->new( -file => "$ARGV[0].parts.fasta",-format=>'fasta'); open(SCAFF,">$ARGV[0].parts.scaff");...2078 days ago
To convert just one specific read group to fastq
# Stop script on error. set -uex # The SRR BioProject number for the sequencing data. PROJECT=PRJNA257197 # The number of datasets to subselect from the project. N...1557 days ago
Installing docker for Bioinformatics on Ubuntu !
jit@jit-HP-Pro-3335-MT:~/Downloads$ sudo apt-get remove docker docker-engine docker.io containerd runc Reading package lists... Done Building dependency tree Read...1528 days ago
Commands to install conda in Ubuntu !
jit@jit-HP-Pro-3335-MT:~/Downloads$ mkdir jittmp jit@jit-HP-Pro-3335-MT:~/Downloads$ cd jittmp/ jit@jit-HP-Pro-3335-MT:~/Downloads/jittmp$ curl -O https://repo.anaconda...1526 days ago
Bash script to extract intronic fragments !
#To obtain introns, we simply need the gene and exonic coordinates; #by subtracting the exonic regions from the genic region, we have the intronic region. gunzip -c ge...1381 days ago
Bash script to get intergenic region from genome files !
#For the intergenic region, we will require the size of the chromosomes. wget http://xxx.chrom.sizes cat xxx.chrom.sizes | sed 's/^chr//' | sed 's/Cp/Pt/' > tmp mv tmp...1381 days ago
Install and set up i-adhore for synteny and wgd analysis ! -- step by step --
#Need to download i-adhore-3.0.01.tar.gz from https://wdiceryfd4rjn74bjhhtz2k5di--bioinformatics-psb-ugent-be.translate.goog/webtools/i-adhore/licensing/ #Follow these...1218 days ago
Commandline to Extract a list of specific read IDs from a bam file
#Save all the IDs in IDs.txt file #Run it on BAM file samtools view file.bam | fgrep -w -f IDs.txt1011 days ago
965 days ago
Perl script for six frame translation !
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; # Path to your input nucleotide sequence file in FASTA format my $input_fasta = 'path/to/your/input.fasta';...109 days ago