Originally posted by vivek_
View Post
as I have read on the documentation you posted, the output should be constituted with 14 columns. But in my case, I have 17 !
I think that I have :
1. Chromosome 1
2. Position 1
3. Orientation 1
4. Chromosome 2
5. Position 2
6. Orientation 2
7. Type of a SV
8. Size of a SV
9. Confidence Score
10. Total number of supporting read pairs
11. Total number of supporting read pairs from each map file
But starting from here I'm not so sure
12. Estimated allele frequency
And this following I don't have:
13. Software version
14. The run parameters
My lines look like this:
Code:
#Chr1 Pos1 Orientation1 Chr2 Pos2 Orientation2 Type Size Score num_Reads num_Reads_lib Sample_Barke_reheadered.bam Sample_Bowman_reheadered.bam Sample_Bulbosum_reheadered.bam Sample_FT11_reheadered.bam Sample_Igri_reheadered.bam Sample_Pubiflorum_reheadered.bam chr1H_part1 13958 241+168- chr1H_part1 14031 241+168- ITX -212 34 3 ../barleyBAMs/Sample_Barke_reheadered.bam|3 NA NA NA NA NA NA chr1H_part1 22659 14+28- chr1H_part1 22826 14+28- ITX -147 89 6 ../barleyBAMs/Sample_Bulbosum_reheadered.bam|6 NA NA NA NA NA NA chr1H_part1 51554 3+2- chr1H_part1 51654 3+2- ITX -157 38 2 ../barleyBAMs/Sample_Pubiflorum_reheadered.bam|2 NA NA NA NA NA NA chr1H_part1 87327 42+24- chr1H_part1 87429 42+24- ITX -179 99 18 ../barleyBAMs/Sample_Barke_reheadered.bam|5:../barleyBAMs/Sample_Bulbosum_reheadered.bam|12:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|1 NA NA NA NA NA NA chr1H_part1 89261 11+5- chr1H_part1 89314 11+5- ITX -176 74 4 ../barleyBAMs/Sample_Bulbosum_reheadered.bam|4 NA NA NA NA NA NA chr1H_part1 94809 25+24- chr1H_part1 94860 25+24- ITX -154 99 15 ../barleyBAMs/Sample_Bulbosum_reheadered.bam|14:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|1 NA NA NA NA NA NA chr1H_part1 95422 4+4- chr1H_part1 95467 4+4- ITX -169 55 3 ../barleyBAMs/Sample_Bulbosum_reheadered.bam|3 NA NA NA NA NA NA chr1H_part1 97915 153+134- chr1H_part1 98001 153+134- ITX -170 99 32 ../barleyBAMs/Sample_Barke_reheadered.bam|10:../barleyBAMs/Sample_Bulbosum_reheadered.bam|15:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|7 NA NA NA NA NA NA chr1H_part1 99599 14+8- chr1H_part1 99668 14+8- ITX -154 41 2 ../barleyBAMs/Sample_Pubiflorum_reheadered.bam|2 NA NA NA NA NA NA chr1H_part1 100374 7+3- chr1H_part1 100449 7+3- ITX -164 39 3 ../barleyBAMs/Sample_Barke_reheadered.bam|1:../barleyBAMs/Sample_Bulbosum_reheadered.bam|2 NA NA NA NA NA NA chr1H_part1 124801 45+81- chr1H_part1 133975 25+2- ITX 8473 99 17 ../barleyBAMs/Sample_Barke_reheadered.bam|11:../barleyBAMs/Sample_Igri_reheadered.bam|6 0.21 0.00 0.03 0.09 0.25 0.02 chr1H_part1 124801 45+64- chr1H_part1 134111 6+2- ITX 5569 45 6 ../barleyBAMs/Sample_Barke_reheadered.bam|3:../barleyBAMs/Sample_Bulbosum_reheadered.bam|2:../barleyBAMs/Sample_Igri_reheadered.bam|1 0.21 0.00 0.03 0.10 0.25 0.02 chr1H_part1 125588 18+2- chr1H_part1 126418 1+15- DEL 890 45 8 ../barleyBAMs/Sample_Barke_reheadered.bam|5:../barleyBAMs/Sample_Igri_reheadered.bam|3 0.21 0.02 0.21 0.06 0.38 0.05 chr1H_part1 127470 8+33- chr1H_part1 133703 9+33- INV 5823 99 22 ../barleyBAMs/Sample_Barke_reheadered.bam|16:../barleyBAMs/Sample_Igri_reheadered.bam|6 0.16 0.00 0.01 0.06 0.18 0.00
Thank you very much, have a nice day
C.
Leave a comment: