| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,763,813 | T→C | E447E (GAA→GAG) | abpA ← | CP4‑57 prophage; anti‑bacteriophage protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,763,813 | 0 | T | C | 97.0% | 106.0 / ‑3.9 | 33 | E447E (GAA→GAG) | abpA | CP4‑57 prophage; anti‑bacteriophage protein |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (9/23); total (10/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.03e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01 | |||||||||||
GAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGC > NC_000913/2763739‑2763886 | gagTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTc < 2:497751/80‑1 (MQ=255) gTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGt < 1:497752/80‑1 (MQ=255) ccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGc > 2:497801/1‑81 (MQ=255) ccACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGc < 2:497753/81‑1 (MQ=255) aCGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGa < 1:497754/81‑1 (MQ=255) aCGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGa < 1:497755/81‑1 (MQ=255) gaacagaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAg > 2:497797/1‑81 (MQ=255) acagaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAg > 1:497803/1‑79 (MQ=255) acagaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCa < 2:497756/81‑1 (MQ=255) cagaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAg < 2:497757/81‑1 (MQ=255) agaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAgg < 2:497758/81‑1 (MQ=255) gaacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGt < 2:497759/81‑1 (MQ=255) aacaAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTa < 2:497760/81‑1 (MQ=255) aaaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTGGtt < 1:497761/81‑1 (MQ=255) aaaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGtt < 1:497762/81‑1 (MQ=255) aaCGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTa > 1:497791/1‑81 (MQ=255) cGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATc < 2:497763/81‑1 (MQ=255) gAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTc > 1:497794/1‑81 (MQ=255) aCGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCttt < 1:497764/81‑1 (MQ=255) gAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGt < 2:497765/80‑1 (MQ=255) gAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTc < 2:497766/81‑1 (MQ=255) aGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCg < 2:497767/81‑1 (MQ=255) tGAGCATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAgg < 1:497768/81‑1 (MQ=255) cATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑ttt > 1:497789/1‑81 (MQ=255) gAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑ttttttt > 1:497811/1‑81 (MQ=255) gAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑tttttt < 1:497769/80‑1 (MQ=255) aGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑tttttttttg > 1:497804/1‑78 (MQ=255) cTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTTtgcgcgcgg > 1:497807/1‑78 (MQ=255) ccATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTAt < 1:497770/81‑1 (MQ=255) aaGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGc > 1:497805/1‑80 (MQ=255) aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt < 2:497771/80‑1 (MQ=255) aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt < 1:497772/80‑1 (MQ=255) aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt < 1:497773/80‑1 (MQ=255) | GAGTCGTAGCCCACGAACAGAACAAACGTGAAACGATCGAGGTGAGCATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGC > NC_000913/2763739‑2763886 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |