Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 512,643 | A→G | 44.1% | 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 | 44.1% | 16.0 / 58.0 | 52 | intergenic (+711/‑6052) | recR/speA_2 | Recombination protein RecR/Arginine decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (14/15); new base G (15/8); total (29/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
AATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA > USA300TCH1516_ALE/512509‑512773 | aaTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAACACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAAATATTAAATTAAATGTTAATcaaa > 2:222585/1‑141 (MQ=255) aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt > 1:186814/1‑141 (MQ=255) tttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGCATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGta < 1:307599/141‑1 (MQ=255) ttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat > 1:244275/1‑141 (MQ=255) gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata < 1:154847‑M1/36‑1 (MQ=255) ttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa > 1:87048/1‑141 (MQ=255) gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt < 1:256610‑M1/46‑1 (MQ=255) tcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTg < 2:64549‑M1/47‑1 (MQ=255) aGCAGTTTGAAAGCCTGTATAGTAATTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg < 1:175461/141‑1 (MQ=255) nttGAAAGCCTGTAGAGTATTTATTTGTTGATGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt < 1:228233/140‑1 (MQ=255) gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg < 1:269830/141‑1 (MQ=255) gtcaaaattacgttacgctataggagaacgaagacaggcagtttttaaatcaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa > 2:269781‑M1/86‑141 (MQ=255) aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa < 1:155826/141‑1 (MQ=255) aGCCTGTATAGTATATATTCGTAGAAGCTAACAAAACAACTAAAGATAGGAAAAGACTTCCAAGAGTGATGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTATCGGAAAGcaa > 2:156226/1‑141 (MQ=255) ataGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaaa > 2:130242/1‑140 (MQ=255) cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa > 1:2668‑M1/76‑141 (MQ=255) tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg < 1:2518/141‑1 (MQ=255) ttatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaa < 2:231636/141‑1 (MQ=255) tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa < 1:169830/141‑1 (MQ=255) ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATCAATTCTTGTCGGTAAGAAAAAAGAACATTCACAACTg > 2:70567/1‑141 (MQ=255) gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa > 1:305784/1‑141 (MQ=255) gtagaatgttgccaggcagttttttaatcaaattgtggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa > 2:310755‑M1/64‑141 (MQ=255) cgttgaaaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa > 1:226305‑M1/58‑141 (MQ=255) aaacaaaacaACTCAACTTAAAAAATAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAata < 2:123708/141‑1 (MQ=255) aaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt < 2:44029/140‑1 (MQ=255) atgcagttttttattagaattttagttaaaaaataatatctacaatataaAAAAAGTTATTGACTTAAAGGTTGATAAAAGGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATGGAAAACTGAATGACAATATGTCa > 1:42548‑M1/51‑141 (MQ=255) aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa > 1:334277/1‑141 (MQ=255) taaaattAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa < 2:219611/136‑1 (MQ=255) atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa > 1:50693‑M1/38‑141 (MQ=255) aacaaattttggataaaaaataaaataaacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAACGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAACACTGAATGACAATATGTCAACGTTAATTCCaa > 2:302090‑M1/38‑141 (MQ=255) aaaaaattttagtaaaaaaaaaaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa > 2:263230‑M1/38‑141 (MQ=255) tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTCATTCCAAAAAACGt > 2:43165‑M1/31‑141 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:310783‑M1/114‑1 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:261837‑M1/114‑1 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:212639‑M1/114‑1 (MQ=255) ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATAAATTGTTGTCGGTAAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 2:105142‑M1/114‑1 (MQ=255) cTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc < 1:169581/141‑1 (MQ=255) aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTTGTTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt > 1:288530/1‑141 (MQ=255) aaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa < 1:125044/141‑1 (MQ=255) aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa > 1:82595‑M1/16‑141 (MQ=255) aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa > 1:64549‑M1/16‑141 (MQ=255) aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa < 1:290225‑M1/126‑1 (MQ=255) aattaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt > 1:278841/1‑141 (MQ=255) aattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt < 1:336264/141‑1 (MQ=255) ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTTATTGATGTCTGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg > 2:154847‑M1/5‑141 (MQ=255) ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg > 2:303541‑M1/5‑141 (MQ=255) aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa > 1:150866/1‑141 (MQ=255) aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTa > 1:332666/1‑141 (MQ=255) gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGCAGTTAAAGGGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa > 1:324772/1‑141 (MQ=18) ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATAAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAATTACAAAAATAATTTAGTATTTATGaa > 2:119956/1‑140 (MQ=255) ttGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa > 2:96717/1‑141 (MQ=255) cTTAAATGTTAATAAACTGCCGAACTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCa < 1:213955/141‑1 (MQ=17) | AATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCA > USA300TCH1516_ALE/512509‑512773 |
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: ‑ |