Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,929,796 | T→C | 38.2% | intergenic (‑879/‑289) | USA300HOU_1794 ← / → USA300HOU_1796 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,929,796 | 0 | T | C | 38.2% | 20.9 / 29.7 | 36 | intergenic (‑879/‑289) | USA300HOU_1794/USA300HOU_1796 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (12/10); new base C (4/10); total (16/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT > CP000730/1929658‑1929936 | ttGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATAAGTAAGATGTGACTGTGCTTGTATATTTGc < 1:331504/141‑1 (MQ=25) cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgat > 2:31861/1‑141 (MQ=255) taataaAATCGCTTACTCTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCCAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatga < 1:3476/141‑1 (MQ=255) ttcTTTCTTAAATTAACAATTTCTATTCTCATAGAATTTGTCCCATTAAGTGTAGACGTTTCGATCTATCAGCTAGAATCATATATCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTagag < 1:44145/139‑1 (MQ=37) tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc < 1:22435/141‑1 (MQ=255) aTTCTCATAGAATTTGACCAATGAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATTACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAtt < 1:171873/141‑1 (MQ=255) tctcATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGAAGCTATAATAGTAGAGTCTCCTATCTCTCAAGCGTCAATTTa > 1:44320/1‑141 (MQ=255) tAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGc < 2:254274/141‑1 (MQ=255) aTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag > 2:88590/1‑141 (MQ=255) ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGCTGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag > 2:182178/1‑141 (MQ=255) gTCCAATTAAGTGTTGACGTTTCAATCGATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga < 2:61848/141‑1 (MQ=255) gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga > 2:126965/1‑141 (MQ=255) cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgt < 2:109113/141‑1 (MQ=255) aTTAAGTGTAGACGATTCAATCTATCAGCTAGCATCATATAACTTATCAGAAGCAAGTGACTATGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgtgt > 1:239493/1‑141 (MQ=255) tAAGTGTAGATGATTCAATCTAACAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGta‑t < 2:218844/141‑1 (MQ=255) tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTAt > 2:56290/1‑141 (MQ=255) aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTGTAATTGCCTATTGTATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTGGACGCAGAGAGGGGGTGTA‑TAA‑GGTGTTGCTTATTTTGGTTCAcgtca > 1:175636/1‑141 (MQ=255) gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg > 2:244606/1‑141 (MQ=255) tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt > 2:184090/14‑141 (MQ=255) tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt > 2:197334/14‑141 (MQ=255) ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt < 1:33254/130‑1 (MQ=255) ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt < 1:33288/130‑1 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACAACATAGCACCAGTCATCAGTGGCtgtg < 2:206528/131‑1 (MQ=255) gaagtagatatgCTTGTATATCCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCCAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg < 2:80923/131‑1 (MQ=255) gaagtagatatgCTTGTACATTCCCCGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg < 1:371/131‑1 (MQ=255) aagtatgtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGATAGCAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc < 2:285818/132‑1 (MQ=255) aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc < 1:14010/132‑1 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg > 1:120736/8‑139 (MQ=255) tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCCa > 1:152367/1‑141 (MQ=255) gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt > 2:164834/1‑141 (MQ=255) gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg > 2:8652/5‑136 (MQ=255) gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg < 1:178286/137‑6 (MQ=255) ttgCTTGTATATCCACAGATGGTATAATAATAGAGTCGCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTTCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGt < 2:103232/140‑1 (MQ=255) cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt > 1:14277/1‑141 (MQ=255) ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta < 1:54762/141‑5 (MQ=255) atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc < 1:285647/141‑1 (MQ=255) tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt > 2:108555/1‑141 (MQ=255) | TTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT > CP000730/1929658‑1929936 |
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: ‑ |