| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 888,014 | 0 | T | G | 23.3% | 90.4 / 10.2 | 43 | noncoding (44/81 nt) | rybB | small regulatory RNA RybB |
| Reads supporting (aligned to +/- strand): ref base T (6/27); new base G (10/0); total (16/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.18e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATCCTGGCC > NC_000913/887877‑888150 | ggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTg < 1:477822/140‑1 (MQ=255) aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 2:1300025/140‑1 (MQ=255) aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 2:139977/140‑1 (MQ=255) aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 2:489317/140‑1 (MQ=255) aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 2:339041/140‑1 (MQ=255) cATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcac > 2:900423/1‑140 (MQ=255) ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAAAggggg > 2:274424/1‑139 (MQ=255) ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCTCCAAACAGGGa > 1:568867/1‑140 (MQ=255) ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCAAAAATggggg > 1:686287/1‑139 (MQ=255) ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGa > 1:1026895/1‑140 (MQ=255) ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAggggg > 2:1216468/1‑139 (MQ=255) cccgTAAGCCGCTATCGGTCTAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 2:1269906/139‑1 (MQ=255) accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 1:900423/140‑1 (MQ=255) accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 2:897713/140‑1 (MQ=255) accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt < 1:994804/140‑1 (MQ=255) gCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg < 1:1216468/140‑1 (MQ=255) tATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCa < 2:1026895/140‑1 (MQ=255) tCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCCCAAAAAAGGGGCAACAAAGAAAAGCAgg > 2:1277238/1‑139 (MQ=255) gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa < 2:37241/140‑1 (MQ=255) gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa < 2:751841/140‑1 (MQ=255) gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa < 1:968718/140‑1 (MQ=255) gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAgggggg > 2:1044924/1‑138 (MQ=255) agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAtt < 2:275973/140‑1 (MQ=255) gATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATg < 2:84420/140‑1 (MQ=255) agGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGaa < 1:208884/140‑1 (MQ=255) agGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAAAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGaa > 2:208884/1‑140 (MQ=255) gTTGCAGGGTATTAGATAAGTTTTAGATAACAAAAAACCCATCAGCCTTGAACCGATATGGCCGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGTAGTGGCGCTAATTAAGACTGATGCCCTGCGGAAAAg < 2:1255227/140‑1 (MQ=255) ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGt > 2:529459/1‑140 (MQ=255) ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGACATCAAAAAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGaaaaat > 2:300426/1‑138 (MQ=255) ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAAAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAgg > 2:826619/1‑139 (MQ=255) gtagtaGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCgg < 1:353409/140‑1 (MQ=255) gtagtaGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCgg > 2:923950/1‑140 (MQ=255) aaGTTTTAGATAACAAAAAACCCATAAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAAGGGACATCAAAAAAAAGCAGGGGGACTAATTAAGACTGATGCCCTGCGGAAAAGGTCTTCGGTTGTGCaaa > 2:824977/1‑140 (MQ=255) gTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaaaa > 1:1297915/1‑140 (MQ=255) ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACCTCAAAGAAAAACCGGGGCACTAATTAAGACCGATGCCCTGCGGGAAAGTTCTGCGGTTGTGGAAAAAAAt > 2:1080911/1‑140 (MQ=255) cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt > 1:276577/1‑140 (MQ=255) aaaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATTGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTc < 2:1250933/140‑1 (MQ=255) aaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAg < 1:529459/140‑1 (MQ=255) aCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGc < 2:276577/140‑1 (MQ=255) aTCAACCTTGAACCGAAATGGCGGGGGTGATTGGGTCCTCCAAAAAGGGAAAAAAAACAAAAACTGTGGCACCAAAGAAAAATTTTGTGGTGCCGAGAAATACTTCGGTTGGGGAAAAAAAATTTATTTTCAgggtgact > 2:907533/1‑135 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAATTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACtt > 1:477616/1‑140 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTc < 1:560712/130‑1 (MQ=255) tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTc > 2:560712/1‑130 (MQ=255) ttGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGttt < 1:826619/140‑1 (MQ=255) tGAACCGAAATGGCGGGGGTGATTAGGCCCCCAAAAAAGGGCCAACCAAGAAAAGAAGTGGCACTAATTTAGACTTAACCCCTCCGGAAAAGTTTTGCGCGTTTTCAAAAAAATTTCATTTTTACGGCAACTTCCAGttt > 2:1107563/1‑140 (MQ=255) aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGttt < 2:311821/137‑1 (MQ=255) aaCCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGttt > 1:311821/1‑137 (MQ=255) ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTTAGGGCAACTT‑CAGTTTTATcc > 2:481505/1‑140 (MQ=255) aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTaa > 2:139749/1‑140 (MQ=255) aaaTGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTaa > 2:1123314/1‑140 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATc < 2:686287/140‑1 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATc < 1:907533/140‑1 (MQ=255) aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATc < 1:139749/140‑1 (MQ=255) ggggTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGGTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATCCTGGcc > 1:651385/1‑140 (MQ=255) | GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTT‑CAGTTTTATCCTAATCCTGGCC > NC_000913/887877‑888150 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |