Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,875,483 | 0 | G | C | 41.1% | 69.1 / 56.6 | 56 | intergenic (+156/‑28) | yeaB/sdaA | predicted NUDIX hydrolase/L‑serine deaminase I |
Reads supporting (aligned to +/- strand): ref base G (13/20); new base C (23/0); total (36/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCAT > REL606/1875410‑1875556 | aTTCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTTGTTCGTTACTGGAACTc > 1:673545/1‑76 (MQ=255) aTTCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTc > 1:214856/1‑76 (MQ=255) aTTCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTc > 1:109246/1‑76 (MQ=255) ttCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTcc < 1:739046/76‑1 (MQ=255) tCATGTGAATCGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCa > 1:1073407/1‑76 (MQ=255) tCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCa > 1:898166/1‑76 (MQ=255) cATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAg > 2:127923/1‑76 (MQ=255) cATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCCCTTACAATATGCGCTGTTATTAGTTTGTTACTGGAACTCCAg > 1:212107/1‑76 (MQ=255) aTGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGt > 1:464475/1‑76 (MQ=255) aTGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGt > 1:88740/1‑76 (MQ=255) tgtgAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACANTATGCGCTGTTATTAGTTCGTNTCTGGAACTCCAGTc > 1:285211/1‑76 (MQ=255) tgtgAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTc < 2:551878/76‑1 (MQ=255) gtgAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCa > 1:858826/1‑76 (MQ=255) gtgAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTAATGGAACTCCAGTCa > 1:861440/1‑76 (MQ=255) tgAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCAc > 1:447199/1‑76 (MQ=255) tgAATAGTTAAGCCAGTCGCCGCGTACCCTCGTACACTAANNGCTGTTATTANNTAGTTACTGGAACTACAGTCAc > 1:1047098/1‑76 (MQ=255) gAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTNCAGTCAcc > 1:1024291/1‑76 (MQ=255) gAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCAcc > 1:214942/1‑76 (MQ=255) gAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCAcc > 1:256432/1‑76 (MQ=255) gAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCAcc > 1:52576/1‑76 (MQ=255) gAATAGTTAAGCCAGTCGCCGCGTTCCCCCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTNCAGTCAcc > 1:1024783/1‑76 (MQ=255) aaTAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCt > 1:669748/1‑76 (MQ=255) aTAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCtt < 2:911625/76‑1 (MQ=255) aTAGTTAAGCCAGTCGCCGCGTTCCATCTTATACTATGCCNNGTTATTAGTTNNTGACTGGAACTCCAGTCACCtt > 1:1054052/1‑76 (MQ=255) tAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTg > 1:596515/1‑76 (MQ=255) tAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTg > 1:236588/1‑76 (MQ=255) aGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGt < 1:165300/76‑1 (MQ=255) aGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGt > 1:975457/1‑76 (MQ=255) aGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGt > 1:922248/1‑76 (MQ=255) aGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGCTACTGGAACGCCAGTCACCTTGt > 1:1004148/1‑76 (MQ=255) aGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTAGGATCTGTTATTANTTCGTTACTGGAACTCCAGTCACCTTGt > 1:96569/1‑76 (MQ=255) gTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTc < 1:55624/76‑1 (MQ=255) ttAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaa > 1:948812/1‑74 (MQ=255) ttAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaa > 1:366923/1‑74 (MQ=255) ttAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaa > 1:857399/1‑74 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAg > 2:53939/1‑76 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaat > 1:469634/1‑73 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaat > 1:299891/1‑73 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaat > 1:352487/1‑73 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaat > 1:525274/1‑73 (MQ=255) tAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAACTCCAGTCACCTTGtaat > 1:709209/1‑73 (MQ=255) aaGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAgg > 2:851765/1‑76 (MQ=255) aaGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCAAGTCACCTTGTCAgg > 2:33905/1‑76 (MQ=255) aGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGa < 1:584203/76‑1 (MQ=255) ccAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGt < 2:363184/76‑1 (MQ=255) cAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGta < 1:863672/76‑1 (MQ=255) tCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGcatta < 2:1150607/76‑1 (MQ=255) ttCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGAtt > 1:675012/1‑76 (MQ=255) ttCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGAtt < 1:448660/76‑1 (MQ=255) ttCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGAtt < 2:391663/76‑1 (MQ=255) ttCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGAtt < 2:297415/76‑1 (MQ=255) ctctTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGta < 2:142462/76‑2 (MQ=255) acTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGa > 1:271860/1‑76 (MQ=255) tATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACa < 2:1024290/76‑1 (MQ=255) aTGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACAt < 1:303526/76‑1 (MQ=255) tGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATg > 1:438281/1‑76 (MQ=255) cTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTa < 2:1017803/76‑1 (MQ=255) tGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTaa < 2:1035031/76‑1 (MQ=255) tGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTaa < 2:401772/76‑1 (MQ=255) tattaGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGt > 2:96171/1‑76 (MQ=255) gTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGa < 2:662995/76‑1 (MQ=255) tACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTc > 1:608449/1‑76 (MQ=255) gAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCa > 1:514995/1‑76 (MQ=255) gaGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCAt < 1:610354/75‑1 (MQ=255) | ATTCATGTGAATAGTTAAGCCAGTCGCCGCGTTCCCTCTTACACTATGCGCTGTTATTAGTTCGTTACTGGAAGTCCAGTCACCTTGTCAGGAGTATTATCGTGATTAGTCTATTCGACATGTTTAAGGTGGGGATTGGTCCCTCAT > REL606/1875410‑1875556 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |