| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,449,251 | A→C | 100% | S18R (AGC→CGC) | PP_5556 → | bacteriophage lysis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,449,251 | 0 | A | C | 80.1% | 10.4 / ‑4.0 | 18 | S18R (AGC→CGC) | PP_5556 | bacteriophage lysis protein |
| Reads supporting (aligned to +/- strand): ref base A (3/0); new base C (13/1); total (16/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGG > NC_002947/3449119‑3449377 | gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgcg > 1:197366/1‑136 (MQ=255) aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCCGGCCGCCGCGGTGTGGCAGGCCacc > 1:147159/1‑134 (MQ=255) tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGgcgccg > 1:201637/1‑134 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCAGGCCCCCGCCGATAAGGCGCAGCTc > 1:221127/1‑136 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTc > 1:14233/1‑136 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCCCCGCCGATAAGGCGCAGcgc > 1:246813/1‑134 (MQ=255) gcgagcgCTCGCGGTGCGGCCGGGCGCGCTGGGGCCACCCGCCGGGGCGCCGGGGGGCTGCGGCGGGCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg < 1:152039/136‑1 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg > 2:10210/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg > 1:6732/1‑136 (MQ=255) gcgggggCCACCCGCCGGGGCGCGGTGGCGCGGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt < 1:214895/132‑1 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt > 1:152829/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCCGCCTGGGGGTGGGAGGCCCACGCCGTTAAGGCGCCGCTGGCGGCTCAGGCCGGAGGCTTCGGCGAGCAgg > 2:6367/1‑135 (MQ=255) ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGGGGCCCCCGCCGATAAGGCGCAGCTCGCGGATCAGGCCGCGGACTTCGGCAAGCAGTTGggggg > 2:232232/1‑135 (MQ=255) ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCAGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCGCGCGGATCAGGCCGATGACTTCGGCAAGCAGTTGGCggg > 2:61094/1‑135 (MQ=255) ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGCCAACGCCGATAAGGCGCAGCGCGCTGCTAAGGCCGCAGACATCGGCGAGCAGGTGGAggg > 2:9818/1‑135 (MQ=255) ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAAGCGCCGCGCGCAGGTCAGGCCGAGGGCTTCGGGAAGCAGGTGGCggg > 2:227432/1‑135 (MQ=255) gTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGAGGTCTTGGGCAAGCAGTTGGCGGGGGAGGATcgcg > 2:51380/1‑136 (MQ=255) ctgGCCGCCCCCGCGGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCAACGCCTATAATGCGCAGCTCGCGGATAAGGCCGAGGACATAGGCCACCACTTGGGGGGGAAAGATAGCGCGCACGGTCGGGAGCggg > 2:264866/1‑136 (MQ=255) | GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGG > NC_002947/3449119‑3449377 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |