smithlabcode / preseq

Software for predicting library complexity and genome coverage in high-throughput sequencing.
https://preseq.readthedocs.io
GNU General Public License v3.0
78 stars 16 forks source link

Question on lc_extrap input #39

Closed stevenlockton-progenity closed 5 years ago

stevenlockton-progenity commented 5 years ago

For the lc_extrap function, the documentation says you can use BED or BAM as inputs. No mention of WALT mr files. But mr files seem to work just fine. So can I confirm this: Are mr files a valid input for preseq lc_extrap?

timydaley commented 5 years ago

From my understanding, mr files are an extension of bed files but with extra columns. Indeed, we have a tool to convert paired end bam files to mr files, so preseq can use them. They should work just fine as long as the first columns match the bed format.