| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,752 | T→G | I966S (ATC→AGC) | 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 | G | 97.9% | 181.1 / ‑4.1 | 47 | I966S (ATC→AGC) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (28/18); total (29/18) | |||||||||||
| 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.43e-01 | |||||||||||
GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG > REL606/4163679‑4163823 | gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGGCCTGTCTGAAGAACTGCAGAGcc > 1:891709/1‑76 (MQ=255)gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGcc > 1:1109067/1‑76 (MQ=255)gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGcc > 1:1013220/1‑76 (MQ=255)gACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGcc > 1:792406/1‑76 (MQ=255) aaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGa < 2:620603/76‑1 (MQ=255) aaCGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGa < 2:1531831/76‑1 (MQ=255) gTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGc > 2:925763/1‑76 (MQ=255) tGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCg > 1:1098325/1‑76 (MQ=255) tGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCt < 2:237529/76‑1 (MQ=255) aaTCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTc < 2:140931/76‑1 (MQ=255) tCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAg > 2:1095089/1‑76 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc > 2:415346/1‑76 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc > 2:1156091/1‑76 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc > 2:1568682/1‑76 (MQ=255) cGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGc > 2:729092/1‑76 (MQ=255) agaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGt < 2:885111/76‑1 (MQ=255) agaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGt > 2:1191983/1‑76 (MQ=255) agaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGt > 2:1343876/1‑76 (MQ=255) gaaATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTa < 2:614604/76‑1 (MQ=255) aTGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATcc < 2:896583/76‑1 (MQ=255) tGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCg < 2:974192/76‑1 (MQ=255) cAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGAATCCgtg < 2:1113294/76‑1 (MQ=255) cAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGGCTGTTCAGCCGTATCCgtg > 2:821268/1‑76 (MQ=255) cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg < 2:1477801/76‑1 (MQ=255) cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg < 2:643144/76‑1 (MQ=255) cTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctg < 2:448831/76‑1 (MQ=255) cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg > 2:1244667/1‑76 (MQ=255) cAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtg < 1:53154/76‑1 (MQ=255) aaCAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCgtgctgtgct > 1:1274322/1‑76 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt > 1:1491410/1‑76 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt < 2:194287/76‑1 (MQ=255) aGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt < 2:299620/76‑1 (MQ=255) ncGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg > 1:1229175/2‑76 (MQ=255) gCGAAGAAAGACCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAg > 2:1402405/1‑76 (MQ=255) gcccTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCg < 2:624317/74‑1 (MQ=255) aCCTGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGt > 1:1529868/1‑76 (MQ=255) tGTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGa > 1:1611472/1‑76 (MQ=255) tGTCTGAAGAACGGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGa < 1:778077/76‑1 (MQ=255) tccTGAAGACCTCCAGCGCCTCGAGGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa < 2:211597/74‑1 (MQ=255) gTCTGAAGAACTGCAGAGCTGCGACGGGGGCCTGTTCGGCCGTCTCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa < 2:863991/76‑1 (MQ=255) gTCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGaa < 2:1492949/76‑1 (MQ=255) tCTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAg < 2:1465336/76‑1 (MQ=255) cTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGc > 1:307257/1‑76 (MQ=255) cTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGc > 1:628978/1‑76 (MQ=255) cTGAAGAACTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGc > 2:1253044/1‑76 (MQ=255) gaagaaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg > 1:549638/1‑76 (MQ=255) gaagaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTg > 2:450439/1‑76 (MQ=255) agaaCTGCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTgag > 2:461199/1‑76 (MQ=255) aCTGCGGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAg < 1:1386518/76‑1 (MQ=255) gCAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTc > 2:1591410/1‑76 (MQ=255) cAGAGCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCg > 2:44378/1‑76 (MQ=255) | GACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCG > REL606/4163679‑4163823 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |