| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,669,508 | 0 | T | A | 30.6% | 145.4 / 18.2 | 62 | I29F (ATT→TTT) | ynfL | putative DNA‑binding transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base T (17/26); new base A (19/0); total (36/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.48e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCA > NC_000913/1669459‑1669557 | aTTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAat < 1:4476100/50‑1 (MQ=255) tttGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAat‑a < 1:4615461/50‑1 (MQ=255) tttGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAat‑a < 1:7260657/50‑1 (MQ=255) tttGTTGCTCCAGCGCCTGAATCTGCTGACTTAGAGGCGGTTGCGAAat‑a < 1:2936807/50‑1 (MQ=255) ttgttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGATGGGAAaa‑at > 1:3848480/1‑50 (MQ=38) tgttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAat‑at < 1:1620046/49‑1 (MQ=255) gttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGGAAA‑Att > 1:2732404/1‑49 (MQ=38) gttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑Att < 1:1340980/49‑1 (MQ=255) tgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGACTAGAAAA‑ATTCAg > 1:7746898/1‑50 (MQ=25) gCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAg < 1:185104/49‑1 (MQ=255) gCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCaag > 1:6568036/1‑48 (MQ=39) cTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGGAAA‑ATTTAGGc > 1:7384877/1‑50 (MQ=38) cTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAgg < 1:6768665/49‑1 (MQ=255) cTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAgg < 1:2123482/49‑1 (MQ=255) cTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGc < 1:7885269/50‑1 (MQ=255) tCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTTCGAAAA‑ATTTAgggg > 1:5960653/1‑48 (MQ=37) tCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGTGGAATAATTCAgcg > 1:7806245/1‑48 (MQ=25) tCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGGAAT‑ATTCAgggg > 1:5480796/1‑48 (MQ=39) tCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGGAAA‑ATTTAgggg > 1:3884335/1‑48 (MQ=37) ccAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTTAGgggg > 1:7357671/1‑50 (MQ=38) ccAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAGGCgg > 1:900924/1‑50 (MQ=255) cAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCgg < 1:2014658/49‑1 (MQ=255) aGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTTAgggggg > 1:3151870/1‑49 (MQ=255) aGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTTAGGCggg > 1:766188/1‑49 (MQ=255) aGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAGGGGGGc > 1:3646180/1‑50 (MQ=255) aGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAGGCggg > 1:3718711/1‑49 (MQ=255) gcgcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCa < 1:1354273/50‑1 (MQ=255) cgcCTGAATCTGCTGACTTAGCGGCGGGTGCGGAAA‑ATTTAGGCGGGGAg > 1:5900301/1‑50 (MQ=255) gcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAgc < 1:2787587/50‑1 (MQ=255) gcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAgc < 1:4756031/50‑1 (MQ=255) gcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAAGCGGGGAgc > 1:3037229/1‑50 (MQ=255) gcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAAGCGGGCAgc > 1:2244950/1‑50 (MQ=255) ccTGAATCTGCTGACTTAGCGGCGGTTGCGGAAA‑ATTTAGGCGGGCAgcg > 1:186860/1‑50 (MQ=255) ccTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAgcg < 1:1465846/50‑1 (MQ=255) cTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAgcgg < 1:3791118/50‑1 (MQ=255) cTGAATCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTCAGGCGGGCCgcgg > 1:8201295/1‑50 (MQ=255) tGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTTAAGCGGGGAgcggc > 1:6531226/1‑50 (MQ=255) aaTCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAgcggc < 1:908361/48‑1 (MQ=255) aaTCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGgcgc < 1:8313562/50‑1 (MQ=255) aaTCTGCTGACTTAGCGGCGGTTGCGAAAT‑AATTAGGCGGGGAGGGgggc > 1:2667977/1‑50 (MQ=255) tCTGCTGACTTAGCGGCGGTTGCGAAAA‑ATTTAGGCGGGCAGCGgcgcgc > 1:3320432/1‑50 (MQ=255) tCTGCTGACTTAGCGGCGGTTGCGAAAA‑AATTAGGCGGGGGGCGgcgcgc > 1:5003014/1‑50 (MQ=255) ctgACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGccc < 1:1491556/48‑1 (MQ=255) ctgACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGcc < 1:7361907/47‑1 (MQ=255) ctgACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCg < 1:4370544/49‑1 (MQ=255) cTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATg < 1:7537209/50‑1 (MQ=255) aGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGc < 1:2985317/48‑1 (MQ=255) gcggcgGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAg < 1:819630/49‑1 (MQ=255) ggcggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGctc > 1:1802950/1‑50 (MQ=255) gcggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGctct < 1:7367464/50‑1 (MQ=255) gcggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGctct > 1:7278339/1‑50 (MQ=255) cggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCtt < 1:595264/50‑1 (MQ=255) cggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCtt > 1:1082102/1‑50 (MQ=255) cggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCtt > 1:6261816/1‑50 (MQ=255) ggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTc < 1:7516911/50‑1 (MQ=255) ggTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTc > 1:2615372/1‑50 (MQ=255) gTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTcc > 1:8204697/1‑50 (MQ=255) gTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTcc > 1:1291810/1‑50 (MQ=255) ttGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCg < 1:367494/50‑1 (MQ=255) tGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCg > 1:7814596/1‑49 (MQ=255) tGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCg > 1:5968983/1‑49 (MQ=255) gCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgc > 1:806619/1‑49 (MQ=255) gAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaac < 1:2844606/50‑1 (MQ=255) gAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaac < 1:6874208/50‑1 (MQ=255) gAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaa > 1:5034715/1‑49 (MQ=255) gAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaa > 1:1237734/1‑49 (MQ=255) aaaT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaaca < 1:7733446/50‑1 (MQ=255) aaaT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaa > 1:1135945/1‑48 (MQ=255) at‑atTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagc < 1:3839679/50‑1 (MQ=255) at‑atTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacag > 1:5375912/1‑49 (MQ=255) t‑atTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagca < 1:7610029/50‑1 (MQ=255) t‑atTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagca < 1:7345048/50‑1 (MQ=255) | ATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAAT‑ATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCA > NC_000913/1669459‑1669557 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |