Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 3,935,823 | Δ1 bp | 100% | intergenic (+35/+1) | ywbE → / ← rlmI | conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 3,935,823 | 0 | T | . | 92.6% | 91.3 / 2.1 | 27 | intergenic (+35/+1) | ywbE/rlmI | conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base . (15/10); minor base A (1/1); total (16/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC > BS168_reference_genome/3935695‑3935952 | tatTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCt > 1:59207/1‑139 (MQ=255) tGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCGACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCGTTTTGA‑TTAttttttggtct > 1:112454/1‑132 (MQ=255) cGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAg > 6:23396/1‑138 (MQ=255) aaaTCAAACTTTCATCCGCTGGCCTTTAATGGCCGTCTTGAAGTGGCCCGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc < 7:185163/139‑1 (MQ=255) cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg > 1:127309/1‑138 (MQ=255) cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg > 5:61348/1‑138 (MQ=255) cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg > 5:25538/1‑138 (MQ=255) cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg > 5:209318/1‑138 (MQ=255) tCTTGAAGATGGCGGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa < 7:31490/139‑1 (MQ=255) tCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCGGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaa > 8:81040/1‑138 (MQ=255) ggCCCCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa < 7:88931/138‑1 (MQ=255) ggcttCGCCCGCTACAAGCAGATCTTTTGTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTT‑AATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGaaa < 4:146683/133‑1 (MQ=255) cATCGGCCGAGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTT‑AATTATTTTTTATGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCCTTTTGAAAtcct > 8:116015/1‑139 (MQ=255) gACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGt < 5:18523/139‑1 (MQ=255) agcagATCGTTTCTACATAAAAAAAAGCTGCCCCCGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt < 7:81040/139‑1 (MQ=255) cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc > 5:110693/1‑139 (MQ=255) gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa < 6:61348/138‑1 (MQ=255) cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg < 6:125522/139‑1 (MQ=255) ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCTTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc > 8:220228/1‑139 (MQ=255) ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc > 1:100522/1‑139 (MQ=255) ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc > 2:171891/1‑139 (MQ=255) ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc > 2:68944/1‑139 (MQ=255) ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAAAGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc > 8:114641/1‑139 (MQ=255) aGCTTTTTCA‑TTATTTTTTCTGCAGCAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc < 7:131241/139‑1 (MQ=255) aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc < 5:23396/139‑1 (MQ=255) cTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc > 7:148130/1‑84 (MQ=255) cTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc < 8:148130/84‑1 (MQ=255) | TATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC > BS168_reference_genome/3935695‑3935952 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |