| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 512,643 | A→G | 48.4% | intergenic (+711/‑6052) | recR → / → speA_2 | Recombination protein RecR/Arginine decarboxylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 512,643 | 0 | A | G | 48.4% | ‑3.3 / 58.2 | 46 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
| Reads supporting (aligned to +/- strand): ref base A (19/5); new base G (7/15); total (26/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.50e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01 | |||||||||||
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAT > USA300TCH1516_ALE/512504‑512771 | ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa < 1:154799/141‑1 (MQ=255) aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAACTAAATTGAATCACTAACGAAAATTAATATTCAAAATTGATTGACTCACATCTAAACAAAAt > 2:270147/1‑141 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggaagtttttttgtatatttgtgcacaaaaaataaaatggtcaagattaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑Gtata < 1:292531‑M1/36‑1 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgcaaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑Gtata < 2:23721‑M1/36‑1 (MQ=255) tagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTAtaattaat < 1:64420‑M1/42‑1 (MQ=255) agcgtcgatgctagctgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggccaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATaattaatt < 2:258430‑M1/43‑1 (MQ=255) aacttacgttccgctagagtagaacgttgtcaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGaaaaa > 1:152326‑M1/82‑141 (MQ=255) cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAAt < 1:246841/141‑1 (MQ=255) gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt > 1:32336/1‑141 (MQ=255) tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg > 1:227049‑M1/73‑141 (MQ=255) tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg < 1:77413‑M1/69‑1 (MQ=255) ataatgtagtagaaagatgcaaggcaggatttgaataaaatattaattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGa > 2:129113‑M1/72‑141 (MQ=255) tttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGCAAAATGAACATTGaa > 2:102224/1‑141 (MQ=255) gctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAAGGAAAATTaaaa < 1:197882‑M1/72‑1 (MQ=255) tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGGAAGAAAAGTGAACATTGaaaa > 1:220955/1‑141 (MQ=255) atttGTTGAGGCAAACAAAACAACTCAACTTAAGAACTAACTTGACTTACTAACGACAACTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTCGTCTGTAAGAAAAATGAACATTGAAAAc > 2:78121/1‑141 (MQ=255) tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTAACATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg < 1:113030/141‑1 (MQ=255) aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAGAACTGAATGAc > 1:78442/1‑141 (MQ=255) ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAATTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGCTAAGAAAAATGAACATTGAAAACTGAATGAcc > 2:30153/1‑140 (MQ=255) gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGAATGACAat < 1:145561‑M1/85‑1 (MQ=255) aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa < 2:125275/141‑1 (MQ=255) acaacTCAACTTAAGAAATCACTTGCATTACTAACGCACATTAATTTTAAAAAGTTATTCACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTCTCTCTAAGAAAACTGAACATTGAAGACTGACTGACAATATGtaacc > 2:70559/1‑137 (MQ=255) acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaaa > 1:243540/1‑140 (MQ=255) acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTAACACAAAAGTTAAGAAAATACCATCA‑CAGTTTTTTTCTGTAAGAAAAATGAACATTGAAAACTGAATAACAAAATGACaaa > 1:56209/1‑140 (MQ=255) agttttgtaatcaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTGAAATGTTGATAAAAT‑GTATAATTAATGCTGGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt > 2:56730‑M1/47‑141 (MQ=255) aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAACAGAATGACAAGATGTcgaagt > 1:104148/1‑136 (MQ=255) aactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAATAATTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt < 2:17630/141‑1 (MQ=255) ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTc < 2:31698‑M1/101‑1 (MQ=255) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTTCaaaa < 1:213441‑M1/106‑1 (MQ=255) caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTATAAGTTAATTCCaaaa < 1:79803‑M1/106‑1 (MQ=255) tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt > 1:258430‑M1/31‑141 (MQ=255) ttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGGAAGAAAGATGGACATTGAAAACTGAATGACAATATGTCAACGTTAATTCGACAAAACGTa < 1:81653‑M1/112‑1 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:171421‑M1/114‑1 (MQ=255) tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt < 1:211328/141‑1 (MQ=255) taaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc < 1:227123‑M1/124‑1 (MQ=255) tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTAAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTaa > 2:170276/1‑140 (MQ=255) aattTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAATTTACAAACATTAttt > 1:21028/1‑141 (MQ=255) ataaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg > 2:142063‑M1/5‑141 (MQ=255) taaAAAAAGTTATTGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTaga < 1:226977‑M1/138‑2 (MQ=255) ttAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGACACAAGGCTGACGATTTGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa > 1:292030/1‑141 (MQ=25) aaAAAAAGTTATTGACTTAAATGTTGGTAAAAT‑GTATAATTAATTCATGTCGGTAAGATACAAAAGCGCTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTa < 2:19078‑M1/139‑1 (MQ=255) tAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt > 1:285722/1‑141 (MQ=255) aaaaGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt > 2:220259/1‑141 (MQ=255) aaaaGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAttt > 2:155498/1‑141 (MQ=255) ttGACTTAAATGTTGATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATATCAACGTTAATTCCAAAAAACGTAAATATAAGTTACAAACATTATTTAGTATTTATGAGCTa > 1:4173/1‑141 (MQ=255) ttGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATATATGAGCTa > 2:289618/1‑141 (MQ=255) gACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt > 1:237719/1‑141 (MQ=255) | GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAT‑GTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAT > USA300TCH1516_ALE/512504‑512771 |
| 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: ‑ |