| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AM260480 | 1,408,076 | 1 | . | T | 33.3% | 29.6 / 25.0 | 18 | coding (282/369 nt) | h16_B1250 | 5‑Carboxymethyl‑2‑hydroxymuconate isomerase |
| Reads supporting (aligned to +/- strand): ref base . (7/5); new base T (3/3); total (10/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 = 8.01e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACTCAGGCCACCACCTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTGGCGCACGGCGGTGTCGCGGCCGTCCAGGATCGACAGGCGCGCGTGGACAAAGGCCCGCTCGACCGCGTCCA > AM260480/1407986‑1408197 | aaCTCAGGCCACCACCTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTgcgc > 1:1064090/1‑132 (MQ=255) ccaccaCCTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCg < 1:238669/110‑1 (MQ=255) ccaccaCCTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCg > 2:238669/1‑110 (MQ=255) accaccTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCTATCACGTCGCAGACGGCCTGCGCCAGCGCCTgg < 2:28996/132‑1 (MQ=255) accaccTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCTATCACGTCGCAGACGGCCTGCGCCAGCGCCTgg > 1:28996/1‑132 (MQ=255) tGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGAAGCCGTGCGCGC‑CGGCtcc > 1:897881/1‑87 (MQ=38) tGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGAAGCCGTGCGCGC‑CGGCtcc < 2:897881/89‑3 (MQ=38) tAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTGGCGCACGGCGGTGTcgcg > 1:716450/1‑134 (MQ=255) cgcCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑gcg < 2:116594/58‑1 (MQ=255) cgcCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑gcg > 1:116594/1‑58 (MQ=255) cgcCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTg < 1:782694/93‑1 (MQ=255) cgcCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTg > 2:782694/1‑93 (MQ=255) ctACCGACACCTGCACGTCCTGGCCGGTGCCGAAGCCGTGCGCGC‑CGGCTCCGGCGATc > 1:585050/1‑59 (MQ=18) ctACCGACACCTGCACGTCCTGGCCGGTGCCGAAGCCGTGCGCGC‑CGGCTCCGGCGATc < 2:585050/59‑1 (MQ=18) cGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTGGCGCACGGCGGTGTCGCGGCCGTCCAGGATCGACAGGCGCGCGTGGACAAAGGCCCGCTc > 2:487996/1‑133 (MQ=255) gCCGGTGCCGAAGCCGTGCGCGC‑CGGCTCCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTgg > 2:178868/1‑67 (MQ=21) gCCGGTGCCGAAGCCGTGCGCGC‑CGGCTCCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTgg < 1:178868/67‑1 (MQ=21) ggTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTGGCGCACGGCGGTGTCGCGGCCGTCCAGGATCGACAGGCGCGCGTGGACAAAGGCCCGCTCGACCGCGTCCa < 1:487996/134‑1 (MQ=255) cGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCGGACGGCCTGCGCCAGCGCCTGgcgaa < 2:591378/63‑3 (MQ=255) cGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTCGCGCa > 1:591378/1‑63 (MQ=255) | AACTCAGGCCACCACCTGCTTGGCATAGGTTGCCCGCGCCATCTCGCAGACCTCTACCGACACCTGCACGTCCTGGCCGGTGCCGGAGCCG‑GCGCGCACGGCTTCGGCGATCACGTCGCAGACGGCCTGCGCCAGCGCCTGGCGCACGGCGGTGTCGCGGCCGTCCAGGATCGACAGGCGCGCGTGGACAAAGGCCCGCTCGACCGCGTCCA > AM260480/1407986‑1408197 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |