correct fastq files for the expression analysis
I am planning to perform expression analysis on RNAseq data but instead of fastq files, I received cram files. I used samtools to convert cram to fastq files. if you have experience, would you please let me know if the resulting fastq would have any differences with the original fastq files (I do not have access to the original fastq files). here is the command I used:
samtools view -b -T test.cram > test.bam
samtools fastq -1 test_R1.fastq.gz -2 test_R2.fastq.gz test.bam
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Questions about a bug when transferring cram file to bam file
written by me •
-
Read Classification from ONT
written by maevalefeuvre •
-
Different number of reads when converting data from FASTQ to BAM and CRAM to FASTQ
written by ManuelDB
-
FASTQ to BAM to CRAM to FASTQ
written by ManuelDB
-
Merging CRAM files
written by Matteo Ungaro •
-
No differentially expressed genes after multiple testing correction in mice
written by Sara •
-
bedtools bamtofastq generates fastq which cannot be gzip correctly?
written by francois •
-
samtools fastq compression
written by bsmith030465
-
bam to fastq conversion, number of reads do not match
written by Anushka •
-
Different number of total reads in different alignments using the same original readset
written by malteherold •
You should name sort (
samtools sort -norsamtools collate) the BAM files before converting to fastq. If your CRAM file did not include unmapped reads then you will not be able to get those. Fastq data you end up with should be identical to information in your CRAM file.Since cram is reference based format, my understanding is that user must know the reference file used in generating cram and the same should be supplied in generating bam.
Since OP has
-Tin their command they must have provided the reference when doing the actual conversion.Two additional caveats:
Original CRAM files could have contained unaligned fastq data and will not need a reference : Is it possible to directly convert fastq to CRAM ?
samtoolsis able to find a reference based on UR field (if that information exists) : Converting CRAM to BAM without reference fastaI would not create cram files without reference (First caveat) nor do I recommend it. If service/cram provider does that, no comments on that.
Second caveat works only if reference file is located as per UR field, on OP machine/network path/public URL. Since it is from a third party, I think samtools may not be able to find the reference. Otherwise, user has to provide by
-T.If OP could post UR field from cram file, it would help better understanding the location of reference file used in cram file.