Question: Counting starting and ending positions of reads in reference genome
0
gravatar for bobbyle0210
9 months ago by
bobbyle021010
bobbyle021010 wrote:

Hi, Currently, I am working with NGS processing data. After getting SAM/BAM file, now, I want to output a text file that indicates the starting (5 end) and ending (3 end) position of each read in the corresponding reference sequence as well as the sequence of the read. The output file may look like:


ref_seq_1 starting_point read_1_sequence ending_point


ref_seq_1 starting_point read_2_sequence ending_point


ref_seq_2 starting_point read_3_sequence ending_point


I guess that I may need to use bedtools but I am still not sure which command should I use to obtain such information. Could you guys give me some suggestions? Thank you so much.

rna-seq sam bam ngs bedtools • 220 views
ADD COMMENTlink modified 9 months ago by Pierre Lindenbaum129k • written 9 months ago by bobbyle021010
0
gravatar for Pierre Lindenbaum
9 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

https://bedtools.readthedocs.io/en/latest/content/tools/bamtobed.html

ADD COMMENTlink written 9 months ago by Pierre Lindenbaum129k
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: 1347 users visited in the last hour