Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 759,433 | 0 | A | C | 57.9% | 3.5 / 15.5 | 21 | intergenic (+42/‑126) | mmuM_1/CJ019_00729 | Homocysteine S‑methyltransferase/L‑methionine gamma‑lyase |
Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (11/0); total (19/2) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.21e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCTTTTGCAGC‑GCCC‑GGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATA > ACB122/759303‑759559 | aCTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGCAACGCAACCTGGACTAA‑GCTCCCCAGTGGCCCGGCGGCTTTTGACGC‑GCCC‑GGCCCGCATccc > 4:113899/1‑134 (MQ=255) gTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCTTTTGCAGC‑GCCC‑GGCCCGCCGCCATTaaa > 6:178307/1‑134 (MQ=255) tGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGGCCACGTGGCCGGGGGGGGTTTGGCGC‑GCCCGGGCCCCCGCCCATTaaa > 2:300338/1‑134 (MQ=255) tGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCCGGCCGCGTTTGCAGC‑GCCC‑GGCCCGCAGCCATTaaaa > 8:134666/1‑134 (MQ=255) ggggTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGGGGCGTGTGCGGCGGCCC‑CGCCGCCCCCCCGTAAAAAAAAc > 7:158319/1‑135 (MQ=255) tCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCGTTTGCCGC‑GCCC‑GGCCCGCCGCCATTAAAAAAAACGCCCGTCAGACCCggc > 4:62639/1‑134 (MQ=255) ggcggTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGACACGCAACCTGGCCTAA‑GCGCCCCCGTGGCCGGGCGGGTGTTTCCGC‑GCCC‑CGCCCGGCGCCATTAAAAAAAACCCCAGGCAGACCCGGCtt > 4:167261/1‑135 (MQ=255) ggcggTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCCGGCGGGGTTTTGGCC‑CCCC‑GCCCCCGCGCCCATAAAAAAAACCCCCCGCAGACCCggcgg > 1:371728/1‑132 (MQ=255) cccGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCTTTTGCAGC‑GCCC‑GGCCCGCCGCCATTAAAAAAAAACCCAGTCAGCCCCGGGTTCCCACGGTTTTAAGgccgcc > 7:191940/1‑135 (MQ=255) ccGGACCACACCCCCGCGGTTAAACG‑AACCCGGCCAAACGCGCCCCGGGGGCGGGGGGGCTTTTGCAGC‑GCCG‑GGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCg < 6:56509/135‑1 (MQ=255) acaTCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCCGGCGGCGTTTGGAGC‑GCCC‑GGCCCGCCGCCAATAAAAAAAAACCCACGCCGCCCCCGCGTCCCACGGTTTTAAGGCCTCCGGCagcga > 6:202184/1‑131 (MQ=255) gCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCTTTTGCCGC‑GCCC‑GGCCCGCAGCCATTAAAAAAAACGCCAGTCCGCCCCGGATTCCATCGTTTTTCAGGCCGCCGTCAtggcccg > 7:256286/1‑127 (MQ=255) gCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCCGGCGGCTTTTGCAGC‑GCCC‑CGCCCGCCGCCAATAAAAAAAAACCCCCTCAGCACCCGATAGCCTCGTTTTTAAGGGCTGCCTGCTCTcccg > 3:158992/1‑132 (MQ=255) cGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCCGGGGGCATTTTGAGC‑GCCC‑CGCCCCCCGCCCATAAAAAAAAACCCACTCCGACCCGGCTTCCCTACGTTTTGAGGCCACCGGGATCGCCCAATAg > 2:292403/1‑135 (MQ=255) tGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCGTTTGCAGC‑GCCC‑GGCCCGCCGCCCTTAAAAAAAACCCCCGTCAGCCCCGGCTTCCCTCGTTTTTGAGGCCTCCGGCCTCTCCCGATAGTCCGc > 6:111587/1‑135 (MQ=255) gcgcCCACGGGGCGGGGCGGCTTTTGACGG‑GCCC‑GGCCGGCTGCGATTAAAAAAATCGCCAGTCAGACCCGGATGCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATc < 8:290331/134‑1 (MQ=255) aCGTGGCCGGGCGGCTTTTGCCGC‑GCCC‑GGCCCGCCGCCAATAAAAAAAACGCCAGTCAGCCCCGGATTCCAAGGTTTTTCAGGCCTCCGGCATCTCCCGATAGTCCGCGGCAGGGCGCTGCTGGAACGGGGTac > 2:58945/1‑135 (MQ=255) c‑ggcccgCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGata > 4:176586/1‑135 (MQ=255) c‑ggcccgCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGata > 3:211013/1‑135 (MQ=255) c‑ggcccgCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGata > 3:184564/1‑135 (MQ=255) c‑ggcccgCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGata > 3:127618/1‑135 (MQ=255) | ACTGGCGTGCGCGGGGTGCCAGCATGGTCGGCGGTTGCTGTGGCATCGGCCCGGAGCACATCGCCGAGCTGAAACGCAACCTGGCCTAA‑GCGCCCACGTGGCCGGGCGGCTTTTGCAGC‑GCCC‑GGCCCGCAGCCATTAAAAAAAACGCCAGTCAGACCCGGATTCCATCGTTTTTCAGGCCTCCGTCATCTCCCTATAGTCCGCTGAAGGCCGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATA > ACB122/759303‑759559 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |