Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 1,087,670 | A→T | 82.6% | L5I (TTA→ATA) | ykoE ← | Putative HMP/thiamine permease protein YkoE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 1,087,670 | 0 | A | T | 82.6% | 48.0 / 6.7 | 23 | L5I (TTA→ATA) | ykoE | Putative HMP/thiamine permease protein YkoE |
Reads supporting (aligned to +/- strand): ref base A (2/2); new base T (12/7); total (14/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01 |
AAACCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCTAGATCAGGTT > USA300TCH1516_ALE/1087538‑1087805 | aaaCCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTg > 2:118848/1‑141 (MQ=255) aaaCCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTg < 1:266790/141‑1 (MQ=255) aCCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCCCACCATTGTAATCAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGAc > 1:263839/1‑141 (MQ=255) aCCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACCCAATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGAc > 1:263863/1‑141 (MQ=255) gTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACAtttt > 2:165969/1‑141 (MQ=255) tGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATaaaaa > 2:118279/1‑141 (MQ=255) aaGGCCTGTAGCTTGCACACCATTGTCAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGACATTTTTACATCCTCCTAATaaaaaaa < 1:169307/141‑1 (MQ=255) cATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACAATCCATAGGaa > 2:178453/1‑141 (MQ=255) aaacaaacCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTAt < 1:118848/141‑1 (MQ=255) aacaaacCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATg < 1:118279/141‑1 (MQ=255) ccaTAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACaa < 1:221653/141‑1 (MQ=255) ccaTAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACaa > 2:169947/1‑141 (MQ=255) tAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGCAATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCAATAGGAAAGTTATGCGTTCACAATATATATTAGTaaa > 2:12439/1‑141 (MQ=255) cTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACa < 1:152873/141‑1 (MQ=255) aTACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACCTCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGATATGCGTTCACAATATATATTAGTAAAACATATGta > 2:218769/1‑141 (MQ=255) aTACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCCCAAAATATATTAGTAAAACattagtt > 2:103769/1‑136 (MQ=255) tACTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAg < 1:165969/141‑1 (MQ=255) aCTAGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGt > 1:197414/1‑141 (MQ=255) caGTATTTCAGATAGCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTc < 1:261812/140‑1 (MQ=255) tttCAGATAGCTTTATACCTTTCGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCt < 1:85705/141‑1 (MQ=255) gCTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCTAGATCAGGt > 2:36165/1‑141 (MQ=255) cTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATGGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCTAGATCAGGtt > 1:233379/1‑141 (MQ=255) cTTTATACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAAAATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCTAGAGTAGGtt > 1:233350/1‑141 (MQ=255) | AAACCCCGTTTGTTAATTGTTCGAGGTGAAGGCCTGTAGCTTGCACACCATTGTAAACAAACCACCATAAATTGTAAATAACTGCGAATACTACTGAAATAAGTACGGTTACTAGTATTTCAGATAGCTTTAAACCTTTTGACATTTTTACATCCTCCTAATAAAAAAACGCACAACCATCCATAGGAAAGTTATGCGTTCACAATATATATTAGTAAAACATATGTATAGTAACACTTTCCTACGCTAGTTCAAGCTAGATCAGGTT > USA300TCH1516_ALE/1087538‑1087805 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |