Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,460,315 | T→C | 13.5% | intergenic (‑347/+53) | USA300HOU_1356 ← / ← ctpA | possible acetyltransferase/S41A family carboxy‑terminal peptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,460,315 | 0 | T | C | 13.5% | 67.0 / 5.9 | 30 | intergenic (‑347/+53) | USA300HOU_1356/ctpA | possible acetyltransferase/S41A family carboxy‑terminal peptidase |
Reads supporting (aligned to +/- strand): ref base T (13/12); new base C (2/2); total (16/14) | |||||||||||
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.20e-01 |
TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA > CP000730/1460168‑1460461 | tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAGTTTa > 2:30537/1‑149 (MQ=255) tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTAGCTTTGCAAATTTa > 1:290512/1‑149 (MQ=255) gCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTag < 1:162944/149‑1 (MQ=255) cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg > 1:221152/1‑149 (MQ=255) cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg > 2:201034/1‑149 (MQ=255) tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg > 1:9423/1‑149 (MQ=255) tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg > 1:168583/1‑149 (MQ=255) atatTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCgtgt > 1:147127/1‑149 (MQ=255) aTGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACata < 1:235409/149‑1 (MQ=255) aaTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataata > 1:6746/1‑149 (MQ=255) aTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATAATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataa > 1:85055/1‑149 (MQ=255) gTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataataa > 2:48555/1‑149 (MQ=255) cGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAttt > 2:66673/1‑149 (MQ=255) gatgatAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAtttt < 1:56179/149‑1 (MQ=255) aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAgtgt < 2:9423/149‑1 (MQ=255) aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAAATTTAGTAgtgt < 2:168583/149‑1 (MQ=255) ttACATTATTTGTAAATGTCTAATTTATATCTTTTAGAGTTATACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTAtt < 2:288991/149‑1 (MQ=255) tcatTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTaaa > 1:27860/1‑149 (MQ=255) aTGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCa < 1:201034/149‑1 (MQ=255) ttATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCga < 1:32911/149‑1 (MQ=255) atCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAc > 2:216511/1‑149 (MQ=255) gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat > 2:239611/1‑149 (MQ=255) gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat > 1:187033/1‑149 (MQ=255) aGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTtta < 2:145290/149‑1 (MQ=255) tCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGc < 2:211313/149‑1 (MQ=255) taaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTGAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCtttt > 1:157093/1‑149 (MQ=255) aaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCttttt < 2:2252/149‑1 (MQ=255) cctAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGt < 2:159953/149‑1 (MQ=255) aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa < 2:264622/149‑1 (MQ=255) aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATAGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa < 1:108746/149‑1 (MQ=255) tttAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTtta < 2:6746/149‑1 (MQ=255) | TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA > CP000730/1460168‑1460461 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |