FASTQ

convert bam to fastq and remap it with bwa mem

Sometimes, you get a processed bam file from other resources, but you want to map it with another aligner. For example, the old BWA aln does not keep the…

rename a bunch of files with bash by regular expression

I am at the MSU NGS 2014 course. My colleague wanted to follow the khmer protocol with his own data, but one of the steps has to use a certain file name…

quality control of your fastq file! my first post after I got my PhD!

I passed my defense on Tuesday, and now I am a PhD! what a long process, it took me six years. I’ve been enjoying the PhD experience though.

qseq to fastq conversion

I was looking at some old ChIP-seq data with raw file in qseq format. I want to convert them to fastq file for mapping with bowtie.

use sed to print out every n lines, split pair-end FASTQ file

#print out the first two lines of every 4 lines. -n flag suppress all of the other lines and only print the line\ you specified. -e option tells sed to…

FASTA/Q file manipulation tool readfq

it is developed by Heng Li in princeton.

Awk Command to Count Total, Unique, and the Most Abundant Read in a FASTQ file

I was reading through a paper on comparative ChIP-Seq when I found this awk gem that lets you get some very basic stats very quickly on next generation…

Cool use of Unix paste with NGS sequences:FASTQ to FASTA

While browsing SeqAnswers.com today I came across a post where Uwe Appelt provided a couple of lines of Unix shell wizadry to solve some problem. What…