Counting starting and ending positions of reads in reference genome
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.
• 1,174 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
bobbyle0210 •
•
written
by
bobbyle0210 •
1 answer
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Strange offsetting in deeptools?
written by Elijah •
-
How to make bam from start pos of read 1 and end pos of read 2 from paired end bam files
written by genomics_student •
-
Get upstream and downstream exon's starting site and ending site of an retained introns
written by Kai_Qi
-
Assembling Short Sequence
written by joseph.landry •
-
How to extract a peptide sequence (interval) from a FASTA file?
written by Peter •
-
bedtools not working properly
written by wstla27 •
-
changing the position of a column from a file to another in perl
written by rijumukherjee79800 •
-
Retrieve the sequence based on the start and end position in the cuffmerged.gtf
written by Chao.wang2 •
-
What to do when reads pairs are almost same,like 90-100 % overlap each other in paired reads?
written by crivenster
-
Inferring Amplicons From Bam/Sam Data
written by Nicholas Mancuso
