Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,767,445 | (C)7→8 | intergenic (+90/‑265) | rnlB → / → yfjP | CP4‑57 prophage; antitoxin RnlB/CP4‑57 prophage; putative GTP‑binding protein YfjP |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,767,438 | 1 | . | C | 88.2% | 37.6 / 0.7 | 17 | intergenic (+83/‑272) | rnlB/yfjP | CP4‑57 prophage; antitoxin RnlB/CP4‑57 prophage; putative GTP‑binding protein YfjP |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (3/12); total (5/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.35e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01 |
TGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTT > NC_000913/2767312‑2767558 | tGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑ccccc > 1:302042/1‑133 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTcaca < 2:302042/133‑1 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATAATAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTcaca < 1:351658/133‑1 (MQ=255) gatgTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGt < 2:185957/133‑1 (MQ=255) atgTACGCTAAGTCCTGTCAATTGATTTTTAAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTc < 1:187422/133‑1 (MQ=255) tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTc < 2:289527/133‑1 (MQ=255) aaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 2:50594/133‑1 (MQ=255) aaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 2:475061/133‑1 (MQ=255) aaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTc < 1:437428/133‑1 (MQ=255) taGAATTCGCCCTTTCGGCATGATCTTATCTTTTCCAATATAATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa < 2:183126/133‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGCCTAAATTTCCCCCCCTTTAAACTTCTCaa > 1:27537/1‑133 (MQ=255) aCTTGTCCTATATGATCACATACCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg < 1:138990/133‑1 (MQ=255) gTCCAATATGATCAAATAGCATTAAA‑C‑ACCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCAACGACTAAATTTACCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTATGGAAAt > 1:398517/1‑133 (MQ=255) tCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg < 1:498838/128‑1 (MQ=255) tCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg > 2:498838/1‑128 (MQ=255) tCCAATATGATCAAATAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAGTTCCCCACCTTAAAACTTCACAAGAATGTCACGAATTATCTGGAAATc > 2:134462/1‑133 (MQ=255) tCCAATATGATCAAATAGCATT‑AA‑CGCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTACTCAATAGTGTGACGAATTTTGTGGAAATc > 1:167577/1‑133 (MQ=255) aaaTAGCATT‑AACC‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt < 1:281847/133‑1 (MQ=255) | TGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑CCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTT > NC_000913/2767312‑2767558 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |