Unusual indel peaks in bam file
0
1
Entering edit mode
4 months ago
boczniak767 ▴ 870

HI,

I'm reanalysing project GSE124340 (RNA-seq).

Quality control of input fastq's with FastQC haven't detected any problems however after alignment I've got indel frequencies in BamQC like this.

indel

I've used STAR aligner: the command used in loop

~/bin/STAR_2.7.10b/Linux_x86_64/STAR --runThreadN 24 --readFilesIn ${srr}_1.fastq.gz ${srr}_2.fastq.gz \
--genomeDir /media/mj/bb71e88b-2122-487c-a602-a2552c680021/star-index --outSAMtype BAM SortedByCoordinate \
--outFileNamePrefix ${srr} --readFilesCommand zcat`

So I've used default settings. Also The data is from the reference genome (B73 reads aligned to B73 genome) and I've not seen such peaks in several other project I've analysed.

Should I try to remove this effect (how?) or it is safe to procceed?

RNA-seq • 249 views
ADD COMMENT

Login before adding your answer.

Traffic: 1927 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