| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 16,981 | G→A | noncoding (30/59 nt) | sokC → | small regulatory RNA antitoxin SokC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 16,981 | 0 | G | A | 94.7% | 40.8 / ‑3.8 | 19 | noncoding (30/59 nt) | sokC | small regulatory RNA antitoxin SokC |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (10/8); total (11/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01 | |||||||||||
TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCG‑GATGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCA > NC_000913/16736‑17176 | ccccccgctgtgacttttcTCGGTTTCGTAACCCGTAAAAACATCAACCACGGTCTGGCCA‑TTCGTGATG‑GCACCCCAAAAGAGTTTTTTGCTCGTTCCCGTCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCCTCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATACAGGGGGCCTCGGTTTGA‑TAGTa < 2:55964/231‑1 (MQ=255) ccgTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCG‑GATGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAACATACCCTTTTT‑CCTCTTTTCTGCTCTCTTCTCCTTA‑CCCTTACAGTCAGTAAACGGCCCCCTACATGTGTTCCACATCTAGGGGGCCTCGGGTTGT‑TAGTAAAATATcc > 2:175364/1‑250 (MQ=255) ttGCCCTTATTCGGATTCGTAACCCTTGAAACCAGCACCCCCCGTCTGCCCAGTTCG‑GATGTGACCCTCACA‑GAGGTTTTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACt < 1:215000/250‑1 (MQ=255) tgccctTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCGGTTGG‑GAGGTGACCCTCACA‑GAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTc < 1:182613/246‑1 (MQ=255) gATTCGTAAGCCGTGCAAACAGCACACTACGTCTGGCCAGTTCG‑GTTGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCCCTACGGCGGTGATAAAGAAGAAAATCAGGGCGGCAATCATCACCCTATG‑CTGCTTTATTGCTCTCTTCTCCTTG‑TCCCTACTGTCCGTAAGAGGCACTCTTCATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTttct > 2:58058/1‑251 (MQ=255) aGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCG‑GATGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGAGATACAGATGACGATCAGGGCGACAAACATCGCCTTATG‑CGGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAAGCACTCTACATGTGTTCAGCATATAGGGGGGCTCGGGTTGG‑TAGTAATAAATCACTCGGGGGTTTTTTCTATCTGCCg > 1:263543/1‑251 (MQ=255) ccAGTTCG‑GAGGTGAACCACCCA‑GGGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑TCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCt < 1:58058/251‑1 (MQ=255) gTTCG‑GATGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCTCa > 1:346756/1‑250 (MQ=255) cggatgtcatcatctcagggtccTTTTCTCGTACCCCGCGCCGCACCTACGGTGGTGATACAGATTACGATCATGGCGACAATCATCGCCTTAAG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGGCAGGAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGTGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATATCCGCAAGCa < 2:345543/229‑1 (MQ=255) tgtgAACCTCACA‑GAGGTTTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCTCAAGCAccc < 1:310384/250‑1 (MQ=255) gtgAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAAAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCTCAAGCACccgc > 1:184312/1‑251 (MQ=255) cAGCGCCGCCACTACGGCGGTGATAAAGATGACGATCCGGGCGACAAACATTGCCTTTTG‑CTGCTTCATTGCCCTATTCTCCTTG‑ACCTTACGGTCCGTAAAAGGCCCTCTACATGTGGTCAGCATATAGGGGGCCTCGGGGTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATTACTGGGACAGATAGCCTCAAGCACCCCCCGCTATTATATCGCTCTCTTTAACCCAttt > 1:383673/1‑251 (MQ=255) gcgggtaacgaattGACGATCAGGGCGAAATTCTTC‑CCTTTTGGTTGCTTCATTGCTCTCTT‑TCCTTGCACCTTGCGGTCAGAAAGATGCTCTCTCCATGTGTAAGGCATATGGGGGGCCT‑GGGTTGAATAGAAAAATACCCCTCGGGGTTTTTCTTAAGCTGC‑GTTCCAGTTAAGGCCTGAGACAAAAACCCACAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTa < 1:234487/238‑1 (MQ=255) gACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGATGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCt < 2:64297/148‑1 (MQ=255) gACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCt > 1:64297/1‑148 (MQ=255) cGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGACAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTCTTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAgccgactc > 1:53027/1‑244 (MQ=255) cAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTATTCCTGACGAAGCCTCGCATTTTTGTgg > 2:178962/1‑251 (MQ=255) ttaagcagctTCATTGCTCTCTACGCATTG‑AACTTACGGTCAGTAAGAGGCACTCTACATGAGGTCAGCATATAGGGGGCCTCGGTTTGA‑TAGTAAAATATCACTCGGGGCTTTTCTCGCTCTGACGTTC‑GGCTAATGACTAATGCAGATAGCATCAAGCACCAGACGCAATTATATCGCTCACGTTAACACATTGTGTTTTATAGATTCTAATCATGAAGACACCCAGCATTTTCGTGGCATAAtatttt > 2:314108/8‑251 (MQ=255) cAGCATATAGGAGGCCTCGGGTTGA‑TGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCa > 2:182476/1‑222 (MQ=255) | TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCG‑GATGTGAACCTCACA‑GAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATG‑CTGCTTCATTGCTCTCTTCTCCTTG‑ACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTC‑AGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCA > NC_000913/16736‑17176 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |