Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,008,240 | 0 | T | G | 25.5% | 96.4 / 11.4 | 47 | T41P (ACC→CCC) | ECB_01907 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (5/30); new base G (12/0); total (17/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGTGTTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAGGA > REL606/2008167‑2008315 | gtgtgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGgtgt < 2:579632/76‑1 (MQ=255) gtgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAg < 2:635271/76‑1 (MQ=255) gtgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAg < 2:179709/76‑1 (MQ=255) tgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGt < 1:120934/76‑1 (MQ=255) gtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGtt < 2:498264/76‑1 (MQ=255) ttCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTa < 2:298986/76‑1 (MQ=255) cAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGc > 1:966913/1‑76 (MQ=255) aaaCCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCa > 2:231513/1‑76 (MQ=255) aaCCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCAt < 1:283952/76‑1 (MQ=255) ccAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATaa < 2:265284/76‑1 (MQ=255) aGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCATAAgg > 1:1077169/1‑76 (MQ=255) gtgtCGTATTGCTGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCtgtg < 2:588491/76‑1 (MQ=255) gtgtCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCtgtg < 2:818729/76‑1 (MQ=255) gtCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCtgtggc > 1:165816/1‑74 (MQ=255) tCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCCGTGccg > 1:707764/1‑74 (MQ=255) cGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGt < 1:702071/76‑1 (MQ=255) gTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTa < 2:239081/76‑1 (MQ=255) aTTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACg < 2:193630/76‑1 (MQ=255) ggATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCCTAAGGTGCGGGCTGTGCAGTACGGgcgc > 1:890168/1‑76 (MQ=255) aTGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGGTAGCATAAGGGGCGGCCTGGGCAGGACGGGCGccc > 2:652619/1‑76 (MQ=255) tGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCg < 1:984992/76‑1 (MQ=255) ggTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCATAAGGGGCGGCCCGTGCAGGACGGGGGccccg > 2:1021513/1‑74 (MQ=255) gTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTc < 2:204312/76‑1 (MQ=255) gTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGGGCGGCCTGTGCAGTACGGGCGCCCCTc > 1:457191/1‑76 (MQ=255) gTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCTGTGCAGTACGGGCGCCCGTc > 2:957529/1‑76 (MQ=255) tATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCTGGGCAGTACGGGCGCCCCTcc > 2:658819/1‑75 (MQ=255) tCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGa < 2:14022/76‑1 (MQ=255) cTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGAc < 1:189820/76‑1 (MQ=255) cTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGGCAGAc > 2:1109950/1‑76 (MQ=255) gTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCTGTGCCGGACGGGCGCCCCTCAGACGGcc > 2:237899/1‑76 (MQ=255) tGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCCGTGGAGTACGGGCGCCCGGCCGACGGCCg > 2:333790/1‑76 (MQ=255) gCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGt < 2:1087574/76‑1 (MQ=255) cACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGtt < 2:170065/76‑1 (MQ=255) cACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGtt < 2:101382/76‑1 (MQ=255) cACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGtt < 1:942893/76‑1 (MQ=255) aCTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTg < 2:311391/76‑1 (MQ=255) tGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGcc < 1:449188/76‑1 (MQ=255) ccACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCCGACGTCCGTTGccct > 1:288433/1‑74 (MQ=255) cacacaAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATc < 2:880413/76‑1 (MQ=255) acacaAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCa < 2:605633/76‑1 (MQ=255) aGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGc < 1:516337/76‑1 (MQ=255) aGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGc < 1:92959/76‑1 (MQ=255) gTGGCGGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCg < 2:88770/76‑1 (MQ=255) gTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCg < 1:932594/76‑1 (MQ=255) tGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGa > 2:963056/1‑76 (MQ=255) tGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGa > 1:585957/1‑76 (MQ=255) ggCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGAt > 2:713028/1‑76 (MQ=255) gCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATa < 1:996812/76‑1 (MQ=255) cAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATaa > 2:138106/1‑76 (MQ=255) cAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATaa < 2:1139231/76‑1 (MQ=255) cAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATaa < 1:132957/76‑1 (MQ=255) tgtAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAGGa > 2:1092617/1‑76 (MQ=255) tgtAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAGGa > 2:895991/1‑76 (MQ=255) tgtAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAGGa > 1:436613/1‑76 (MQ=255) | GTGTGTTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAGGA > REL606/2008167‑2008315 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |