Closed damientully closed 7 years ago
Hi Damien,
That's definitely my bad; it sounds like read quality information is missing for your read(s), I didn't know that was possible. Are you able to make a new bam or sam file containing just one read (or a few reads) to send to me for investigation? Running
samtools view -h NameOfDamiensBam.bam | head > JustAFewReadsForChris.sam
from the command line should do the trick. NB that will produce mapped reads in sam format (uncompressed bam), which you can open with a plain text editor (and e.g. change the bases in the sequence, if the sequence itself is sensitive information). You can copy-paste the contents of that sam file here, or email it to me...
Cheers
Just emailed you a sample file! Thanks, D
Pull the code and try again :)
Seems to be working! Thanks.
Hi Chris,
I am having this problem running a few bam files and get the following error message:
Any suggestions on what the source of this error could be?
Many thanks, Damien