Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • HLA identification

    Hello,

    I would like to identify the HLA alleles of some samples sequenced with Solid (50bp single end). My first idea was to extract a sam file containing the reads mapped in HLA locus (using samtools) then convert it in fastq (using picard) and finally use seq2HLA. The problem is that seq2HLA need paired end read. I saw that HLAminer could also do the job. The problem is that I can manage to type a correct command line : I tried : perl HLAminer.pl -a myfile.sam ; but it get this :
    Usage: HLAminer.pl [v1.0.5]
    Derivation of HLA class I and II predictions from Shotgun sequence datasets
    --------------------------------------------------------------------
    HPTASR (HLA Predictions by Targeted Assembly of Shotgun Reads):
    -b blastn alignments.........................<tig_vs_hla-ncbi.coord>
    -r reciprocal blastn.........................<hla_vs_tig-ncbi.coord>
    -c contig fasta file.........................<TASRhla200.contigs>
    -z minimum contig size.......................<200>
    ------------------------------- OR ---------------------------------
    HPRA (HLA Predictions by Read Alignment):
    -a sam alignments............................<ngs_vs_hla.sam>
    --------------------------------------------------------------------
    -h hla fasta file............................<HLA_ABC_CDS.fasta>
    -p P-designation file........................<hla_nom_p.txt>
    -i minimum % sequence identity...............<99>
    -q minimum log10 (phred-like) expect value...<30>
    -s minimum score.............................<1000>
    -n consider null alleles (1=yes/0=no)........<0>
    -l label (run name) -optional-

    Could you help me ?

    Thanks in advance

  • #2
    you have to specify the hla reference sequence file (-h) at a minimum..

    -h hla fasta file............................<HLA_ABC_CDS.fasta>

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Exploring the Dynamics of the Tumor Microenvironment
      by seqadmin




      The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
      07-08-2024, 03:19 PM
    • seqadmin
      Exploring Human Diversity Through Large-Scale Omics
      by seqadmin


      In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
      06-25-2024, 06:43 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 07:20 AM
    0 responses
    23 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-16-2024, 05:49 AM
    0 responses
    38 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-15-2024, 06:53 AM
    0 responses
    43 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-10-2024, 07:30 AM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Working...
    X