Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,455,028 | 0 | C | T | 28.1% | 44.4 / 18.3 | 32 | M194I (ATG→ATA) | fbp | fructose‑1,6‑bisphosphatase 1 |
Reads supporting (aligned to +/- strand): ref base C (17/6); new base T (7/2); total (24/8) | |||||||||||
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.07e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTA > NC_000913/4454950‑4455105 | gTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCAt > 2:921305/1‑81 (MQ=255) cACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTg > 2:921299/1‑81 (MQ=255) aCCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTgg < 2:921244/81‑1 (MQ=255) ccGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGcac < 1:921245/81‑1 (MQ=255) gTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAg > 2:921283/1‑81 (MQ=255) cGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAgg > 2:921295/1‑79 (MQ=255) gAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAg > 2:921294/1‑80 (MQ=255) aaaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGa > 2:921304/1‑80 (MQ=255) aaaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGa > 2:921314/1‑80 (MQ=255) aaaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGa > 2:921302/1‑80 (MQ=255) aaCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGaa > 2:921263/1‑80 (MQ=255) aCTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGaaaa > 2:921303/1‑81 (MQ=255) cTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAAc < 2:921246/81‑1 (MQ=255) tAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGc > 2:921284/1‑81 (MQ=255) aaTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGcc > 1:921291/1‑81 (MQ=255) aaTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGcc > 1:921271/1‑81 (MQ=255) tGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGa < 1:921247/81‑1 (MQ=255) aGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGc > 1:921320/1‑80 (MQ=255) aGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGc > 2:921281/1‑80 (MQ=255) ttCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAg > 1:921288/1‑81 (MQ=255) tCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCg > 1:921296/1‑81 (MQ=255) cGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATc < 2:921248/79‑1 (MQ=255) gATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCgtaag > 1:921285/1‑80 (MQ=255) tGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCgtaagta > 2:921311/1‑80 (MQ=255) aGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAgg < 2:921249/81‑1 (MQ=255) tAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCg > 1:921326/1‑81 (MQ=255) ttttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGaa < 2:921250/81‑1 (MQ=255) ttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAaca > 2:921298/1‑81 (MQ=255) ttCTCCGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCAt > 1:921318/1‑81 (MQ=255) ccGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg > 2:921319/1‑81 (MQ=255) cGGGAAGCGTATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggt > 2:921310/1‑81 (MQ=255) gggAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtg < 2:921251/81‑1 (MQ=255) gcgtATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGgtgt > 1:921308/1‑80 (MQ=255) gcgcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTa < 1:921252/81‑1 (MQ=255) cgcATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTa > 1:921329/1‑80 (MQ=255) | GTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTA > NC_000913/4454950‑4455105 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |