Entering edit mode
4.2 years ago
ss20g16
•
0
How would I convert a file that is fq.gz (fastq.gz) to BED file format? I'm working on the Iridis command line (MobaXterm; Linux command line) - what code can I use to achieve this? thanks
May I ask why you are doing that? (by the way, do you mean fastq.gz?).
You can get something like:
Sequence_name\t0\tlength_of_read
You should use bowtie2 or bwa to map fq.gz to some reference genome (human:hg38/hg19) first to get the sam file and use bedtools to convert sam file to bed file.