Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,313,176 | (T)8→7 | noncoding (93/93 nt) | micF → | small regulatory RNA MicF |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,313,169 | 0 | T | . | 75.0% | 30.2 / 2.9 | 21 | noncoding (86/93 nt) | micF | small regulatory RNA MicF |
Reads supporting (aligned to +/- strand): ref base T (4/1); new base . (6/9); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01 |
AAATGCAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAAT‑AAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑TTT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTG‑CAATTATTTTTGTTGGGCTTTTTG‑TAGGTTATTTGTACAGCAAAATGGCGC‑TTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGG > NC_000913/2312926‑2313386 | caaaggggaaacttttttaAAGATTATTGCGGAATGGCGAAATAAAGCACCTAACATCAAGCAATAATAATTC‑AGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTCAAAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTTCCTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACt‑a < 1:1263962/286‑1 (MQ=255) cAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAAT‑AAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAG‑AAAACT‑ATAACGCTTGTTAAc < 1:1263961/306‑1 (MQ=255) cGAAATAAAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGaaa < 1:1263966/257‑1 (MQ=255) aaGCAATAATAATTCAAGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCa < 1:1263968/293‑1 (MQ=255) caatataatTCAAGGTTTAAAAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAACAAAAACCTTCACTCGCAACTAGAATAACTCCCGCTATC‑TCATT‑ACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTT‑CCCCTATTTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCTTTACACAC‑tt > 1:415792/5‑209 (MQ=255) taatTCAAGGTT‑‑‑AAAATCAATTAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGtt < 1:1263969/244‑1 (MQ=255) taaaTCAA‑TAACTTATTC‑TAAGTATTTGACAGCACTGAATGTC‑AAAAC‑‑AAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGAC‑TG‑ACTTCGATT‑‑G‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGtt < 1:1263970/278‑1 (MQ=255) tCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑TTT‑TTTAT‑CCCTTCTTTACACAC‑TTTTCATTATTCt > 1:415795/1‑197 (MQ=255) ttAAGTA‑TTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑‑GA‑‑TG‑CCTCGCATTCGG‑ttt > 1:1263972/1‑148 (MQ=255) aTGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATaa > 1:1263974/1‑225 (MQ=255) aaaaacaaaaaCCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCtgtg > 1:415798/2‑162 (MQ=255) aacCTTCACTCGCAACTAGAAT‑ACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGAGCTGACTTCGGATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCt < 1:415796/156‑1 (MQ=255) gCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑‑T‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGaaaactat < 1:1263973/158‑5 (MQ=255) tAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑TTT‑TTTTATCCCTTCTTTACACAC‑TTTTCa > 1:415800/1‑123 (MQ=255) tAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑AT‑ACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTG‑CAATTATTTTTGTTGGGCTTTTTG‑TAGGTTATTTGTACAGCAAAATGGCGCTTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACgg > 1:1263976/1‑310 (MQ=255) ttatttatTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACTAATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTG‑CAATTATTTTTGTTGGGCTTTTTG‑TAGGTTATTTGTACAGCAAAATGGCGC‑TTGTACATCTATTTCCCCCAATGCAGgatga > 1:1263978/1‑275 (MQ=255) attaCCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCG‑‑‑TT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAG‑AAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCa < 1:1263977/166‑1 (MQ=255) tcattT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACCGGGA‑‑TG‑CCTCGCATTCGG‑TTTGTTTTCACCCTTCTTTACACAC‑TTTTCA‑TATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCttgt < 1:415802/169‑3 (MQ=255) tGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑TTT‑TTTTATCCCTTCTTTACACAC‑TTTTCa > 1:415803/1‑67 (MQ=255) ttACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTTCG‑GTT‑TTTTTACCCTTCTTTACACAC‑ttt > 1:415804/1‑61 (MQ=255) tCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑‑TT‑TTTTTACCCTTCtttt > 1:415805/1‑40 (MQ=255) aaaCC‑GGA‑‑TG‑CCTCGCATTCGGTTTT‑TTTTTACCCTTCTTTACACACTTTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTGTCAATTATTTTTGTTGGGCTTTTTGTTAGGTTAtttt > 1:1263980/2‑184 (MQ=255) | AAATGCAAGGGAACTTTTTAAGATTATTGCGGAATGGCGAAAT‑AAGCACCTAACATCAAGCAATAATAATTCAAGGTT‑‑‑AAAATCAA‑TAACTTATTCTTAAGTATTTGACAGCACTGAATGTC‑AAAAC‑AAAACCTTCACTCGCAACTAGAATAACTCCCGCTATCATCATTAACTTTATTTATTACCGTCATTCATTT‑CTGAATGTCTGTTTACCCCTA‑TTTCAACC‑GGA‑‑TG‑CCTCGCATTCGG‑TTT‑TTTTTACCCTTCTTTACACAC‑TTTTCATTATTCTGTGCTACCACAGAAAAACT‑ATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAATCTTTGTTTATTTG‑CAATTATTTTTGTTGGGCTTTTTG‑TAGGTTATTTGTACAGCAAAATGGCGC‑TTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGG > NC_000913/2312926‑2313386 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |