Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,765,794 | (C)7→8 | intergenic (‑18/‑124) | yfjM ← / → rnlA | CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,765,787 | 1 | . | C | 98.2% | 194.0 / ‑3.4 | 57 | intergenic (‑11/‑131) | yfjM/rnlA | CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (20/36); total (21/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.68e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01 |
GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT > NC_000913/2765650‑2765919 | ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg > 2:271592/1‑138 (MQ=255) ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTg < 1:271592/138‑1 (MQ=255) tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTTTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑cc‑cc < 1:831152/139‑1 (MQ=255) ctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑cc‑ccc < 1:297176/139‑1 (MQ=255) ctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑cc‑ccc < 2:426108/139‑1 (MQ=255) tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑cc‑cccc < 1:370249/139‑1 (MQ=255) tCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGAc > 2:629862/1‑139 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 1:194436/139‑1 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 1:715794/139‑1 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 1:868433/139‑1 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 2:865616/139‑1 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 1:477368/139‑1 (MQ=255) tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 2:276932/139‑1 (MQ=255) cGTAGAATTTTCTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGATTGCCC‑CCCCCGGTAGACGATCCTGccc < 1:682767/139‑1 (MQ=255) tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCta > 2:228948/1‑139 (MQ=255) tAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCta > 1:876866/1‑139 (MQ=255) ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 2:819507/139‑1 (MQ=255) ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 1:361520/139‑1 (MQ=255) ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 2:1034378/139‑1 (MQ=255) ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 1:169046/139‑1 (MQ=255) ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 2:464814/139‑1 (MQ=255) aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 1:106414/99‑1 (MQ=255) aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc > 2:106414/1‑99 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 2:646074/139‑1 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 1:851589/139‑1 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg > 2:851589/1‑139 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 1:228948/139‑1 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg > 1:825746/1‑139 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 2:178972/139‑1 (MQ=255) aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 2:825746/139‑1 (MQ=255) tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTg < 1:833167/94‑1 (MQ=255) tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTg > 2:833167/1‑94 (MQ=255) cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt < 2:423748/139‑1 (MQ=255) cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCTCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt > 1:539603/1‑139 (MQ=255) ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca > 2:676945/1‑139 (MQ=255) ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca > 2:410264/1‑139 (MQ=255) tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg < 1:319082/139‑1 (MQ=255) tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg < 1:676945/139‑1 (MQ=255) gaTCAGGAAGAGGTTTCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg < 2:833538/139‑1 (MQ=255) aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa > 1:543260/1‑122 (MQ=255) aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa < 2:543260/122‑1 (MQ=255) ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg > 1:493126/1‑125 (MQ=255) ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg < 2:493126/125‑1 (MQ=255) ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag > 1:662168/1‑139 (MQ=255) ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag > 2:135248/1‑139 (MQ=255) tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt < 2:662168/139‑1 (MQ=255) ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 2:200250/73‑1 (MQ=255) ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg > 1:200250/1‑73 (MQ=255) aCGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa < 1:171966/139‑1 (MQ=255) cGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc > 1:109802/1‑50 (MQ=255) cGATGGACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATc < 2:109802/50‑1 (MQ=255) ggACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa > 2:553032/1‑100 (MQ=255) ggACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa < 1:553032/100‑1 (MQ=255) ggACATCATCCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt < 2:539603/139‑1 (MQ=255) atcCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCATATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg < 2:876866/139‑1 (MQ=255) tcCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg < 1:922191/55‑1 (MQ=255) tcCCCCATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAg > 2:922191/1‑55 (MQ=255) ccATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGccc > 2:696978/1‑40 (MQ=255) ccATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGccc < 1:696978/40‑1 (MQ=255) ccATTTTGTGGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc < 2:330483/139‑1 (MQ=255) gtgGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg < 1:935345/71‑1 (MQ=255) gtgGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg > 2:935345/1‑71 (MQ=255) tgGACTGCCC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCCGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTCTTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt > 2:170889/1‑139 (MQ=255) | GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CC‑CCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT > NC_000913/2765650‑2765919 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |