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")...2061 days ago
Generate simulated polyploid genome !
...rozygous3percent.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheterozygous3percent.fasta #generated 50X of Illumina paired-end reads sim_reads --depth 50...2051 days ago
2002 days ago
1603 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 s...1601 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...1557 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 % Rec...1553 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....1540 days ago
1523 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... Do...1509 days ago