| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 4,163,698 | T→G | I948S (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,698 | 0 | T | G | 98.2% | 218.1 / ‑7.3 | 54 | I948S (ATC→AGC) | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (33/20); total (33/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTT > REL606/4163625‑4163770 | gtATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCg < 1:62693/76‑1 (MQ=255)gtATCCGGTACGGTTATCGACGATCNGCTCTTTACTAGCGATGGCGTNGAAGAAGACAAACGTGCGCTGGAAAGCg < 1:258561/76‑1 (MQ=255) tATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCga > 2:1359797/1‑76 (MQ=255) tATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCga > 2:1640469/1‑76 (MQ=255) tATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCga > 2:791470/1‑76 (MQ=255) ccGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAANAGACAAACGTGCGCTGGAAAGCgaaga < 2:231/76‑1 (MQ=255) gTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAAt > 2:984242/1‑76 (MQ=255) aCGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGc > 2:440743/1‑76 (MQ=255) cGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCa > 2:1007830/1‑76 (MQ=255) ggTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAg > 1:888279/1‑76 (MQ=255) ttATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCt < 1:907600/76‑1 (MQ=255) tATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTc > 1:1209111/1‑76 (MQ=255) aTCGACGTTCAGGTCTTTACTCGCGATCGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCa > 1:1564537/1‑76 (MQ=255) tCGACGTTCAGGTCTTTACTCGCGATGGCGTAGCAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCaa > 1:1619677/1‑76 (MQ=255) aCGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGACGAAATGCAGCTCAAACa < 1:1758077/76‑1 (MQ=255) aCGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACa > 1:240479/1‑76 (MQ=255) aCGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACa < 1:1290176/76‑1 (MQ=255) ttCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGc < 1:802086/76‑1 (MQ=255) tCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCg > 1:825673/1‑76 (MQ=255) cAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCgatga > 2:143016/1‑56 (MQ=255) cAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCgaaga < 1:143016/56‑1 (MQ=255) cAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCga > 2:1367635/1‑76 (MQ=255) cAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCga > 2:480928/1‑76 (MQ=255) ggNCTTTAATCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCgnag > 2:1442299/1‑76 (MQ=255) gTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCgaaga < 2:1483960/76‑1 (MQ=255) cTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCTCTCGAAATCGANGAAATGCAGCTCAAACAGGCGaagaaa < 1:1543892/76‑1 (MQ=255) cTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGCAAAGCGAAGAAATGTAGCTCAAACAGGCGaagaaa > 2:1403536/1‑76 (MQ=255) ttACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGaagaaaga > 1:114927/1‑76 (MQ=255) ttACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGaagaaaga > 2:1543697/1‑76 (MQ=255) aCTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGAcc > 2:1599401/1‑76 (MQ=255) tCGCGATGGCGTAGAAAAAGACAACCGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTg < 2:654944/76‑1 (MQ=255) tCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAATAAATGCAGCTCAAACAGGCGAAGAAAGACCTg < 2:1190164/76‑1 (MQ=255) gcgATGGCGTAGAAAAAGTCAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTc < 1:534339/76‑1 (MQ=255) gATGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTg < 2:100427/76‑1 (MQ=255) tGGCGTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTgaa < 2:1256401/76‑1 (MQ=255) gTAGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAAc > 1:446338/1‑76 (MQ=255) aGAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTg > 2:204138/1‑76 (MQ=255) gAAAAAGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACTTGTCTGAAGAACTGc < 2:154345/76‑1 (MQ=255) aaaaGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAg > 2:1315415/1‑76 (MQ=255) aaaGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGa > 1:1027453/1‑76 (MQ=255) aaaGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGa > 1:1156637/1‑76 (MQ=255) aaaGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGa > 2:704646/1‑76 (MQ=255) aaGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGAt < 2:1533226/76‑1 (MQ=255) aGACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATc < 1:1407182/76‑1 (MQ=255) gACAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATcc > 1:1244627/1‑76 (MQ=255) cAAACGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTc > 2:412238/1‑76 (MQ=255) aaaCGTGCGCTGGACAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGCAGAACTTCAGATCCTCg < 1:140184/76‑1 (MQ=255) aaaCGTGCGCTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCg > 2:544844/1‑76 (MQ=255) aCGTGCGCTTGAAAGCGAAGNAATGCAGCTAAAACANGCGAAGAAAGACCCGACTGAAGAACTGCAGATCCTCGaa < 2:1755808/76‑1 (MQ=255) tGCGCTGGAAAGCGAAGAAAAGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCg > 1:937527/1‑76 (MQ=255) cgcTGGAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAAATGCAGATCCTCGAAGCggg > 2:701302/1‑76 (MQ=255) tGGAAAGCGAAGAAATGCAGCTCAAACAGGTGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCt < 2:1543858/76‑1 (MQ=255) ggAAAGCGTAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTg > 2:56799/1‑76 (MQ=255) ggAAAGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTg > 2:472656/1‑76 (MQ=255) aaaGCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGtt > 2:867821/1‑76 (MQ=255) | GTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTT > REL606/4163625‑4163770 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |