Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,777,009 | 0 | C | G | 28.8% | 135.3 / 21.3 | 52 | R3P (CGC→CCC) | rplT | 50S ribosomal protein L20 |
Reads supporting (aligned to +/- strand): ref base C (35/2); new base G (0/15); total (35/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.03e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTTCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAAC > REL606/1776959‑1777059 | gTTTCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTAcgc < 1:974863/51‑1 (MQ=255) gTTTCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTAcgc > 1:3623231/1‑51 (MQ=255) tttCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTAcgcg < 1:6525545/51‑1 (MQ=255) ttCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGa > 1:6457028/1‑51 (MQ=255) ttCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGa > 1:6117076/1‑51 (MQ=255) ttCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGa > 1:3199660/1‑51 (MQ=255) ttCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGa > 1:613546/1‑51 (MQ=255) cAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGc < 1:2912906/51‑1 (MQ=255) cAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGc > 1:221893/1‑51 (MQ=255) cAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGc < 1:7160275/51‑1 (MQ=255) anaTTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCa > 1:5560450/1‑51 (MQ=255) aaTTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCat > 1:6060836/1‑51 (MQ=255) aTTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCata > 1:2495806/1‑51 (MQ=255) aTTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCata > 1:188408/1‑51 (MQ=255) ttCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATAtgt > 1:4565932/1‑51 (MQ=255) ttCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATAtgt > 1:5450025/1‑51 (MQ=255) ttCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATAtgt > 1:286643/1‑51 (MQ=255) ttCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATAtgt > 1:4019706/1‑51 (MQ=255) ttCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATAtgt > 1:4217484/1‑51 (MQ=255) cTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGc > 1:1832420/1‑51 (MQ=255) ttGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGct > 1:2993534/1‑51 (MQ=255) ttGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGct > 1:3727784/1‑51 (MQ=255) tgtgACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGctc > 1:4149764/1‑51 (MQ=255) tgtgACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGctc > 1:6910002/1‑51 (MQ=255) cgtgcacgtgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTg > 1:2991130/1‑51 (MQ=255) cgtgcacgtgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTg > 1:5066792/1‑51 (MQ=255) tgcacgtgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTa > 1:5123672/1‑51 (MQ=255) gcacgtgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTAt > 1:4058453/1‑51 (MQ=255) gtgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCtat > 1:2983974/1‑51 (MQ=255) gtgaaATAACCCCACGTTTTGCGGGAGCCATATGTGCTCTCCTGTATCtat < 1:6200983/51‑1 (MQ=255) tgcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCtata > 1:6301139/1‑51 (MQ=255) gcaATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTatat > 1:3149730/1‑51 (MQ=255) aaTACCACCACTTTTTCCGGGAGCCATATGTGCTCTCCTGTATCta < 1:2835653/46‑1 (MQ=255) aaCACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTaa > 1:1548551/1‑51 (MQ=255) aaCACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTaa > 1:3354572/1‑51 (MQ=255) acacCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAAt > 1:6643259/1‑51 (MQ=255) ccccacTTTTTACGGGAGCCATATGTGCTCTCCTGTATCta < 1:4675179/39‑1 (MQ=37) ccccacGTTTTCCGGGAGCCATATGTGCTCTCCTGTATCTatat < 1:4745554/42‑1 (MQ=38) ccccacGTTTTACGGGAGCCATATGTGCTCTCCTGTATCtata < 1:4888485/41‑1 (MQ=39) accacGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTa > 1:239805/1‑51 (MQ=255) acGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTaaaa > 1:2774881/1‑51 (MQ=255) acGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTaaaa > 1:5647304/1‑51 (MQ=255) ttttGCGGGAGCCATATGTGCTCTCCTGTATCTATATTCt < 1:1675296/40‑1 (MQ=255) ttttCCGGGAGCCATATGTGCTCTCCTGTATCtat < 1:3661880/35‑1 (MQ=37) ttttCCGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTaaaa < 1:2961479/48‑1 (MQ=255) ttttCCGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTa < 1:5762536/45‑1 (MQ=255) ttttCCGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAg < 1:2387853/50‑1 (MQ=255) ttttACGGGAGCCATATGTGCTCTCCTGTATCTATATTc < 1:6574320/39‑1 (MQ=255) ttttACGGGAGCCATATGTGCTCTCCTGTATCTATATTCt < 1:1029111/40‑1 (MQ=255) ttttACGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTaaaaa < 1:7051297/49‑1 (MQ=255) ttttACGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTaaa < 1:917507/47‑1 (MQ=255) ttttACGGGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTa < 1:6516713/45‑1 (MQ=255) ttttACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGt > 1:5112462/1‑51 (MQ=255) tttACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGtt > 1:5269726/1‑51 (MQ=255) aCGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTaaa > 1:3767302/1‑51 (MQ=255) cgcgAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTaaaa > 1:6943431/1‑51 (MQ=255) gcgAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTaaaaa > 1:454256/1‑51 (MQ=255) cgAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAAc > 1:6643168/1‑51 (MQ=255) cgAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAAc > 1:2875490/1‑51 (MQ=255) cgAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAAc > 1:6197954/1‑51 (MQ=255) | GTTTCAAAATTTTCTTGTGACGTGCACGTGCAATAACACCACGTTTTACGCGAGCCATATGTGCTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAAC > REL606/1776959‑1777059 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |