| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,648,148 | 0 | A | C | 25.0% | 51.2 / 12.0 | 24 | G1380G (GGT→GGG) | yfhM | alpha2‑macroglobulin |
| Reads supporting (aligned to +/- strand): ref base A (17/1); new base C (0/6); total (17/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.20e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTC > NC_000913/2648084‑2648210 | gCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCa < 1:3738576/65‑1 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 1:3332867/1‑65 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc < 1:1218379/65‑1 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 1:1791438/1‑65 (MQ=255) aGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCAcc > 1:1762742/1‑65 (MQ=255) tttCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCat > 1:2689607/1‑65 (MQ=255) ttcttcCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCata > 1:520079/1‑65 (MQ=255) ccAGCAAGGAGAGCATTACCGCGTTGTCGCGCAGTAAACTCCCGAAACCCCCCAGCCATATCCGc < 1:3264879/65‑1 (MQ=255) gaaaGGAGAGCATTACCGCGTTGTCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCtca < 1:2738357/63‑1 (MQ=255) gaaaGGAGAGCATGAACGCGTTGTCGCGCAGTGAACTCCCGTAATCCCCCAGCCATATCCGCtca < 1:2458861/63‑1 (MQ=255) aagaagAGCATTACCGCGTTGCCGCGCAGTAAACTACCTTAACCACCCAGCCATATCCGCtcatc < 1:1603363/65‑1 (MQ=255) aaGGGGAGCATTACCGCGTTGTCGCGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCtcatc < 1:528315/65‑1 (MQ=255) ggAGAGCATTACCGGTTTGTCGCGCAGTGAACTCCCGTAACCCCCCAGCCATATCCGCTCATCAc < 1:338742/65‑1 (MQ=255) gaggGCATAAACCGGTTGGCGCGAGGTAAACTACCCAATTCCCCCAGCCATATCCGCTCATCACt < 1:103719/65‑1 (MQ=255) gagaGTATAACCGCGTTGGCCCCCAGAAAATTTCGGTAACCACCCAGCCATATCCGCTCATCACt < 1:280314/65‑1 (MQ=255) gagaGCTTTAACGCGTTGGCGCGCAGTAAACTACCGTAATCCCCCAGCCATATCCGCTCATCACt < 1:3134211/65‑1 (MQ=255) gagaGCATTATGGCGTTGCGGCGTACTGATCTTCCATAACCCCCCAGCCATATCCGCTCATCACt < 1:1124550/65‑1 (MQ=255) gagaGCATTACCGGGTTGCCGCGCAGTAAACTACCGTAACCCCCCAGCCATATCCGCTCATCACt < 1:2028463/65‑1 (MQ=255) gagaGCATTACCGCGTTGGCGCGGAGTAACTTCCCTTAATCACCCAGCCATATCCGCTCATCACt < 1:2125935/65‑1 (MQ=255) gagaGCATTAAGGGGTTGCCGCGCAGTGACCTCCCGTAACCCCCCAGCCATATCCGCTCATCACt < 1:1381038/65‑1 (MQ=255) gagaGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACt > 1:3991566/1‑65 (MQ=255) gagaGCATAACCGCGTTGTGGCGCGGAAAACTCCCATAATCCCCCAGCCATATCCGCTCATCACt < 1:953820/65‑1 (MQ=255) aTTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcg > 1:544998/1‑65 (MQ=255) aTTAACGCGTTGTCGCCCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTAcg > 1:1463906/1‑65 (MQ=255) cgcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACgcggcg > 1:1499935/1‑65 (MQ=255) gcgTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGt > 1:3858711/1‑65 (MQ=255) ttGTCGCCCGTTAAACTCCCGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGtttt < 1:516189/65‑1 (MQ=255) gcgcAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcg > 1:2588057/1‑65 (MQ=255) gcAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgcc > 1:1842443/1‑65 (MQ=255) aGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccag > 1:224501/1‑65 (MQ=255) aGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTcagcgccag > 1:2680102/1‑65 (MQ=255) cTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAAt > 1:2614480/1‑65 (MQ=255) aCCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCg > 1:718122/1‑65 (MQ=255) cttAATCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct < 1:2867147/63‑1 (MQ=255) cGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct > 1:3449825/1‑65 (MQ=255) cGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct > 1:3285021/1‑65 (MQ=255) cGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct > 1:3234693/1‑65 (MQ=255) cGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct < 1:1585975/65‑1 (MQ=255) cGTAACCCCCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGct < 1:2400801/65‑1 (MQ=255) aaTCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTCTTCAGCGCCAGCGCAATCGcttct > 1:1335788/1‑65 (MQ=255) tCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCttcttc > 1:1928283/1‑65 (MQ=255) | GCAGTTTATTTTCTTCCAGCAAGGAGAGCATTAACGCGTTGTCGCGCAGTGAACTACCGTAATCACCCAGCCATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTC > NC_000913/2648084‑2648210 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |