Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_017635 | 4,733,046 | 0 | A | G | 29.3% | 103.1 / 10.8 | 58 | E189G (GAA→GGA) | bdcR | TetR/AcrR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (21/20); new base G (17/0); total (38/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.76e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGGCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGAAAGG > NC_017635/4732999‑4733095 | ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAg > 2:1125127/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAg > 2:1266840/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAg > 2:272255/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAg > 2:209367/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAg > 1:908371/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAg > 1:236477/1‑50 (MQ=255) ggggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAg > 1:1703934/1‑50 (MQ=255) gggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGc > 1:1315253/1‑50 (MQ=255) gggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGc > 1:1985707/1‑50 (MQ=255) ggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAgag > 2:1814977/1‑48 (MQ=255) ggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGc > 1:924155/1‑49 (MQ=255) ggCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGc > 1:678221/1‑49 (MQ=255) cACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGa > 1:1041786/1‑49 (MQ=255) cACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGa > 1:43885/1‑49 (MQ=255) cACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGa > 1:1102883/1‑49 (MQ=255) cACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGa > 1:725907/1‑49 (MQ=255) cACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGa > 1:564481/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGAt > 1:885367/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGAt > 1:665319/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGAt > 1:645016/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGAt > 1:637411/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGAt > 1:556948/1‑49 (MQ=255) aCTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGATa < 2:440797/50‑1 (MQ=255) cAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGATaaa > 1:953566/1‑49 (MQ=255) cAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGATaaa > 1:1865194/1‑49 (MQ=255) cAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGATaaa > 1:1007942/1‑49 (MQ=255) aTCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGATAAAAAc > 1:362608/1‑50 (MQ=255) cGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGGAGAGGAAAAAAatc > 2:29531/1‑50 (MQ=25) caacTCTGTGCAACAGCTGCACTGGCGGGGGGAGCGATAAAAAAtcttct > 1:1952654/1‑50 (MQ=39) aacTCTGTGCAACAGCTGCACTGGCGGGGGGGACGATAAAAAAtattct > 2:220733/1‑49 (MQ=25) aacTCTGTGCAACAGCTGCACTGGCGGGGGGAAGGATAAAAACTCTtctc > 1:1935259/1‑50 (MQ=38) ctctGTGCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCaa < 1:1787248/50‑1 (MQ=255) tgtgCAACAGCTGAACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAg < 2:1356922/48‑1 (MQ=255) gtgCAACAGCTGCACTGGCGGGGGGGGGGGAAAAAAAACCTCTCAAgg > 2:821773/1‑48 (MQ=12) gtgCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAg > 1:1434438/1‑50 (MQ=255) tgCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAgg < 2:1025119/47‑1 (MQ=255) tgCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAgg < 1:1738049/47‑1 (MQ=255) tgCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAgg < 1:1916329/47‑1 (MQ=255) tgCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGa < 2:859655/48‑1 (MQ=255) gCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGGGt > 1:2043864/1‑49 (MQ=255) cAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTg < 2:124708/49‑1 (MQ=255) aaCAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGAt < 1:926023/50‑1 (MQ=255) aaCAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGAt < 1:1617690/50‑1 (MQ=255) aaCAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGAt < 1:1798091/50‑1 (MQ=255) cAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGc < 1:561317/50‑1 (MQ=255) cAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGc < 1:1202711/50‑1 (MQ=255) cAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGc < 2:1424671/50‑1 (MQ=255) aGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTTTCAAGGAGGGATGCt > 1:333763/1‑50 (MQ=255) gCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTg < 1:1809431/50‑1 (MQ=255) cTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTgg < 1:616361/50‑1 (MQ=255) cACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCt > 1:963355/1‑50 (MQ=255) cACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCt < 1:336011/50‑1 (MQ=255) gCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGa < 1:1605062/48‑1 (MQ=255) gCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGACTTGAt < 2:476717/49‑1 (MQ=255) ggggAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCg > 1:1658316/1‑49 (MQ=255) gggAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGa > 2:299593/1‑49 (MQ=255) ggAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGaaa < 2:144145/50‑1 (MQ=255) ggAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGa > 1:1139786/1‑48 (MQ=255) aaGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGAAAgg < 2:1774566/50‑1 (MQ=255) | GGGGCACTCAATCGAACAACTCTGTGCAACAGCTGCACTGGCGGGGGAAGCGATAAAAACTCTTCTCAAGGAGTGATGCTGGCCTTGATCCGAAAGG > NC_017635/4732999‑4733095 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |