Find and replace in multifasta or fasta header with perl onliner
You have a fasta file and you want to replace: "|" You are told to replace that by "_" perl -i -p -e "s/\|/_/g" genome.fasta -i = inplace editing -p = loop over lines and print each line (after processing) -e = command line script1668 days ago
Samtools commands for bioinformatician !
## count mapped reads samtools view -c -F 260 mapping_file.bam ### converting sam file into fasta samtools fasta reads_mapped.sam > reads.fasta ### converting sa...1655 days ago
Perl script to remove duplicated lines !
#!/usr/bin/perl use strict; use warnings; { $_ = ; my $next_line; while( $next_line = ) { #print "current line: $_ -- next line: $next_line$/";...1612 days ago
Bash script to download SRA file !
#We can use the sratoolkit to directly pull the sequence data (in paired FASTQ format) from the archive. fastq-dump is in the SRA toolkit. It allows directly downloading d...1611 days ago
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...1611 days ago
Bash commandline to install Anaconda !
#The line begins with $ are the commands $ mkdir tmp $ cd tmp/ $ curl -O https://repo.anaconda.com/archive/Anaconda3-2019.03-Linux-x86_64.sh % Total % Rece...1607 days ago
Bash command to install Miniconda !
(base) [wsu29@bladeamd-2 tmp]$ wget https://repo.anaconda.com/miniconda/Miniconda3-latest-Linux-x86_64.sh --2020-02-01 13:25:25-- https://repo.anaconda.com/miniconda/Min...1607 days ago
Bash command to install bwa, samtools, picard !
(base) [wsu29@bladeamd-2 tmp]$ conda install bwa Collecting package metadata (repodata.json): done Solving environment: - The environment is inconsistent, please check...1607 days ago
Bash command to install GATK, Bedtools and SnpEff !
(my_GATK) [wsu29@bladeamd-2 tmp]$ conda install gatk Collecting package metadata (repodata.json): done Solving environment: done ## Package Plan ## environment l...1607 days ago
Installing ggplot2 and its dependencies on Ubuntu !
jit@jit-HP-Pro-3335-MT:~/Downloads/MitoHunter/minidot/bin$ sudo R [sudo] password for jit: R version 3.4.4 (2018-03-15) -- "Someone to Lean On" Copyright (C) 2018 Th...1602 days ago