Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,261,359 | T→A | 15.4% | intergenic (+53/‑310) | yjbM → / → dusA | uncharacterized protein/tRNA‑dihydrouridine synthase A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,261,359 | 0 | T | A | 15.4% | 64.2 / 6.8 | 26 | intergenic (+53/‑310) | yjbM/dusA | uncharacterized protein/tRNA‑dihydrouridine synthase A |
Reads supporting (aligned to +/- strand): ref base T (10/12); new base A (2/2); total (12/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01 |
ATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACAT > NC_000913/4261175‑4261627 | aTGCGAAATGCTTAACCAAAATTAATTAGNTTCGCGATGTCGTTTGACCCAACCTNAAATGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGNGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg < 3:90291/275‑1 (MQ=255) aTGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa > 1:90291/1‑264 (MQ=255) ctgACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTCGGTTCTTTTTTGTCTATCAATACTATTAGCAGGgaga < 7:958854/146‑1 (MQ=255) ctcTTCTAAATTAATTTATGATTCGTTCTGGGTTAAATCGGAGAAATAATATGTTAAGTTCTTCGATCGTATGGAAGGGATTATATAAAAGGACCCAATATTTAGTGGGTTATTTTTTCTCTATCAATACTATTAGCAGGGAGataag > 8:958854/3‑146 (MQ=255) cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta > 8:823002/1‑102 (MQ=255) cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta < 7:823002/102‑1 (MQ=255) tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat > 8:736396/1‑84 (MQ=255) tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat < 7:736396/84‑1 (MQ=255) tAAATCGGAGAAATAATATGTGCAGTTCTTCGATAGTATGGAAGGCAGTATAAATAAGTGCCCATGATTTATTAGGTTCTTATATCTCTCTCAATAGTATTAGCGGGGTGGTATATCAGGAGAGATTAAGCTGTGATTTtttattttt < 8:106051/148‑3 (MQ=255) gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc < 7:587407/97‑1 (MQ=255) gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc > 8:587407/1‑97 (MQ=255) ataaaaTGTGAAGTTTTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTa < 2:491665/274‑1 (MQ=255) tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCCCCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTCATAAGATTATTTCCATCACTTCGTTAGAGCTTCATGCCTTGAGAGGATCTCATTTTTATTTTGCAATGAGTCAGGCGCTTCCTGTTGTTATGGTATAGTTCTCCGATATTGACCCTCCTGCATtgtg > 4:486487/1‑275 (MQ=255) tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa < 7:1114376/81‑1 (MQ=255) tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa > 8:1114376/1‑81 (MQ=255) ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt < 7:304949/91‑1 (MQ=255) ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt > 8:304949/1‑91 (MQ=255) aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt > 7:451406/1‑111 (MQ=255) aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt < 8:451406/111‑1 (MQ=255) aaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg < 8:830833/91‑1 (MQ=255) aaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg > 7:830833/2‑92 (MQ=255) aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc < 7:147092/58‑1 (MQ=255) aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc > 8:147092/1‑58 (MQ=255) cccAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt < 6:1063305/79‑1 (MQ=255) cccAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt > 5:1063305/1‑79 (MQ=255) ccatatTTATTGGGTTCTTTTTTCTCAATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTAAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat < 2:486488/273‑1 (MQ=255) cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc < 5:713886/87‑1 (MQ=255) cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc > 6:713886/1‑87 (MQ=255) cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTa < 5:381507/118‑1 (MQ=255) cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTa > 6:381507/1‑118 (MQ=255) | ATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACAT > NC_000913/4261175‑4261627 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |