Yesterday’s attempt at producing a bedgraph was a failure and a prodcuct of a major brain fart. The worst part is that I was questioning what I was doing the entire time, but still went through with the process! Yeesh!
The problem was that I tried to take our Trinity-assembled transcriptome and somehow align that to our genome. This can’t work because each of those assemblies don’t know the coordinates used by the other. So, as was the case, you end up with a bedgraph that shows zero coverage for all genome contigs.
Anyway, here’s the correct procedure!
Used HISAT2 on our HPC Mox node to align our RNAseq reads to our Olurida_v081 genome assembly:
SBATCH script files:
PERFORM GENOME INDEXING & ALIGNMENT
– 20180925_oly_RNAseq_genome_hisat2.sh
SORT & INDEX ALIGNMENT OUTPUT
– 20180925_oly_RNAseq_genome_sort_bam.sh
RESULTS
Output folder:
Sorted BAM file (58GB):
Will get the sorted BAM file converted to a bedgraph showing genome coverage for use in IGV.