Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,795,650 | 0 | A | C | 29.5% | 99.2 / 10.9 | 44 | V39G (GTG→GGG) | celF | cryptic phospho‑beta‑glucosidase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base A (27/4); new base C (0/13); total (27/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGC > REL606/1795575‑1795722 | aCGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATcca > 2:1381638/1‑76 (MQ=255) gCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc > 1:86486/1‑76 (MQ=255) gCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc > 1:666673/1‑76 (MQ=255) cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg > 1:1034973/1‑76 (MQ=255) aGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGc < 1:1012170/76‑1 (MQ=255) cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCa < 2:1310720/76‑1 (MQ=255) cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCa < 1:1523780/76‑1 (MQ=255) ttATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATa > 2:988071/1‑76 (MQ=255) ttATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATa > 2:111587/1‑76 (MQ=255) tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCCCCAGCCATaa < 2:375827/76‑1 (MQ=255) tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa > 1:404638/1‑76 (MQ=255) cGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCg > 1:1043236/1‑76 (MQ=255) gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCgg < 1:1622031/76‑1 (MQ=255) ttGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGc > 1:1431525/1‑76 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACACCCACCAGCCATAATTCGCTGACCGGCa < 1:1172694/76‑1 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa > 2:628391/1‑76 (MQ=255) cAGAGATCGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAAtt < 2:1009458/76‑1 (MQ=255) tCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGAt > 2:972906/1‑76 (MQ=255) cGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa < 2:80752/76‑1 (MQ=255) aaaaTAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAAc > 2:1985687/1‑76 (MQ=255) aaaTAATCTCCAGTTTCGTTTTACCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg < 2:188759/76‑1 (MQ=255) aaaTAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg > 1:1120315/1‑76 (MQ=255) taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt > 2:1037047/1‑76 (MQ=255) aataTCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta < 1:2016639/76‑1 (MQ=255) atatCCAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa < 1:1738732/76‑1 (MQ=255) atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa > 2:1075365/1‑76 (MQ=255) atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata > 1:1484723/1‑76 (MQ=255) atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata > 1:1293359/1‑76 (MQ=255) ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa > 1:1264135/1‑76 (MQ=255) gTTTCGCTTTACCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATca < 2:199309/76‑2 (MQ=255) tttCGCTTTCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCt < 1:987160/76‑1 (MQ=255) ttCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTTGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt < 1:793424/76‑1 (MQ=255) ttCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt < 1:988071/76‑1 (MQ=255) ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt > 2:1444268/1‑76 (MQ=255) ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt > 1:167738/1‑76 (MQ=255) tCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTc < 1:1313328/76‑1 (MQ=255) cGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc < 1:1978469/76‑1 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 2:39776/1‑76 (MQ=255) cTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAg > 1:1506200/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 1:93302/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 1:23973/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 2:735082/1‑76 (MQ=255) accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc > 1:42982/1‑76 (MQ=255) ccTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCg < 1:206808/76‑1 (MQ=255) ccTCCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCg < 1:179055/76‑1 (MQ=255) ngACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGt > 1:213738/2‑76 (MQ=255) gACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGta < 1:1829301/76‑1 (MQ=255) aCATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtat > 2:1411872/1‑76 (MQ=255) aTCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAg > 1:1242002/1‑76 (MQ=255) tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc > 1:1271177/1‑76 (MQ=255) | ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGC > REL606/1795575‑1795722 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |