|
| 1 | +##fileformat=VCFv4.2 |
| 2 | +##FILTER=<ID=PASS,Description="All filters passed"> |
| 3 | +##bcftoolsVersion=1.9+htslib-1.9 |
| 4 | +##bcftoolsCommand=mpileup -B -C 50 -d 250 --fasta-ref covid19-refseq.fasta --threads 10 -Ou SRR10903402.fq.sam.sorted.bam.rmdup |
| 5 | +##reference=file://covid19-refseq.fasta |
| 6 | +##contig=<ID=NC_045512.2,length=29903> |
| 7 | +##ALT=<ID=*,Description="Represents allele(s) other than observed."> |
| 8 | +##INFO=<ID=INDEL,Number=0,Type=Flag,Description="Indicates that the variant is an INDEL."> |
| 9 | +##INFO=<ID=IDV,Number=1,Type=Integer,Description="Maximum number of reads supporting an indel"> |
| 10 | +##INFO=<ID=IMF,Number=1,Type=Float,Description="Maximum fraction of reads supporting an indel"> |
| 11 | +##INFO=<ID=DP,Number=1,Type=Integer,Description="Raw read depth"> |
| 12 | +##INFO=<ID=VDB,Number=1,Type=Float,Description="Variant Distance Bias for filtering splice-site artefacts in RNA-seq data (bigger is better)",Version="3"> |
| 13 | +##INFO=<ID=RPB,Number=1,Type=Float,Description="Mann-Whitney U test of Read Position Bias (bigger is better)"> |
| 14 | +##INFO=<ID=MQB,Number=1,Type=Float,Description="Mann-Whitney U test of Mapping Quality Bias (bigger is better)"> |
| 15 | +##INFO=<ID=BQB,Number=1,Type=Float,Description="Mann-Whitney U test of Base Quality Bias (bigger is better)"> |
| 16 | +##INFO=<ID=MQSB,Number=1,Type=Float,Description="Mann-Whitney U test of Mapping Quality vs Strand Bias (bigger is better)"> |
| 17 | +##INFO=<ID=SGB,Number=1,Type=Float,Description="Segregation based metric."> |
| 18 | +##INFO=<ID=MQ0F,Number=1,Type=Float,Description="Fraction of MQ0 reads (smaller is better)"> |
| 19 | +##FORMAT=<ID=PL,Number=G,Type=Integer,Description="List of Phred-scaled genotype likelihoods"> |
| 20 | +##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> |
| 21 | +##INFO=<ID=ICB,Number=1,Type=Float,Description="Inbreeding Coefficient Binomial test (bigger is better)"> |
| 22 | +##INFO=<ID=HOB,Number=1,Type=Float,Description="Bias in the number of HOMs number (smaller is better)"> |
| 23 | +##INFO=<ID=AC,Number=A,Type=Integer,Description="Allele count in genotypes for each ALT allele, in the same order as listed"> |
| 24 | +##INFO=<ID=AN,Number=1,Type=Integer,Description="Total number of alleles in called genotypes"> |
| 25 | +##INFO=<ID=DP4,Number=4,Type=Integer,Description="Number of high-quality ref-forward , ref-reverse, alt-forward and alt-reverse bases"> |
| 26 | +##INFO=<ID=MQ,Number=1,Type=Integer,Description="Average mapping quality"> |
| 27 | +##bcftools_callVersion=1.9+htslib-1.9 |
| 28 | +##bcftools_callCommand=call -mv -Ov -o SRR10903402.fq.sam.sorted.bam.rmdup.vcf; Date=Thu Mar 26 16:34:45 2020 |
| 29 | +##bcftools_filterVersion=1.9+htslib-1.9 |
| 30 | +##bcftools_filterCommand=filter '-e%QUAL<10 ||(RPB<0.1 && %QUAL<15) || (AC<2 && %QUAL<15) || (DP4[0]+DP4[1])/(DP4[2]+DP4[3]) > 0.3' SRR10903402.fq.sam.sorted.bam.rmdup.vcf; Date=Thu Mar 26 16:52:08 2020 |
| 31 | +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SRR10903402.fq.sam.sorted.bam.rmdup |
| 32 | +NC_045512.2 1526 . GCCTATTGGGTTCC GCC 10.7919 PASS INDEL;IDV=1;IMF=1;DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=40 GT:PL 1/1:40,3,0 |
| 33 | +NC_045512.2 12643 . TA TACGCA 10.7919 PASS INDEL;IDV=1;IMF=1;DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=40 GT:PL 1/1:40,3,0 |
| 34 | +NC_045512.2 26062 . GGTGTTG GGTGTTGTGTTG 12.6563 PASS INDEL;IDV=1;IMF=1;DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=42 GT:PL 1/1:42,3,0 |
0 commit comments