I have a chipseq data where I am looking for mitochondrial genes differential binding. The samples are for human cell lines. My question is if I align reads with hg19 genome which as such will have mitochondrial genome will that give me a correct representation or I have to specifically take mitochondrial genome as reference for using any aligner
Thanks
1 answer
EIther way, you will get the information.
Map all reads and filter the Bed/Bam file for chrM or just subset the reference fasta index before mapping for chrM and then map and proceed with downstream analysis. The former will take more time, depending on the file size.
Oh, I never thought about that, have a no idea. People at BWA/Bowtie support could answer.
There are papers for circular genome mappings though.
http://www.dcs.kcl.ac.uk/staff/csi/publications/ITAB_2010.pdf
Hmm, being discussed here, 2 strategies
- Either add the beginning of the genome/chromosome to the end for the length of the reads, so you don't loose anything.
- With Bowtie, use
-M 1so that the read is assigned randomly. (quite not clear to me) - Use NovoAlign
http://www.researchgate.net/post/NGS_sequencing_issue_When_the_regerence_genome_is_circular
•
written
by
Sukhi Singh
Log in to answer this question.
-
Generate Read counts from bam file
written by singhankit973 •
-
bowtie2; using for multiple fastq files, linux loop code
written by Farzaneh •
-
Where can I find a fasta for mitochondrial genes?
written by blur
-
microRNA- sequencing analysis
written by kanwarjag
-
question about homologous binding motifs
written by kanwarjag
-
Chipseq overlap of peaks
written by kanwarjag
-
why location of binding sites are different while human sequence is the same?
written by na.cna30 •
-
Chip-Seq For Mitochondrial Genes
written by kanwarjag
-
Overlapping Pe Reads Alignments
written by kanwarjag
-
Chipseq Data On Prostate Cancer
written by Bioinfosm
