Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 253,973 | C→T | 11.8% | intergenic (+186/‑2) | USA300HOU_0226 → / → USA300HOU_0227 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 253,973 | 0 | C | T | 11.8% | 84.6 / 4.1 | 34 | intergenic (+186/‑2) | USA300HOU_0226/USA300HOU_0227 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (16/14); new base T (2/2); total (18/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01 |
TTGAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGA > CP000730/253840‑254109 | ttGAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAATAATCCTCACCATTGaaa > 1:163675/1‑141 (MQ=255) aaTAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAaga > 1:174666/1‑141 (MQ=255) aaaTGCATATGTCGCCAAGCCATTGACAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta < 2:262168/141‑1 (MQ=255) aaTGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtat > 2:219221/1‑141 (MQ=255) aTGCATATGTCGCCAAGCCATTGCAAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATAACCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata < 2:181163/141‑1 (MQ=255) tGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat > 1:58500/1‑141 (MQ=255) gCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata < 2:60479/141‑1 (MQ=255) cGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc < 2:327523/141‑1 (MQ=255) aTTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTCTACTTTATATAAGTATACAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGGGAg < 2:285329/141‑1 (MQ=255) aaaaTTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtggt < 1:277153/131‑4 (MQ=38) aaaaTTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtggt > 2:277153/1‑128 (MQ=38) acacCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATAATAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa > 2:94732/1‑141 (MQ=255) cacCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaa > 1:145967/1‑141 (MQ=255) acCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaa < 1:51780/141‑1 (MQ=255) aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGaaaa < 1:11942/87‑1 (MQ=255) aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGaaaa > 2:11942/1‑87 (MQ=255) cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTca > 1:227690/1‑141 (MQ=255) aagaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAATCCTGCTGTCACAATGGTCATa > 1:234249/1‑141 (MQ=255) gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc > 1:190017/1‑141 (MQ=255) tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa < 1:171628/141‑1 (MQ=255) tataATTATATATCGATGGTCATACCCATTCACAGTAACAATCCTCACCATGGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa < 1:135577/141‑1 (MQ=255) tataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAACAAAAATCCTGTTGTCACAATGGTCATAGACACGACATACTTTa < 1:158121/141‑1 (MQ=255) tataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa < 1:123859/141‑1 (MQ=255) aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg > 1:93079/1‑141 (MQ=255) gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt < 2:203312/141‑1 (MQ=255) gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa > 2:271831/1‑141 (MQ=255) aTGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatat > 2:1922/1‑141 (MQ=255) ggTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataa < 1:262295/141‑1 (MQ=255) tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTAACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataata > 2:40829/1‑141 (MQ=255) aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGaa < 2:227690/141‑1 (MQ=255) gcaGTAACAATCATCACCATTCAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCGTGAAATCACAATGGTCATAGAAACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCa > 2:303771/2‑141 (MQ=255) aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAg > 1:289683/1‑141 (MQ=255) cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATgg < 1:327325/141‑1 (MQ=255) tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGa > 2:108249/1‑141 (MQ=255) | TTGAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGA > CP000730/253840‑254109 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |