Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,907 | A→G | T231A (ACC→GCC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,907 | 0 | A | G | 90.5% | 62.4 / ‑1.0 | 21 | T231A (ACC→GCC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (9/10); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01 |
GTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGA > REL606/3928832‑3928980 | gTTCTGGTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTa > 1:123532/1‑76 (MQ=255) ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCt > 1:165692/1‑76 (MQ=255) ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCt > 1:763241/1‑76 (MQ=255) gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCttt < 2:595778/76‑1 (MQ=255) ctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTg < 2:445362/76‑1 (MQ=255) tgAACGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCAGGTAAANGGTGAAGTTGTTGCTTCTGCCTTGGa < 2:680757/76‑1 (MQ=255) aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACg > 2:765241/1‑76 (MQ=255) aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACg < 2:725073/76‑1 (MQ=255) aaCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGc > 2:279127/1‑76 (MQ=255) cGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCAt < 2:316481/76‑1 (MQ=255) gtagTAACCGAGATGCAGCGTCCGGNAAAAGGAGAAGTTGTTGATTCTGCCTNTGACGAACCNGCATCTCGNCAcn > 1:675443/3‑75 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTc < 1:244579/76‑1 (MQ=255) aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCgg > 2:408381/1‑76 (MQ=255) ggCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGa < 2:439257/75‑1 (MQ=255) gCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTCACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATg < 1:525184/76‑1 (MQ=255) cGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATgg < 1:301113/76‑1 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTGCCTTTAGCGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCga > 1:776808/1‑76 (MQ=255) aaGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGaa < 2:767478/76‑1 (MQ=255) gttgttGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc > 1:284055/1‑76 (MQ=255) tgttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc > 1:644473/1‑76 (MQ=255) tgttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc > 2:81472/1‑76 (MQ=255) ttCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt < 2:511364/76‑1 (MQ=255) cTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAGACGCCTGGTTGa < 2:123532/76‑1 (MQ=255) | GTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGA > REL606/3928832‑3928980 |
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: ‑ |