| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,213,803 | 0 | T | C | 40.3% | 45.2 / 89.2 | 67 | S253P (TCT→CCT) | rpoD | RNA polymerase, sigma 70 (sigma D) factor |
| Reads supporting (aligned to +/- strand): ref base T (22/18); new base C (13/14); total (35/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.25e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCT > NC_000913/3213756‑3213852 | cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCt > 1:3469775/1‑50 (MQ=255)cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt > 1:2401323/1‑50 (MQ=255) gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCt > 1:3422720/1‑49 (MQ=255) gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTg < 1:1465271/50‑1 (MQ=255) aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg > 1:4029339/1‑49 (MQ=255) ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGt < 1:1724866/50‑1 (MQ=255) ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg > 1:2496611/1‑49 (MQ=255) ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt < 1:1063137/50‑1 (MQ=255) gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa < 1:3696032/50‑1 (MQ=255) tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTAt > 1:2326955/1‑50 (MQ=255) tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt < 1:838113/50‑1 (MQ=255) gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTc < 1:3607729/50‑1 (MQ=255) cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCa < 1:3072691/50‑1 (MQ=255) cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa < 1:584526/50‑1 (MQ=255) gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCaaa < 1:1715426/50‑1 (MQ=255) tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAAc < 1:3017485/50‑1 (MQ=255) cACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCaaa > 1:641781/1‑48 (MQ=255) cACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCaaa > 1:262842/1‑48 (MQ=255) cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa > 1:697186/1‑48 (MQ=255) aCGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAg < 1:2464318/50‑1 (MQ=255) aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg < 1:2841214/50‑1 (MQ=255) cGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGt < 1:2047840/50‑1 (MQ=255) cGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGt < 1:581553/50‑1 (MQ=255) cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt > 1:1025677/1‑50 (MQ=255) gCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGtt < 1:815583/50‑1 (MQ=255) gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt > 1:911825/1‑50 (MQ=255) gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt < 1:1481523/50‑1 (MQ=255) gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt > 1:2045325/1‑50 (MQ=255) cTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTc > 1:2712720/1‑50 (MQ=255) cTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTc > 1:269968/1‑50 (MQ=255) cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc > 1:355423/1‑50 (MQ=255) cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc < 1:157799/50‑1 (MQ=255) ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc < 1:3694235/48‑1 (MQ=255) cGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGcc > 1:4303019/1‑50 (MQ=255) cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc < 1:2060144/50‑1 (MQ=255) cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc > 1:1928071/1‑50 (MQ=255) gCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCt < 1:265423/50‑1 (MQ=255) gCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCt < 1:122961/50‑1 (MQ=255) gCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCt < 1:2555910/50‑1 (MQ=255) aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTg < 1:2665959/50‑1 (MQ=255) aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTg > 1:3847893/1‑50 (MQ=255) aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTg > 1:3644702/1‑50 (MQ=255) aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTg < 1:2568638/50‑1 (MQ=255) aGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTg < 1:4255058/50‑1 (MQ=255) aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt < 1:2119841/49‑1 (MQ=255) aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt < 1:3192788/49‑1 (MQ=255) aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg < 1:228237/50‑1 (MQ=255) ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc < 1:4225767/50‑1 (MQ=255) gAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGcc < 1:1879335/50‑1 (MQ=255) agagATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg > 1:129596/1‑49 (MQ=255) gagaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg > 1:3148131/1‑48 (MQ=255) gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg > 1:3652011/1‑48 (MQ=255) gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGaa < 1:602996/50‑1 (MQ=255) agaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg > 1:283899/1‑50 (MQ=255) aTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg > 1:1120560/1‑48 (MQ=255) aTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg > 1:4152952/1‑48 (MQ=255) tCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg > 1:518466/1‑50 (MQ=255) tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg > 1:116530/1‑50 (MQ=255) ccTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt > 1:2268827/1‑50 (MQ=255) gAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt > 1:1673763/1‑49 (MQ=255) gAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt > 1:1435803/1‑49 (MQ=255) gAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt > 1:745550/1‑49 (MQ=255) gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt > 1:3265691/1‑49 (MQ=255) aCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc > 1:555241/1‑49 (MQ=255) aCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACt < 1:1596409/50‑1 (MQ=255) aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc > 1:2139716/1‑49 (MQ=255) gTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTAcc > 1:3546510/1‑50 (MQ=255) tCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGGAGCAGTTTGACTACCt < 1:3457691/50‑1 (MQ=255) tCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTAc < 1:2509119/48‑1 (MQ=255) tCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCt > 1:509528/1‑50 (MQ=255) | CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCT > NC_000913/3213756‑3213852 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |