Hello everybody,
I have two fastq files (reverse and forward) with an average of 150pb per read from a bacterial comunity (Illumina HiSeq). With those files I've assembled with metaVelvet. In order to get more restrictive information about the coverage of some extracted genomes from a metagenome, I want to map the reads to the assembly but only those which are 95% or higher similar, i.g keep only the reads that have: 0.95 >= match / (mismatch + match). I've read this: https://zombieprocess.wordpress.com/2013/05/21/calculating-percent-identity-from-sam-files/ but since I have a lot of reads ( more than 30 millions) I was hopping that there was a more automatic way to do this.
I appreciate if you could help me.
E.O.
0 answers
No answers yet.
Log in to answer this question.
-
Difficulties to assess bulk RNA-seq results
written by charles.chapus •
-
merging of two assemblies from different short-reads techniques
written by Andrew Yu •
-
filtering SAM/BAM to remove hits spanning short combined alignment lengths and low counts
written by Bryan •
-
BBmap params for accurate MAG relative abundance estimation?
written by Bryan •
-
Differing number of reads for read 1 and 2 in fastq's from a subset bam file
written by michael •
-
bcl2fastq with xGen Dual Index UMI Adapters to produce 3 read and 2 index fastqs
written by segrossk •
-
Forum: Calling Single-Barcode Samples from Mixed Runs as Dual-Barcode Samples | Possible Illumina R…
written by Charles Warden
-
How to segregate metagenomic sequences in prokaryotic and eukaryotic groups
written by dlpp •
-
ChIP-Seq alignment with random site for duplicated
written by kuragari.ch •
-
What's the best way to map CAGE-seq data to the genome?
written by dustar1986
BBMap can do this with the flag "idfilter=0.95". It can also tag the output reads with their identity using the flag "idtag" so you can manually filter them later.