2040 days ago
1739 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 sam...1617 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 sor...1573 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 The...1564 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...1225 days ago
Create random 1000 INDEL in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuSNP.simseq.genome.fa -indel_count 1000 -prefix simuINDEL [Sun Jan 10 16:14:00 2021] Starting simuG .....1225 days ago
Create random 1000 CNVs in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuINDEL.simseq.genome.fa -cnv_count 100 -prefix simuCNV [Sun Jan 10 16:24:20 2021] Starting simuG .. [Su...1225 days ago
Create random 5 inversions in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuCNV.simseq.genome.fa -inversion_count 5 -prefix simuINV [Sun Jan 10 16:30:40 2021] Starting simuG .. [...1225 days ago
Create random 2 translocations in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuINV.simseq.genome.fa -translocation_count 2 -prefix simuTRANS [Sun Jan 10 17:12:58 2021] Starting simuG ....1225 days ago