Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,431,130 | C→T | 26.6% | D29D (GAC→GAT) | PP_3042 → | terminase large subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,431,130 | 0 | C | T | 26.6% | 22.0 / 5.0 | 15 | D29D (GAC→GAT) | PP_3042 | terminase large subunit |
Reads supporting (aligned to +/- strand): ref base C (5/6); new base T (2/2); total (7/8) | |||||||||||
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 |
ATTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGT > NC_002947/3431007‑3431239 | aTTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGAt > 2:193811/1‑136 (MQ=255) tCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCAGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATc < 2:122203/136‑1 (MQ=255) atAAAGAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGATGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTTTGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCcgc < 1:317460/135‑1 (MQ=255) cccACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGc > 2:435988/1‑136 (MQ=255) cccACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGc > 2:438534/1‑136 (MQ=255) ggTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTaccgcaccg > 2:333788/1‑136 (MQ=255) aCCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCa > 1:366691/1‑62 (MQ=255) aCCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCa < 2:366691/62‑1 (MQ=255) cgcgAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt > 1:377590/1‑105 (MQ=255) cgcgAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt < 2:377590/105‑1 (MQ=255) gAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGccc < 1:435988/136‑1 (MQ=255) gAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGccc < 1:438534/136‑1 (MQ=255) cGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt < 1:96651/90‑1 (MQ=255) cGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt > 2:96651/1‑90 (MQ=255) aCGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGt < 1:193811/136‑1 (MQ=255) | ATTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGT > NC_002947/3431007‑3431239 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |