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

iCount xlsites--many reads with second-start not falling on segmentation #209

Open
ryanmarina opened this issue May 27, 2022 · 0 comments
Open

Comments

@ryanmarina
Copy link

I am attempting to determine cross-linking site with mapped reads, but certain libraries seem to give the same warning about several second starts not falling along segmentation borders when no custom segmentation file is provided (default):

Detecting cross-links...
All records in BAM file: 310737
Reads not mapped: 0
Mapped reads records (hits): 310737
Hits ignored because of low MAPQ: 0
Records used for quantification: 310737
Records with invalid randomer info in header: 310737
Records with no randomer info: 0
Ten most frequent randomers:
: 310737
There are 175499 reads with second-start not falling on segmentation. They are reported in file: /home/rmarina/scratch/iclip_data_download/beds/iCLIP_TDP-43_Brain_High_mm_beds/iCLIP_TDP-43_Brain_High_mm_skipped.bam
Saved to BED file (uniquely mapped reads): /home/rmarina/scratch/iclip_data_download/beds/iCLIP_TDP-43_Brain_High_mm_beds/iCLIP_TDP-43_Brain_High_mm_unique.bed
Saved to BED file (multi-mapped reads): /home/rmarina/scratch/iclip_data_download/beds/iCLIP_TDP-43_Brain_High_mm_beds/iCLIP_TDP-43_Brain_High_mm_multiple.bed

Can you explain what this means with no segmentation file provided, and if providing a custom segmentation reference will help recover some of these reads?

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