Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 498,343 | G→A | 12.9% | L192L (CTG→CTA) | dctS → | two‑component sensor histidine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 498,343 | 0 | G | A | 12.9% | 83.7 / 3.4 | 31 | L192L (CTG→CTA) | dctS | two‑component sensor histidine kinase |
Reads supporting (aligned to +/- strand): ref base G (12/15); new base A (2/2); total (14/17) | |||||||||||
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 |
AACCAAATAGGCGTGGTGCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGATCGCTATTGACAATCGTCTC > BS168_reference_genome/498212‑498466 | aaCCAAATAGGCGTGGTGCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGc > 7:143714/1‑138 (MQ=255) aaTAGGCGTGGTGCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTAGCCCGCCatat > 1:87129/1‑138 (MQ=255) aaTAGGCGTGGTGCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTAGCCCGCCatat > 1:90036/1‑138 (MQ=255) tggtgCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAAc > 3:78664/1‑138 (MQ=255) tggggCTTGTCGCAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTAGGCCTTGCTGGTTCCTATTTGCTGGCCCGCCATATCAAAAAACa < 1:84425/139‑1 (MQ=255) gtgCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACaaa > 6:90710/1‑139 (MQ=255) gtgCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACaaa > 3:90902/1‑139 (MQ=255) gtgCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTCGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACaaa > 4:185443/1‑139 (MQ=255) cTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGt > 6:74235/1‑139 (MQ=255) gTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGtt < 3:129658/137‑1 (MQ=255) gAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCt < 5:124652/139‑1 (MQ=255) aaaCGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGa < 3:212475/139‑1 (MQ=255) tCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGt < 1:5401/139‑1 (MQ=255) gCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTc < 8:143714/138‑1 (MQ=255) cTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCg < 8:168526/139‑1 (MQ=255) ccTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTAt < 8:85608/139‑1 (MQ=255) aTTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAAc < 4:78664/138‑1 (MQ=255) tAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGa < 5:74235/138‑1 (MQ=255) tGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTgg < 5:127941/62‑1 (MQ=255) tGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTgg > 6:127941/1‑62 (MQ=255) tGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGtt > 2:91357/1‑87 (MQ=255) tGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGtt < 1:91357/87‑1 (MQ=255) ccTTCATCGTTGTGCTGCCACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGACCCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTc > 1:115569/1‑139 (MQ=255) cATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGat < 4:90902/138‑1 (MQ=255) cATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTAGCCCGCCATATTAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGat < 2:87129/138‑1 (MQ=255) cATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTAGCCCGCCATATTAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATCAAGAACGGACAGCGACGTTTCACTCGat < 2:90036/138‑1 (MQ=255) aTCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGatg > 3:176408/1‑138 (MQ=255) gtgCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTc < 7:151577/70‑1 (MQ=255) gtgCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTc > 8:151577/1‑70 (MQ=255) acacTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGa > 7:210/1‑139 (MQ=255) gCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGATCGCTATTGACAAtcgt < 7:76201/138‑1 (MQ=255) gCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGATCGCTATTGACAAtcgt > 8:76201/1‑138 (MQ=255) gTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTGCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGATCGCTATTGACAAtcgtctc < 2:11848/137‑1 (MQ=255) | AACCAAATAGGCGTGGTGCTTGTCGGAAAAACGCTGCCGGGGATCGCTGACATCCTGCTTCATTTAAAGCGTGATATCGCCTTCATCGTTGTGCTGACACTCGGCTTTGGCCTTGCTGGTTCCTTTTTGCTGGCCCGCCATATCAAAAAACAAATGTTTCAGCTTGAACCCCATGAGATTGTCCGAATGTATGAAGAACGGACAGCGACGTTTCACTCGATGAATGAAGGGGTGATCGCTATTGACAATCGTCTC > BS168_reference_genome/498212‑498466 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |