Short coding test for Genomics Aotearoa Genome Graphs project position.
Generate a reproducible workflow (e.g., a shell script, R markdown file, or Jupyter notebook) that accomplishes the following:
-
Download BAM file of aligned Oxford Nanopore sequence data: BAM file
-
Output the total number of reads in the file
-
Produce a histogram of read lengths
-
Count the number of reads that are longer than 10kb
-
Calculate the average Phred quality score for each read, and output these averages for the first ten reads.
Once you have completed these tasks, please email an appropriate file (e.g., .sh, .Rmd, .ipynb) containing the required commands to carry out each step, and the output that is produced (in html or pdf format) to: [email protected]