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 | 90.9% | 96.2 / 2.4 | 33 | 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 . (3/0); new base C (13/17); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGT > NC_000913/2767303‑2767572 | ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cct > 2:813949/1‑138 (MQ=255) ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑ccc < 1:897594/139‑1 (MQ=255) tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cccc‑c > 1:710743/1‑139 (MQ=255) tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑cacc‑acc > 1:899724/1‑139 (MQ=255) tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑accc‑ccc > 1:108772/1‑139 (MQ=255) gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑caca < 2:710743/139‑1 (MQ=255) ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTc < 2:838477/139‑1 (MQ=255) atgatgTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCAA‑CCCTACAAAAAGTCCTg > 1:950627/1‑139 (MQ=255) aGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGATTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAggg > 1:190592/1‑133 (MQ=255) aGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAggg < 2:190592/133‑1 (MQ=255) tttctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTc < 2:899724/139‑1 (MQ=255) gtataagagacaggTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc < 1:50575/126‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑caca > 2:825944/1‑62 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑caca < 1:825944/62‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGt > 2:818578/1‑66 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGt < 1:818578/66‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTc < 1:848730/67‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTc > 2:848730/1‑67 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CAAAACGTCCTGCATAGGGAACACGTTTTTCCATGTGAACCCAAGACAAAATTTACCCCACTTTAAAcctgtatcttataa > 2:50575/1‑126 (MQ=255) atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg < 2:950627/139‑1 (MQ=255) tAACTTGTCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt < 1:813949/139‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTACTTTTCCCCCCCTTTAAACTTCTCAATAATGTCAcc > 2:603093/1‑118 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCAcc < 1:603093/118‑1 (MQ=255) tCCAATATGATCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga > 2:177233/1‑139 (MQ=255) tCCAATATGATCAAATAGCATTAA‑CCCCACCAT‑CACAACGTCCTGCATATTGAACACGTATTCCCCTGTGCACCCACGACTAAACTTCCCCCCCTTTAAACTTCTCAACAATGTAACCAATTTTCTGGAAATCCTataa > 1:62906/1‑137 (MQ=255) aTCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCCCCCACGACTAAATTTccc > 2:935898/1‑81 (MQ=255) aTCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc < 1:935898/81‑1 (MQ=255) tCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTcccc > 2:704984/1‑56 (MQ=255) tCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTcccc < 1:704984/56‑1 (MQ=255) tCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg < 1:92045/131‑1 (MQ=255) tCAAATAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTCTTCACCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg > 2:92045/1‑131 (MQ=255) aTAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa > 2:586168/1‑139 (MQ=255) aTAGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGACACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa > 2:338357/1‑139 (MQ=255) aGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa < 1:586168/139‑1 (MQ=255) aGCATTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa < 1:338357/139‑1 (MQ=255) aTTAACCCCC‑CCCT‑CACAACGTCCTGGATAGGCAACACGTTTTCCCCTGTGCAcc > 2:510937/1‑55 (MQ=255) aTTAACCCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 1:510937/55‑1 (MQ=255) ttAA‑CCCCACCCT‑CAAAACGTCCTGCATAGGGAACACGATTTCCACAGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTACCAAAAAATTCAACCAATTTTAAGCAAACCGTATGAGACATTTCCTGGTAAATgga > 2:531426/1‑138 (MQ=255) | CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCC‑CCCT‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGT > NC_000913/2767303‑2767572 |
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: ‑ |