Sorry, title is wrong - make that SRMA...
I am trying to perform local realignment on a couple of BAM alignment files. I got stuck with GATK and tried SRMA - seems more intuitive and easier to use. However, the output BAM files I have gotten appear to be 'broken' - I think because of quality scores?
Here's what I did:
java -Xmx2048m -jar /usr/local/bin/srma/srma-0.1.15.jar I=G502sorted.bam O=G502realigned.bam R=LASV096/Mosaik/reference.fasta
samtools index G502realigned.bam
The programs complete successfully, but the resultant .bam file can't be read by FastQC and when I try to import it into Geneious I get a 'phred score above 102' error. Did SMRA somehow change my quality scores? Any way I can correct for this? I can open up the file in Tablet and it has been realigned, but I can see quality scores, etc. in Tablet.
Thanks
I am trying to perform local realignment on a couple of BAM alignment files. I got stuck with GATK and tried SRMA - seems more intuitive and easier to use. However, the output BAM files I have gotten appear to be 'broken' - I think because of quality scores?
Here's what I did:
java -Xmx2048m -jar /usr/local/bin/srma/srma-0.1.15.jar I=G502sorted.bam O=G502realigned.bam R=LASV096/Mosaik/reference.fasta
samtools index G502realigned.bam
The programs complete successfully, but the resultant .bam file can't be read by FastQC and when I try to import it into Geneious I get a 'phred score above 102' error. Did SMRA somehow change my quality scores? Any way I can correct for this? I can open up the file in Tablet and it has been realigned, but I can see quality scores, etc. in Tablet.
Thanks
Comment