Failed to replicate results SRX175350
parveshbarak opened this issue · 0 comments
parveshbarak commented
I downloaded SRR641640 run of SRX175350 as a fastq file from NCBI website.
I then used the below set of commands to convert the fastq file to pat format:
- Index the hg19.fa file
bwa index -a bwtsw hg19.fa
- Generate a .sam file using reference and input fasta/fastq files:
bwa mem reference_genome.fasta input.fastq > output.sam
- Convert the output.sam file to sorted bam file:
samtools sort -o sorted_output.bam output.sam
- Index the sorted_output.bam file:
samtools index sorted_output.bam
- Convert the .bam file to .pat.gz format:
wgbstools bam2pat sorted_output.bam
The pat file I got should ideally contain all the entries of the pat file provided in the tutorial section of UXM tools for SRX175350. But that was not the case and therefore I also got different results on deconvoluting the pat file.
Attaching the given pat file and output (files starting with name Given_) in the tutorial for SRX175350 and the ones I got (files starting with MyExp_) with SRR641640 run of SRX175350
Given_Lung_STL002.pat.gz
Given_Lung_STL002.pdf
MyExp_Lung_STLOO2.pdf
MyExp_Lung_STLOO2.pat.gz