Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,795,650 | 0 | A | C | 32.5% | 78.8 / 11.4 | 40 | V39G (GTG→GGG) | celF | cryptic phospho‑beta‑glucosidase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base A (26/1); new base C (0/13); total (26/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTAT > REL606/1795579‑1795719 | cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg > 1:224464/1‑76 (MQ=255) cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCa > 1:575514/1‑76 (MQ=255) ttATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATa > 1:561789/1‑76 (MQ=255) atcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAt > 1:1381636/1‑76 (MQ=255) caatctTCCGTTGGTAGAGAGCAAAAATTATATCCAGTTCCTCTTTTCCCCCTTCGACACCCACCAGCCATAATTc < 2:825111/76‑1 (MQ=255) catcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTg > 2:1353437/1‑76 (MQ=255) ccGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGAcc < 2:1328794/76‑1 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa > 2:1250511/1‑76 (MQ=255) cAGAAATCAAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAAtt < 2:46899/76‑1 (MQ=255) agagATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTc > 1:894263/1‑76 (MQ=255) agagATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTc > 2:1273561/1‑76 (MQ=255) agagATCAAAAATAATACCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTc < 1:347121/76‑1 (MQ=255) gagaTCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCg > 2:521089/1‑76 (MQ=255) tcaaaaaTAATATCCATTTTCCCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGAt < 1:1128677/73‑1 (MQ=255) tCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGAt < 1:154008/76‑1 (MQ=255) caaaaaTAATATCCATTTTCTCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa < 2:621435/74‑1 (MQ=255) cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa > 1:1150855/1‑76 (MQ=255) aaaaaaTAATACCCATTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa < 2:498921/74‑1 (MQ=255) aaaaTAATATCCATTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAAc < 1:240289/76‑1 (MQ=255) taataTCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt < 2:76020/76‑1 (MQ=255) taataCCCAGTTCCGCTTTACCCCCTCCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt < 2:666408/76‑1 (MQ=255) aataTCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta < 1:92281/76‑1 (MQ=255) aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta > 1:1401386/1‑76 (MQ=255) atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa > 2:216198/1‑76 (MQ=255) tatTCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat > 2:1074300/1‑76 (MQ=255) ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa > 2:951279/1‑76 (MQ=255) ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa > 2:951264/1‑76 (MQ=255) ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa > 2:854108/1‑76 (MQ=255) cAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAAt > 1:1292419/1‑76 (MQ=255) tCGCTTTACTACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTc > 1:802163/1‑76 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 2:1094628/1‑76 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 1:1306440/1‑76 (MQ=255) tttCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt < 1:579482/76‑1 (MQ=255) ttACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa < 1:710831/76‑1 (MQ=255) ttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa > 2:851901/1‑76 (MQ=255) tACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa > 2:666274/1‑76 (MQ=255) caccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTc > 2:844553/1‑76 (MQ=255) caccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTc > 2:399160/1‑76 (MQ=255) ccTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCg > 2:856278/1‑76 (MQ=255) gACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGta < 1:1402627/76‑1 (MQ=255) acacCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtat < 2:905314/76‑1 (MQ=255) aCATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtat > 1:1543255/1‑76 (MQ=255) | CAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTAT > REL606/1795579‑1795719 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |