Bash script to alignment of short reads against reference genome !
bwa mem -t 40 -R '@RG\tID:K12\tSM:K12' \ E.coli_K12_MG1655.fa SRR1770413_1.fastq.gz SRR1770413_2.fastq.gz \ | samtools view -b - >SRR1770413.raw.bam sambamba so...1583 days ago
Pack a perl program with their dependencies on Ubuntu !
#Follow steps to create your own executable ./web jit@jit-HP-Pro-3335-MT:~/Downloads/autoConTAMPR/bin$ sudo apt install libpar-packer-perl Reading package lists... Don...1535 days ago
Install Ragout genome assembler
$ conda install -c bioconda ragout Collecting package metadata (repodata.json): done Solving environment: done ## Package Plan ## environment location: /home/ana...1488 days ago
Onliner to split the multifasta to singlefasta files !
#Split the multifasta to singlefasta # Multi fasta #Single fasta awk '$0 ~ "^>" { match($1, /^>([^:]+)/, id); filename=id[1]} {print >> filename".fa"}' sequence.fasta1424 days ago
get GC across the entire CDS !
#look at GC across the entire CDS. gffread -x - -g | \ seqtk comp - | \ awk -v OFS="\t" '{ print $1, "0", $2, ($4 + $5) / $2 }'1417 days ago
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.gz1389 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 tm...1389 days ago
Bash script to handle Multifasta files
#Convert all lowercase residues to uppercase in a FASTA sequence file $ awk 'BEGIN{FS=" "}{if(!/>/){print toupper($0)}else{print $1}}' input.fasta > output.fasta #Re...1381 days ago
1239 days ago
Create random 10000 SNPs in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq SGDref.R64-2-1.dups.fa -snp_count 10000 -prefix simuSNP [Sun Jan 10 16:05:57 2021] Starting simuG .. [Sun...1235 days ago