| 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, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,767,438 | 1 | . | C | 77.8% | 60.2 / 13.7 | 27 | intergenic (+83/‑272) | rnlB/yfjP | CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein |
| Reads supporting (aligned to +/- strand): ref base . (4/2); new base C (9/12); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 | |||||||||||
GTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑A‑C‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATG > NC_000913/2767313‑2767570 | gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑caca < 1:466483/139‑1 (MQ=255) tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACaa > 2:556270/1‑139 (MQ=255) cccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATCAACTTGTCCAATATGATCAAATAGCATTGA‑C‑C‑CC‑CCC‑CT‑CACAACGt < 1:819116/139‑1 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACC‑GTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCt > 1:869741/1‑139 (MQ=255) gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑A‑CTCA‑CCC‑TT‑CACAACGTCATGAATAAGGAaccaat > 2:100072/1‑135 (MQ=255) tGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTc < 1:556270/139‑1 (MQ=255) gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑A‑C‑CCACCA‑CT‑CACACTGTCCTGCATGTAGAACACGTATTcc > 1:20993/1‑139 (MQ=255) gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACC‑GTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATAATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTccc < 2:869741/139‑1 (MQ=255) tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑caca > 2:283810/1‑101 (MQ=255) tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑caca < 1:283810/101‑1 (MQ=255) ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CCCAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑c‑c‑cc < 2:1158622/139‑1 (MQ=255) cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTaaa < 1:821573/139‑1 (MQ=255) tCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑A‑A‑AC‑CCC‑CTATAAAACCTCCTTCATAGGGTAAACTTTTTTACATTTGCACCCAACTCCAAAT‑T‑T‑CTC‑ATCCCTGTAAACCTCTCcataa > 2:4820/1‑139 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑A‑A‑CC‑CCC‑CT‑AACAACGTCCTGCATAGAGAACACGTTTGCCACTGTGCAACATAGACTAAAT‑T‑T‑C‑CTCCCCCTTTAAAATTCTCAATAATg > 2:814009/1‑139 (MQ=255) aCTTGTCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATc < 1:814009/139‑1 (MQ=255) tCCAATATGATCAAATAGCATT‑ACC‑C‑CC‑CCC‑CT‑caca < 2:118381/37‑1 (MQ=255) tCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑caca > 1:118381/1‑37 (MQ=255) tCCAATATGATCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACTTCCTGCATAGGGTACACGTTTTCACCTGTGCACCCACGACTAAATAT‑C‑C‑C‑CCCCCTTTAAACTTCTCATTATTTTCACCAATTTTTTGGACATCTTatg > 2:799778/1‑139 (MQ=255) tCCAATATGATCAAATAGCATT‑A‑A‑C‑CCACGC‑CT‑CACAACGCCCTGCCTAGAGAACACGTGTTCCCCTGCGGACGCACGTCTAAAT‑T‑TCC‑C‑CCCCCTTTCAGCTTCTCAATAATGTCACCAATGTGCTGGAAATCCTatg > 1:68180/1‑139 (MQ=255) tCCAATATGATCAAATAGCATT‑A‑A‑C‑CC‑CCCGCT‑CACAACATCCTGCAAAGGGTACACGATCTCCCCTGTGCACCCACGACTAAAT‑TAA‑A‑C‑ACTCCGTTAAACATCTCAATAATGTCACCAAACTTCTggccagcccatg > 2:1094166/1‑129 (MQ=255) aTCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc < 1:100072/139‑1 (MQ=255) tCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTTCATAGGGAACACGTTTTCCACTCTGCACCCACGACTCAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGACAACCTATGACGCCttttc > 2:781852/1‑137 (MQ=255) tCAAATAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCCATTTTCTGGAAATCCTATGATGCCTTTcc > 2:116170/1‑139 (MQ=255) aaaTAGCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg < 1:4820/139‑1 (MQ=255) aGCATT‑A‑ACC‑CC‑CCC‑CT‑CACGTCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCTAAGTCTACAT‑T‑T‑C‑CGCCCCCTTTAATCTTCCTAATAAGGTCAGCAATTTTCTGGAAATCCTCTGATGACTTTCCCGGTaa < 2:20993/139‑1 (MQ=255) gCATT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGGTTTCCCCTGTGCACCCACGACTAAATGT‑C‑C‑C‑CCCCCTTTAAACCTCTCAATGATGGCCCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa > 2:902331/1‑139 (MQ=255) aTT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCAcc > 2:1276239/1‑100 (MQ=255) aTT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCAcc < 1:1276239/100‑1 (MQ=255) aTT‑A‑ACC‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCC‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg < 1:1067383/139‑1 (MQ=255) | GTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑A‑A‑C‑CC‑CCC‑CT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑C‑C‑CCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATG > NC_000913/2767313‑2767570 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |