BBSplit: Read Binning Tool for Metagenomes and Contaminated Libraries
BBSplit internally uses BBMap to map reads to multiple genomes at once,...uo; flag. The normal use of BBSplit is like Seal, either quantifying how many reads go to each reference, or spli...lmonella, you could do this:bbsplit.sh in=reads.fq ref=ecoli.fa,s...2346 days ago
2336 days ago
PerlOneLiner for Bioinformatician
FILE SPACING------------ # Double space a fileperl -pe '$\="\n"'perl -pe 'BEGIN { $\="..."@list"' # Create an array from a string@months = split ' ', "Jan Feb Mar Apr May Jun...in the stringperl -le 'print join ", ", map { ord } split...2333 days ago
List of visualization tools for genome alignments
Genome browsers are useful not only for showing fina...or.htm JVM Java Visual Mapping tool for NGS reads http://www.springer.com/cda/...gnments of huge-scale single-end and pair-end short reads http://omictools.com/mapview...2315 days ago
Awk for Bioinformatician and computational biologist
...ions} /pattern2/ {Actions}' file The working of Awk is as followsAwk reads the input files one line at a...TATATAAwk has number of builtin variables. For each record i.e line, it splits the record delimited by whit...2311 days ago
Some useful Bioinformatics links
...apping of nanopore sequencing reads with GraphMap | Nature Commun...s/ncomms11307 lumpy-sv/extractSplitReads_BwaMem at master ·...sv/blob/master/scripts/extractSplitReads_BwaMem jts/nanocorrect:...ols.html How To Filter Mapped Reads With Samtools https://www.bi...2273 days ago
2249 days ago
Gap filling or Contigs extensions tools !
There are many tools to perform gap...to fill the gap within paired reads" or "Toward almost closed gen...local re-assemblies using 454 reads. We used GAPresolution but it...of interest, and the set S of reads those contigs contain. Let C...ctor, and generate paired-end reads for this library; then try st...2196 days ago
Frequent Paired-end reads (PE 2x100) mapping command lines
bowtie2 -x hs37m -X 650 -q -1 r1.fq -2 r2.fq -S r12.bowtie2.sam bwa aln hs37m.fa r1.fq > r1.sai && bwa aln hs37m.fa r2.fq > r2.sai \ &&am...2213 days ago
LINKS scaffolder bloomfilter setting !
➜ bin git:(master) ✗ ls -l total 68 drwxrwxr-x 3 urbe...-s file-of-filenames, full path to long sequence reads or MPET pairs [see below] (Multi-FASTA/fastq format, required) -m MPET reads (default -m 1 = yes, default...2182 days ago