User contributions for WikiSysop

A user with 79 edits. Account created on 19 March 2024.
Jump to navigation Jump to search
Search for contributionsExpandCollapse
⧼contribs-top⧽
⧼contribs-date⧽
(newest | oldest) View (newer 50 | ) (20 | 50 | 100 | 250 | 500)

19 March 2024

  • 18:1018:10, 19 March 2024 diff hist +11,253 N Program 2019Created page with "'''REMEMBER TO BRING A LAPTOP COMPUTER FOR EXERCISES''' Lectures and exercises will take place at the Section for Bioinformatics at the Technical University of Denmark, '''Building 208, Room 062''' in Lyngby ([https://www.google.com/maps/place/DTU+bygning+208%2F+DTU+Building+208/@55.7863894,12.520095,357m/data=!3m1!1e3!4m5!3m4!1s0x46524e62f280d429:0x1bbb0474519984b5!8m2!3d55.787537!4d12.5182208 map]). The courses has two main parts, the first half is lectures and exerc..." current
  • 18:0918:09, 19 March 2024 diff hist +2,502 N Metagenomic assembly solutionCreated page with "Q1. We cant really find the bell shaped distributions in our samples - except for in MH0032 where there are two very small bell shaped coverage distributions. Q2. This is because we have many organisms with relative low abundance. This makes it very hard to distinguish them using coverage information. Q3. It will probably perform much like SOAPdenovo. Q4. There arent really any major differences in the coverage distributions between the assemblers, the metagenome asse..." current
  • 18:0918:09, 19 March 2024 diff hist +23,886 N Metagenomic assembly exerciseCreated page with "<H2>Overview</H2> <p>In this exercise we will try <i>de novo</i> assemble a metagenomic dataset of Illumina paired end reads from a stool sample. The data is part of the MetaHit project and was published [http://www.nature.com/nature/journal/v464/n7285/full/nature08821.html here]. These are just two samples of 396, in the next exercise (tomorrow) you will analyze data from 124 samples. Today, you will try to: <OL> <LI>Analyze k-mer distributions <LI>Assemble using SOA..." current
  • 18:0718:07, 19 March 2024 diff hist +11,508 N Program 2020Created page with " '''REMEMBER TO BRING A LAPTOP COMPUTER FOR EXERCISES''' Lectures and exercises will take place at the Section for Bioinformatics at the Technical University of Denmark, '''Building 210, Room H162''' in Lyngby ([https://goo.gl/maps/hyifvFrZ5cPLroCv8 map]). The courses has two main parts, the first half is lectures and exercises and the last half is project work ending with the exam on Friday 24th of January. This term we will be using Piazza for class discussion. The..." current
  • 18:0618:06, 19 March 2024 diff hist +12,550 N Program 2021Created page with " <!--- '''REMEMBER TO BRING A LAPTOP COMPUTER FOR EXERCISES''' --> Lectures and exercises will take place on Discord (https://discord.gg/JxV3pHyHgV). Please register with your full name. Will use Discord for online classes and collaboration with your project partners. The course has two main parts, the first half is lectures and exercises and the last half is project work ending with the exam on Friday 22th of January. This term, we will be using Piazza for class dis..." current
  • 18:0518:05, 19 March 2024 diff hist +1,343 N CfDNA exercise answersCreated page with "'''Q1''' Using: <pre> samtools stat /data/shared/exercises/cfdna/patient_1.bam > patient1.stat samtools stat /data/shared/exercises/cfdna/patient_2.bam > patient2.stat plot-bamstats -p patient1 patient1.stat plot-bamstats -p patient2 patient2.stat firefox patient1.stat firefox patient2.stat </pre> Both peak at 168bp however patient 1) clearly has the greatest variance and has an overabundance of short DNA fragments compared to patient 2. '''Q2''' We first run: <p..." current
  • 18:0518:05, 19 March 2024 diff hist +3,363 N CfDNA exerciseCreated page with "<H2>Overview</H2> First: <OL> <LI>Navigate to your home directory: <LI>Create a directory called "cfdna" <LI>Navigate to the directory you just created. </OL> Blood was drawn from 2 patients. One of those is a healthy patient, the other one has lung squamous cell carcinoma. It was previously [https://www.cell.com/fulltext/S0092-8674(15)01569-X reported] that tumors can leave a signature in the blood plasma via cell-free DNA, namely: # The distribution of fragment l..." current
  • 18:0418:04, 19 March 2024 diff hist +4,574 N Logging on to Cloud systemCreated page with " <HR> <H2>Overview</H2> In this exercise, we will prepare our computers to log on to our very own cloud system on Computerome. This is done via a Virtual Private Network (VPN). <HR> <H3>Student accounts</H3> <p>To be able to perform the exercises and project you will need an account to log in to the cloud computers. Open the [https://docs.google.com/spreadsheets/d/1TAj59JQurp6ene6kOlrug9mRJSHpdsIS8U6KKvkcLN0/edit?usp=sharing google doc], find student id and determin..." current
  • 18:0318:03, 19 March 2024 diff hist +12,968 N Program 2022Created page with " '''REMEMBER TO BRING A LAPTOP COMPUTER FOR EXERCISES''' Lectures and exercises will take place on Discord (https://discord.gg/FBb2edFW). Please register with your full name. Will use Discord for online classes and collaboration with your project partners. The course has two main parts, the first half is lectures and exercises and the last half is project work ending with the exam on '''Friday 21st of January'''. <!-- This term, we will be using Piazza for class disc..." current
  • 18:0218:02, 19 March 2024 diff hist +12,699 N Program 2023Created page with "'''REMEMBER TO BRING A LAPTOP COMPUTER FOR EXERCISES''' Lectures will be in person in building [https://goo.gl/maps/k4wYkMjTJ2HLHuyN8 303A] in auditorium 45. Offline discussions will take place on Discord (https://discord.gg/4yVB6vMG2n). Please register with your '''full name'''. Will use Discord for online classes and collaboration with your project partners. <!-- Lectures and exercises will take place on Discord (https://discord.gg/FBb2edFW). Please register with..." current
  • 18:0118:01, 19 March 2024 diff hist +132 N MediaWiki:SidebarCreated page with " * navigation ** https://teaching.healthtech.dtu.dk/|Course List ** https://teaching.healthtech.dtu.dk/22126/|Course 22126 * TOOLBOX" current
  • 18:0018:00, 19 March 2024 diff hist +4,071 N Cancerseq exercise answersCreated page with "'''Q1''' We run: <pre> gatk Mutect2 -R /home/databases/references/human/GRCh38_full_analysis_set_plus_decoy_hla.fa -I /home/projects/22126_NGS/exercises/cancer_seq/TCRBOA2-N-WEX_recaled.bam -I /home/projects/22126_NGS/exercises/cancer_seq/TCRBOA2-T-WEX_recaled.bam -normal TCRBOA2-N-WEX -L chr10:3100000-5100000 --germline-resource /home/databases/databases/GRCh38/somatic-hg38_af-only-gnomad.hg38.vcf.gz -O TCRBOA2.vcf.gz </pre> Then either: <pre>..." current
  • 17:5917:59, 19 March 2024 diff hist −8 Cancerseq exerciseNo edit summary current
  • 17:5817:58, 19 March 2024 diff hist +9,763 N Cancerseq exerciseCreated page with " <H2>Overview</H2> Adapted from an original exercise by Marcin Krzystanek and Aron Eklund. First: <OL> <LI>Navigate to your home directory: <LI>Create a directory called "cancerseq" <LI>Navigate to the directory you just created. </OL> In this exercise, we learn the differences between standard genotyping and genotyping for cancer-seq: <OL> <LI>Somatic mutation calling <LI>Interpretation of the resulting somatic mutations </OL> <H2>Somatic mutation calling</..."
  • 17:5617:56, 19 March 2024 diff hist +4,682 N Genomic epidemiology solutionCreated page with "Q1. We can use the fastx tools to get this information <pre> fastx_readlength.py --i reads_R1.fastq.gz --gz fastx_readlength.py --i reads_R2.fastq.gz --gz </pre> This should give you the answer. Reads: 2*1408847 = 2,817,694; Bases: 284587094 Task1 - can be accomplished using a quality control tool like fastqc. <pre> mkdir fastqc fastqc -o fastqc reads_R*.fastq.gz firefox fastqc/reads_R1_fastqc.html fastqc/reads_R2_fastqc.html & </pre> Q2. There are no overrepresented..." current
  • 17:5617:56, 19 March 2024 diff hist +15,177 N Genomic epidemiology exerciseCreated page with "<H3>Overview</H3> You are a bioinformatician working with epidemiology in Copenhagen. In the last couple of days more and more people have been admitted to the Hospitals in the region with severe bacteremia (bacteria in the bloodstream) and several patients, especially infants and young children, have had fatal outcomes. Luckily the probable causative infectious bacteria has been isolated and sequenced on a desktop sequencer and now it is your job to find out: <OL> <LI..." current
  • 17:5417:54, 19 March 2024 diff hist 0 N File:Pca.pngNo edit summary current
  • 17:5417:54, 19 March 2024 diff hist 0 N File:Bc heatmap.pngNo edit summary current
  • 17:5317:53, 19 March 2024 diff hist 0 N File:Wilcoxon significant.pngNo edit summary current
  • 17:5317:53, 19 March 2024 diff hist 0 N File:Sample depth.pngNo edit summary current
  • 17:5317:53, 19 March 2024 diff hist +2,619 N QuantitativeMetagenomicsSolutionCreated page with "'''Q1. How many samples do we have and how many genes?''' <pre>> str(Counts) int [1:251436, 1:401] </pre> 251436 genes and 401 individuals '''Q2. What's the sample depth range?''' <pre>> range(sampleDepth) [1] 1533776 41391478 </pre> 1533776 to 41391478 File:Sample_depth.png The figure shows the number of samples (persons) on the y-axis containing the displayed number of reads on the x-axis '''Q3. How many species are there in total?''' <pre> str(taxCounts)..." current
  • 17:5217:52, 19 March 2024 diff hist 0 N File:Log fold diff sign.pngNo edit summary current
  • 17:5217:52, 19 March 2024 diff hist 0 N File:Bmi class.pngNo edit summary current
  • 17:5117:51, 19 March 2024 diff hist 0 N File:Comparing species abundance.pngNo edit summary current
  • 17:5117:51, 19 March 2024 diff hist 0 N File:Comparing shannon.pngNo edit summary current
  • 17:5017:50, 19 March 2024 diff hist 0 N File:Comparing richness.pngNo edit summary current
  • 17:5017:50, 19 March 2024 diff hist 0 N File:Downsized shannon.pngNo edit summary current
  • 17:4917:49, 19 March 2024 diff hist 0 N File:Downsized richness.pngNo edit summary current
  • 17:4817:48, 19 March 2024 diff hist 0 N File:Downsized sampledepth.pngNo edit summary current
  • 17:4817:48, 19 March 2024 diff hist 0 N File:Raw sampledepth.pngNo edit summary current
  • 17:4717:47, 19 March 2024 diff hist 0 N File:Raw richnessVSshannonZoom.pngNo edit summary current
  • 17:4717:47, 19 March 2024 diff hist 0 N File:Raw shannon.pngNo edit summary current
  • 17:4617:46, 19 March 2024 diff hist 0 N File:Raw richness.pngNo edit summary current
  • 17:4517:45, 19 March 2024 diff hist +10,393 N QuantitativeMetagenomicsCreated page with " <H3>Overview</H3> If you need to use metagenomics for your final project, we have a more thorough workflow that you will need to use https://teaching.healthtech.dtu.dk/22136/index.php/22136:Course_plan_autumn_2020 here. Since metagenomics data is often very large, it requires a lot of computational resources and time, we have cheated a little bit and prepared some data for you in advance! In this exercise we have done the assembly and counting across a cohort of..." current
  • 17:4417:44, 19 March 2024 diff hist 0 N File:DDseq 350OTU.pngNo edit summary current
  • 17:4417:44, 19 March 2024 diff hist 0 N File:DDseq 100OTU.pngNo edit summary current
  • 17:4417:44, 19 March 2024 diff hist 0 N File:PrincipalComponents.pngNo edit summary current
  • 17:4317:43, 19 March 2024 diff hist 0 N File:Rawkaijubar.pngNo edit summary current
  • 17:4317:43, 19 March 2024 diff hist +3,019 N Kaiju solutionCreated page with "<b> Q1: What is nr_euk? And how do the choice of database influence the results of Kaiju?</b> nr stands for non-redundant and indicates, that each entry is only found once within the database. BLASTS' nr database contains bacteria, archea and fungi. euk indicates that microbial eukaryotes and fungi are included. The database should be intended as the target of the search, so it is of importance that it contains the organisms you are searching for. <b> Q2: Explain the..." current
  • 17:4217:42, 19 March 2024 diff hist +21,585 N Kaiju exerciseCreated page with "<!-- =22136 Metagenomics and Microbiome Analysis Kaiju exercise= --> ==Introduction== [http://kaiju.binf.ku.dk/ Kaiju] is a protein-based sensitive taxonomic classification of high-throughput sequencing reads from metagenomic whole genome sequencing or metatranscriptomics experiments. <br> Kaiju translates metagenomic sequencing reads into the six possible reading frames and searches for maximum exact matches (MEMs) of amino acid sequences in a given database of annotat..." current
  • 17:4117:41, 19 March 2024 diff hist +2,972 N Longread exercise answersCreated page with "'''Q1''' Counting all the lines minus the header gives us: <pre> zcat BGI_hg38_chr20.vcf.gz |grep -v "^#"|wc -l </pre> 1878 variants '''Q2''' We can try the following: <pre> zcat BGI_hg38_chr20.vcf.gz |grep -v "^#" |cut -f 10 |sed "s/:.*//g"|sort | uniq -c |sort -n </pre> <OL> <LI>grep -v "^#" grep -v: Inverts the match, i.e., selects lines that do not match the given pattern. "^#": The pattern to match lines starting with a hash (#). These lines are usually he..." current
  • 17:4117:41, 19 March 2024 diff hist +7,233 N Longread exerciseCreated page with "<H2>Overview</H2> First: <OL> <LI>Navigate to your home directory: <LI>Create a directory called "longread" <LI>Navigate to the directory you just created. </OL> We will phase some variants using [https://www.biorxiv.org/content/10.1101/085050v2 WhatsHap] (no not the messaging app). First, what is phasing? Phasing means that we determine which base is on the same chromosome as another base for neighboring variants. Let's consider a small example with just two varia..." current
  • 17:3917:39, 19 March 2024 diff hist 0 N File:Rnaseq fig3.pngNo edit summary current
  • 17:3917:39, 19 March 2024 diff hist 0 N File:Rnaseq fig2.pngNo edit summary current
  • 17:3917:39, 19 March 2024 diff hist 0 N File:Rnaseq fig1.pngNo edit summary current
  • 17:3717:37, 19 March 2024 diff hist +18,721 N Rnaseq exercise answersCreated page with " <div class="page-content has-page-title"> <div id="overview-and-background" class="section level1"> <h1>Overview and background</h1> <div id="groups" class="section level2"> <h2>Groups</h2> <p>Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.</p> </div> <div id="assignment-notes" class="section level2"> <h2>Assignment notes</h2> <p>While some question..." current
  • 17:3717:37, 19 March 2024 diff hist +11,080 N Rnaseq exerciseCreated page with " <div class="page-content has-page-title"> <div id="overview-and-background" class="section level1"> <h1>Overview and background</h1> <div id="groups" class="section level2"> <h2>Groups</h2> <p>Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.</p> </div> <div id="assignment-notes" class="section level2"> <h2>Assignment notes</h2> <p>While some question..." current
  • 17:3617:36, 19 March 2024 diff hist +2,833 N Ancient DNA exercise answersCreated page with "'''Q1''' the read length is about 100bp but the actual insert size is unknown. '''Q2''' very low, less than 1% '''Q3''' About 40bp. '''Q4''' About 25%. '''Q5''' As and Gs '''Q6''' The sample indeed looks ancient. If we did not see DNA fragmentation or damage it could be indicative of present-day human contamination. '''Q7''' <pre> wc -l world.fam wc -l world.bim </pre> 297 samples and 587772 SNPs. '''Q8''' <pre> cut -f2 world.sampleInfo.txt | tail -n +2..." current
  • 17:3517:35, 19 March 2024 diff hist +11,971 N Ancient DNA exerciseCreated page with "<H2>Overview</H2> Adapted from Martin Sikora. First: <OL> <LI>Navigate to your home directory: <LI>Create a directory called "adna" <LI>Navigate to the directory you just created. </OL> We will try to # Authenticate ancient DNA # do some basic population genetics <h2> Data authentication</h2> Authentication involves making sure that the DNA that you have extracted from my fossil and sequenced is indeed from the fossil and not some modern contaminant. A big differe..." current
  • 17:3517:35, 19 March 2024 diff hist +1,480 N Denovo solutionCreated page with "Q1. Illumina Q1A. discarded contains reads that are too short, pair1 and pair2 files contain the read pairs were both passed trimming and singleton are reads were one of the two pairs were discarded. Q2. Around 84 Q3. N = (M*L)/(L-K+1) = (84*99)/(99-15+1) = 97.84 Genome_size = T/N = (213721367+212523694)/97.84 = 4.35Mb Q4. Mean = 259 ; SD = 11 Q5. It is lower, this means that the actual kmer peak we found (unless you found one higher than 84) is higher (this would g..." current
(newest | oldest) View (newer 50 | ) (20 | 50 | 100 | 250 | 500)