Hi everyone !
This is my first time using Cufflinks.
I'm trying to get a fasta corresponding to the RNA-seq from different tissues from D. suzukii. I first used hisat2 in order to align my reads against a reference, here the command I used in Hisat2 :
hisat2 -p 25 --max-intronlen 50000 -x /media/DATAPART1/Documents/suzukii_assembly/annotation/indexes/suz_canu3 -1 '/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_ant.R1.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_ovi.R1.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_pro.R1.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_tar.R1.fastq.gz' -2 '/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_ant.R2.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_ovi.R2.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_pro.R2.fastq.gz','/media/DATAPART1/Documents/suzukii_assembly/annotation/EST_evidences/suz_tar.R2.fastq.gz' | samtools view -Sbo suzukii_rnaseq_ibdm.bam -
Then I sorted the bam file using this command :
samtools sort /media/DATAPART1/Documents/suzukii_assembly/annotation/suzukii_rnaseq_ibdm.bam /media/DATAPART1/Documents/suzukii_assembly/annotation/suzukii_rnaseq_ibdm_sorted.bam
Something that bothered me at this step is, first it took a long time but well, and second, the sorted file was shorter than the original file (37.2G -> 26.2G), and at this point I don't know if this was normal, but I tried anyway to launch Cufflinks using this sorted file.
At first, I used this command :
cufflinks --library-type fr-firststrand -o /media/DATAPART1/Documents/suzukii_assembly/annotation/cufflinks/hisat2 -p 25 /media/DATAPART1/Documents/suzukii_assembly/annotation/suzukii_rnaseq_ibdm_sorted.bam
But I get strange errors like this :
BAM record error: found spliced alignment without XS attribute
I read a post on Biostar (don't remember which one), and the problem could be the --library-type, so I tried with this option on :
And then it seems to work for a moment, but it froze at 14% of the sequence (I tried 2 times) :
You are using Cufflinks v2.2.1, which is the most recent release. [17:08:27] Inspecting reads and determining fragment length distribution. Processing Locus tig00000383:187849-193055 [*** ] 14%
Cufflinks didn't stop, it didn't throw any error, it is still doing stuff apparently (I'm watching it with htop), but it run all nigh, and I think that something went wrong...
So, is there something wrong with this particular locus ? Is the sorting step a problem ? Or I've just used the wrong library type ?
Thanks for your help !