Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,929,372 | A→G | I386V (ATC→GTC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,929,372 | 0 | A | G | 98.5% | 285.2 / ‑3.9 | 68 | I386V (ATC→GTC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (38/29); total (38/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.41e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01 |
CGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATG > REL606/3929300‑3929444 | cGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCa > 1:1450244/1‑76 (MQ=255) gTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCAt < 2:544166/76‑1 (MQ=255) ttCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCAtt < 1:228579/76‑1 (MQ=255) ttCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCAtt < 1:346746/76‑1 (MQ=255) tACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCg > 1:1297766/1‑76 (MQ=255) tACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACCGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCg > 2:1078537/1‑76 (MQ=255) cccGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGAt > 1:912978/1‑76 (MQ=255) cccGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGAt < 1:1430446/76‑1 (MQ=255) cccGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGAt < 2:1809194/76‑1 (MQ=255) cGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATgg < 2:1541856/76‑1 (MQ=255) tAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGc > 2:1347444/1‑76 (MQ=255) tAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGc < 1:380419/76‑1 (MQ=255) tAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGc > 2:1508867/1‑76 (MQ=255) aaaGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCg < 1:581632/76‑1 (MQ=255) aagaagAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGa > 2:203785/1‑76 (MQ=255) aagaagAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGa > 2:1068274/1‑76 (MQ=255) aagaagAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGa > 2:459875/1‑76 (MQ=255) aagaagAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGa < 1:567056/76‑1 (MQ=255) agaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGGTCCTGCGCAAAGTCATTCACCCGATGGGCGaa < 1:1182067/76‑1 (MQ=255) gaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGaaa > 1:1037424/1‑76 (MQ=255) gaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGaaa > 1:1037448/1‑76 (MQ=255) gaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGaaa > 1:1037450/1‑76 (MQ=255) aagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAAt > 2:1485120/1‑76 (MQ=255) agagCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATc < 2:860943/76‑1 (MQ=255) agCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGa > 1:1648196/1‑76 (MQ=255) ctgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATg > 2:1627095/1‑76 (MQ=255) cACGTCTCAGGAAGAACTGCAGAAAATGTTGATCCTGCGCCAAGTCATTCCCCCGATGGGCGAAATCGATGCAATg < 1:516504/76‑1 (MQ=255) cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATg < 2:17530/76‑1 (MQ=255) cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATg > 1:379062/1‑76 (MQ=255) cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATg < 2:1848731/76‑1 (MQ=255) cGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGa > 2:1278517/1‑76 (MQ=255) tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTc > 1:959743/1‑76 (MQ=255) tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTc > 1:1066786/1‑76 (MQ=255) cAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTcc < 1:64734/76‑1 (MQ=255) aGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCt > 2:229751/1‑76 (MQ=255) ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTc < 2:1438807/76‑1 (MQ=255) aagaaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCAt > 1:124009/1‑76 (MQ=255) aagaaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCAt > 1:98784/1‑76 (MQ=255) agaaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCAtt > 1:1319059/1‑76 (MQ=255) agaaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGACGGGCGAAATCGATGCAATGGAATTCCTCAtt > 2:1143832/1‑76 (MQ=255) aaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaa > 1:1430483/1‑76 (MQ=255) aaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaa > 1:1700791/1‑76 (MQ=255) aaCTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaa > 1:825515/1‑76 (MQ=255) cTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaata < 1:985849/76‑1 (MQ=255) cTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaata < 1:1035192/76‑1 (MQ=255) cTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaata < 1:944719/76‑1 (MQ=255) cTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaata < 1:754033/76‑1 (MQ=255) cTGCAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaata < 1:391158/76‑1 (MQ=255) cAGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAAc > 2:571534/1‑76 (MQ=255) aGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACt < 1:895664/76‑1 (MQ=255) aGAAAATGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACt > 2:853414/1‑76 (MQ=255) aaaaTGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTgg < 1:1916228/76‑1 (MQ=255) aaaaTGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTgg > 2:1358162/1‑76 (MQ=255) aaaaTGTGGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTgg < 1:1935388/76‑1 (MQ=255) tgtgGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt < 2:1426243/76‑1 (MQ=255) tgtgGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt < 1:515165/76‑1 (MQ=255) tgtgGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt > 1:1910166/1‑76 (MQ=255) tgtgGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt > 1:1292453/1‑76 (MQ=255) tgGATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATga < 2:1197738/76‑1 (MQ=255) tgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATga < 2:91966/76‑1 (MQ=255) ggATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgac < 2:1346494/76‑1 (MQ=255) gATCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgacc < 2:1173339/76‑1 (MQ=255) aTCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgacca > 1:982192/1‑76 (MQ=255) tCCTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaa > 1:648019/1‑76 (MQ=255) cTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaaga > 1:272352/1‑76 (MQ=255) cTGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaaga < 1:1720287/76‑1 (MQ=255) tGCGCAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaagac > 1:1834015/1‑76 (MQ=255) cAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAt > 1:785214/1‑76 (MQ=255) cAAAGTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAt > 1:765739/1‑76 (MQ=255) aaagTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATg < 1:1134375/76‑1 (MQ=255) | CGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATG > REL606/3929300‑3929444 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |