| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 1,542,428 | G→T | 11.5% | L181L (CTG→CTT) | ylaA → | conserved protein of unknown function |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 1,542,428 | 0 | G | T | 11.5% | 93.0 / 5.2 | 35 | L181L (CTG→CTT) | ylaA | conserved protein of unknown function |
| Reads supporting (aligned to +/- strand): ref base G (12/19); new base T (2/2); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01 | |||||||||||
GTACAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTGCAGTAATGCCATATACGAT > BS168_reference_genome/1542300‑1542556 | gTATAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGtaga < 4:162453/138‑1 (MQ=255)gTACAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGaa < 6:42656/139‑1 (MQ=255)gTACAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGaa < 5:304369/139‑1 (MQ=255)gTACAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGACTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGaa < 7:71450/139‑1 (MQ=255) aaTGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTTTTGGAGAACGGAt < 8:274297/139‑1 (MQ=255) tttttAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGAt < 5:47854/139‑1 (MQ=255) gTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTAc < 8:211167/138‑1 (MQ=255) tttCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGAtggtac > 5:163124/1‑132 (MQ=255) acagggaaTACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATAc < 5:239583/135‑1 (MQ=255) aaaggaaTACACTGATGGGTTGGAGCCGTTTGGATTAGTGTTTCTGTTTTTCCTCATCACTCTTCCCATTGCGTTGAATGTAGTGGTTTTAAAGTTGTATTGGAGAACGGATAAGATCTTACTATGATTTTTTTATACg < 2:205245/138‑1 (MQ=255) ggatggAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGc < 5:104246/138‑1 (MQ=255) gCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTtg < 6:228780/139‑1 (MQ=255) agtGCATTGGTGTTTCTGGTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGGGTTTTTGAAGCGGTTTTGGAGAACGGATAAGATCTTACAAGGCTTATTTTATACGCGGCTTATCGGCAGGCTTGTCTtgctgc > 6:131773/3‑139 (MQ=255) aTTGCATTCGTGTTTCTGTTTTTCCTCATCACTCGTCCCATTGCGGTGAATGTAGTGTTTTTGAAGCTGTATTGGAGAACGGATAAGATCATACTAGGCTTATTTTATACGCTGATTATCGGCCGGCTTGTCTttctgc > 8:147436/1‑139 (MQ=255) aTTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTGAAGCTGTATTTGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTtgctgc > 1:113427/1‑139 (MQ=255) aTTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTGAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATAGGACAGATTGTCTtgctgc > 5:55745/1‑139 (MQ=255) aTTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGAt > 8:245690/1‑84 (MQ=255) aTTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGAt < 7:245690/84‑1 (MQ=255) tCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTTTATTGGAGAACGGAt > 3:174323/1‑70 (MQ=255) tCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTTTATTGGAGAACGGAt < 4:174323/70‑1 (MQ=255) tCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTTTATTGGAGAACGGAt > 3:170220/1‑70 (MQ=255) tCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTTTATTGGAGAACGGAt < 4:170220/70‑1 (MQ=255) ttCATCATCACTCTTCCCATTGCGTTGATTGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCaa > 5:248608/1‑138 (MQ=255) atcaCTCTTCCCATTGCGTTGAATGTAGTGGTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCa < 4:89069/138‑1 (MQ=255) cccATTGCGTTGAATGTAGTGTTTTTCTAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTa < 7:3428/137‑1 (MQ=255) ccATTGCGTTGAATGTAGTGGTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATgg < 4:252443/139‑1 (MQ=255) gCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTa > 3:76665/1‑131 (MQ=255) gCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTa > 5:49771/1‑131 (MQ=255) gCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTa < 4:76665/131‑1 (MQ=255) gCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTa < 6:49771/131‑1 (MQ=255) tAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTGCAGTAATGCCata < 7:147436/139‑1 (MQ=255) ttttCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGCCTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTGCAGTAATGCCATATACGag > 2:18272/1‑138 (MQ=255) ttttCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTGCAGTAATGCCATATACGat > 2:308459/1‑139 (MQ=255) cAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTg < 8:215864/116‑1 (MQ=255) cAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTg > 7:215864/1‑116 (MQ=255) | GTACAAATGATTTTTAGTTTCTTTTCGTATTGGAAGGAATACCCGGATGGGATGGAGCCGATTGCATTCGTGTTTCTGTTTTTCATCATCACTCTTCCCATTGCGTTGAATGTAGTGTTTTTCAAGCTGTATTGGAGAACGGATAAGATCTTACTAGGCTTATTTTATACGCTGATTATCGGACTGCTTGTCTTGCTGCTTGCGCCATATCCCAATCAAATCTCAAAAGATTATGGTGCAGTAATGCCATATACGAT > BS168_reference_genome/1542300‑1542556 |
| 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: ‑ |