| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | USA300TCH1516_ALE | 2,372,753 | +T | coding (1023/1335 nt) | pbuG ← | Guanine/hypoxanthine permease PbuG |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 2,372,752 | 1 | . | T | 95.7% | 62.4 / ‑3.4 | 23 | coding (1024/1335 nt) | pbuG | Guanine/hypoxanthine permease PbuG |
| Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (10/12); total (11/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.78e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAA‑TAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGA > USA300TCH1516_ALE/2372616‑2372887 | gCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAA‑ttag > 2:63966/1‑138 (MQ=255) cTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAACGAACAGCGCTAA‑ttagg > 1:83235/1‑137 (MQ=255) ggCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAACGTTTGCAGCCATTACTCCACCTACCCCCACTAAGGCTGGTCTCGTAACTCCACTTGTTACCACAGCGATTAACAGGCTAAAGAACAGGGCTAA‑TTGGgaa > 2:225758/1‑141 (MQ=255) cAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTAGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAGGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGc > 1:158185/1‑141 (MQ=255) ctttATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCaa < 2:222100/140‑1 (MQ=255) ctcAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTg > 2:11041/1‑141 (MQ=255) aGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTc > 2:212786/1‑141 (MQ=255) aGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTc > 2:212731/1‑141 (MQ=255) gCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCt < 1:170189/141‑1 (MQ=255) gCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCt < 1:182123/141‑1 (MQ=255) gttGCAGCCATTAATACAGATAGTACCACTAAGGGTGGGGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACCGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCAcc < 2:126498/140‑1 (MQ=255) cATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCaa < 2:97310/141‑1 (MQ=255) tACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTCGCACCTCCAGCGACACCTg > 1:28665/1‑141 (MQ=255) tACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTg > 2:195943/1‑141 (MQ=255) cTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATaagaa < 1:48264/141‑1 (MQ=255) tGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATaagaag < 2:40250/141‑1 (MQ=255) tgtCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAgttg < 1:197190/141‑1 (MQ=255) tCGGAACTGCAGTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAGCAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAgttgtt < 1:224256/141‑1 (MQ=255) cTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTCTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAgtt < 1:120490/141‑1 (MQ=255) ttAACGGACTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACa > 2:113073/1‑141 (MQ=255) aCTAAAGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATACTACAATCGTAg > 1:179794/1‑141 (MQ=255) aaGAACAGCGCTAATTAGGAAACAACAACCTGTTACGATCCTCGCAAATCCTGATTTAGACCATGCCCGAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTaa > 1:113231/1‑141 (MQ=255) aaNAGCGCTAATTAGGAAACAACACACTCTGGTCTTGCGATCAAATCCTGTTATGGGACCGACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGa < 1:150124/141‑1 (MQ=255) aaCAGCGCTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGa < 1:220558/141‑1 (MQ=255) gcTAATTAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTga > 1:167107/1‑141 (MQ=255) | GCTGGCACTGCAACTTCAAAACTTTTCCAGTTTATCTCAGCAAAGTTTGCAGCCATTAATACACCTACCACCACTAAGGCTGGTGTCGTAACTGCACTTGTTACAACAGCGATTAACGGACTAAAGAACAGCGCTAA‑TAGGAAACAACAACCTGTTACGATACTAGCAAATCCTGTTCTAGCACCTACAGCAACACCTGAAGTTGATTCGATATAAGAAGTTGTTGTAGTTGTTCCAAAAACTGATCCTACAATCGTAGCTAAAGAATCTGA > USA300TCH1516_ALE/2372616‑2372887 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |