Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

about to-mr function #72

Open
SYSUMSD opened this issue Apr 7, 2024 · 1 comment
Open

about to-mr function #72

SYSUMSD opened this issue Apr 7, 2024 · 1 comment

Comments

@SYSUMSD
Copy link

SYSUMSD commented Apr 7, 2024

Hi. I am using gc_extrap to process some BAM files, and my version of Preseq is 3.2.0. According to your usage examples, BAM files need to be converted to MR files. However, I found that Preseq version 3.2.0 does not include a function to perform this conversion. While searching other websites, I discovered that Preseq previously had functions such as BAM2MR or to-mr, but these are not present in version 3.2.0. Currently, Preseq includes functions like c_curve, lc_extrap, gc_extrap, bound_pop, and pop_size. I am using Bedtools to convert these BAM files to BED files and then processing the BED files with the command 'preseq gc_extrap -B -o test1.txt test.bed' instead of using MR files. Can I run gc_extrap like this?

@SYSUMSD
Copy link
Author

SYSUMSD commented Apr 7, 2024

I have run the example data in your data directory using samtools and gc_extrap.
This is my code :
samtools sort -@ 30 -o SRR1106616_5M_subset.sort.bam SRR1106616_5M_subset.bam
bedtools bamtobed -i SRR1106616_5M_subset.sort.bam > SRR1106616_5M_subset.sort.bed
preseq gc_extrap -B -o testbed.txt SRR1106616_5M_subset.sort.bed.
But I got a error:
ERROR: reads unsorted in: SRR1106616_5M_subset.sort.bed
prev = chr1 9807080 9807090 SRR1106616.105083640 255 +
curr = chr1 9790490 9790500 SRR1106616.160558132 255 +
Increase seg_len if in paired end mode

Is there a mistake in my process?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant