| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | USA300TCH1516_ALE | 37,558 | G→A | intergenic (+582/‑255) | USA300HOU_RS14915 → / → ugpQ_1 | hypothetical protein/Glycerophosphodiester phosphodiesterase, cytoplasmic |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 37,558 | 0 | G | A | 93.1% | 47.5 / ‑3.5 | 29 | intergenic (+582/‑255) | USA300HOU_RS14915/ugpQ_1 | hypothetical protein/Glycerophosphodiester phosphodiesterase, cytoplasmic |
| Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (11/16); total (13/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.92e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01 | |||||||||||
ATGCAAAATTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTAAGTTAAAAGCAATCCTAAGTAAAATTGCAGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCATAA > USA300TCH1516_ALE/37419‑37688 | aTGCAAAATTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAggg > 1:341420/1‑141 (MQ=255)aTGCAAAATTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGAAGTTCTCAGTAAAATTGCTGATAAggg > 1:344558/1‑141 (MQ=255) aaaTTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAACAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTaca > 1:379172/1‑139 (MQ=255) aaTTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTacag > 2:143501/1‑138 (MQ=255) gATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGAttt < 2:218573/141‑1 (MQ=11) caaGTGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTGGATTACAAAATGGAGCTTGGGACATAAATGAttttt < 1:296984/140‑1 (MQ=14) agagGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTa < 2:822566/141‑1 (MQ=14) agagGTAAGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTa < 2:822641/141‑1 (MQ=11) agGTACGTTAAAAGCAGTTCAAGGCAAAGTTGCAGATAAGAGGTGCGTTAAAAGGAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTaaa > 1:640516/1‑141 (MQ=14) gcaacgctAACAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTaaaa < 1:963906/134‑1 (MQ=14) acGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAatg > 1:621201/2‑141 (MQ=17) cGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAatga > 1:907092/1‑141 (MQ=14) ttAAAAGCAGTTCTAGGCAACGTTGGCGCTAAGAGGTGCGTTAAAAGCAGCTCGCAGCAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAACTTGGGACATAAATGATTTTTTAAAAatgaga > 1:391889/1‑141 (MQ=12) ttAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAACATGGAGCTTGGGACATAAATGATTTTTTAAAAatgaga < 2:308637/141‑1 (MQ=17) aaaGCAGTTCTAGGCAAAATTGCAGATAAGAGCTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTATGCTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAatgagatga < 1:660832/141‑1 (MQ=14) aaGTAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAatgagatgag < 1:198634/141‑1 (MQ=17) aaGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAatgagatgag < 1:388163/141‑1 (MQ=18) tctaggcaaaaTTGCAGATAAGAGGTGCGTTAAAAGTAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGa < 1:252364/134‑1 (MQ=18) ctaggcaaaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGAt > 2:756856/7‑141 (MQ=18) taggcaaaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGAtt < 1:851400/136‑1 (MQ=18) aggcaaaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTa < 1:879526/137‑1 (MQ=21) aggcaaaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTa > 2:293827/5‑141 (MQ=18) aggcaaaaTTGCAGATAAGAGGTGCGATAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTGTGGACATAAATGATTTGTTAAAAATGAGATAAGACGAAGATTa > 2:263529/5‑141 (MQ=17) gcaaaaTTGCAGATAAGAGGAGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAAc < 1:991022/139‑1 (MQ=21) gacaaaCTGCAGATAAGAGGGGAGGTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAAc < 2:153299/138‑1 (MQ=25) aaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTcc > 1:251116/1‑141 (MQ=21) aaaTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTcc > 1:646615/1‑141 (MQ=21) aTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCAt < 2:707206/141‑1 (MQ=25) aTTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCAt > 1:379725/1‑141 (MQ=21) tgcagATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCCTaa > 2:83769/5‑141 (MQ=21) tgcagATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCATaa < 1:97304/137‑1 (MQ=25) | ATGCAAAATTGCTGATAAGGGGTAAGTTAAAAGCAGTTCTCAGTAAAATTGCAGATAAGAGGTACGTTAAAAGCAGTTCTAGGCAAAATTGCAGATAAGAGGTGCGTTAAAAGCAGTTCTCAGTAAAATTGCTGATAAGGGGTAAGTTAAAAGCAATCCTAAGTAAAATTGCAGATAAGGGGTACAGAAAAACTAGACTTGATTACAAAATGGAGCTTGGGACATAAATGATTTTTTAAAAATGAGATGAGACGTAGATTAACTCCATAA > USA300TCH1516_ALE/37419‑37688 |
| 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: ‑ |