Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,561,536 | A→T | 22.2% | L189Q (CTG→CAG) | PP_4047 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,561,536 | 0 | A | T | 22.2% | 29.2 / 5.6 | 18 | L189Q (CTG→CAG) | PP_4047 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (6/8); new base T (2/2); total (8/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01 |
ATCGTCGACGCCCCACAGCACAGTGACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCAATTCATGATAAGGGCGCTCATGAATCAACGTGCTCCCGGGGACAGTCACCCCAACACGAGAAAACAG > NC_002947/4561406‑4561660 | aTCGTCGACGCCCCACAGCACAGTGACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTaa > 2:490581/1‑136 (MQ=255) gtcgACGCCCCACAGCACAGTGACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTa > 2:338712/1‑136 (MQ=255) tcgACGCCCCACAGCACAGTGACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTAc < 1:490581/136‑1 (MQ=255) gACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCa > 1:100960/1‑136 (MQ=255) ccAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCTGGTAACTAc > 1:3346/1‑114 (MQ=255) ccAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCTGGTAACTAc < 2:5783/114‑1 (MQ=255) ccAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCTGGTAACTAc > 1:5783/1‑114 (MQ=255) ccAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCTGGTAACTAc < 2:3346/114‑1 (MQ=255) ccAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATc > 2:227823/1‑136 (MQ=255) aaTCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCAGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCt > 1:374920/1‑136 (MQ=255) tCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTtc < 1:227823/136‑1 (MQ=255) tCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTtc < 1:384195/136‑1 (MQ=255) cTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTtctc < 2:100382/136‑1 (MQ=255) aGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCAtt < 2:100960/136‑1 (MQ=255) ccAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCa < 1:338712/136‑1 (MQ=255) gCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCAATTCATGATAAGGGCGCTCATGaa < 2:394417/136‑1 (MQ=255) tCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCAATTCATGATAAGGGCGCTCATGAATCAACGTGCTCCCGGGGACAGTCACCCCAACACGAGaa > 1:410675/1‑136 (MQ=255) aCCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCAATTCATGATAAGGGCGCTCATGAATCAACGTGCTCCCGGGGACAGTCACCCCAACACGAGAAAACAg < 2:292000/136‑1 (MQ=255) | ATCGTCGACGCCCCACAGCACAGTGACCAATCCTGCCTCAAATGCCCAGTAGCCAAAGAAACCCGCATCGCTTCGTGCTTTGTGCCGATCATGCCAATAACAGTTGCGCATGCCTCGATACCACTGCTTCAGGTAACTACGCACTGCCTGTTGCTGCTCATCCTTCTCAGTATTTAAAGCATTGAGCAATTCATGATAAGGGCGCTCATGAATCAACGTGCTCCCGGGGACAGTCACCCCAACACGAGAAAACAG > NC_002947/4561406‑4561660 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |