| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 76,295 | T→C | 16.2% | intergenic (+353/‑781) | USA300HOU_0072 → / → USA300HOU_0075 | universal stress protein/lysophospholipase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 76,295 | 0 | T | C | 16.2% | 65.1 / 10.6 | 39 | intergenic (+353/‑781) | USA300HOU_0072/USA300HOU_0075 | universal stress protein/lysophospholipase |
| Reads supporting (aligned to +/- strand): ref base T (19/13); new base C (4/2); total (24/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01 | |||||||||||
AGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCT > CP000730/76158‑76420 | aGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTGACAAGGACGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTAt > 1:341056/1‑141 (MQ=25) ggTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGACTACTATCTAAGATATGCATTGAGTTATc > 1:30020/1‑141 (MQ=21) tttacAAAGTTGAATCTATAATAAAAATAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata < 2:423550/137‑1 (MQ=11) gTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata < 1:378014/141‑1 (MQ=35) cgttcAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata < 2:225883/137‑1 (MQ=25) cAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCTTTATAGATc > 2:393187/1‑141 (MQ=255) aaaGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCt > 2:436742/1‑141 (MQ=35) aaGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTg < 1:277130/141‑1 (MQ=35) aTCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATGTGAAATAtt > 1:391156/1‑141 (MQ=21) ttataatataaTTATAACAAAAAGGAGTCTACTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTaa < 2:293384/140‑1 (MQ=255) taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGACAAATTCAggg > 1:378573/1‑141 (MQ=12) gTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTAATCGTGATATATCTGAAATATTAAGGGaa > 1:98799/1‑141 (MQ=17) gTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGAATGTTATCACTGTAGCCGTTGGCTACTAACTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGaa > 1:98814/1‑141 (MQ=17) aaaaTTTTAACAAAAAGGAGTCTTCTGTAAGAACTAATACAAATAAAAACAATTAAACAAGGATGTTATCACTGTAGCCGTTGGCTACAATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACg < 2:8243/139‑1 (MQ=11) taATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACgt > 1:181296/1‑141 (MQ=17) aTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgg < 2:431278/141‑1 (MQ=32) tcttctGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAACCGTTCATCAtt > 1:448916/1‑141 (MQ=17) aCTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTGCTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTAcc < 1:258712/141‑1 (MQ=17) ttCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTgg < 1:436782/141‑1 (MQ=32) cAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATAAACATATATGTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGt > 1:107080/1‑141 (MQ=11) gATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAAATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCCACGATCTACCGTAGGGTTc > 1:206051/1‑141 (MQ=11) aaaCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGaa > 2:284445/1‑141 (MQ=32) aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTCCGATATGATTTGGAATATCGTGATATATCTGAAATATGAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTTTACCGTTGGCTTTAAGAATATGCACCaa > 1:271973/1‑141 (MQ=11) aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCTCCaa > 1:285521/1‑141 (MQ=17) aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCACCaa > 2:390706/1‑141 (MQ=17) cAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAAt > 2:323664/1‑141 (MQ=32) cAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAAt > 1:269993/1‑141 (MQ=32) aGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATAACTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACTGTCTACCGTTGGGTTCAAGAATATGCCCCAAttt > 1:25148/1‑141 (MQ=17) gTTATCGCTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTatat < 2:417626/141‑1 (MQ=25) tATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCa > 2:38302/1‑141 (MQ=35) tATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCa < 1:428579/141‑1 (MQ=35) gTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTgg < 1:331602/141‑1 (MQ=25) ccGTTGGCTACTATCTAAGATTTGCAATGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaaga < 1:241744/141‑1 (MQ=21) gTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaa > 1:365246/1‑141 (MQ=255) aacggCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaa > 1:463146/4‑141 (MQ=255) tGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaaaa > 1:34335/1‑141 (MQ=35) ggCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAg < 1:436742/141‑1 (MQ=35) gctaacaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGc < 1:15377/135‑1 (MQ=25) aaGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCACATTTGGAAGAAAAAGCATAAAAAGCt > 1:156301/1‑141 (MQ=38) | AGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCT > CP000730/76158‑76420 |
| 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: ‑ |