Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP009273 | 4,173,770 | T→G | D866E (GAT→GAG) | rpoB → | RNA polymerase, beta subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,173,770 | 0 | T | G | 97.8% | 144.1 / ‑4.1 | 45 | D866E (GAT→GAG) | rpoB | RNA polymerase, beta subunit |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (18/26); total (19/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.22e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01 |
TGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT > CP009273/4173706‑4173832 | tGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGaa > 1:3773983/1‑68 (MQ=255) gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt < 1:103901/68‑1 (MQ=255) gggTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAt < 1:2769837/68‑1 (MQ=255) ccGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCgg > 1:1117613/1‑68 (MQ=255) ccGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTCGAGGAATCCgg < 1:3264586/68‑1 (MQ=255) aaGAGATCACCGCTGAGATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATc < 1:2976632/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAGTGGAGGAATCCGGTATCg < 1:2200259/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:248086/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:3487082/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAAACCGGTATCg > 1:5618898/1‑68 (MQ=255) gagaTCACCGCTGACATCCCTAACGTGTGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGt < 1:5495891/68‑1 (MQ=255) gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt < 1:319848/68‑1 (MQ=255) gaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGttt < 1:2953839/68‑1 (MQ=255) aTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTa > 1:3318245/1‑68 (MQ=255) tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc < 1:2051163/68‑1 (MQ=255) tCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTAc > 1:1514885/1‑68 (MQ=255) cACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACa < 1:5194898/68‑1 (MQ=255) aCCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACAt > 1:3930812/1‑68 (MQ=255) ccGCTGACATCCCGAACGTGGGTGAAGCTGCGCTGTCCAAACTGGAGGAATCCTGTCTCTTATACaca > 1:829288/1‑66 (MQ=255) gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg > 1:2762522/1‑68 (MQ=255) tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg > 1:666624/1‑68 (MQ=255) aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGCATCGTTTACATTGGTGCGGa < 1:1031040/68‑1 (MQ=255) cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg > 1:2133322/1‑68 (MQ=255) cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg > 1:1413315/1‑68 (MQ=255) cccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAg < 1:1175522/68‑1 (MQ=255) ccGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGt < 1:830413/68‑1 (MQ=255) gAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGa < 1:4439/68‑1 (MQ=255) cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg > 1:89875/1‑68 (MQ=255) cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg > 1:5633471/1‑68 (MQ=255) cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg > 1:2296559/1‑68 (MQ=255) gTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg < 1:5625353/68‑1 (MQ=255) gTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc < 1:5272949/68‑1 (MQ=255) gAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATAGTTTACTTTGGTGCTGATGTGCCCGGTGGCGa > 1:4188391/1‑68 (MQ=255) aaGCTGCGCTCTCCAAACTGGAGGAAACCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc < 1:1889990/68‑1 (MQ=255) gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt < 1:4973068/68‑1 (MQ=255) cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt > 1:5418139/1‑68 (MQ=255) cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt < 1:3318420/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:662239/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:2501615/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:3111780/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGGGACCGGTGGCGACATTCtggttg < 1:5118337/68‑1 (MQ=255) cTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt > 1:5104871/1‑68 (MQ=255) cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggttaggt > 1:1614745/1‑68 (MQ=255) cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt > 1:1351355/1‑68 (MQ=255) cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt < 1:102865/68‑1 (MQ=255) | TGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT > CP009273/4173706‑4173832 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |