Does Anybody know how to convert SAM files to equivalent SRA ? I downloaded a SRA file and used sam-dump to convert it to sam format. Thereafter, I used bam-load to convert sam file to SRA. The reference fasta and configuration files are specified for bam-load. The output file is valid. however when I use vdb-diff to compare it with the original SRA, it shows they are different
Here is an example: 1) prefetch SRR11433881 -o SRR11433881.sra 2) sam-dump -u SRR11433881.sra > SRR11433881.sam 3) vdb-dump SRR11433881.sra -T REF -C NAME,SEQ_ID -f tab | uniq > config.txt 4) dump-ref-fasta SRR11433881.sra > ref.fasta 5) bam-load -r ref.fasta -k config.txt -o SRR11433881_out.sra SRR11433881.sam 6) vdb-diff SRR11433881_out.sra SRR11433881.sra
src[ 1 ] : SRR11433881.sam src[ 2 ] : SRR11433881_out.sra
- rows : all
- progress : hide
- intersect: no
- max err : 1
- col-by-col: no
2021-10-30T01:08:35 vdb-diff.2.11.2 int: file system table of contents incomplete while parsing archive file within file system module - accession #1 cannot be opened as table or database 0 differences discovered ( rc = 808124559 )
0 answers
No answers yet.
Log in to answer this question.