Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 315,244 | G→A | intergenic (+16/‑47) | eaeH → / → insE1 | pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 315,244 | 0 | G | A | 92.9% | 106.9 / 0.1 | 43 | intergenic (+16/‑47) | eaeH/insE1 | pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A |
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (18/21); total (21/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 |
GGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAG > NC_000913/315113‑315365 | ggCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac < 1:1099536/139‑1 (MQ=255) gcgcAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcg < 2:1237433/139‑1 (MQ=255) gCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac > 2:435320/1‑132 (MQ=255) gCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac < 1:435320/132‑1 (MQ=255) tATGAACAGTATTCTGGCGATGAAGTCGGGCTGTTTGGTAAAGATCAGAGCCCGAAAGACCGGCATTCGATTTCTGCCGGGGTGACCTATACGCCATTGCCTCTTACCAAGAAATAATGGACACGCGGCTAAGTGAGTa < 2:466239/139‑1 (MQ=255) tATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa < 2:101940/139‑1 (MQ=255) aCAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc < 1:280737/139‑1 (MQ=255) cAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct < 1:385716/139‑1 (MQ=255) cAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct < 1:66520/139‑1 (MQ=255) aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTa > 1:260169/1‑124 (MQ=255) aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTa < 2:260169/124‑1 (MQ=255) tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag > 2:375459/1‑139 (MQ=255) aTGACGTCGGGCTGTTTGGTTAAGATAAGCGCCAGAAAGACCCGCATGCTAGTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTGACCCAGCAATATTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGCGa > 2:76783/1‑139 (MQ=37) aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt > 2:977032/1‑130 (MQ=255) aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt < 1:977032/130‑1 (MQ=255) tAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc > 2:824242/1‑112 (MQ=255) tAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtc < 1:824242/112‑1 (MQ=255) aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGACAGAGGTGACTCACATGACAAAAACAGTAACaa > 2:337323/1‑139 (MQ=255) gcgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa < 1:1182695/126‑1 (MQ=255) gcgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa > 2:1182695/1‑126 (MQ=255) cgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt < 1:688511/139‑1 (MQ=255) aGATAGACCCGCATGCTATTACTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac > 1:913799/1‑111 (MQ=255) aGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac < 2:913799/111‑1 (MQ=255) aGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtcag > 2:396645/1‑100 (MQ=255) aGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtcag < 1:396645/100‑1 (MQ=255) gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc > 1:238271/1‑91 (MQ=255) gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc < 2:238271/91‑1 (MQ=255) gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa > 1:432880/1‑120 (MQ=255) gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa < 2:432880/120‑1 (MQ=255) ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg < 1:580618/51‑1 (MQ=255) ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg > 2:580618/1‑51 (MQ=255) ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtc < 2:275932/74‑1 (MQ=255) ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtc > 1:275932/1‑74 (MQ=255) ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg > 1:177556/1‑139 (MQ=255) ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac > 2:876206/1‑83 (MQ=255) ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac < 1:876206/83‑1 (MQ=255) taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTca > 1:775654/1‑73 (MQ=255) taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTca < 2:775654/73‑1 (MQ=255) gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTg > 1:775037/1‑139 (MQ=18) tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa > 2:68864/1‑139 (MQ=18) tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa > 2:362210/1‑139 (MQ=18) ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAg > 2:1277026/1‑139 (MQ=14) ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAg < 1:68864/139‑1 (MQ=14) | GGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAG > NC_000913/315113‑315365 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |