Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,223,605 | C→T | 22.3% | Y202Y (TAC→TAT) | priA → | primosomal protein N prime |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,223,605 | 0 | C | T | 22.3% | 30.8 / 5.6 | 18 | Y202Y (TAC→TAT) | priA | primosomal protein N prime |
Reads supporting (aligned to +/- strand): ref base C (7/7); new base T (2/2); total (9/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01 |
AAAGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGA > CP000730/1223471‑1223730 | aaaGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTg < 1:256161/141‑1 (MQ=255) tgatATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGc < 2:117489/141‑1 (MQ=255) gaaAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTgaa < 2:151200/141‑1 (MQ=255) cGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGaataaa > 1:208993/1‑140 (MQ=255) aCTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACt < 1:11940/141‑1 (MQ=255) ccAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg < 2:58787/141‑1 (MQ=255) ccAAGCGTACTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg < 2:59206/141‑1 (MQ=255) ccAAGCGTACTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg < 2:59208/141‑1 (MQ=255) gtgtCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCaa < 2:26201/141‑1 (MQ=255) gtgtCATTGAAGGGTATCATCCTGATGAAGTATTAGATAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCaa < 2:55443/141‑1 (MQ=255) atcatcCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATg > 1:184862/1‑141 (MQ=255) atcCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGAc > 2:235804/1‑141 (MQ=255) aGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTg > 1:254800/1‑141 (MQ=255) aGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTg > 2:53440/1‑141 (MQ=255) gagaAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGaaaaa > 2:69841/1‑141 (MQ=255) agaAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAat > 1:199055/1‑141 (MQ=255) aGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGa > 1:203390/1‑141 (MQ=255) aGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGa > 1:203335/1‑141 (MQ=255) | AAAGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGA > CP000730/1223471‑1223730 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |