Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 76,889 | G→A | 22.8% | intergenic (+947/‑187) | 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,889 | 0 | G | A | 22.8% | 73.2 / 11.5 | 35 | intergenic (+947/‑187) | USA300HOU_0072/USA300HOU_0075 | universal stress protein/lysophospholipase |
Reads supporting (aligned to +/- strand): ref base G (12/15); new base A (2/6); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.31e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 |
TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA > CP000730/76751‑77027 | tATCAAAATATCAATACGGCAAAGAATACTTTAAAAGGTATTGATTGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACg < 2:262835/141‑1 (MQ=11) atatCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTg > 1:230292/1‑141 (MQ=11) aacggCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTTTAAATGGTATGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACTCTGACATGAt < 1:288281/137‑1 (MQ=1) tacaaGGGATTGAATGTACTTACGCTATTTATAAAAAGAAACGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGACACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCta < 1:119266/138‑1 (MQ=1) ggtcttgggtggaTTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttt > 2:344883/13‑141 (MQ=1) agtgcaTTTACGCTCTATATAAAAAGAACCGCACGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttttttactt < 1:20127/136‑5 (MQ=1) ttgaCGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc > 2:111397/4‑141 (MQ=1) tcgatataAAAATAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcca < 2:225532/138‑2 (MQ=1) tccatctaAAAAGAACCGCAGGTCTCTTCACATCTTCGGATTTTGGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAAcc < 2:217049/135‑1 (MQ=1) cctacagaaaaaGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcca < 1:218019/134‑2 (MQ=1) tatataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCAGCCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCta > 2:162075/1‑141 (MQ=11) atataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCtat > 1:145919/1‑141 (MQ=31) ataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaa < 2:21387/141‑1 (MQ=17) aaaaaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa < 1:237099/141‑1 (MQ=34) aaaaaGAACCGCAGGTCTCATCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa > 1:238324/1‑141 (MQ=21) aaaaGAAACGCAGCTCTTTTCAGATCTACGGATGTTCGCTATGCCACGAAATTAGCACCATGATAGCAAGTTAAGCGAACACTGACATTATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaaa < 2:86919/141‑1 (MQ=11) aCCGCAGGCCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATgg < 2:347105/141‑1 (MQ=255) cGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGAGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAACCATGGGc > 1:303927/1‑141 (MQ=255) gCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAAAACCTACAAAAAACATGGGCt < 1:213476/141‑1 (MQ=255) gCAGGTCTCTTCAGATCTACGGACTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCt < 2:107346/141‑1 (MQ=255) gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaaa > 2:115110/1‑141 (MQ=255) ttCAGATCTACGGATTTCTGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt < 2:18655/141‑1 (MQ=255) ctCAGAGCTACGGATTTTCGCCATGCCACGCAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATCTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt < 2:256252/140‑1 (MQ=255) gcgtgttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat < 2:195641/136‑1 (MQ=255) aCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat > 1:122281/1‑141 (MQ=255) cGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatata < 1:163463/141‑1 (MQ=255) aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa > 1:171136/1‑141 (MQ=255) aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa > 2:321990/1‑141 (MQ=255) ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTAAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAt < 2:50317/141‑1 (MQ=255) cGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACt > 2:398822/1‑141 (MQ=255) aCGAAAGTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAtt < 1:96047/141‑1 (MQ=255) atcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa > 1:334940/1‑141 (MQ=255) tcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAGTGATTTGGGTCaa < 2:279530/141‑1 (MQ=255) tAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATACAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGAGGTAtt > 2:1362/1‑141 (MQ=255) aaGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa < 2:122281/141‑1 (MQ=255) | TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA > CP000730/76751‑77027 |
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: ‑ |