Header Leaderboard Ad

Collapse

BBMap (aligner for DNA/RNAseq) is now open-source and available for download.

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Franz
    replied
    How is coverage calculated in bbmap?

    Hi Brian and everyone else,

    I am using bbmap for my metagenomes. Currently, I am comparing the coverage of different genes across metagenomes. BBmap has the flag -rpkm, which calculates fold coverage, RPKM and FPKM. But I couldnĀ“t find the respective formulas, how exactly it is calculated. When I tried to re-calculate these parameters, I came to different results.
    Thank you very much in advance! I hope this thread is at the right position in SEQanswers.

    Best,
    Franz

    Leave a comment:


  • turnersd
    replied
    Originally posted by HESmith View Post
    Those genotype calls are based on the read data, which are summarized in the INFO fields of the VCF. If you change the genotypes, those field data will be inconsistent.

    It would be useful if you explained more clearly what you're trying to accomplish.
    Yes, completely correct about GT not matching up with INFO fields. I'll probably end up removing all the info fields altogether because GT is really the only thing I'm interested in.

    I've used a published simulation framework to simulate pedigrees having genotype data, and what I'm trying to do is specify a rate at which I want to mutate genotypes, i.e., to mimic allele drop-out or drop-in, each at specified rates. Mimicking a het to hom-ref would be easy if variant alleles were all I cared about (ie just deleting random rows from the VCF), but I _do_ care about genotypes at invariant sites.

    My plan was to write some custom bash/python/bcftoolsy things to get out genotypes and mutate at a specified rate then stick them back into the VCF. I was just wondering if there was some existing tool to do something like this in bbmap (or elsewhere). I can't seem to find anything.

    Leave a comment:


  • HESmith
    replied
    Originally posted by turnersd View Post
    Not exactly. I literally want to replace the 0/1s or 0/0s or 1/1s in the VCF, not generate reads.
    Those genotype calls are based on the read data, which are summarized in the INFO fields of the VCF. If you change the genotypes, those field data will be inconsistent.

    It would be useful if you explained more clearly what you're trying to accomplish.

    Leave a comment:


  • turnersd
    replied
    Originally posted by HESmith View Post
    You can do it in two steps: 1) use mutate.sh on the reference genome; 2) use randomreads.sh of the mutated reference to generate the data (I assume that's what you mean by sample).
    Not exactly. I literally want to replace the 0/1s or 0/0s or 1/1s in the VCF, not generate reads.

    Leave a comment:


  • HESmith
    replied
    You can do it in two steps: 1) use mutate.sh on the reference genome; 2) use randomreads.sh of the mutated reference to generate the data (I assume that's what you mean by sample).

    Leave a comment:


  • turnersd
    replied
    Insert mutations in a VCF

    I see mutate.sh takes a vcf and inserts mutations in a reference genome. I need something a little different. I want to specify a mutation rate and potentially a sample, and insert mutations at that specified rate in the specified sample. Anything in bbmap that might help with this?

    Leave a comment:


  • reliscu
    replied
    All "N" reads in paired reads not being filtered if other read is "good"

    Version 38.94

    This thread details the "bug":
    https://www.biostars.org/p/9493307/#9493339

    Leave a comment:


  • silask
    replied
    bbmap ignores minid

    bbmap ignores minid parameter

    This is for version 38.93 and 38.92


    https://sourceforge.net/projects/bbmap/
    Last edited by silask; 09-28-2021, 11:44 PM.

    Leave a comment:


  • GenoMax
    replied
    Originally posted by lmusgrove View Post
    Hi Brian,

    I'm running BBMap with the rpkm output option and would like to know how to see mapped rate for each read file. I run multiple files consecutively with nohup. Here's my code:

    bbmap.sh ref=data/Assembly.fasta \
    in1=data/clean/A_1.clean.fq.gz \
    in2=data/clean/A_2.clean.fq.gz \
    rpkm=data/fpkm/A.fpkm \
    t=5 &

    bbmap.sh ref=data/Assembly.fasta \
    in1=data/clean/B_1.clean.fq.gz \
    in2=data/clean/B_2.clean.fq.gz \
    rpkm=data/fpkm/B.fpkm \
    t=5

    etc etc

    So far it's worked really well. However, while the stdout file shows the mapped rates it doesn't tell me which read files relate to which stats. It just has a number of repeated --- Results 1 ---- records and I don't know which is which.

    Is there a flag I need to add to ensure I can see which stats are for which files?

    Thanks!

    Lisa
    You should be able to use
    Code:
    covstats=<file>         Per-scaffold coverage info.
    for each of your commands. You can also capture the stderr/out to a file to get statistics. https://linuxize.com/post/bash-redirect-stderr-stdout/

    Leave a comment:


  • lmusgrove
    replied
    Finding mapped rate for rpkm output

    Hi Brian,

    I'm running BBMap with the rpkm output option and would like to know how to see mapped rate for each read file. I run multiple files consecutively with nohup. Here's my code:

    bbmap.sh ref=data/Assembly.fasta \
    in1=data/clean/A_1.clean.fq.gz \
    in2=data/clean/A_2.clean.fq.gz \
    rpkm=data/fpkm/A.fpkm \
    t=5 &

    bbmap.sh ref=data/Assembly.fasta \
    in1=data/clean/B_1.clean.fq.gz \
    in2=data/clean/B_2.clean.fq.gz \
    rpkm=data/fpkm/B.fpkm \
    t=5

    etc etc

    So far it's worked really well. However, while the stdout file shows the mapped rates it doesn't tell me which read files relate to which stats. It just has a number of repeated --- Results 1 ---- records and I don't know which is which.

    Is there a flag I need to add to ensure I can see which stats are for which files?

    Thanks!

    Lisa

    Leave a comment:


  • mewu3
    replied
    pileup.sh explication

    Hello,

    Can some please kindly explain the output file of pileup.sh ?
    • basecov
    • bincove
    • covstat

    How the coverage is calculated ?

    Leave a comment:


  • mewu3
    replied
    Originally posted by HESmith View Post
    "Exception in thread "main" java.lang.AssertionError: Attempting to output paired reads to different sam files."

    Typically, BBMap tools keep paired reads together. You're attempting to write aligned and unaligned reads to separate files, which violates that function.
    Thank you !

    Leave a comment:


  • GenoMax
    replied
    @mewu3: Since paired-end reads are aligned together you should use a single "out=output.sam". If you wanted to capture unmapped reads into separate files then you would want to do that as "outu1=R1.unmapped.fq outu2=R2.unmapped.fq".

    You may be able to write them out as an unmapped sam file "outu=unmapped.sam" but then again you should use only one output for that. This is untested.

    Leave a comment:


  • kesmarl
    replied
    Thank you very so much !

    I downloaded it

    Leave a comment:


  • HESmith
    replied
    "Exception in thread "main" java.lang.AssertionError: Attempting to output paired reads to different sam files."

    Typically, BBMap tools keep paired reads together. You're attempting to write aligned and unaligned reads to separate files, which violates that function.

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    Improved Targeted Sequencing: A Comprehensive Guide to Amplicon Sequencing
    by seqadmin



    Amplicon sequencing is a targeted approach that allows researchers to investigate specific regions of the genome. This technique is routinely used in applications such as variant identification, clinical research, and infectious disease surveillance. The amplicon sequencing process begins by designing primers that flank the regions of interest. The DNA sequences are then amplified through PCR (typically multiplex PCR) to produce amplicons complementary to the targets. RNA targets...
    03-21-2023, 01:49 PM
  • seqadmin
    Targeted Sequencing: Choosing Between Hybridization Capture and Amplicon Sequencing
    by seqadmin




    Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...
    03-10-2023, 05:31 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 12:26 PM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-17-2023, 12:32 PM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-15-2023, 12:42 PM
0 responses
21 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-09-2023, 10:17 AM
0 responses
68 views
1 like
Last Post seqadmin  
Working...
X