Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,947 | A→G | E244G (GAA→GGA) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,947 | 0 | A | G | 98.4% | 254.0 / ‑4.8 | 61 | E244G (GAA→GGA) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (34/26); total (35/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
CGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTC > REL606/3928877‑3929021 | cGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATg < 2:1352176/75‑1 (MQ=255) gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATgg < 2:492567/75‑1 (MQ=255) gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATgg > 1:1151549/1‑75 (MQ=255) gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAAGTTGCGGGAATgg < 1:1186414/75‑1 (MQ=255) tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGt < 2:35059/75‑1 (MQ=255) tGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGa > 1:1076914/1‑75 (MQ=255) ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGAt < 1:1013202/75‑1 (MQ=255) ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGAt > 2:358329/1‑75 (MQ=255) ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGAt < 1:744895/75‑1 (MQ=255) gTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATc < 2:1322232/75‑1 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCg > 1:1188062/1‑75 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCga < 1:785636/75‑1 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCgag > 1:350805/1‑75 (MQ=255) aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGaa > 1:1127847/1‑75 (MQ=255) ggTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAg > 2:1164789/1‑75 (MQ=255) tGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGc > 1:815750/1‑75 (MQ=255) gAAGTTGTTGCTTCTACGTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGcc > 2:33478/1‑75 (MQ=255) gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGcc > 2:163208/1‑75 (MQ=255) gttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGc > 2:147488/1‑75 (MQ=255) gttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGc > 1:1521254/1‑75 (MQ=255) tgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCt > 2:214234/1‑75 (MQ=255) cTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTgg > 2:656235/1‑75 (MQ=255) ttCTAGCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGt > 1:1241230/1‑75 (MQ=255) ttCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGt > 2:22103/1‑75 (MQ=255) tCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGtt > 2:1239731/1‑75 (MQ=255) tCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGtt > 1:487406/1‑75 (MQ=255) tACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa > 1:514916/1‑75 (MQ=255) tACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa > 1:183649/1‑75 (MQ=255) tttGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac < 1:331027/75‑1 (MQ=255) ttGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca < 1:968923/75‑1 (MQ=255) ttGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca > 1:698122/1‑75 (MQ=255) gACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag > 1:623480/1‑75 (MQ=255) gACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag > 2:478937/1‑75 (MQ=255) gACGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag > 2:1496011/1‑75 (MQ=255) aCGAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga < 2:336813/75‑1 (MQ=255) gAACCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaa < 1:1150813/75‑1 (MQ=255) aaCCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag < 2:599338/75‑1 (MQ=255) aaCCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag < 1:691320/75‑1 (MQ=255) aCCCGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga > 2:705130/1‑75 (MQ=255) ccGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg < 1:1578410/75‑1 (MQ=255) ccGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg > 2:1383386/1‑75 (MQ=255) ccGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg > 2:1018111/1‑75 (MQ=255) cGCATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGt > 2:413084/1‑75 (MQ=255) cATCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTa < 2:1084788/75‑1 (MQ=255) aTCTCGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTAt < 2:104809/75‑1 (MQ=255) tctcGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc < 1:1598392/75‑1 (MQ=255) tctcGCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc < 1:793078/75‑1 (MQ=255) gCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc > 2:503827/1‑75 (MQ=255) gCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc < 2:156699/75‑1 (MQ=255) gCCACGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc < 2:381693/75‑1 (MQ=255) aCGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGc > 1:553068/1‑75 (MQ=255) aCGTTCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGc > 1:469664/1‑75 (MQ=255) gtgcAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc < 2:494293/72‑1 (MQ=255) ttCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg < 1:147488/75‑1 (MQ=255) tCAGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGa > 2:969827/1‑75 (MQ=255) aGGTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACt > 2:537421/1‑75 (MQ=255) ggTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTc < 2:694748/75‑1 (MQ=255) gTTGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc > 1:161276/1‑75 (MQ=255) ttGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa < 2:692051/75‑1 (MQ=255) ttGCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa < 2:904298/75‑1 (MQ=255) gCGGGAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc > 1:795292/1‑75 (MQ=255) aaaTGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc > 2:502953/1‑75 (MQ=255) | CGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTC > REL606/3928877‑3929021 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |