Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,515,909 | 0 | A | C | 48.9% | 25.5 / 34.3 | 45 | V260G (GTG→GGG) | ypfI | predicted hydrolase |
Reads supporting (aligned to +/- strand): ref base A (14/9); new base C (0/22); total (14/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.81e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAA > REL606/2515837‑2515981 | aacaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGAccacca < 2:3696786/76‑1 (MQ=255) acaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCAccag < 1:2610133/76‑1 (MQ=255) aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc < 1:2433761/76‑1 (MQ=255) aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GCCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc < 2:1562581/76‑1 (MQ=255) acGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca < 1:2166013/76‑1 (MQ=255) acGTTCGAGGAAAACGCGATACCAGTTGAGGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca < 2:3992043/76‑1 (MQ=255) cGTTCGAGTAAAACGCGATCCCAGTTGATGCAACAAT‑GGGGCAGGTATGGCTGCGGCTCCATCGACCAccagccag < 2:2704345/76‑1 (MQ=255) cGTTCGAGGAAACCGCGATCCCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCAccagccag < 1:2073816/76‑1 (MQ=255) agaaaacCGCAACACCAGTTGATGCAAAAAT‑GGCGCAGGTATGGCGGCGCCTTCATCGACCCCCAGCCAGTCGG‑Ct < 2:2600680/73‑1 (MQ=255) ggAAACCCCGATCCCAGTTGATGCAACAAG‑GGCGCAGGTTTGGCGGCGGCTTCATCGCCCCCCAGCCAGTCGG‑Ctt < 2:1358715/76‑1 (MQ=255) ggAAAACGCTATCCCAGTTGATGCAACAAG‑GGCGCAGGTAGGCCTGCGCCTTCACCGCCCCCCAGCCAGTCGG‑Ctt < 2:2002979/76‑1 (MQ=255) aaaCGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTTCTCTACCACCCGCCAGTGGG‑CGTgtt > 1:2751897/1‑74 (MQ=255) cgcgACACCATTTGATGCAACAAG‑GGCGCGGGTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCAt < 1:4440026/76‑1 (MQ=255) atcccATTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCt < 2:2858372/73‑1 (MQ=255) aTACCAGTGGATGCAAAAAG‑GGCGCAGGTATGGCTGCGCCTTCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCt < 2:3050233/76‑1 (MQ=255) tACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTg < 1:3022381/76‑1 (MQ=255) tACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCCTTATTGACCACCAACCAGGCGGCCTTGCACTTCGCt > 1:4943090/1‑76 (MQ=255) tACCAGTTGATGCAACAAG‑GGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTg < 1:3704639/76‑1 (MQ=255) acccAGTTGATCCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTg < 1:2347648/74‑1 (MQ=255) acccAGTTAATGCAACAAT‑GGCGCGGGTAGGGCTGCGTCTTCATCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTg < 1:4905397/74‑1 (MQ=255) ccAGTTGA‑GCCACAATGGGCGCAGGTTTGCCGGCGTCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGc < 1:4177546/76‑1 (MQ=255) cAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCt < 2:4675110/76‑1 (MQ=255) tgatgaAACAAG‑GGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑a < 2:718631/76‑1 (MQ=255) gaggCAACAAG‑GGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑aa < 1:1508493/73‑1 (MQ=255) gATGCAAAAAG‑GGCGCGGGTAGGGCTGCGGCTTCATCGCCCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑aa < 1:3246139/76‑1 (MQ=255) gATCCAACAGC‑GGCGCAGGTGTGGTGGCGCCTTCACCGACCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑aa < 2:1819830/76‑1 (MQ=255) ggCAACAGG‑GGCGCGGTTGTGGCTGCGGCTTCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAgg < 1:1361733/75‑1 (MQ=255) gggCGCAGGTAGGGCGGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGgc < 2:2148673/75‑1 (MQ=255) cgcgGGTATGGCGGCGCCTTCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCa < 1:3451779/76‑1 (MQ=255) gggTGTGGTTGCGCCTTCCCCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATa < 1:914762/75‑1 (MQ=255) gggTAGGGCGGCGGCTTCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATa < 1:2024234/75‑1 (MQ=255) ttATGGCGGCGCTTTCAACCCCCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATaaa < 1:270247/75‑1 (MQ=255) gTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGG‑CTTTCTCATCGCTGGCTAAC‑AAAGGATCCCGCGCAATaaa > 1:4436648/1‑76 (MQ=255) gTATCGCTGCGGCTTCGTCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATaaa < 1:305283/76‑1 (MQ=255) tttGGCTGCGCCTTCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAg < 2:1637176/74‑1 (MQ=255) ttgggCTGCGCCTCCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAg < 1:95727/73‑1 (MQ=255) tATGGCTGCGGCTTCATCGACCACCAGCCAGTCGG‑CTTTCTCCTCGCTGGGTTAA‑AAAGCATCCCGCGCCATaaaa > 1:2017659/1‑75 (MQ=255) tATGGCGGGGGCTTCCTCGCCCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAg < 2:1064041/76‑1 (MQ=255) tGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCg < 2:4899165/76‑1 (MQ=255) ggCTGCGGCTTCATCCACCACCAGCCAGTCGG‑CTTGCTCATCGCCGGCGAACAAAGGCACCCGG‑GCAATAAAGCGa > 2:1329825/1‑76 (MQ=255) gcggCGGCTTCACCCCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGaa < 1:4614099/73‑1 (MQ=255) gCTGCGCCTTCAGCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGaa < 1:3774394/76‑1 (MQ=255) tGCGGCTTCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGGATCCGGCGCCATAAAGCGGAAc > 1:1839335/1‑76 (MQ=255) tGCGGCTTCATCCACCACCAGCCAGTCGG‑CTTGATCATCGCTGGCTAAC‑AAGGGATCCGGCCGCATAAAACGaaaa > 1:3514069/1‑75 (MQ=255) tGCGGCTTCACCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAAc < 1:2856372/76‑1 (MQ=255) gCTTCGCCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTtct < 2:3304998/76‑1 (MQ=255) gCTTCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGTTAAC‑AAGGCATCCGGCGCAATAAAGCGGAACTtct > 1:2925931/1‑76 (MQ=255) ttCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAAATTCTCg > 1:2875563/1‑76 (MQ=255) ttCATCGACCACCAGCCAGTCAG‑CTTGGTCATCGATGGCTAAC‑GACGGATCCCGAGCAAAAAAGTGAAACTtctct > 2:4902746/1‑75 (MQ=255) cATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGcc > 1:1018821/1‑76 (MQ=255) aTCGACCACCAGCCAGTCGG‑CTTGTCCATGGCTGGCTAAC‑AAGGCATCCGGGGCAATAGAGCGAAACTTTTCGccc > 2:794515/1‑76 (MQ=255) tCGCCCCCCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCg < 2:1961667/76‑1 (MQ=255) cGACCACCCGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑ACAGCATCCGGCTCAAATAAGAGAAACTTCTCGccctc > 2:4975555/1‑74 (MQ=255) cGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc > 1:2706907/1‑76 (MQ=255) cGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc > 1:20901/1‑76 (MQ=255) cGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc > 2:524353/1‑76 (MQ=255) accaccAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCaa > 1:3211739/1‑76 (MQ=255) | AACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAAT‑GGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGG‑CTTGCTCATCGCTGGCTAAC‑AAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAA > REL606/2515837‑2515981 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |