Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,929,804 | A→G | 43.2% | intergenic (‑887/‑281) | 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,804 | 0 | A | G | 43.2% | 5.7 / 35.4 | 35 | intergenic (‑887/‑281) | USA300HOU_1794/USA300HOU_1796 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (11/8); new base G (5/10); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01 |
TAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT > CP000730/1929664‑1929944 | 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‑GG‑TGATGCTTAt > 2:56290/1‑141 (MQ=255) aGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTGTAATTGCCTATTGTATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTGGACGCAGAGAGGGGGTGTA‑TAA‑GG‑TGTTGCTTATTTTGGTTCAcgtca > 1:175636/1‑141 (MQ=255) gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg > 2:244606/1‑141 (MQ=255) tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt > 2:197334/14‑141 (MQ=255) tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt > 2:184090/14‑141 (MQ=255) ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt < 1:33254/130‑1 (MQ=255) ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt < 1:33288/130‑1 (MQ=255) gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACAACATAGCACCAGTCATCAGTGGCtgtg < 2:206528/131‑1 (MQ=255) gaagtagatatgCTTGTATATCCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCCAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg < 2:80923/131‑1 (MQ=255) gaagtagatatgCTTGTACATTCCCCGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg < 1:371/131‑1 (MQ=255) aagtatgtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGATAGCAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc < 2:285818/132‑1 (MQ=255) aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc < 1:14010/132‑1 (MQ=255) gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg > 1:120736/8‑139 (MQ=255) tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCCa > 1:152367/1‑141 (MQ=255) gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt > 2:164834/1‑141 (MQ=255) gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg > 2:8652/5‑136 (MQ=255) gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg < 1:178286/137‑6 (MQ=255) ttgCTTGTATATCCACAGATGGTATAATAATAGAGTCGCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTTCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGt < 2:103232/140‑1 (MQ=255) cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt > 1:14277/1‑141 (MQ=255) ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta < 1:54762/141‑5 (MQ=255) atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc < 1:285647/141‑1 (MQ=255) tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt > 1:122592/9‑141 (MQ=255) tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt > 1:91212/6‑141 (MQ=255) tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt > 2:108555/1‑141 (MQ=255) cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg > 1:152294/5‑141 (MQ=255) cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg > 2:42977/5‑141 (MQ=255) gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGAAACTTATTTTTGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTg > 1:77367/1‑141 (MQ=255) cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc > 2:329878/3‑141 (MQ=255) cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc < 2:139902/139‑1 (MQ=255) agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt > 1:40252/2‑141 (MQ=255) atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt > 1:289981/1‑141 (MQ=255) | TAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT > CP000730/1929664‑1929944 |
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: ‑ |