| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,728 | A→G | K958R (AAA→AGA) | 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,728 | 0 | A | G | 97.3% | 142.9 / ‑4.3 | 37 | K958R (AAA→AGA) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (20/16); total (21/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
TTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGG > REL606/4163656‑4163800 | ttACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGaagagaga > 2:1330925/1‑76 (MQ=255)ttACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGaagagaga > 1:524684/1‑76 (MQ=255) aCTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGAcc > 1:140393/1‑76 (MQ=255) cTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCt > 1:851609/1‑76 (MQ=255) cgcgATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGt > 1:1334781/1‑76 (MQ=255) cgATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCt > 2:2059839/1‑76 (MQ=255) gATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTg < 2:777214/76‑1 (MQ=255) ggCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTgaag > 1:361507/1‑76 (MQ=255) cGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTgaagaa < 2:508855/76‑1 (MQ=255) tAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACt < 1:1636940/76‑1 (MQ=255) aGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTg < 2:146330/76‑1 (MQ=255) aGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTg > 1:2016163/1‑76 (MQ=255) aaaaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCa > 1:261601/1‑76 (MQ=255) aaaaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAg < 1:707046/76‑1 (MQ=255) aaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGAt < 2:1172208/76‑1 (MQ=255) aCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGaa < 2:1250765/76‑1 (MQ=255) cGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAg > 2:509234/1‑76 (MQ=255) tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCt > 1:814290/1‑76 (MQ=255) aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTc > 2:787447/1‑76 (MQ=255) tCGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAg < 2:775902/76‑1 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGc > 1:1150702/1‑76 (MQ=255) gaagaaATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGcc < 2:1673271/76‑1 (MQ=255) gaagaaATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGcc > 2:1474037/1‑76 (MQ=255) gaagaaATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGcc > 1:1963351/1‑76 (MQ=255) agaaATGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGt > 1:1259587/1‑76 (MQ=255) aTGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATcc > 1:1639071/1‑76 (MQ=255) tGCAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCg > 1:659090/1‑76 (MQ=255) cAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtg > 1:195392/1‑76 (MQ=255) cAGCTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtg < 2:341320/76‑1 (MQ=255) cTCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg < 1:597855/76‑1 (MQ=255) tCAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgt < 1:1366485/76‑1 (MQ=255) cAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg > 2:1017046/1‑76 (MQ=255) cAAACAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg < 2:1854855/76‑1 (MQ=255) cAAACAGGCAAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg < 1:983226/76‑1 (MQ=255) cAGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg > 1:92981/1‑76 (MQ=255) aGGCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt > 1:1587221/1‑76 (MQ=255) ggCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGCGCTGGTa > 2:900526/1‑76 (MQ=255) gCGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg < 2:1009897/76‑1 (MQ=255) cGAAGAGAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGc < 1:1304809/76‑1 (MQ=255) gaagagAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc < 2:1715687/76‑1 (MQ=255) aagagagaCCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg > 1:566179/1‑76 (MQ=255) agagagaCCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCgg > 2:232383/1‑76 (MQ=255) | TTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGG > REL606/4163656‑4163800 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |