Question: Ways To Detect Bias In Dna Sampling For Genomic Sequencing
1
gravatar for Lee Katz
6.1 years ago by
Lee Katz2.9k
Atlanta, GA
Lee Katz2.9k wrote:

Hi all, I need to detect DNA sampling bias. We extracted DNA using two different library preps and then multiplexed the samples.

I think that sampling bias will basically be something where each read will be more likely to start with a certain GC content or sub sequence. Therefore I think the best way to analyze any bias will be to assemble with an assembly output file (in my case, I'll do velvet, with the output afg file). After that, I am not exactly sure what to do. I could look places with low depth, e.g., at contig breaks, and see what the first 10 bp are at each place.

Another way might be to look at the first 10 bp at each read and see if there is any bias too. Gibbs sampling, %GC, some kind of kmer analysis?

But anyway, what is a good standard or common-sense way to do this? Thank you BioStar community!

genome sequencing • 3.0k views
ADD COMMENTlink modified 6.1 years ago by Zev.Kronenberg11k • written 6.1 years ago by Lee Katz2.9k
5
gravatar for Matt Shirley
6.1 years ago by
Matt Shirley8.7k
Cambridge, MA
Matt Shirley8.7k wrote:

FastQC will give you over-represented subsequences, as well as kmers and GC content.

ADD COMMENTlink written 6.1 years ago by Matt Shirley8.7k

as well as nucleotide distributions along each position of the reads in your sample ...

ADD REPLYlink written 6.1 years ago by Steve Lianoglou4.9k
1
gravatar for Zev.Kronenberg
6.1 years ago by
United States
Zev.Kronenberg11k wrote:

I would sugest trying MEME.

You can convert your fastq to fasta and sample some reads (for speed).

If you have a motif (rather than just elevated gc content) then you will get a nice motif plot.

http://tools.genouest.org/tools/meme/doc/examples/meme_example_output_files/meme.html

ADD COMMENTlink modified 6.1 years ago • written 6.1 years ago by Zev.Kronenberg11k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1169 users visited in the last hour