Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,735,608 | 0 | G | C | 52.9% | 6.3 / 21.6 | 17 | A409P (GCC→CCC) | ascF | fused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component |
Reads supporting (aligned to +/- strand): ref base G (7/1); new base C (0/9); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGCATCG > REL606/2735533‑2735683 | cGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg < 2:5014279/76‑1 (MQ=255) gCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑c < 2:5185554/76‑1 (MQ=255) ccTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑cc < 2:5553564/76‑1 (MQ=255) ctaaaaCCTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑cca < 2:4567455/73‑1 (MQ=255) cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑cca < 1:1087741/76‑1 (MQ=255) cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑cca < 1:1953740/76‑1 (MQ=255) tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑ccag < 2:3125856/76‑1 (MQ=255) aaaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑ccagcc < 2:2551098/76‑1 (MQ=255) aaaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑ccagcc < 1:2081515/76‑1 (MQ=255) aaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg‑ccagcca < 1:4124092/76‑1 (MQ=255) aCGTCCTCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTTTGGCGGGGCTTC‑CCAGCCAc < 2:6175801/76‑1 (MQ=255) aCGTCCTCTTATCCCCAGTCTTATCAGCGGTTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTG‑CCAGCCAc < 2:3346272/76‑1 (MQ=255) cccGCTTGTCGCCACGCTTCTCGGCGTCTTTTTTGGCGGGGCGGTTGCCGGTATGGGGGGGTTTC‑CCAGCCACTCa < 2:445707/75‑1 (MQ=255) cGCTTGTCGCCATTCTTACCAGCGGCTTTTTTTGGGCGGCGATGGACGGTGTGGCGGGGCTTC‑CCAGCCACTCAAt < 2:3099384/76‑1 (MQ=255) gtttATCGCCGTTTTTATCAGCGGCTTTTTTTGCGGCGCGGTTGCGGGTATGGGGGGGCTTC‑CCAGCCACTCAATg < 2:2278808/74‑1 (MQ=255) cattCTTATCAGCGGCTTTATTTGCGGCGCGGTGGCCGGTAGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCg < 1:2105578/73‑1 (MQ=255) tgTCTTAGCAGCGGCTTTATGAGCGGGGCGGTCGCGGGTCGGGCGGGATTTG‑CCAGCCACTCAATGGCAGCGCCgg < 2:934261/75‑1 (MQ=255) tgTCTTACCGGCGGTTTTATTGGCGGCGCGTTTGCGGGTAGGGGGGGGCTCC‑CCAGCCACTCAATGGCAGCGCCgg < 2:2213712/75‑1 (MQ=255) aGTTTTATCAGGGGTTTTTTTTGGGGCGCGGTGCCCGGTAGGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCgg < 2:2352718/76‑1 (MQ=255) aGTCTTATCGGCGGTTTTATTTGCGGCCGGTTTGCCGGTAGGGCGGGCCTTC‑CCAGCCACTCAATGGCAGCGCCgg < 2:1316216/76‑1 (MQ=255) aGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCgg < 2:3100706/76‑1 (MQ=255) gttttATCGGCGGTTTTTTTGGCGGCGCGGTTGCCGGTATGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCggg < 1:1666980/73‑1 (MQ=255) gTCTTACCGGCGCTTTTATTTGCGGGGGGTTTGCCGGTAGGGGGGGGCTGC‑CCAGCCACTCAATGGCAGCGCCggg < 2:3890420/76‑1 (MQ=255) tCTTATCAGGGGTTTTTTTTGCGGCGCGTTGGCCGGTAGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCgggg < 1:5334170/76‑1 (MQ=255) tttACCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTAGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGc < 1:4505537/75‑1 (MQ=255) taaTCAGCGGCTTTCTTTGCGGCGCGGTGGCCGGTGGGGCGGGGTCTC‑CCAGCCCCTCAATGGCAGCGCCGGGGCt < 1:4737705/74‑1 (MQ=255) ttgcAGCGGTTTTATTGGCGGCGGGGTGGCGGGTATGGGGGGCCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTa < 1:1344469/73‑1 (MQ=255) tATCAGCGTCTTTATTTGCGGCGGGGTTGCGGGTTGGGCGGGCCTCG‑CCAGCCACTCAATGGCAGCGCCGGGGCTa < 2:4561824/76‑1 (MQ=255) aTCAGCGGTTTTATTGGCGGCGCGGTTGCCGGTAGGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTAt < 1:5356120/76‑1 (MQ=255) tcggCGGTTTTATTTGCGGCGCTGTCGCCGTTAGGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTAtt < 2:6155099/73‑1 (MQ=255) tCAGCGGTTTTATTTGCGGCGCGGTTGCGGGTATGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTAtt < 2:750449/76‑1 (MQ=255) cAGCGGCTTTATTTGCGGCGGGGTTGCGGGTAGGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTAttt < 1:1397922/76‑1 (MQ=255) tATTTGGGGCGCGGTTTCGGGTAGGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTg < 2:3123684/76‑1 (MQ=255) ttggcggcgGGTTGGCGGGTAGGGCGGGCTTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑a < 1:2017301/73‑1 (MQ=255) gggcgTGGTGGCCGGTTGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AAttc < 1:3673147/75‑1 (MQ=255) gggcgCGGTTGCGGGTAGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AAttc < 1:31558/75‑1 (MQ=255) cggcgCGGTTCCGGATATGGCGGGCCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AAttc < 1:5573638/76‑1 (MQ=255) gcgcGGTGCCCGTTAGGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AAttctt < 2:450213/76‑1 (MQ=255) gcccgGTGGCCGGTATGGGGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AAttctt < 2:5626107/73‑1 (MQ=255) gggTTGCCGGTAGGGCGGGCCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGa < 1:1112960/75‑1 (MQ=255) ttGCCGGTATGGCGGGGCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATcc < 2:4645079/76‑1 (MQ=255) tGCGGGTAGGGCGGGGCTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCg < 1:4050695/76‑1 (MQ=255) gggCGGGTCTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATcc < 2:4407047/75‑1 (MQ=255) gggggggCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCa < 1:2784602/73‑1 (MQ=255) gggggggCTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCa < 1:1208107/73‑1 (MQ=255) ggCGGGGCTTG‑CCAGCCACTCAATGGCAGGGCCGGGGGTATTTTACAACGTGCAAATTTTTCGTTCCGCGGAATca > 2:2706160/1‑75 (MQ=255) ggCGGGGCTTG‑CCAGCCACTCAATGGCAGGGCCGGGGCTATTTACCAGCGTGG‑AAATCTTCCATCCCGCGAAACCa > 2:2240595/1‑76 (MQ=255) ggCGGGGCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTTTTTGATCCGGCGGATccc > 2:1184668/1‑75 (MQ=255) ggCGGGGCTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCa < 1:1534283/76‑1 (MQ=255) ggggggCCTGCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCa < 2:378937/74‑1 (MQ=255) gggtCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATg < 1:2139572/76‑1 (MQ=255) ggggCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATg > 1:366900/1‑76 (MQ=255) ggggCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCATCGTGC‑AATTCTTCGATCCGGCGAATCCAATg < 2:6168669/76‑1 (MQ=255) ggggCTGC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATg < 1:6066234/76‑1 (MQ=255) ggggCTCC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATg < 2:6031893/76‑1 (MQ=255) gctgCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATg < 2:1965515/73‑1 (MQ=255) gggCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGa < 1:6189702/76‑1 (MQ=255) gggCTTC‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGa < 1:4642109/76‑1 (MQ=255) gggCTTC‑CAAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGa < 1:388963/76‑1 (MQ=255) ggCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAg < 1:3152048/76‑1 (MQ=255) gCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGc > 2:2272733/1‑76 (MQ=255) gCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCCGCGAATCCAATGAGc > 2:367742/1‑76 (MQ=255) ttG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGCAt > 2:4508876/1‑76 (MQ=255) cg‑ccagccaCTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGCATc < 2:1271093/75‑1 (MQ=255) g‑ccagccaCTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGCATCg > 1:3989906/1‑76 (MQ=255) | CGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTG‑CCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGC‑AATTCTTCGATCCGGCGAATCCAATGAGCATCG > REL606/2735533‑2735683 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |