| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 3,684,870 | C→T | 100% | pseudogene (189/582 nt) pseudogene (1252/1645 nt) |
ilvG ← ilvG ← |
ECK3760:JW3741:b3768; acetolactate synthase II, large subunit, C‑ter fragment ECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 3,684,870 | 0 | C | T | 96.9% | 87.3 / ‑3.7 | 32 | pseudogene (189/582 nt) pseudogene (1252/1645 nt) | ilvG ilvG | ECK3760:JW3741:b3768; acetolactate synthase II, large subunit, C‑ter fragment ECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (9/22); total (9/23) | |||||||||||
| 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.98e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCCAGCTCTTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACAACGGTATCGTTCGGTCGCGCGACTTGTGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGTACCTAAACCGCTGGAGG > AP009048/3684808‑3684937 | cccAGCTCTTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTcgcgc > 1:3725123/1‑71 (MQ=255) gCTCTTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACt < 1:4553787/71‑1 (MQ=255) tctTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTcgcg > 1:476074/1‑64 (MQ=255) tctTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTcgc > 1:778626/1‑63 (MQ=255) gCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTtgc < 1:3684443/68‑2 (MQ=255) tcatcaTGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCaa < 1:3440394/69‑1 (MQ=255) catcatGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTgcgc > 1:4901506/1‑63 (MQ=255) atcatGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTtgcg < 1:1818965/61‑3 (MQ=255) tcatTAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCaa < 1:1063935/66‑1 (MQ=255) tcatGAAAGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTc > 1:2267931/1‑48 (MQ=39) aaGAGCCGTCACCGGAGATACAGACGACGGTATCGTTTGGTcgcgc < 1:2343200/46‑1 (MQ=255) ccGTCACCGGAGATACAGACGACGGTATCGTTTGGTcg < 1:2142063/38‑1 (MQ=38) ccGTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTtgc < 1:1902698/48‑2 (MQ=255) cGTCACCGGATATACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAAcc < 1:1095742/70‑1 (MQ=255) gTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTGcgcg < 1:3699167/49‑1 (MQ=255) gTCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCaa < 1:1436267/71‑1 (MQ=255) tCACCGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCaaa < 1:1657701/71‑1 (MQ=255) aCCGGAGATACAGACGACGGTACCGTTTGGTCGCGCGACt < 1:2498334/40‑1 (MQ=37) ccGGAGATACAGACGACGGTATCGTTTGGTCGCGCGACTtg < 1:3965343/41‑1 (MQ=38) ggAGAAACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTa < 1:3291936/59‑1 (MQ=255) agaTACAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCAAAACCCATg > 1:3030120/1‑71 (MQ=39) tACAGACGACGGTATCGTTTGGTCGCGCGACTtgc < 1:3210038/35‑2 (MQ=37) aCAGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGtgc > 1:1484733/1‑69 (MQ=38) aGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCaa < 1:4585561/39‑1 (MQ=37) aGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAAccgac < 1:5143910/44‑3 (MQ=37) aGACGACGGTATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCa < 1:2509921/56‑1 (MQ=255) tATCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCgg < 1:13402/40‑1 (MQ=38) tCGTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGTGCCTAAGCCGCTgg < 1:3369233/70‑1 (MQ=38) tCGTTCGGTCGCGCGACTTGTGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGTACCTAAACCGCTgg < 1:1044510/70‑1 (MQ=255) gTTTGGTCGCGCGACTTGCGCGCCAACCGCTGCCGGTAAACCAAAACCCATGGTGCCTAAGCCGCTGGAgg > 1:2708162/1‑71 (MQ=38) gTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCaaaa > 1:2640317/1‑46 (MQ=39) gTTTGGTCGCGCGACTTGCGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGTGCCTAAGCCGCTGGAg < 1:1227030/70‑1 (MQ=38) | CCCAGCTCTTGCACATTCATCATGAAAGAGCCGTCACCGGAGATACAGACAACGGTATCGTTCGGTCGCGCGACTTGTGCGCCAACCGCCGCCGGTAAACCAAAACCCATGGTACCTAAACCGCTGGAGG > AP009048/3684808‑3684937 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |