This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Where can I get fastq sequence files of length 150 and more.

I want to perform alignment of the human sequence reads with the reference genome. I need reads length 150 and more (500) to test some algorithm. Where can I find such type of reads, both single and paired-end. I got reads from 1000 genome project around 100 length, however I want reads of length more.

alignment sequence
• 2,495 views • link
updated by lakhujanivijay lakhujanivijay • written by AHW • AHW

2 answers

Use Sequence Read Archive advanced search and provide read length as 150 and species as Homo sapiens. You can add more filters if you want. The query should look something like this:

(150[ReadLength]) AND "Homo sapiens"[orgn:__txid9606]

And here is a result for paired-end sequencing data with read length 150.

• 0 views • link
written by Gungor Budak Gungor Budak

This is a very basic requirement and a lot of tools are available to simulate artificial reads from the genome under question. More interestingly, you can define the number, length and quality of reads also. One such well documented program ( ArtificialFastqGenerator )is here

Here is the publication link

Another very sophisticated tool is ART (courtesy: Biostars Handbook) which can mimic the sequencing platforms very well.

• 0 views • link
written by lakhujanivijay lakhujanivijay

Log in to answer this question.

More posts like this