Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,161,105 | G→A | E84K (GAA→AAA) | 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,161,105 | 0 | G | A | 97.7% | 163.9 / ‑3.9 | 45 | E84K (GAA→AAA) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (20/23); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01 |
TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC > REL606/4161032‑4161176 | tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa > 1:1052425/1‑76 (MQ=255) tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa < 1:40210/76‑1 (MQ=255) agagCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATg > 1:222732/1‑76 (MQ=255) agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc > 2:368864/1‑76 (MQ=255) agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc < 1:713823/76‑1 (MQ=255) agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc > 1:1807076/1‑76 (MQ=255) tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGCAATGTCaaa > 1:1282670/1‑76 (MQ=255) tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCaaa < 1:613074/76‑1 (MQ=255) aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt < 2:925780/76‑1 (MQ=255) aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt > 2:1690714/1‑76 (MQ=255) aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt < 2:825590/76‑1 (MQ=255) tAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGc > 2:524633/1‑76 (MQ=255) aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCg > 2:75262/1‑76 (MQ=255) aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAACCCGTGGCg > 2:75245/1‑76 (MQ=255) ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg < 2:1468207/76‑1 (MQ=255) ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg < 2:1856855/76‑1 (MQ=255) tCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGa < 2:583249/76‑1 (MQ=255) ccGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGAc > 1:149582/1‑76 (MQ=255) gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt > 2:573768/1‑76 (MQ=255) gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt > 2:956096/1‑76 (MQ=255) aGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTa > 2:164592/1‑76 (MQ=255) cTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTAtt > 2:1508739/1‑76 (MQ=255) tCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCt < 2:509168/76‑1 (MQ=255) aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc > 1:333893/1‑76 (MQ=255) aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc > 2:1469339/1‑76 (MQ=255) ccGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGtt < 2:1117745/76‑1 (MQ=255) gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa > 1:1201179/1‑76 (MQ=255) gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa < 2:563303/76‑1 (MQ=255) gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa > 2:376307/1‑76 (MQ=255) ttGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACt < 1:1000082/76‑1 (MQ=255) gCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCg > 1:124138/1‑76 (MQ=255) aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt < 1:791103/76‑1 (MQ=255) aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt > 1:1497435/1‑76 (MQ=255) aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCATGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt > 1:1345740/1‑76 (MQ=255) ccGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTgg < 2:1151066/76‑1 (MQ=255) cGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGt < 2:1509742/76‑1 (MQ=255) ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 1:402648/76‑1 (MQ=255) ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg < 2:376228/76‑1 (MQ=255) tgtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGAt < 1:1615596/76‑1 (MQ=255) gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc < 1:376307/76‑1 (MQ=255) gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc < 1:1071087/76‑1 (MQ=255) gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc > 1:1131320/1‑76 (MQ=255) tGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTAt < 2:814515/76‑1 (MQ=255) gACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATg < 2:963818/76‑1 (MQ=255) cGTGCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg < 1:387146/76‑1 (MQ=255) cGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg < 1:1500552/76‑1 (MQ=255) tCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcg < 1:111121/76‑1 (MQ=255) ccAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcgc < 1:166771/76‑1 (MQ=255) | TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC > REL606/4161032‑4161176 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |