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 | 88.0% | 70.3 / 3.5 | 25 | 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 (8/14); total (11/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.17e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA > NC_000913/2767301‑2767552 | aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑aaa > 1:1422389/1‑138 (MQ=255) ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑a‑cc > 2:743919/1‑139 (MQ=255) tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c < 1:109423/139‑1 (MQ=255) tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c > 1:1000212/1‑139 (MQ=255) tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGGCAATGGATTTATGAAACTCTTCTTTCGAAACGAACCAGTATAGAATTCGCACTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c > 1:976492/1‑139 (MQ=255) gTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑cccc < 2:322691/139‑2 (MQ=255) gTCTCCCGATGATGTACGCCAAGTCCTGTCGATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTcaca < 1:613161/139‑1 (MQ=255) gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑cccctcgca < 2:819473/139‑7 (MQ=255) gTCTCCCGATGATGTACGCCAAGTCCAGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTcaca < 2:976492/139‑1 (MQ=255) ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTc < 2:1000212/139‑1 (MQ=255) ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAac < 2:1516529/139‑1 (MQ=255) tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTTACAACGTCCTGCATAGGGAACACGTTTTccc < 1:1407764/139‑1 (MQ=255) tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c > 1:261953/1‑95 (MQ=255) tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c < 2:261953/95‑1 (MQ=255) tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc > 2:651362/1‑92 (MQ=255) tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc < 1:651362/92‑1 (MQ=255) aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc > 1:1251569/1‑86 (MQ=255) aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc < 2:1251569/86‑1 (MQ=255) tcttctTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTaa < 2:924612/139‑1 (MQ=255) aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa < 2:1422389/139‑1 (MQ=255) tgtataagagacagATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt < 1:925705/126‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCC‑CCC‑CCCTCACAACGTCCTGCATAGGGAac < 2:616195/81‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCC‑CCC‑CCCTCACAACGTCCTCCATAGGGaaa > 1:616195/1‑80 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc < 2:641780/102‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc > 1:641780/1‑102 (MQ=255) tGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTGTCAATAATGTCACCAATTTTCTGGAAATCCTat > 1:524146/1‑139 (MQ=255) gTCCAATATGATCAAATAGCATTAAA‑‑C‑CCC‑GCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAACTTCACCAATATTCTGGAAATCCTatg > 1:1279786/1‑139 (MQ=255) tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACTTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga > 1:625888/1‑139 (MQ=255) tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATa < 2:1087858/49‑1 (MQ=255) tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATa > 1:1087858/1‑49 (MQ=255) tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa < 1:464737/100‑1 (MQ=255) tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTGTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATcccctgt > 2:194899/1‑134 (MQ=255) tCCAATATGATCAAATAGCATT‑AA‑‑CACCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa > 2:464737/1‑100 (MQ=255) tCCAATATGATCAAATAGCATT‑AA‑‑C‑CCCTCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTTCCCCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCCGGAAATCCTatga > 2:370148/1‑139 (MQ=255) tttgatcaTAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc > 2:46497/8‑73 (MQ=255) tttgatcaTAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc < 1:46497/66‑1 (MQ=255) cATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc < 2:811556/56‑1 (MQ=255) cATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc > 1:811556/1‑56 (MQ=255) | AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA > NC_000913/2767301‑2767552 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |