Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,404,055 | T→A | I21N (ATC→AAC) | rpsR → | 30S ribosomal protein S18 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,404,055 | 0 | T | A | 98.8% | 326.8 / ‑4.4 | 80 | I21N (ATC→AAC) | rpsR | 30S ribosomal protein S18 |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base A (50/29); total (50/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
GACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCAC > REL606/4403984‑4404127 | gaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGAc > 1:966137/1‑76 (MQ=255) aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACt < 1:1088396/76‑1 (MQ=255) aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACt < 2:744878/76‑1 (MQ=255) aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACt < 2:1512571/76‑1 (MQ=255) cTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACta > 1:1493858/1‑76 (MQ=255) cTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACta > 2:587159/1‑76 (MQ=255) tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtat > 2:1306556/1‑76 (MQ=255) tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtat > 1:768561/1‑76 (MQ=255) aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtata > 2:1175140/1‑76 (MQ=255) aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtata > 2:1216517/1‑76 (MQ=255) aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtata > 1:234149/1‑76 (MQ=255) ccATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCTGAAGGCGTTCAAGAGAACGACTATaaa > 1:1184536/1‑76 (MQ=255) ccATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATaaa < 2:829196/76‑1 (MQ=255) atatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGa > 2:1504643/1‑76 (MQ=255) atatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGa > 1:367969/1‑76 (MQ=255) atatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGa < 1:210208/76‑1 (MQ=255) atatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGa > 2:58310/1‑76 (MQ=255) tatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGat > 1:1580663/1‑76 (MQ=255) tatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGat > 1:53686/1‑76 (MQ=255) tatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGat > 2:56001/1‑76 (MQ=255) atGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGata > 1:1491188/1‑76 (MQ=255) tGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGatat > 1:847059/1‑76 (MQ=255) gCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCg < 1:79119/76‑1 (MQ=255) cACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGc > 2:766609/1‑76 (MQ=255) aCGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCt < 2:1333674/76‑1 (MQ=255) cGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTa < 2:1214086/76‑1 (MQ=255) gTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTAc > 1:727826/1‑76 (MQ=255) tATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGc > 2:719430/1‑76 (MQ=255) tATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGc > 1:817449/1‑76 (MQ=255) aTTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCt > 2:1190972/1‑76 (MQ=255) aTTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCt < 1:1561142/76‑1 (MQ=255) aTTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCt < 1:31557/76‑1 (MQ=255) tttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTg > 1:1379418/1‑76 (MQ=255) ttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGa > 1:1328626/1‑76 (MQ=255) ttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGa > 1:303645/1‑76 (MQ=255) ttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGa > 2:1123328/1‑76 (MQ=255) ttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGa > 2:1653459/1‑76 (MQ=255) ttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGa > 1:526684/1‑76 (MQ=255) tCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaa > 2:811646/1‑76 (MQ=255) ccGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaa > 1:1176761/1‑76 (MQ=255) ccGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaa < 2:84813/76‑1 (MQ=255) ccGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaa > 1:238715/1‑76 (MQ=255) ccGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaa > 2:314950/1‑76 (MQ=255) tcgtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAAc < 1:684388/76‑1 (MQ=255) tcgtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAAc < 1:560882/76‑1 (MQ=255) tcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTAc > 2:831234/1‑76 (MQ=255) tcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTAc < 2:1578465/76‑1 (MQ=255) gCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACAt < 1:1412906/76‑1 (MQ=255) gTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCAcc > 2:898697/1‑76 (MQ=255) ttCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCg < 2:262445/76‑1 (MQ=255) tCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGa > 1:675076/1‑76 (MQ=255) tGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaaa < 2:1495185/76‑1 (MQ=255) gCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAg > 1:457728/1‑76 (MQ=255) gCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAg < 1:196496/76‑1 (MQ=255) cGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg < 1:578033/76‑1 (MQ=255) cGTTACACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg < 1:171133/76‑1 (MQ=255) ttCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTa < 1:1088483/76‑1 (MQ=255) tCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTaa > 1:539691/1‑76 (MQ=255) cACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAACGCGGTAAg > 2:44891/1‑76 (MQ=255) aCCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGa > 2:923544/1‑76 (MQ=255) aCCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGa > 1:469823/1‑76 (MQ=255) cgcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGAtt < 2:1132838/76‑1 (MQ=255) gcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg < 2:359929/76‑1 (MQ=255) gcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg < 1:451783/76‑1 (MQ=255) gcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg < 1:1657607/76‑1 (MQ=255) gAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTcc < 2:938448/76‑1 (MQ=255) aaGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTccc < 1:456901/76‑1 (MQ=255) ggCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCaa < 1:1518040/76‑1 (MQ=255) gTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGcc > 1:14273/1‑76 (MQ=255) ttCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg < 1:947187/76‑1 (MQ=255) ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg > 1:210258/1‑76 (MQ=255) ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg < 1:1216517/76‑1 (MQ=255) ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg > 2:1133076/1‑76 (MQ=255) ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg > 1:737345/1‑76 (MQ=255) tCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGt > 1:1574415/1‑76 (MQ=255) tCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGt > 2:989930/1‑76 (MQ=255) cAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTa > 1:1310715/1‑76 (MQ=255) agagAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATc < 1:134729/76‑1 (MQ=255) gagaTCGACTATAAAGATATCCCTACGCTGAAAAACTACATCCCCAAAGGCGGTAAGATTTTCCCAAGCCGTATCa < 2:543110/76‑1 (MQ=255) gagaACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCa > 2:415518/1‑76 (MQ=255) gagaACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCa < 1:1174484/76‑1 (MQ=255) agaACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCAc > 2:827471/1‑76 (MQ=255) agaACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCAc > 1:980839/1‑76 (MQ=255) | GACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCAC > REL606/4403984‑4404127 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |