| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 1,307,499 | G→T | 16.9% | K251N (AAG→AAT) | USA300HOU_RS06455 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 1,307,499 | 0 | G | T | 16.9% | 37.3 / 4.4 | 24 | K251N (AAG→AAT) | USA300HOU_RS06455 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (6/14); new base T (2/2); total (8/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.78e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-02 | |||||||||||
ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA > USA300TCH1516_ALE/1307362‑1307637 | aTGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATgaagaa > 2:23506/1‑140 (MQ=255) gTTACCATTAAACAAAGTCGTTCAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCTAATTGAGAGAAAATGAAAATCAAGAAATACATTACCA‑AGAAATGAAGAAacaga < 2:9500/140‑1 (MQ=255) gTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACCAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAACAAATACATTACGA‑AGGAATGAAGAAtcaga < 2:176756/140‑1 (MQ=255) ccATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAATTGAAGAAACAGACAGt < 1:145141/140‑1 (MQ=255) cATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTc > 2:85786/1‑140 (MQ=255) caaAGTCGTTGAACATCCAGATTTCTTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTa < 2:67182/140‑1 (MQ=255) cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg > 2:3555/1‑140 (MQ=255) cAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTg > 2:236233/1‑140 (MQ=255) atttatttGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGga > 1:203144/1‑140 (MQ=255) aaGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAACTCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAATCCAAAGTAAGAATAATTTAGATTTTGGAGATCAGTTtggtggtg > 1:111810/1‑140 (MQ=255) tGTACAATGGAAGATTTTGAGTATATTATCGAGTTGAGAGAAAATGAAAATCAAGAAATTCATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTGAGATTTTGGAGATCCGTTTGGTGGtggttttgg < 1:105824/140‑1 (MQ=255) tGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAgtgt < 2:211742/140‑1 (MQ=255) attatCGAATGCAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAaggagg < 2:163341/140‑1 (MQ=255) tatCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCt < 1:34244/140‑1 (MQ=255) atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa < 2:284409/140‑1 (MQ=255) atCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGGCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTa < 1:23506/140‑1 (MQ=255) cGAATTGAGAGAAAATGAAAATAAAGAAATAAAATACGATATGATTAAAGAAA‑TGACAGAAGGAAATAATTATAACATAAATGTGGGATATCCGTTGGGTGATGGTTTTTGTTTCTTCGAGTAAGTGTAAGGAGGCTATg < 2:68560/140‑1 (MQ=255) gaaaatgaaaatCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGAAAGAATAATTTAGATTTTGGAGATCCGGTTGGTGGAGGTGTTGGTTTCTTCGAGCAAGTGTAAGGAGGCAATGAATAGAacac > 2:275712/1‑140 (MQ=255) cAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTCGGTGGTGGTTTTGGTTTCATCGAGTATGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAGgc > 2:158912/1‑140 (MQ=255) aTACATTACAA‑AGAAAAGAAGAAACAGACAGTCCAAAGTAAGAATAATGTAGATATTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAGGTGTAAGAAGGCTATGAATTGAAAACGCTCGGTGAAGCGCTaaa > 2:61812/1‑140 (MQ=255) cagtACGA‑AGAAATGAAGAAACAGACTGTCCAAAGTAATAATGATTTAGATGTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTAGAGTAAGTGTAAGGAGGCTATGTATTGAAAACGGTCGGTGAAGCGCTAAAAgg < 2:177005/137‑1 (MQ=255) aTTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGGGTAATTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGt < 1:184029/140‑1 (MQ=255) aTTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa < 1:3555/140‑1 (MQ=255) aTTACGA‑AGAAATG‑ATAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTa < 1:236233/140‑1 (MQ=255) aaTGAAGAAACCGCCAGTCCAAAGTAATAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAGCGGTCGGTGAAGCGCTAAAAGGTAGACGTGaaa < 1:158912/140‑1 (MQ=255) agaaACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTa > 1:49128/1‑140 (MQ=255) | ATGACGTTACCATCAATCAAAGTCGTTGAACATCCAGATTTATTTGTTGAAACAAGTACTTGTACAATGGAAGATTTTGACTATATTATCGAATTGAGAGAAAATGAAAATCAAGAAATACATTACGA‑AGAAATGAAGAAACAGACAGTCCAAAGTAAGAATAATTTAGATTTTGGAGATCCGTTTGGTGGTGGTTTTGGTTTCTTCGAGTAAGTGTAAGGAGGCTATGAATTGAAAACGGTCGGTGAAGCGCTAAAAGGTAGACGTGAAAGGTTA > USA300TCH1516_ALE/1307362‑1307637 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |