Question: generate fastq with bed file
0
gravatar for sacha
3.1 years ago by
sacha1.3k
France
sacha1.3k 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 • 1.7k views
ADD COMMENTlink modified 3.1 years ago • written 3.1 years ago by sacha1.3k

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 3.1 years ago by Pierre Lindenbaum111k
1
gravatar for sacha
3.1 years ago by
sacha1.3k
France
sacha1.3k 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 3.1 years ago • written 3.1 years ago by sacha1.3k

That generates a fasta, not a fastq :)

ADD REPLYlink written 16 months ago by endrebak650

fasta to fastq is easy !

ADD REPLYlink written 15 months ago by sacha1.3k
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: 637 users visited in the last hour