maiziezhoulab / RegionIndel

Large indel detection in region-based phased diploid assemblies from linked-reads
MIT License
0 stars 0 forks source link

human data only? #3

Open pdimens opened 3 months ago

pdimens commented 3 months ago

Hello,

I see that the chr_num parameter excepts only an integer and that the code then expects the contigs to be named chr_. Could this not be generalized to all species if chr_num was just chr as a string and the code then used >chr (etc.) to pull out the relevant contig name?

maiziex commented 2 months ago

Yes, you can. Ensure that the chromosome names in the reference genome FASTA file follow this format.

maiziex commented 2 months ago

Yes, you can. Ensure that the chromosome names in the reference genome FASTA file follow this format.

- Maizie (Xin) Zhou, Ph.D. Assistant Professor Biomedical Engineering, Computer Science, and Data Science Institute Vanderbilt University

5919 Stevenson Center

Nashville, TN 37235

Cell: 336-473-3540

https://lab.vanderbilt.edu/maizie-zhou-lab/


From: Pavel V. Dimens @.> Sent: Friday, July 5, 2024 8:27 AM To: maiziezhoulab/RegionIndel @.> Cc: Subscribed @.***> Subject: [maiziezhoulab/RegionIndel] human data only? (Issue #3)

Hello,

I see that the chrnum parameter excepts only an integer and that the code then expects the contigs to be named chr. Could this not be generalized to all species if chr_num was just chr as a string and the code then used >chr (etc.) to pull out the relevant contig name?

— Reply to this email directly, view it on GitHubhttps://github.com/maiziezhoulab/RegionIndel/issues/3, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABMOIQ2NMYQ6G7UG452LBV3ZK2NK3AVCNFSM6AAAAABKNGXBA2VHI2DSMVQWIX3LMV43ASLTON2WKOZSGM4TENRRGE2DAMI. You are receiving this because you are subscribed to this thread.Message ID: @.***>

pdimens commented 2 months ago

Would it not be possible to modify the logic such that the user inputs the full contig/chromosome name rather than forcing Contig names to be chr_ ?

maiziex commented 2 months ago

Unfortunately, some parts of the code logic need to be modified to accommodate this. The easiest approach for now is to modify the contig name. We plan to make this applicable to a general case in a future version.

- Maizie (Xin) Zhou, Ph.D. Assistant Professor Biomedical Engineering, Computer Science, and Data Science Institute Vanderbilt University

5919 Stevenson Center

Nashville, TN 37235

Cell: 336-473-3540

https://lab.vanderbilt.edu/maizie-zhou-lab/


From: Pavel V. Dimens @.> Sent: Tuesday, July 9, 2024 6:48 AM To: maiziezhoulab/RegionIndel @.> Cc: Zhou, Maizie @.>; Comment @.> Subject: Re: [maiziezhoulab/RegionIndel] human data only? (Issue #3)

Would it not be possible to modify the logic such that the user inputs the full contig/chromosome name rather than forcing Contig names to be chr_ ?

— Reply to this email directly, view it on GitHubhttps://github.com/maiziezhoulab/RegionIndel/issues/3#issuecomment-2217436967, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABMOIQ7XLCI6CPL6PSB6CTTZLPEZBAVCNFSM6AAAAABKNGXBA2VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJXGQZTMOJWG4. You are receiving this because you commented.Message ID: @.***>