Home Project List Latest Recipes
Login

Bioinformatics Data Analysis

Project 0  Data 21  Recipes 26 Results
View Result
Results for How to generate variant calls from sequencing data: AF086833, SRR1553500, 100000

Results generated by running the recipe.

Completed
• Runtime 2 minutes • updated 5.1 years ago by Istvan Albert
Recipe Code Rerun Recipe Edit Result Delete Result

Results generated by running the recipe.

Parameters used during the run:

  • Genome accession number:AF086833
  • The SRR run number:SRR1553500
  • The number of reads:100000
Output Messages
Messages printed to the standard output stream:


        

            
        

    


    

        
Other Messages

        
Messages printed to the standard error stream:

        
+ ACC=AF086833
+ SRR=SRR1553500
+ mkdir -p refs
+ REF=refs/AF086833.fa
+ BAM=SRR1553500.bam
+ efetch -db=nuccore -format=fasta -id=AF086833
+ seqret -filter -sid AF086833
+ bwa index refs/AF086833.fa
+ samtools faidx refs/AF086833.fa
+ fastq-dump -X 100000 --split-files SRR1553500
+ R1=SRR1553500_1.fastq
+ R2=SRR1553500_2.fastq
+ TAG='@RG\tID:SRR1553500\tSM:SRR1553500\tLB:SRR1553500'
+ bwa mem -R '@RG\tID:SRR1553500\tSM:SRR1553500\tLB:SRR1553500' refs/AF086833.fa SRR1553500_1.fastq SRR1553500_2.fastq
+ samtools sort
+ samtools index SRR1553500.bam
+ bcftools mpileup -Ovu -f refs/AF086833.fa SRR1553500.bam
[mpileup] 1 samples in 1 input files
+ bcftools call -vm -Ov genotypes.vcf
+ freebayes --ploidy 1 -f refs/AF086833.fa SRR1553500.bam


    





        

    






    


    

        
            

                Powered by the 
                  release 2.3.6