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 | 96.7% | 97.1 / ‑3.6 | 30 | 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 . (1/0); new base C (11/18); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATG > NC_000913/2767305‑2767554 | tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cccc‑c < 1:859652/139‑1 (MQ=255) tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cccc‑c > 2:1362087/1‑139 (MQ=255) ctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAAc < 1:901915/139‑1 (MQ=255) gatgatGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCt < 1:1362087/139‑1 (MQ=255) tgTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATa < 1:986185/139‑1 (MQ=255) aGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACg > 2:656970/1‑139 (MQ=255) aaTGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCt < 2:404602/139‑1 (MQ=255) ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 1:656970/139‑1 (MQ=255) aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc > 2:758286/1‑114 (MQ=255) aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 1:758286/114‑1 (MQ=255) aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑ccc > 2:1372852/1‑67 (MQ=255) aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑ccc < 1:1372852/67‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTcaca > 1:629902/1‑62 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTcaca < 2:629902/62‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cccc > 2:1198215/1‑53 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cccc < 1:1198215/53‑1 (MQ=255) ttAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa < 1:342155/139‑1 (MQ=255) ttAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa < 1:290058/139‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCAt < 1:88222/48‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCAt > 2:88222/1‑48 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAac < 1:1072894/55‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAac > 2:1072894/1‑55 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAac < 1:247890/55‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAac > 2:247890/1‑55 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTccc > 1:511324/1‑65 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTccc < 2:511324/65‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgt > 1:610022/1‑69 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgt < 2:610022/69‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccc < 1:628385/91‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCCCACAACGTCCTGCATAGGGAAAACGGTCTCCCCTGTGCACCCACGACTCAATTTCTCGCCCTTTAAACTTCCACATAATGTCAACAATTTTCTGGACATCCTatga > 2:764769/1‑139 (MQ=255) tCCAATATGATCAAATAGCATTAA‑CCCCACCCTCACAACGTCCTGCATAGGGAACAAGATTTCCCCTGTGCACCCACGACTAAATTTcccc > 2:628385/1‑91 (MQ=255) cAATATGATCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg < 1:92740/139‑1 (MQ=255) aTCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTccc > 2:1473381/1‑56 (MQ=255) aTCAAATAGCATTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTccc < 1:1473381/56‑1 (MQ=255) aTTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc > 2:1188845/1‑55 (MQ=255) aTTAACCCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 1:1188845/55‑1 (MQ=255) | TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATG > NC_000913/2767305‑2767554 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |