low unique mapped reads from STAR
0
0
Entering edit mode
3.0 years ago
dorentinah • 0

Hello,

I am having low unique mapped read. I am a beginner and I might have done any mistake.

So first I downloaded a reference genome file and annotation file from ensemble.

https://drive.google.com/drive/folders/1goBTntkXcH3sKGAkrK-piSVII5FdRcUa

then I did genome indexing with the following code

STAR --runThreadN 64 --runMode genomeGenerate --genomeDir STAR_folder --genomeFastaFiles Mus_musculus.GRCm39.dna.toplevel.fa --sjdbGTFfile Mus_musculus.GRCm39.103.gtf

Next the alignment was done for this sample SRR1552444.fastq.gz

With this code:

STAR --genomeDir STAR_folder/ --readFilesIn SRR1552444.fastq.gz --readFilesCommand zcat --outFileNamePrefix STAR_Align/444_ --outSAMtype BAM SortedByCoordinate --runThreadN 16 --sjdbGTFfile Mus_musculus.GRCm39.103.gtf

The problem was that the unique mapped reads were only 15% and from literature it is a least 75 %

Do you know where the mistake might be?

Best Regards,

sequencing rnaseq Alignment star • 550 views
ADD COMMENT

Login before adding your answer.

Traffic: 2861 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6