Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,162,470 | A→C | T539P (ACC→CCC) | 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,162,470 | 0 | A | C | 96.9% | 257.9 / ‑1.8 | 65 | T539P (ACC→CCC) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (29/34); total (30/35) | |||||||||||
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.32e-01 |
TATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA > REL606/4162397‑4162542 | tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc < 2:1065663/76‑1 (MQ=255) tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc > 1:1069451/1‑76 (MQ=255) tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc > 1:1363552/1‑76 (MQ=255) tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg < 2:1131630/76‑1 (MQ=255) tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg < 1:1791782/76‑1 (MQ=255) aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGa > 2:473841/1‑76 (MQ=255) ccAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGaa > 1:458615/1‑76 (MQ=255) gAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGt < 2:1160405/76‑1 (MQ=255) aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTg > 2:339497/1‑76 (MQ=255) acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGc < 1:77014/76‑1 (MQ=255) caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCa > 2:1246604/1‑76 (MQ=255) aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCan < 1:472763/76‑2 (MQ=255) aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg > 2:1919884/1‑76 (MQ=255) cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGc < 1:956506/76‑1 (MQ=255) ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCt > 1:113455/1‑76 (MQ=255) cGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGAACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCtt < 1:770955/76‑1 (MQ=255) tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa > 2:185116/1‑76 (MQ=255) tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa > 2:136414/1‑76 (MQ=255) gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa < 2:1589251/76‑1 (MQ=255) tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAg < 1:639839/76‑1 (MQ=255) cTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGCGCAGGCTTCGAAGt > 1:559577/1‑76 (MQ=255) tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGtt < 2:206602/76‑1 (MQ=255) tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGtt > 1:1021588/1‑76 (MQ=255) gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTc > 1:1393243/1‑76 (MQ=255) agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg < 2:815298/76‑1 (MQ=255) agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg < 2:1945026/76‑1 (MQ=255) aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag < 2:997069/76‑1 (MQ=255) aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag < 1:1189235/76‑1 (MQ=255) aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag < 1:1448253/76‑1 (MQ=255) ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga < 2:1759479/76‑1 (MQ=255) ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga > 1:1041907/1‑76 (MQ=255) ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga < 1:488063/76‑1 (MQ=255) ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGGTCgaga > 1:601761/1‑76 (MQ=255) tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc > 1:1164345/1‑76 (MQ=255) aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg < 2:247888/76‑1 (MQ=255) aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg > 1:352645/1‑76 (MQ=255) cGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGt < 1:502227/76‑1 (MQ=255) gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa < 2:1072322/76‑1 (MQ=255) cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac > 2:487282/1‑76 (MQ=255) caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac < 1:330853/76‑1 (MQ=255) caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac < 1:474624/76‑1 (MQ=255) aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc > 1:862877/1‑76 (MQ=255) gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa < 1:1066264/76‑1 (MQ=255) gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa < 2:311399/76‑1 (MQ=255) gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa > 2:1674918/1‑76 (MQ=255) tATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc > 1:1457528/1‑76 (MQ=255) aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt > 2:1965647/1‑76 (MQ=255) aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt < 2:514945/76‑1 (MQ=255) aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt > 1:1514872/1‑76 (MQ=255) aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt < 2:1457528/76‑1 (MQ=255) aTCTCCGCACTCGGCCCAGGCGGTCTGCCACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt > 2:764013/1‑76 (MQ=255) tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg > 2:213096/1‑76 (MQ=255) tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg > 2:74884/1‑76 (MQ=255) cgcgCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGt < 1:642087/76‑1 (MQ=255) cACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg < 2:1867104/76‑1 (MQ=255) aCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc > 2:993607/1‑76 (MQ=255) aCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc > 1:1548414/1‑76 (MQ=255) cGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTa < 2:1598799/76‑1 (MQ=255) ggCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAt < 2:842642/76‑1 (MQ=255) gCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg > 2:1266268/1‑76 (MQ=255) gCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg < 1:137097/76‑1 (MQ=255) gCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt < 1:551511/76‑1 (MQ=255) ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCga < 1:970943/36‑1 (MQ=255) ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCga > 2:970943/1‑36 (MQ=255) ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg < 2:779662/76‑1 (MQ=255) ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg > 1:26704/1‑76 (MQ=255) gTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGa > 1:1604476/1‑76 (MQ=255) tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa > 1:536939/1‑76 (MQ=255) tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa < 2:1671039/76‑1 (MQ=255) cTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa < 1:878548/76‑1 (MQ=255) | TATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA > REL606/4162397‑4162542 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |