| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 1,306,605 | G→T | 33.3% | G223V (GGA→GTA) | phaB → | Acetoacetyl‑CoA reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 1,306,605 | 0 | G | T | 33.3% | 16.8 / 5.3 | 12 | G223V (GGA→GTA) | phaB | Acetoacetyl‑CoA reductase |
| Reads supporting (aligned to +/- strand): ref base G (6/2); new base T (2/2); total (8/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.47e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATGGCAAGTGAGTGGCAAGAAGATGAACTTCAAGCAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGTTAATTAAACAAGGGGTGAGGT > USA300TCH1516_ALE/1306472‑1306741 | atGGCAAGTGAGTGGCAAGAAGATGAACTTCAAGCAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACtata < 2:48980/141‑1 (MQ=255) cAAGAAGATGAACTTCAAGCAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGt > 2:124204/1‑141 (MQ=255) gCAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAat > 2:98775/1‑141 (MQ=255) cAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGTAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAata > 2:99896/1‑141 (MQ=255) cAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGTAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAata > 2:99904/1‑141 (MQ=255) gCGATTGGTTTTACCGAGTGAGGTTGCTCATACAAGCGCCTATTTATATCACCCAAAAGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAAAGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGt > 2:7922/1‑141 (MQ=255) aGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGtt > 2:118417/1‑141 (MQ=255) aTGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGttaattaa > 1:78202/1‑141 (MQ=255) aTGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGttaattaa > 1:78210/1‑141 (MQ=255) aCTGTAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGTTAATTAAACAAGGGGTGAGGt < 1:99896/141‑1 (MQ=255) aCTGTAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGTTAATTAAACAAGGGGTGAGGt < 1:99904/141‑1 (MQ=255) aCTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGTTAATTAAACAAGGGGTGAGGt < 1:96272/141‑1 (MQ=255) | ATGGCAAGTGAGTGGCAAGAAGATGAACTTCAAGCAATGATAACTGAATTACCACAACAGCGATTGGTTTTACCGAGTGAGGTTGCTCATACATGCGCCTATTTATATCACCCAAATGCTAGAAGTGTTACTGGAACTATACAGAAAGTTAATGGTGCTTGGTATATTTAAAATATAGAACTATTTTTATGAGCAAATGTTTTTATGATTGATGTTATAATATGAATATTACAATGAATCATAATTTGTTAATTAAACAAGGGGTGAGGT > USA300TCH1516_ALE/1306472‑1306741 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |