bed file for each read

[m_elena_bioinfo]

Good morning SEQanswer community,
I have to create a bed file containing one read for lane and relative coverage, like this:

chr2 1 100 1
chr2 2 101 1
chr2 3 102 1
...
Someone can help me?!?
Thanks a lot!
ME

[dariober]

Hello,

It's bit vague what you want to achieve... See if this points you in the right direction...

If you have a SAM file you can convert it to BED (one read per line) with samtools and bedtools like this:

## Convert SAM file myaln.sam to BAM:
samtools view -b -S -o myaln.bam myaln.sam

## Convert BAM to BED:
bamToBed -i myaln.bam > myaln.bed

If you want to count how many reads map on each gene or transcript of a GTF file:

## Coverage of GTF features (for each gene/transcript, how many reads overlap it?):
coverageBed -abam myaln.bam -b myspecies.gtf > myaln.cov

See also samtools pileup to have the read count at each genome position (not sure if this is what you mean).

Good luck!
Dario

[m_elena_bioinfo]

Thank you very much Dario, it works fine!!!
It's perfect!!