Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 28.6% | 142.4 / 17.5 | 63 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (8/37); new base G (18/0); total (26/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.03e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCAC > REL606/3983650‑3983750 | gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 1:5387237/51‑1 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 1:2548537/51‑1 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 1:4530029/51‑1 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 1:5660611/51‑1 (MQ=255) ggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGcc < 1:4431805/51‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTcc < 1:1639381/50‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTccc > 1:595024/1‑51 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCca < 1:4678023/51‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTccc > 1:5728362/1‑50 (MQ=255) gcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTccc > 1:390818/1‑48 (MQ=39) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaaa < 1:552206/51‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAt < 1:3815158/51‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAt < 1:3452146/51‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAt < 1:1668709/51‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAt < 1:6147423/51‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACAAAt > 1:4988738/1‑51 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAtt < 1:5282075/51‑1 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAtt < 1:5023098/51‑1 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACaaa > 1:6375966/1‑49 (MQ=39) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCcac < 1:1164872/45‑1 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAttt < 1:1003215/51‑1 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAttt < 1:6074968/51‑1 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGAc < 1:3072758/51‑1 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCaaa > 1:5341811/1‑45 (MQ=38) cTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACaaa > 1:6531300/1‑44 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGa < 1:4188478/51‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGa < 1:5671557/51‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAttt > 1:3698897/1‑46 (MQ=37) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCaaa > 1:6131162/1‑43 (MQ=37) aCGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 1:983294/51‑1 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 1:5494652/51‑1 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACaaa > 1:258922/1‑42 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATc < 1:1144787/51‑1 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATc < 1:2542702/51‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCg < 1:6293347/51‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCg < 1:2569542/51‑1 (MQ=255) gCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAttt > 1:5927119/1‑42 (MQ=255) gCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGc < 1:578048/51‑1 (MQ=255) cTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCaaa > 1:2320717/1‑38 (MQ=25) cTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCGCt > 1:809025/1‑51 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTg > 1:5394214/1‑51 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAttt > 1:3893952/1‑39 (MQ=37) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAttt > 1:4353616/1‑39 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGga < 1:1333081/51‑1 (MQ=255) tATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGc > 1:5441490/1‑46 (MQ=255) ttGGTGGCATCGGTAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt > 1:5333339/1‑51 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:4847831/51‑1 (MQ=255) tggtggCATCGGCAGGGGGGCCTCCCACAAAttt > 1:4128794/1‑34 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 1:4388151/51‑1 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 1:2460923/51‑1 (MQ=255) gtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGAc < 1:6013426/51‑1 (MQ=255) tggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACg < 1:6675011/51‑1 (MQ=255) tggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACg < 1:2851213/51‑1 (MQ=255) ggCATCGGCAGGGGGGCCTCCCACAAAttt > 1:2443253/1‑30 (MQ=37) cATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCa < 1:1391192/51‑1 (MQ=255) cATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCa > 1:6192472/1‑51 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 1:4538732/51‑1 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgc < 1:2422610/51‑1 (MQ=255) gCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg > 1:6524267/1‑51 (MQ=255) gCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg > 1:7239785/1‑51 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa > 1:6423592/1‑51 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa > 1:7155497/1‑51 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa > 1:2983936/1‑51 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 1:6738334/51‑1 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac < 1:5122600/51‑1 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac > 1:3741974/1‑51 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac > 1:6992311/1‑51 (MQ=255) | GAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCAC > REL606/3983650‑3983750 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |