| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,272 | C→T | P806L (CCG→CTG) | rpoB → | DNA‑directed RNA polymerase subunit beta |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,163,272 | 0 | C | T | 97.6% | 150.4 / ‑3.2 | 41 | P806L (CCG→CTG) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (18/22); total (19/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.63e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01 | |||||||||||
GCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCAC > REL606/4163201‑4163344 | gCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTgg < 1:459697/76‑1 (MQ=255) cTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGa > 2:932938/1‑76 (MQ=255) tGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGaa < 2:624538/76‑1 (MQ=255) cAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATgg > 2:767271/1‑76 (MQ=255) cGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTAc > 1:108927/1‑76 (MQ=255) ggTCCGTCCACCGACCTCGGTGAACTGGCGCTTNGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACa > 2:1364741/1‑76 (MQ=255) ggTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACa < 1:1004606/76‑1 (MQ=255) gtccgtccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACaa < 2:572809/76‑1 (MQ=255) gtccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTc < 2:634603/76‑1 (MQ=255) tccACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCg < 2:89851/76‑1 (MQ=255) cACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGaa > 1:1285608/1‑76 (MQ=255) gACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACt < 1:892909/76‑1 (MQ=255) cTCGNTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGNTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCa < 1:197111/76‑1 (MQ=255) cTCGGTGAACTGGCGCTTGGTCAGAACATGCGCNTAGCGTTCATGCTGTGGAATGGTTACAACTTCCAAGACTCCa > 2:1375897/1‑76 (MQ=255) cGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATc > 1:1377893/1‑76 (MQ=255) gAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCg > 1:1091813/1‑76 (MQ=255) aaCTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGt > 2:1488271/1‑76 (MQ=255) cTGNCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGANTGGTTACAACTTCGAAGACTCCATCCTCGTAt < 2:1512355/76‑1 (MQ=255) tGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATc > 1:1157113/1‑76 (MQ=255) gcgcTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCg < 2:870850/76‑1 (MQ=255) gcgcTTGGTCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCg < 2:509876/76‑1 (MQ=255) ttGGTCAGAACATGCGCGTAGCGTTCATGCTGTTGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCg > 1:1320909/1‑76 (MQ=255) tCAGAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGtgtt > 1:76249/1‑76 (MQ=255) aGAACATGCGCGTAGCGCTCATGCTGTGGAATGGTTACAACTTCGAAGACGCCATCCTCGTATCAGAGCGtgtngt > 2:1529568/1‑76 (MQ=255) gAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCGATCCTCGTATCCGAGCGTgttgtt < 2:1172122/76‑1 (MQ=255) gAACATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTgttgtt < 1:966892/76‑1 (MQ=255) aaCATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTc < 1:825045/76‑1 (MQ=255) aCATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCa > 2:128445/1‑76 (MQ=255) cATGCGCGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAg < 2:488263/76‑1 (MQ=255) gcgcGTAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGaa > 2:674113/1‑76 (MQ=255) cgTAGCGTTCATGCTGTGCAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGAc < 2:823848/76‑1 (MQ=255) tAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCg < 2:309585/76‑1 (MQ=255) tAGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCg < 1:548122/76‑1 (MQ=255) aGCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGt > 2:183859/1‑76 (MQ=255) aGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGt > 2:320340/1‑76 (MQ=255) gCGTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGtt < 2:401363/76‑1 (MQ=255) gTTCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTc < 1:1264743/76‑1 (MQ=255) tCATGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcac < 1:1179962/76‑1 (MQ=255) aTGCTGTGGAATGGTTACANCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcacca > 1:1516126/1‑76 (MQ=255) tGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcaccac < 2:1091905/76‑1 (MQ=255) tGCTGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTcaccac > 1:631220/1‑76 (MQ=255) | GCTGGCAGACGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCAC > REL606/4163201‑4163344 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |