Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Identical insert size and SD #948

Open
jamespblloyd-uwa opened this issue Jul 30, 2024 · 0 comments
Open

Identical insert size and SD #948

jamespblloyd-uwa opened this issue Jul 30, 2024 · 0 comments

Comments

@jamespblloyd-uwa
Copy link

Is the bug primarily related to salmon (bulk mode) or alevin (single-cell mode)?
salmon (bulk mode)
Describe the bug
I did a lot of independent RNAseq libs recently and all came out with the identical frag_length mean and sd (see below). This surprised me greatly as while the libraries looked consistent on a tapestation, I would not expect them to be this consistent down to >7 decimal places. Also it is unclear how I can access fld.gz to see the raw data. MultiQC only shows me one distribution (for

frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_0h_rep2A_RL5383_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_0h_rep3A_RL5403_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_37C_rep1_RL5385_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_37C_rep2_RL5386_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_37C_rep3_RL5406_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_4C_rep1_RL5387_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_4C_rep2_RL5388_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_4C_rep3_RL5407_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_RT_rep1_RL5384_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_RT_rep2_RL5404_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_3h_RT_rep3_RL5405_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_37C_rep1_RL5390_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_37C_rep2_RL5391_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_37C_rep3_RL5410_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_4C_rep1_RL5392_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_4C_rep2_RL5393_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_4C_rep3_RL5411_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_RT_rep1_RL5389_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_RT_rep2_RL5408_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427

Sample: Gam_6h_RT_rep3_RL5409_quant
frag_length_mean: 250.0000097391447
frag_length_sd: 25.001769513739427```

**To Reproduce**
Steps and data to reproduce the behavior:

Specifically, please provide at least the following information:

  * Which version of salmon was used?
  v1.10.0
  * How was salmon installed (compiled, downloaded executable, through bioconda)?
  Binary
  * Which reference (e.g. transcriptome) was used?
  Moss (Physco) v6.1 that I had generated 
  * Which read files were used?
  
  * Which which program options were used?
`-l A -p 16 --validateMappings --numBootstraps 100 --seqBias --gcBias`

**Expected behavior**
A clear and concise description of what you expected to happen.

**Desktop (please complete the following information):**
 - OS: Linux server: Rocky Linux release 8.5 (Green Obsidian)
 - Version Rocky Linux release 8.5 (Green Obsidian)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant