| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 888,014 | 0 | T | G | 31.6% | 222.5 / 23.3 | 95 | noncoding (44/81 nt) | rybB | small regulatory RNA RybB |
| Reads supporting (aligned to +/- strand): ref base T (19/46); new base G (30/0); total (49/46) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.20e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAT > NC_000913/887965‑888063 | aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGt > 1:6798713/1‑50 (MQ=255) gTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTg < 1:6957973/50‑1 (MQ=255) tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt > 1:2018816/1‑50 (MQ=255) tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGAt > 1:4428267/1‑50 (MQ=255) tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGAt > 1:344192/1‑50 (MQ=255) ttAGATAACAAAAATCCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 1:7165628/49‑1 (MQ=255) ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATg > 1:1269952/1‑50 (MQ=255) ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATg > 1:6713857/1‑50 (MQ=255) ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATg > 1:1834561/1‑50 (MQ=255) tAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg > 1:366554/1‑50 (MQ=255) tAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATg < 1:7041301/49‑1 (MQ=255) tAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGAAgg > 1:549413/1‑50 (MQ=39) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATggg > 1:6570731/1‑50 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg < 1:3715350/49‑1 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg < 1:822551/49‑1 (MQ=255) aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATggg > 1:6383181/1‑50 (MQ=255) gATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATggg < 1:4823393/49‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCt < 1:4776005/50‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCt < 1:2496826/50‑1 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATgggtc > 1:402825/1‑48 (MQ=255) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATTGGCt > 1:3102512/1‑50 (MQ=39) aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCt > 1:7437727/1‑50 (MQ=255) tAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTTATGGGccc > 1:7242971/1‑48 (MQ=38) tAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTc > 1:2881040/1‑50 (MQ=255) aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCca < 1:4889537/50‑1 (MQ=255) aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTTca > 1:6965780/1‑50 (MQ=255) aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCca > 1:7333101/1‑50 (MQ=255) cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcac < 1:3224085/50‑1 (MQ=255) cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcac < 1:7319055/50‑1 (MQ=255) cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCca > 1:3354131/1‑49 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca < 1:5228746/50‑1 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca < 1:685402/50‑1 (MQ=255) aaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTCGATGGGCTCcaca > 1:6137214/1‑49 (MQ=255) aaaaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGTTGGGCCCccca > 1:5485464/1‑49 (MQ=255) aaaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACaa > 1:3494709/1‑49 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGTTTATGGGCTCCACaaa > 1:1757685/1‑49 (MQ=255) aaaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACaaaa > 1:662579/1‑50 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAt < 1:6505503/50‑1 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAt > 1:738827/1‑50 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAt < 1:2523509/50‑1 (MQ=255) aaCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCTTAAAAt > 1:1977748/1‑50 (MQ=255) cccATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCCCCCCAAAATgg > 1:640371/1‑50 (MQ=255) ccATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg > 1:2790810/1‑50 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg < 1:2377437/50‑1 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg < 1:2409812/50‑1 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg < 1:2894921/50‑1 (MQ=255) cATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg > 1:3853690/1‑49 (MQ=255) aTCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGa < 1:6635914/50‑1 (MQ=255) aTCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGa > 1:2945247/1‑50 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGAc < 1:1724899/50‑1 (MQ=255) cAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACa > 1:2953838/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 1:5687193/50‑1 (MQ=255) aaCCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACAt > 1:7245591/1‑50 (MQ=255) aaCCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACAt > 1:2196622/1‑50 (MQ=255) aCCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATc < 1:1219136/50‑1 (MQ=255) aCCTTGAACCGAAATGGCGGGGGTGAAGGGGCCCACAAAATGGGGAcca > 1:5176311/1‑47 (MQ=255) ccTTGAACCGAAATGGCGGGGTTGATGGGGTCCACAAAATGGGGACATc > 1:6624714/1‑49 (MQ=255) ccTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCa < 1:4494085/50‑1 (MQ=255) ccTTGAACCGAAATGGCGGGGTTGAGGGGCTCCACAAAATGGGGACATc > 1:6474331/1‑49 (MQ=255) cTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCaa < 1:4066673/50‑1 (MQ=255) cTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCa > 1:139182/1‑49 (MQ=255) cTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCa > 1:3400072/1‑49 (MQ=255) cTTGAACCGAAATGGCGGGGGTGATGGGCCCCACAAAATGGGGACATCa > 1:1279880/1‑49 (MQ=255) tGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCa > 1:199581/1‑47 (MQ=255) tGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCa > 1:2321537/1‑47 (MQ=255) gAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAg < 1:876357/49‑1 (MQ=255) gAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAg < 1:3816145/49‑1 (MQ=255) gAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGACATCAAAg > 1:7305220/1‑49 (MQ=39) aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaa < 1:2338710/50‑1 (MQ=255) aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaa < 1:789744/50‑1 (MQ=255) aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaa < 1:1635258/50‑1 (MQ=255) ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaaaa < 1:3543809/50‑1 (MQ=255) ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATAAAAGaaaa < 1:6458888/50‑1 (MQ=255) cGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg > 1:2267150/1‑50 (MQ=255) gAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGc < 1:3912135/50‑1 (MQ=255) aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg > 1:5270252/1‑48 (MQ=255) aaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAg < 1:6375537/50‑1 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGt < 1:1872875/50‑1 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGt < 1:5591805/50‑1 (MQ=255) ggCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTgg > 1:995139/1‑50 (MQ=255) gCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTgg < 1:6545527/49‑1 (MQ=255) gCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTgg < 1:4172745/49‑1 (MQ=255) gCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTgg < 1:1451623/49‑1 (MQ=255) ggggTTGATGGGCTCCACAAGATGGGGACATCAAAGAAAAGCAGTGGCAc < 1:4789955/50‑1 (MQ=255) ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc > 1:6822663/1‑50 (MQ=255) ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc < 1:4904002/50‑1 (MQ=255) ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc > 1:6153172/1‑50 (MQ=255) gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACt < 1:5667694/50‑1 (MQ=255) gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACt < 1:3364612/50‑1 (MQ=255) gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACt < 1:3552154/50‑1 (MQ=255) ggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTa < 1:4666912/50‑1 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa < 1:3093348/50‑1 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa > 1:3766337/1‑50 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa < 1:6779986/50‑1 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa > 1:2477971/1‑50 (MQ=255) gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTaa > 1:2884759/1‑50 (MQ=255) ttGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAt > 1:5122401/1‑50 (MQ=255) | AAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAT > NC_000913/887965‑888063 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |