Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,015,758 | G→C | 11.1% | I21M (ATC→ATG) | PP_5266 ← | acetyl‑CoA hydrolase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,015,758 | 0 | G | C | 11.1% | 106.1 / 5.8 | 36 | I21M (ATC→ATG) | PP_5266 | acetyl‑CoA hydrolase family protein |
Reads supporting (aligned to +/- strand): ref base G (15/17); new base C (2/2); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01 |
GGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAAC > NC_002947/6015619‑6015902 | ggCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTg < 2:417359/150‑1 (MQ=255) cATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCa > 1:396464/1‑150 (MQ=255) aGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGAc < 2:547635/140‑1 (MQ=255) aGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGAc > 1:547635/1‑140 (MQ=255) ggAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGc < 1:594326/150‑1 (MQ=255) gTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATg < 2:704177/150‑1 (MQ=255) gTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATg < 1:778752/150‑1 (MQ=255) aTGGTCAGTCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCatatgg > 2:595575/1‑96 (MQ=255) aTGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCatatgg < 1:595575/96‑1 (MQ=255) ggCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCa > 1:99249/1‑150 (MQ=255) gcgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCAt < 2:99249/150‑1 (MQ=255) ccGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCg > 1:144677/1‑102 (MQ=255) ccGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCg < 2:144677/102‑1 (MQ=255) cGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATgggg > 1:527529/1‑150 (MQ=255) gTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGc < 1:450070/150‑1 (MQ=255) cGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGcctc < 1:164486/150‑1 (MQ=255) ggCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAATAGCGGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGTGcctcct > 2:419431/1‑150 (MQ=255) aGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGaa > 1:30941/1‑150 (MQ=255) tAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAggg < 1:728233/150‑1 (MQ=255) aCGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAg > 1:672457/1‑150 (MQ=255) gAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACAt > 1:139947/1‑145 (MQ=255) gAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACAt < 2:139947/145‑1 (MQ=255) gAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCt > 1:228631/1‑150 (MQ=255) gAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCt > 1:617711/1‑150 (MQ=255) gggTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGggt < 2:527529/150‑1 (MQ=255) cccAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGc > 1:224967/1‑150 (MQ=255) gCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAgg < 2:30941/150‑1 (MQ=255) cccAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAggg < 2:672457/150‑1 (MQ=255) cAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCa > 2:558319/1‑150 (MQ=255) aaGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAg > 2:512209/1‑144 (MQ=255) aaGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAg < 1:512209/144‑1 (MQ=255) aggcaggcCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAAc < 1:558319/150‑1 (MQ=255) ccagcCCATATCCATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGaa < 1:581871/85‑1 (MQ=255) ccagcCCATATCCATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGaa > 2:581871/4‑88 (MQ=255) cagcCCATATCCATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGt < 2:3424/68‑1 (MQ=38) cagcCCATATCCATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGt > 1:3424/3‑70 (MQ=38) aggcCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCCACCTCTg < 2:396464/150‑1 (MQ=255) ccATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCa > 1:433517/1‑150 (MQ=255) atatggatatggGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAGCATCCGTGGGTTCAGGATGTGTACATTAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAAc > 2:89539/1‑150 (MQ=255) | GGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAAC > NC_002947/6015619‑6015902 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |