| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,752 | T→A | I966N (ATC→AAC) | 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,752 | 0 | T | A | 96.4% | 93.0 / ‑4.4 | 28 | I966N (ATC→AAC) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base A (19/8); total (19/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.21e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01 | |||||||||||
AAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGA > REL606/4163677‑4163824 | aaGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAt > 2:922119/1‑76 (MQ=255) gcTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGt > 1:322611/1‑76 (MQ=255) aagaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCg < 1:402683/76‑1 (MQ=255) agaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGt > 1:198313/1‑76 (MQ=255) aaaTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTAt < 2:966453/76‑1 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATcc > 2:901964/1‑76 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATcc > 1:529029/1‑76 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATcc < 2:688433/76‑1 (MQ=255) aGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgc > 1:216890/1‑76 (MQ=255) cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCGCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg < 2:958939/76‑1 (MQ=255) cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg > 2:762934/1‑76 (MQ=255) cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg > 2:476943/1‑76 (MQ=255) aaaCAGGCGAAGAAAGACCTGTCTGAAGAACTACAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgc > 2:920195/1‑76 (MQ=255) cAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTgg < 1:257230/76‑1 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt > 2:483414/1‑76 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt > 2:239001/1‑76 (MQ=255) ggCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTa < 2:17304/76‑1 (MQ=255) gCGAAGAAAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg < 2:1041513/76‑1 (MQ=255) gaagaaAGACCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGcc > 1:272695/1‑76 (MQ=255) aagaaagaCCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCg < 1:501355/76‑1 (MQ=255) aaagaCCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg > 2:583832/1‑76 (MQ=255) aaagaCCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTg > 2:558025/1‑76 (MQ=255) agaCCTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGc < 2:321203/76‑1 (MQ=255) cTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg > 1:555373/1‑76 (MQ=255) cTGTCTGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTg > 1:902232/1‑76 (MQ=255) gtgtGAAGAACTGGAGACCCTCGAGGCGGGTCTGTCCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa < 1:920397/73‑1 (MQ=255) tGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCt > 2:380342/1‑76 (MQ=255) tGAAGAACTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGGGGCGTTGAAGCt > 2:88955/1‑76 (MQ=255) gaagaaCTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg > 2:550468/1‑76 (MQ=255) aagaaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTga < 1:918875/76‑1 (MQ=255) agacCTGCAGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgag < 1:726470/76‑1 (MQ=255) aGACCCTGGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGa < 1:840271/76‑1 (MQ=255) aGAACCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGGGGCGTTGTAGCTGAGCAGCTCGa > 2:331532/1‑76 (MQ=255) | AAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGA > REL606/4163677‑4163824 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |