mflamand / Bullseye

Bullseye analysis pipeline for DART-seq analysis
MIT License
12 stars 5 forks source link

How to get cell-sites matrix in scDART-seq? #15

Open MopingXu opened 3 months ago

MopingXu commented 3 months ago

Hi Mflamand, I'm running the scDART-seq data using Bullseye's pipeline, and I have some confusion: How do I get the matrix of cells and sites when there is no barcode information in the final result of single cell analysis? Or I missed some analysis? Looking forward to your reply. Best, Mope

mflamand commented 3 months ago

Hello MopingXu,

It is difficult to see what is the problem without knowing what code you are using.

The options to use will depend on the type of sequencing used and if you have individual bam files for each cell or a merged bam containing all cells.

For the latter, you would typically use use the following options "-mode SingleCell" and "-Cell_ID_pattern X" where X is "10X" (CB:Z:), "SMART" (RG:Z) or the actual string used in the BAM file for the cell ID.

If you have distinct bam files for each cell, I would invite you to look at this repo (@tegowski) for an example: https://github.com/tegowski/scDARTHEKcells

Please let me know if you have more questions.

mflamand commented 3 months ago

Hi,

The columns of the matrix file outputted by parseBAM.pl are : “chromosome, position, #A, #T, #C, #G, #N, #Total, Cell_ID”,

Best,

On Mar 27, 2024, at 11:13 PM, mope @.***> wrote:

You don't often get email from @.*** Learn why this is important https://aka.ms/LearnAboutSenderIdentification
Avertissement automatisé : Ce courriel provient de l'extérieur de votre organisation. Ne cliquez pas sur les liens et les pièces jointes si vous ne reconnaissez pas l'expéditeur.

Hi Mflamand,

Thank you for reply. I don't understand the result of “parseBAM.pl”, what exactly does each column mean? Best

https://can01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fmflamand%2FBullseye%2Fassets%2F136780089%2Fdf6c4210-8fe1-4de7-99ff-49a614a3eb7c&data=05%7C02%7Cmathieu.flamand%40crchudequebec.ulaval.ca%7C521faa6ba9114d5c3ce608dc4ed50870%7Ccb779487921b4201825c1ed13e6f1b40%7C0%7C0%7C638471924073253195%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=rdK1UfAk0GQ6EVLbW3DDZXqjxzgqW4xrRF%2B7nxBELt0%3D&reserved=0 / PNG Image · 7 KB

— Reply to this email directly, view it on GitHub https://can01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fmflamand%2FBullseye%2Fissues%2F15%23issuecomment-2024320037&data=05%7C02%7Cmathieu.flamand%40crchudequebec.ulaval.ca%7C521faa6ba9114d5c3ce608dc4ed50870%7Ccb779487921b4201825c1ed13e6f1b40%7C0%7C0%7C638471924073263433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=eEgrMXflGFVPKqbwZzUir0omKEGX%2B0YqUZ0U7gv0Cc8%3D&reserved=0, or unsubscribe https://can01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FALHDBWDEE6FOCZQJ4BXRHGTY2ODFLAVCNFSM6AAAAABFKNTBQWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAMRUGMZDAMBTG4&data=05%7C02%7Cmathieu.flamand%40crchudequebec.ulaval.ca%7C521faa6ba9114d5c3ce608dc4ed50870%7Ccb779487921b4201825c1ed13e6f1b40%7C0%7C0%7C638471924073271703%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=gDpDQqPGqSwI0Kj0lz5zeIb%2FkgHXyrX6YeGWOuRDPkQ%3D&reserved=0. You are receiving this because you commented.