Hi
I'm trying to recover whole chloroplast sequence data from paired-end Illumina reads from a pool-seq project (25 individuals per pool). We have access to a close reference but not a species-specific chloroplast reference. Following alignment and filtering of reads from one lane (basically following the Popoolation pipeline), I'm getting coverage of about 10% of chloroplast genome at a read depth of 10. This seems very low to me...I expected that because there are many more copies of the cp genome than the nuclear genome, I would get higher cp coverage (our nuclear coverage is at about 70% at a read depth of 25-500). Does anyone have any numbers that I can compare to (e.g. cp coverage from paired-end Illumina data, not necessarily pooled)?
Thanks!
1 answer
I just want to make the note that the chloroplast sequence is also much shorter than that of the nuclear genome. So it could have hundreds or thousands of copies and still be much shorter than the nuclear DNA.
•
written
by
Istvan Albert
Log in to answer this question.
-
Tools for Subsampling Paired-End Reads to 25X Coverage ?
written by Sony Nguyen •
-
De novo assembly of chloroplast genome
written by Eisuan •
-
ATAC-seq - remove mt reads using KMC
written by maxrwjones •
-
Merge single end and paired end illumina raw sequences for SNP calling
written by prashantwaiker •
-
How to extract only the alignments to the chloroplast reference genome of a sam file
written by macielrodriguez2 •
-
Genome coverage with pool of individuals
written by jomaco
-
One sample pooled in a whole lane vs multiple sample in each lane
written by bob •
-
How can I use BLAST to extract chloroplast sequences from DNA reads?
written by AcademicDialysis
-
Multiple sequence alignment and editing following consensus sequence generation in Samtools
written by jal •
-
Best way to remove contaminants to get nuclear genome
written by williamjohn360
