Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 542,820 | C→T | 37.8% | A151T (GCC→ACC) | ylbA ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 542,820 | 0 | C | T | 37.8% | 28.0 / 35.8 | 37 | A151T (GCC→ACC) | ylbA | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (10/13); new base T (6/8); total (16/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 = 8.33e-01 |
TTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTG > W3110S.gb/542752‑542889 | tttGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGc < 1:1435223/69‑1 (MQ=255) aaaTCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAc > 1:903125/1‑71 (MQ=255) aaaTCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAc < 1:1843399/71‑1 (MQ=255) tCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAg < 1:2660452/71‑1 (MQ=255) agcagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa < 1:2439829/71‑1 (MQ=255) gcagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCAc < 1:3649338/71‑1 (MQ=255) agAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACgg > 1:1203271/1‑71 (MQ=255) gAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGc < 1:3663510/49‑1 (MQ=255) gAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACgg > 1:270445/1‑70 (MQ=255) tAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCa > 1:644559/1‑47 (MQ=255) tAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAAcc < 1:302299/59‑1 (MQ=255) tCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCAc < 1:1394654/60‑1 (MQ=255) tCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACg < 1:4616714/61‑1 (MQ=255) cGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGa < 1:553855/49‑1 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCtt > 1:4127261/1‑71 (MQ=255) cATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCt < 1:1405264/70‑1 (MQ=255) cATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa > 1:4229663/1‑54 (MQ=255) gCCTTCATAATGAATGCGTTCCAGTTCGCTGGTAt > 1:2057090/1‑35 (MQ=255) cTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGa < 1:4272108/71‑1 (MQ=255) ttCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGa < 1:1639118/39‑1 (MQ=255) ttCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGt < 1:3022160/52‑1 (MQ=255) tCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAAc < 1:1714735/41‑1 (MQ=255) tCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAAc < 1:4206777/41‑1 (MQ=255) tCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCa < 1:4364957/47‑1 (MQ=255) tCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa > 1:547830/1‑64 (MQ=255) cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGaaa > 1:762415/1‑39 (MQ=255) cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGaaa > 1:1879565/1‑39 (MQ=255) cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa > 1:627842/1‑63 (MQ=255) aTAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACAt < 1:781388/71‑1 (MQ=255) aTAATGAATGCGTTCCAGTTCGCTGGCATTGCCATAAACCAGCCACGGTGCATAGCCTTCTACCGGAACa < 1:558517/70‑1 (MQ=255) aTAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACa > 1:4339035/1‑70 (MQ=255) atgaatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa > 1:2517252/1‑59 (MQ=255) tgaatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATa > 1:937225/1‑50 (MQ=255) tatgCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTAc < 1:984020/56‑1 (MQ=255) tgCGTTCCAGTTCGCTGGCATTGCCATAAACCAGCCACGTTGCATAGCCTTCTACCGGAACATAGCGgcgc < 1:674065/71‑1 (MQ=255) cGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGc > 1:3064999/1‑36 (MQ=255) cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATaaa > 1:4056693/1‑68 (MQ=255) tGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACAt < 1:2196164/48‑1 (MQ=255) gCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTg < 1:2204043/71‑1 (MQ=255) | TTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTG > W3110S.gb/542752‑542889 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |