Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | 124 | 13,836 | 0 | T | C | 65.2% | 47.8 / 39.3 | 46 | intergenic (–/–) | –/– | –/– |
Reads supporting (aligned to +/- strand): ref base T (8/8); new base C (15/15); total (23/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCACAATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > 124/13736‑13852 | gAGCACAATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAAt > 1:406599/1‑101 (MQ=255) gAGCACAATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAAt > 2:1632493/1‑101 (MQ=255) aGCACAATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAAtt < 1:1430964/101‑1 (MQ=255) acaATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCa > 2:1528973/1‑101 (MQ=255) caATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCAc < 1:1607497/101‑1 (MQ=255) aaTGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCAcc > 1:496929/1‑101 (MQ=255) tGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCa > 2:201891/1‑101 (MQ=255) tGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGCAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcgga > 1:163802‑M1/1‑97 (MQ=255) gCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggat < 2:1340785‑M1/101‑6 (MQ=255) gCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggat < 1:593008‑M1/101‑6 (MQ=255) cTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatg < 2:980198‑M1/101‑7 (MQ=255) tGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgc < 1:350100‑M1/101‑8 (MQ=255) tGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTAAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgc < 1:1625450‑M1/101‑8 (MQ=255) ggCGCATCTTATCAGGCCTGCGGGGTGTGCAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggaaggg > 2:606623‑M1/1‑93 (MQ=255) gaccatcTTGTTAGGCCTGCGGGTTGTGTTGAGGGTGTAGGCCGGATAAGTTGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcgg < 2:818846‑M1/98‑10 (MQ=255) cgcATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcggc > 1:783685‑M1/1‑91 (MQ=255) cATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcggcga > 1:917244‑M1/1‑89 (MQ=255) aTCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcggcgat > 1:1428458‑M1/1‑88 (MQ=255) cTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcggcgatta < 1:1303854‑M1/101‑16 (MQ=255) aGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGGcggatgcggcgattaatccag < 1:474786‑M1/101‑22 (MQ=255) aatcctcctggccggataaaACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcg < 1:1485807‑M1/92‑40 (MQ=255) aatcctcctggccggataaaACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcg > 2:1163358‑M1/10‑62 (MQ=255) aatcctcctggccggataaaACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcg < 2:107967‑M1/92‑40 (MQ=255) aatcctcctggccggataaaACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcg > 2:660680‑M1/10‑62 (MQ=255) tcctggccggATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttccc > 1:1386664‑M1/5‑57 (MQ=255) ctggccGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgt > 1:136534‑M1/3‑55 (MQ=255) ggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtag < 2:799310‑M1/101‑49 (MQ=255) ggCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtag > 1:1207585‑M1/1‑53 (MQ=255) gCCGGATAAAACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtagg < 1:682173‑M1/101‑50 (MQ=255) ggATAACACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcg < 2:825653‑M1/101‑53 (MQ=255) taaaacgcggcaagCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggat > 1:388203‑M1/6‑46 (MQ=255) aaaacgcggcaaGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggata > 2:1490661‑M1/5‑45 (MQ=255) aaacgcggcaAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataa < 2:165866‑M1/98‑58 (MQ=255) gacgcgacAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatccccccggccggataaa > 2:67942‑M1/3‑43 (MQ=255) aacgcggcAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataag < 1:1703037‑M1/99‑59 (MQ=255) acgcgaCAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaa < 1:69504‑M1/100‑60 (MQ=255) cgcgACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaac < 1:1317084‑M1/101‑61 (MQ=255) gacaaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacg < 1:989530‑M1/99‑64 (MQ=255) gacaaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacg < 2:814304‑M1/99‑64 (MQ=255) gcaAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataagaagcga > 2:1584319‑M1/2‑37 (MQ=255) acaAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacgg > 1:1182973‑M1/2‑37 (MQ=255) aaGCGTTGCATCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataagacgcgaca < 2:1144464‑M1/101‑67 (MQ=255) aGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaatcgttata < 2:338302‑M1/101‑68 (MQ=255) aGCGTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacggcaa > 2:625096‑M1/1‑34 (MQ=255) gTCGCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacggcaagcg > 1:271624‑M1/1‑31 (MQ=255) gCATCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacggcaaacgtcg > 2:799966‑M1/1‑28 (MQ=255) aTCCGGCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataagacgcgacaagcgtcgca > 1:1614749‑M1/1‑26 (MQ=255) tCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacggcaagcgtctcat < 1:1410598‑M1/101‑77 (MQ=255) tCCGGCATCTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcaatcccccaggccggataaaacacggcaagcgtcgcat < 1:659946‑M1/101‑77 (MQ=255) ggCAACTGCACCGCGCCACTGGcggatgcggcgtgaacgccttatccgccctacatgtgcgttcccgtaggtcggataagacgcgacaagcgtcgcctccg > 2:481358‑M1/1‑22 (MQ=255) | GAGCACAATGCCTGATGCGATGCTGGCGCATCTTATCAGGCCTGCGGGGTGTGTAGCGGGTGTAGGCCGGATAAGGCGCGTCAAGCGTCGCATCCGGCAATTGCACCGCGCCACTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > 124/13736‑13852 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |