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

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 2.1 years ago by Pierre Lindenbaum96k
1
gravatar for sacha
2.1 years ago by
sacha770
France
sacha770 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 2.1 years ago • written 2.1 years ago by sacha770

That generates a fasta, not a fastq :)

ADD REPLYlink written 3 months ago by endrebak570

fasta to fastq is easy !

ADD REPLYlink written 12 weeks ago by sacha770
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: 871 users visited in the last hour