ฉันมี vcf-file (file.vcf) และไฟล์ประชากรหลายไฟล์(HA1.txt, HA1.txt, HA2.txt,...,HA28.txt)
ที่มี ID แต่ละรายการ
ไฟล์ประชากรแต่ละไฟล์มีลักษณะดังนี้ (แต่ละบรรทัดเป็นหนึ่งไฟล์):
$ cat HA1.txt
QQ48_SD1A-37
ED19_SD1A40-3_357
TT335_SD1A-20
HH356_SD1A-7
Q029B_SD1A38
HT73_SD1A-28
HT288_SD1A-24
Q004B_SD1A-1
Q027_SD1A-4
Q096_SD1A-40
ฉันต้องการแยกความลึกของการครอบคลุมสำหรับแต่ละบุคคลในแต่ละประชากร ฉันจะรันคำสั่งนี้ได้อย่างไร: grep -v "^#" file.vcf | cut -f 10 | cut -d ':' -f2
สำหรับแต่ละคนในไฟล์ประชากรและบันทึกผลลัพธ์ของประชากรแต่ละคนในไฟล์แยกกัน
ผลลัพธ์ที่ฉันต้องการสำหรับประชากรแต่ละคนมีลักษณะดังนี้:
. 6 4 6 . 5 . 10 . 7 .
. 9 16 8 3 8 . 16 9 22 .
. 8 11 8 . 8 . 16 8 18 .
6 20 12 20 12 28 3 24 4 26 14
6 25 15 24 13 32 3 25 3 25 15
แต่ละคอลัมน์ในหนึ่งบุคคล!
ใช่ฉันเห็นด้วยและพยายามถามคำถามของฉันในทางที่ดีขึ้น!
—
Anna1364
โปรดแสดงเวอร์ชันที่เป็นตัวอย่างของคุณ
—
αғsнιη
file.vcf
โปรดที่นี่ฉันไม่สามารถจำแนกได้ว่าผลลัพธ์ที่คุณคาดหวังจะสร้างด้วยอินพุตตัวอย่าง HA1.txt
ซึ่งไม่เกี่ยวข้องกับเอาต์พุตและคำสั่งที่คุณใช้อยู่
Chrom1 1308 A 6622.86 T 0448 GT: DP: DPR: RO: QR: AO: QA: GL 0/0: 8: 8,0: 8: 292: 0: 0: 0, -2.40824, -26.6303 0/0: 6: 6,0 : 6: 247: 0: 0: 0, -1.80618, -21.0485 0/1: 18: 18,6: 12: 438: 6: 211: -13.9004,0, -32.5893
—
Anna1364