Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • zorph
    Member
    • May 2010
    • 40

    #1

    Using GATK to call SNPs in RNA-seq data

    has anyone ever used SNP calling from RNA-seq data?
    I am currently trying to do this but GATK is calling the reference as the variant??? I have checked my reference file and it's perfectly correct.

    I don't know if this a result from me using RNA-seq or just my command.

    I am using an older version of gatk because it doesn't require the extensive file modifications that is required of the newer versions.

    My command is below:

    GenomeAnalysisTK.jar
    -T UnifiedGenotyper
    -l INFO
    -I file.bam
    -R Homo_sapiens_assembly18.fasta
    -U ALLOW_UNINDEXED_BAM
    --assume_single_sample_reads sample_name
    --platform SOLEXA
    -dcov 300
    -varout file.vcf
    -nt 1
    -bad_mates

    any advice?

Latest Articles

Collapse

  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2




    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
    Yesterday, 06:46 AM
  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
26 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-28-2026, 09:51 AM
0 responses
17 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
71 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Working...