| 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 | 89.3% | 81.2 / 2.7 | 28 | 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 (10/15); total (13/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.73e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑C‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAA > NC_000913/2767303‑2767568 | ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑c‑c‑c < 1:971868/139‑1 (MQ=255) tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑c‑c‑ccccc > 2:643341/1‑139 (MQ=255) tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑c‑c‑ccacc > 2:301647/1‑139 (MQ=255) ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTc < 1:564551/139‑1 (MQ=255) ccGATGAAGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTGCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTc < 1:1286118/139‑1 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTcc < 2:906686/139‑1 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTcc > 1:906686/1‑139 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑c‑c‑ccccccctc > 2:63422/1‑127 (MQ=255) cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑c‑c‑ccccccctc < 1:63422/131‑5 (MQ=255) tGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCGATATGATCAAATAGCATTAA‑C‑CACCCACTAACATCGTCCTCAATAggacagacatattc > 2:1004608/1‑127 (MQ=255) ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg < 1:643341/139‑1 (MQ=255) tctttctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc < 1:246606/118‑1 (MQ=255) tctttctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc > 2:246606/1‑118 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTcaca > 2:458538/1‑62 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTcaca < 1:458538/62‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCATAACGTCCTGCATAGGGAac < 1:1102918/80‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATCGGGAac > 2:1102918/1‑80 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc < 2:1551998/90‑1 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc > 1:1551998/1‑90 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAACTTTCCCCCCCTTTAAACTTCTCAATAATGt > 1:702881/1‑139 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt > 2:101880/1‑139 (MQ=255) gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑C‑CACCCACTCACAACGTCCTGCATACGGAACACGTTTTACACTGTGAAGCCACGGCTAAAATTCCCCACCTTTATACTTCTCaataattt > 1:1123525/1‑137 (MQ=255) tGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa > 2:239117/1‑111 (MQ=255) tGATATTAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa < 1:239117/111‑1 (MQ=255) tAACTTGTCCAATATGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt < 1:101880/139‑1 (MQ=255) tCCAATATGATCAAATAGCATTAA‑CAC‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga > 2:112113/1‑139 (MQ=255) tGATCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt < 1:483866/139‑1 (MQ=255) tCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc > 2:913994/1‑80 (MQ=255) tCAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc < 1:913994/80‑1 (MQ=255) cAAATAGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg < 1:301647/139‑1 (MQ=255) aGCATTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa < 1:112113/139‑1 (MQ=255) aTTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc > 2:674672/1‑130 (MQ=255) aTTAACC‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc < 1:674672/130‑1 (MQ=255) | CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑C‑C‑CCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAA > NC_000913/2767303‑2767568 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |