Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in
  • vidjil vidjil
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 1,739
    • Issues 1,739
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 105
    • Merge requests 105
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • vidjilvidjil
  • vidjilvidjil
  • Issues
  • #1304
Closed
Open
Issue created Nov 29, 2016 by Vidjil Team@vidjilteamMaintainer

Réarrangements incomplets

Patient MAQ : en lançant Vidjil en ligne de commande ./vidjil -V germline/IGK-INTRON.fa -J germline/IGK-KDE.fa -c clones data.fastq on trouve un clone en Intron-KDE qui ressort fortement (35987 reads)


Pour le patient MAQ en DD2/DD3 : rien. Mais ça vient des germlines. Le DD2 fait 9nt et on utilise une graine qui s'étend sur 11. Il faudrait récupérer la séquence en amont du DD2 et en aval du DD3.


|13+0=13| |rev-compl| actgggggatacgcacagtgctacaaaacctacagagacctgtac

En utilisant ces séquences-là (en local) on a un clone à 44584 reads.


|


|9+0=9| |rev-compl| tatactgatgtgtttcattgtgccttcctac et cette séquence-là pour le DD3 :

M22152|TRDD3*01|Homo sapiens|F|D-REGION|214..226|13 nt|1|


|


Sur bioinfo-inria on a cette séquence-là pour le DD2 :

M22153|TRDD2*01|Homo sapiens|F|D-REGION|34..42|9 nt|1|


Pour récupérer automatiquement les séquences amonts on peut peut-être utiliser l'API du NCBI (pour le faire automatiquement). Sinon on peut les stocker quelque part :)


Pour TRDD201 : http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nuccore&id=M22153&rettype=fasta&retmode=text&to=42 Pour TRDD301 : http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nuccore&id=M22152&rettype=fasta&retmode=text&from=214&to=258


b43f8672 pour récupérer les fichiers avec les séquences amont/aval et cffcd742, 7aaa26fd pour gérer les réarrangements incomplets (merci Mathieu)


@magiraud

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information
Assignee
Assign to
Time tracking