| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 2,732,937 | C→G | 9.3% | R311R (CGC→CGG) | USA300HOU_2590 → | possible S15 family peptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,732,937 | 0 | C | G | 9.3% | 133.0 / 6.7 | 43 | R311R (CGC→CGG) | USA300HOU_2590 | possible S15 family peptidase |
| Reads supporting (aligned to +/- strand): ref base C (18/21); new base G (2/2); total (20/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01 | |||||||||||
TGCTAGTTGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTACCTAACGC > CP000730/2732799‑2733074 | tgCTAGTTGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGcaa > 1:61182/1‑139 (MQ=255) aGTTGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCaaaaa > 1:279931/1‑141 (MQ=255) tGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGGCAAAAATCa > 2:119921/1‑141 (MQ=255) tGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGGCAAAAATCa > 2:119947/1‑141 (MQ=255) tACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCaaagaa < 2:71616/141‑6 (MQ=255) tACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCaaagaa < 2:257127/141‑6 (MQ=255) cacaAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTg > 1:64035/1‑141 (MQ=255) aaGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCAAAAAAAAAAAAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGAtt < 2:13844/141‑1 (MQ=255) ttGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTAc > 2:331874/1‑141 (MQ=255) tGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTAcc > 1:139042/1‑141 (MQ=255) tGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTAcc > 1:67060/1‑141 (MQ=255) tGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTCACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTAcc < 1:256300/141‑1 (MQ=255) ccGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAaaga < 2:131852/141‑1 (MQ=255) ggCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAgaagaa > 2:2249/1‑141 (MQ=255) gCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGaaa < 1:291876/141‑1 (MQ=255) tttGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAAc < 2:89775/141‑1 (MQ=255) tGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGGCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGa < 1:119947/141‑1 (MQ=255) tGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGGCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGa < 1:119921/141‑1 (MQ=255) aaGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGAtt > 1:51085/1‑141 (MQ=255) ttAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAg < 2:279931/141‑1 (MQ=255) cAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACg > 1:96898/1‑141 (MQ=255) gCATCTGAAGAAAAATGGCTATATGTACATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCAt > 2:125676/1‑141 (MQ=255) cTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCa < 2:210944/141‑1 (MQ=255) aaaTGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGaa > 2:50321/1‑141 (MQ=255) aaaTGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGaa > 2:25455/1‑141 (MQ=255) tatGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGAt > 1:66452/1‑141 (MQ=255) gtgCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCaa < 2:171985/141‑1 (MQ=255) tgCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAAt > 1:325779/1‑141 (MQ=255) cATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAAttt < 1:222577/141‑1 (MQ=255) cATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTGGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAATTAGAGATCAAttt < 2:215680/141‑1 (MQ=255) ggACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTtat < 1:114827/141‑1 (MQ=255) aaaGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGc < 1:193937/141‑1 (MQ=255) aaaGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGc > 2:51942/1‑141 (MQ=255) aaaGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGc < 1:95159/141‑1 (MQ=255) gggAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCa < 2:86843/141‑1 (MQ=255) ttACTACGCTAGAGAAAATCTCGAAACCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAAtcagc < 2:276570/141‑1 (MQ=255) actacGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAAtcagcgt < 1:50321/141‑1 (MQ=255) acGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAAtcagcgtcag > 1:175561/1‑141 (MQ=255) aaaaTCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCtt < 1:55445/141‑1 (MQ=255) aaaTCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCttt > 1:242628/1‑141 (MQ=255) tctcGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTAcc > 1:8186/1‑141 (MQ=255) tctcGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTAcc < 2:87170/141‑1 (MQ=255) ctcGAACGCAAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTACCt < 2:61182/141‑1 (MQ=255) aCGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTACCTAACGc > 2:157212/1‑141 (MQ=255) | TGCTAGTTGGTCTACACAAGGTTTGCACAACCGTGGCTCTTTTGAAGGATTTAAACAAGCTGCATCTGAAGAAAAATGGCTATATGTGCATGGACGTAAAGAGTGGGAAAGTTACTACGCTAGAGAAAATCTCGAACGCCAAAAATCATTCTTTGATTTTTACCTTAAAGAAGAAAATAACGATTGGAAAGATACGCCTCATGTCATTTATGAAGTTAGAGATCAATTTTATAAAGGCGAATTCAAATCAGCGTCAGCGTTCCCTTTACCTAACGC > CP000730/2732799‑2733074 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |