# How to convert a .fastq chip-seq file into a wig/bigwig track and how to extract chip-seq peaks?

I have a reference .fasta file and a raw .fastq file with chip-seq data. I am trying to create a bigwig track from the and .fastq and .fasta ref file of the raw signal. Then I would like to do some peak analysis for the chip-seq track.

Can someone help me with googling the write things, I feel really lost while trying to google it and not sure where to start from. I was thinking to do .fastq -> .bam -> wig -> bigwig and then somehow to do peak detection analysis.

– benn
Aug 1 '19 at 8:38

Before you proceed, it is usual to do quality control on the .fastq data. The FASTQC program is a good place to start.

Once you have checked the quality, trimming is the typical next step to remove low-quality data. There are options such as Trimomatic, trimGalore for this.

Once you are sure you have clean data, your first step is alignment - align the sequences in the .fastq file to the reference sequence in the .fasta file. bowtie2 is one option for this process. This gets you to the .bam stage.

These steps will take a little time, some research, and may be somewhat frustrating, but they are the foundation for what you want to do. One option for somewhat automating these steps is to use the Galaxy platform (usegalaxy.org), otherwise these are mostly command-line tools.

The bam to wig and wig to bigwig steps are relatively simple conversions, and Biostars has answers for these.

• What about peak calling once I have a bigwig file?
– 0x90
Aug 1 '19 at 14:12
• I have heard of the R package BayesPeak, but ChIP-Seq is not really my thing. Aug 1 '19 at 15:38
• macs2 should do the work.
– 0x90
Aug 5 '19 at 2:48