Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,164,822 | T→G | F1323V (TTC→GTC) | 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,164,822 | 0 | T | G | 95.3% | 262.6 / 0.2 | 64 | F1323V (TTC→GTC) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base T (3/0); new base G (43/18); total (46/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.53e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01 |
GAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAA > REL606/4164749‑4164895 | gAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTc > 1:582869/1‑76 (MQ=255) aaCGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCa > 2:1085451/1‑76 (MQ=255) aaCGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGACATGCCAGAATCCGTCa < 2:399473/76‑1 (MQ=255) cGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAAc > 1:332722/1‑76 (MQ=255) gTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGt > 1:839412/1‑76 (MQ=255) tCGTACCAAGTTGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCATAATCCGTCAACGTa > 1:157838/1‑76 (MQ=255) cGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTAt > 1:1665961/1‑76 (MQ=255) gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTAtt < 1:59428/76‑1 (MQ=255) tACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTAttg > 2:830259/1‑76 (MQ=255) aCCAAGATGTGTAAAAAGATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATACGTCAACGTAttgt > 1:157901/1‑76 (MQ=255) aCCAAGATGTGTAAAAAGATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATACGTCAACGTAttgt > 1:157918/1‑76 (MQ=255) ccAAGATGTATAAAAACATCGTGGACGGNTACCATCAGATGGAGCCGGTCNTGCCAGAATCCGTCAACGTAttgtt > 1:276275/1‑76 (MQ=255) ccAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTAttgtt < 2:973432/76‑1 (MQ=255) cAAGATGTATAAAAATATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTAttgttg > 1:174746/1‑76 (MQ=255) aaGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGa < 2:1002949/76‑1 (MQ=255) aaGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGa > 2:1275845/1‑76 (MQ=255) aGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGaa < 2:951154/76‑1 (MQ=255) gATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGaaa > 1:937266/1‑76 (MQ=255) aTGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAag > 1:720476/1‑76 (MQ=255) tGTATAAAAACATCGTGGACGGCAACCATCAGATGGACCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAaga > 1:1537154/1‑76 (MQ=255) gTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGCAagag < 1:113948/76‑1 (MQ=255) tataAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAgaga > 1:681854/1‑76 (MQ=255) tataAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAgaga > 1:384640/1‑76 (MQ=255) taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGAtt < 1:1266144/76‑1 (MQ=255) taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGAtt < 1:1052688/76‑1 (MQ=255) naaaaaCATCGTGGANGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTNGAAAGAGAtt < 2:1237705/75‑1 (MQ=255) cATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGAttcgttcg < 2:764697/76‑1 (MQ=255) aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGc > 1:1007097/1‑76 (MQ=255) tCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCt > 2:1313250/1‑76 (MQ=255) gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTgg < 1:598249/76‑1 (MQ=255) tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg > 2:1548543/1‑76 (MQ=255) tGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg > 2:1196858/1‑76 (MQ=255) ggACGGCAACCATCAGATGGAGCCGGGCATGTCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGt > 1:1081763/1‑76 (MQ=255) ggACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGt > 2:377274/1‑76 (MQ=255) gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTa > 2:1298516/1‑76 (MQ=255) gACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTa > 1:1043861/1‑76 (MQ=255) cGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATc > 2:1190931/1‑76 (MQ=255) cGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATc < 1:604469/76‑1 (MQ=255) ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCa < 1:1046244/76‑1 (MQ=255) cAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAAc > 1:255550/1‑76 (MQ=255) cAACCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAAc > 2:1007038/1‑76 (MQ=255) aCCATCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACAt < 1:1295240/76‑1 (MQ=255) aTCAGATGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGa > 2:1357043/1‑76 (MQ=255) tCAGATGGAGCCGGGCATGCCAGAATCCGTCAACCTANTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGaa > 1:1785101/1‑76 (MQ=255) cAGATGGAGCCGGGCATGCCAGAATCCNTCAACGTNNTGTTGGAAGAGATTCGTTCGATGGGTATCAACATCGAAc > 1:1796107/1‑76 (MQ=255) aTGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg > 2:1057309/1‑76 (MQ=255) aTGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg > 1:423127/1‑76 (MQ=255) tGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGa > 1:1010470/1‑76 (MQ=255) tGGAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGa > 1:1470239/1‑76 (MQ=255) ggAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGaa > 2:476085/1‑76 (MQ=255) gAGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAg > 1:724918/1‑76 (MQ=255) aGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGACGa > 1:334310/1‑76 (MQ=255) aGCCGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGa > 1:450306/1‑76 (MQ=255) ccGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACg < 1:1054682/76‑1 (MQ=255) cGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGa > 1:344736/1‑76 (MQ=255) cGGGCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGa > 2:911830/1‑76 (MQ=255) ggCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGt < 1:1022058/76‑1 (MQ=255) gCATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTa > 1:707461/1‑76 (MQ=255) cATGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTaa < 1:369566/76‑1 (MQ=255) tGCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAtt < 1:1190932/76‑1 (MQ=255) gCCAGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCTCTGGGTATCAACATCGAACTAGACGAAGAGTAATtc > 1:21129/1‑76 (MQ=255) aGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCg > 2:608684/1‑76 (MQ=255) aGAATCCGTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCg > 1:578786/1‑76 (MQ=255) aTCCGTCAACGTATTGTTGCAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTc > 1:1546316/1‑76 (MQ=255) ccTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCaa > 2:206414/1‑76 (MQ=255) | GAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAA > REL606/4164749‑4164895 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |