Question: generate fastq with bed file
0
gravatar for sacha
22 months ago by
sacha720
France
sacha720 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 • 847 views
ADD COMMENTlink modified 22 months ago • written 22 months ago by sacha720

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 22 months ago by Pierre Lindenbaum92k
0
gravatar for sacha
22 months ago by
sacha720
France
sacha720 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 22 months ago • written 22 months ago by sacha720

That generates a fasta, not a fastq :)

ADD REPLYlink written 4 days ago by endrebak560
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: 1438 users visited in the last hour