This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to identify every read in a fastq file?

Hi, Everybody is looking for their target reads in a fastq file, and I am just sitting here and can not find a good program to identify the remaining (junk) reads. I am working on ancient DNA (currently horse) illumina reads, and I want to identify the exact organisms (possibly pathogens, human or other contamination, etc) of the remaining reads besides the horse sequences (approx 20-30% of the data contains horse DNA actually). So, can anyone recommend a good program for this task? Thanks in advance :)

genome next-gen blast
• 3,846 views • link
written by gerberd1990 • gerberd1990

WOW, thanks for everybody, I see many valuable information here :)

• 1 views • link
written by gerberd1990 • gerberd1990

4 answers

First you need to align your reads to the expected reference that it may be contaminated with using FastQ Screen Or DeconSeq, then below post to remove it

http://seqanswers.com/forums/showpost.php?p=109308&postcount=6

• 1 views • link
written by Medhat Medhat

This would seem like a good use case for something like Kraken: https://ccb.jhu.edu/software/kraken/

But your ability to assign every read in your experiment to an organism of origin will depend entirely on the completeness of your database.

• 1 views • link
written by Simon Cockell Simon Cockell

This might be worth trying: https://github.com/smangul1/rop, see paper: http://biorxiv.org/content/early/2016/05/13/053041

• 1 views • link
written by WouterDeCoster WouterDeCoster

I think it is specialized in:

discover the source of all reads, which originate from complex RNA molecules, recombinant antibodies and microbial communities.

• 1 views • link
written by Medhat Medhat

Yes, it is not meant for this application but might give some ideas.

• 1 views • link
written by WouterDeCoster WouterDeCoster

There are some suggestions in these previous threads:

• 1 views • link
written by igor igor

Log in to answer this question.

More posts like this