Braker2 with filtered BAM files

Hi, I was thinking after HiSAT2 to run MarkDuplicates and remove aligned reads with lower than 40 quaily score. Then to run BRAKER2 with the filtered reads.

Did anyone try it and did it improve the annotation?

Thank you in advance,

assembly gene sequencing rna-seq

• 766 views • link

written 6.8 years ago by Ric

0 answers

No answers yet.

More posts like this

rRNA contamination after filtering
written 1.8 years ago by sleepystudent •

Hi everyone ! I am re-analyzing a dataset from [GEO][1]. I have noticed that there was a strange GC in some samples, and after pseudo-counting …
Help with proteomics downstream analysis
written 2.2 years ago by Álvaro •

Hello everyone. I have a nLC-MS/MS proteomics dataset with protein abundances for every sample (n = 250). The problem I have is that the dataset …
no unique mapped paired-end reads in BAM file
written 6.5 years ago by Ric

Hi, I would like to only unique mapped paired-end reads. Reads which are not anymore in pair I would like to remove. After the the …
only unique mapped paired-end reads
written 6.6 years ago by Ric

Hi, I would like to only unique mapped paired-end reads. Reads which are not anymore in pair I would like to remove. Are the below …
intersecting two GFF3 missing data from the second file
written 6.7 years ago by Ric

Hi, I have both StringTie (*stringtie_merged.gtf.fasta.transdecoder.genome.gff*3) and BRAKER2 (*augustus.hints_utr.gff3*) GFF3 files which I would like to intersect with bedtools. Whenever StringTie covers the genome region …
Picard MarkDuplicates flag vs remove issue
written 8.2 years ago by s1469060 •

Hi all I have "a quick question" about Picard MarkDuplicates. I have ATAC-seq data, already filtered for mitochondrial and unmapped reads. The initial file has …
map reads to a allopolyploid species reference using STAR or HISAT2: low percentage of Uniquely m…
written 8.7 years ago by clingyun •

Hello, I am mapping rna-seq reads to a reference genome. Very sure, the reads and the genome reference are generated from same variety. They are …
Picard vs samtools rmdup
written 9.3 years ago by xd_d

Hey all, I want to remove duplicates from my bam file. I use picard MarkDuplicates to remove the duplicates. (REMOVE_DUPLICATES=true) After I run picard to …
HISAT2 output statistics does not match the actual output
written 9.3 years ago by Tobias

I have been running HISAT2 for some of my FASTQ files via hisat2 -x genome --sp 1000,1000 -U /home/a.fastq.gz -S /home/output.sam It worked and I …
Is there a correct order in which to filter BAM files?
written 11.3 years ago by James Ashmore

<p>Say I have a BAM file which I want to filter, I want to apply the following processing steps: Remove reads aligned to blacklist regions, …