1849 days ago
682 days ago
splitbam: splits a BAM by chromosomes
...file/bam pattern) REQUIRED. MUST contain __CHROM__ and end with .bam -s assume input is sorted. -x | --index create index. -t | --tmp (dir) tmp file directory -G (file) chr...2636 days ago
Find certain files/documents in Linux OS
.... Find all Empty FilesTo file all empty files under certain path.# find /tmp -type f -empty20. Find all Empty DirectoriesTo file all empty directories under certain path.# find /tmp -type d -empty21. File all Hi...3695 days ago
Linux Sort Commands for Bioinformatics
...ed on its first column, then the second if the first is the same, and so on:sort input.txt sort a huge file (GNU sort ONLY):sort -S 1500M -t $HOME/tmp input.txt > sorted.txt so...3640 days ago
2364 days ago
2203 days ago
Insert the sequence at desire location in multi-fasta file with Perl
#!/usr/bin/perl use warnings; use strict; use Bio::SeqIO; use Bio::Seq; use File::Copy; #ARGV[0] should be in following format --- Keep the coordinate sorte...2327 days ago
Perl script to read multi fasta sequence one by one
#!/usr/bin/env perl use strict; use warnings; #USAGE #perl rohanRun.pl seq.fa my $outfile='tmp.fa'; my $fastaSeq_ref = readfasta ("$ARGV[0]"); my %fastaSeq = %$fastaSeq_r...2234 days ago
Comment on "Nanopolis: polish a genome assembly"
...raft.fa # Align the basecalled reads to the draft sequence bwa mem -x ont2d -t 8 draft.fa reads.fa | samtools sort -o reads.sorted.bam -T reads.tmp - samtools index reads.sorte...2112 days ago
2812 days ago