АвторТема: WGSExtract, FASTQ и разные референсные геномы  (Прочитано 3776 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн Georg

  • Сообщений: 949
  • Страна: ru
  • Рейтинг +308/-6
  • Племянник vk511(Y),vk160(аДНК)
  • Y-ДНК: I1a1b1a1e2 (Y353312)
  • мтДНК: Я U5b2a1b, дети T1a1ct и H5b*
Re: WGSExtract, FASTQ и разные референсные геномы
« Ответ #15 : 14 Февраль 2025, 17:37:45 »
Забрал в Биотехе свой геном в формате fastq.
Хотел в t2t сам перегнать.
В wgse стабильной июньской сборки никак чтото не вижу, куда засунуть эти файлики (4 папки с файлами).

С samtools в командной строке чтото не хочется сталкиваться 
Возможно переоценил свои силы, надо было bam брать)

Оффлайн Daemon2017Автор темы

  • Сообщений: 2553
  • Страна: ru
  • Рейтинг +1315/-19
  • Y-ДНК: R1a-Y35174
  • мтДНК: V7-a2a2a2b*
Re: WGSExtract, FASTQ и разные референсные геномы
« Ответ #16 : 14 Февраль 2025, 17:44:06 »
Забрал в Биотехе свой геном в формате fastq.
Хотел в t2t сам перегнать.
В wgse стабильной июньской сборки никак чтото не вижу, куда засунуть эти файлики (4 папки с файлами).

С samtools в командной строке чтото не хочется сталкиваться 
Возможно переоценил свои силы, надо было bam брать)

Вот тут манул https://madelinemiller.dev/blog/whole-genome-sequencing-wgsextract/
WGSe должен быть, кажется, не древнее 2022 года, чтобы работать с FASTQ, а дальше все легко: импортировали, проанализировали, выровняли, записали в BAM с нужным опорником (hg38, T2T)

Оффлайн Georg

  • Сообщений: 949
  • Страна: ru
  • Рейтинг +308/-6
  • Племянник vk511(Y),vk160(аДНК)
  • Y-ДНК: I1a1b1a1e2 (Y353312)
  • мтДНК: Я U5b2a1b, дети T1a1ct и H5b*
Re: WGSExtract, FASTQ и разные референсные геномы
« Ответ #17 : 14 Февраль 2025, 18:59:43 »
Забрал в Биотехе свой геном в формате fastq.
Хотел в t2t сам перегнать.
В wgse стабильной июньской сборки никак чтото не вижу, куда засунуть эти файлики (4 папки с файлами).

С samtools в командной строке чтото не хочется сталкиваться 
Возможно переоценил свои силы, надо было bam брать)

Вот тут манул https://madelinemiller.dev/blog/whole-genome-sequencing-wgsextract/
WGSe должен быть, кажется, не древнее 2022 года, чтобы работать с FASTQ, а дальше все легко: импортировали, проанализировали, выровняли, записали в BAM с нужным опорником (hg38, T2T)
а какой t2t выбирать?
chm13 v2.0  25 SN (шимпанзе, что ли?)
hg002xy v2.7 25 SN
hg01243 v3 89 SN (начал пока с этим)
hg002xy v2 25 SN

---------------------------
Reference Genome File
---------------------------
The THG1243v3 Reference Genome (hg01243_v3.fna.gz) is required and still missing. The current command is cancelled.
не скачивает референсы и не устанавливает. была когда-то такая ошибка

« Последнее редактирование: 14 Февраль 2025, 19:23:21 от Georg »

Оффлайн Semargl

  • "El sueño de la razón produce monstruos" ©
  • Модератор
  • *****
  • Сообщений: 6318
  • Страна: ru
  • Рейтинг +4665/-5
  • www.semargl.me
    • www.semargl.me
  • Y-ДНК: R1a [CTS3402+]
  • мтДНК: U4a2g
Re: WGSExtract, FASTQ и разные референсные геномы
« Ответ #18 : 14 Февраль 2025, 21:01:56 »
а какой t2t выбирать?
chm13 v2.0  25 SN (шимпанзе, что ли?)
hg002xy v2.7 25 SN
hg01243 v3 89 SN (начал пока с этим)
hg002xy v2 25 SN

chm13 v2.0  25 SN

Оффлайн Capo

  • Сообщений: 26
  • Страна: ru
  • Рейтинг +7/-0
Re: WGSExtract, FASTQ и разные референсные геномы
« Ответ #19 : 14 Февраль 2025, 21:09:07 »
Добрый вечер,
Я тоже забирал в fastq, yfull за 2500₽ сделал мне bam  :)

 

© 2007 Молекулярная Генеалогия (МолГен)

Внимание! Все сообщения отражают только мнения их авторов.
Все права на материалы принадлежат их авторам (владельцам) и сетевым изданиям, с которых они взяты.