This is a test version of Biostars. For the public version, visit https://www.biostars.org.
converting SAM file to SRA

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 )

sra sam

0 answers

No answers yet.

Log in to answer this question.