Question: generate fastq with bed file
0
gravatar for sacha
23 months ago by
sacha740
France
sacha740 wrote:

I want to generate fastq file for testing.. I already check FastQArtificialGenerator, wsgim and Simseq. But it seems I cannot specify a region from a bed file. 

For exemple, from hg19, I only want a fastq where reads map my bed file ( exome) . 

generator fastq bed • 912 views
ADD COMMENTlink modified 23 months ago • written 23 months ago by sacha740

please comment or mark as 'answered' your previous questions: How to get exon number from position from script? ; genome browser with genomic coordinate ; samtools : print single columns at position ; bcftools extract data

ADD REPLYlink written 23 months ago by Pierre Lindenbaum94k
1
gravatar for sacha
23 months ago by
sacha740
France
sacha740 wrote:

I solve it .. 

I just need to create a new reference containing my bed coordinate !

bedtools getfasta -fi ref.fa -bed target.bed -fo newref.fa 
ADD COMMENTlink modified 23 months ago • written 23 months ago by sacha740

That generates a fasta, not a fastq :)

ADD REPLYlink written 5 weeks ago by endrebak570

fasta to fastq is easy !

ADD REPLYlink written 4 weeks ago by sacha740
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: 943 users visited in the last hour