Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 564,385 | Δ1 bp | 9.1% | intergenic (+34/‑139) | rplQ → / → katA | 50S ribosomal protein L17/catalase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 564,385 | 0 | T | . | 9.1% | 127.2 / 7.7 | 44 | intergenic (+34/‑139) | rplQ/katA | 50S ribosomal protein L17/catalase |
Reads supporting (aligned to +/- strand): ref base T (17/22); new base . (2/2); total (20/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01 |
GCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGA > NC_002947/564239‑564533 | gcggcTCACCGCGCATCCGGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGcc < 2:476347/150‑1 (MQ=255) cTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAGGATTGGGGGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCggg > 1:618685/1‑149 (MQ=255) cTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGt < 1:445401/150‑1 (MQ=255) cTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCCGGCGACAACGCGCCTATGGCGTGTTTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAG‑GCCCGGt < 1:12096/150‑1 (MQ=255) tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTGCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGtt < 2:712142/150‑1 (MQ=255) tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGtt < 2:1014280/150‑1 (MQ=255) ccTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGtttt < 2:217881/150‑1 (MQ=255) cTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGttttt > 2:1054787/1‑150 (MQ=255) gcATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGtt < 2:912548/150‑1 (MQ=255) cATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAG‑GCCCGGTTTTCTGtag < 1:278053/150‑3 (MQ=255) ccTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGTTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCCGGC‑CTAGTGCCCGGTTTTTTGTTtata < 2:363010/150‑1 (MQ=255) cGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTg < 1:977259/150‑1 (MQ=255) ggTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTCCAGTCGACAAAGCTGGACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGt > 2:308603/1‑150 (MQ=255) ccGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCTGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAtt < 2:348305/150‑1 (MQ=255) aaCGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTaa < 1:1054787/150‑1 (MQ=255) gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAt < 2:16892/109‑1 (MQ=255) gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAt > 1:16892/1‑109 (MQ=255) gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAAt > 2:1084102/1‑150 (MQ=255) tCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACCGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc > 2:651027/1‑150 (MQ=255) gTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTc > 2:994728/1‑150 (MQ=255) gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCaa > 1:782925/1‑150 (MQ=255) tGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATCCAGGCGAGTAATCTTGCTCGTCAATActctct < 1:308603/150‑1 (MQ=255) gCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc < 1:651027/150‑1 (MQ=255) gctgctGAGTAAGACGTTCTGTCTGCTACAGAAAACGGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGa < 2:359270/150‑1 (MQ=255) gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTa > 2:668219/1‑150 (MQ=255) gACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATg > 1:593117/1‑87 (MQ=255) gACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATg < 2:593117/87‑1 (MQ=255) cGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTg < 2:898415/150‑1 (MQ=255) ttCTGTCTGCTACAGAAAACCGGGC‑CTAGGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTTATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTTTCTCTACATCCGATTAGCCGGCCCTTGaa > 2:209031/1‑150 (MQ=255) tgtctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAc < 2:19310/150‑1 (MQ=255) tctgCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCg > 2:655572/1‑150 (MQ=255) gCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGaaa < 1:282708/150‑1 (MQ=255) cAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc > 2:65655/1‑113 (MQ=255) cAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctc < 1:65655/113‑1 (MQ=255) gAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgaga > 1:324956/1‑150 (MQ=255) aaaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGc < 2:782925/150‑1 (MQ=255) aaCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc > 2:824727/1‑150 (MQ=255) aaCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc > 2:935464/1‑150 (MQ=255) aCCGGGCACTAG‑GCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGccc > 2:751092/1‑150 (MQ=255) aCCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGccc > 1:1072005/1‑150 (MQ=255) cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt < 2:1072005/150‑1 (MQ=255) cGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt > 2:360316/1‑150 (MQ=255) tAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCaa < 1:444565/150‑1 (MQ=255) gTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTCGACGTCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGa > 1:10404/1‑150 (MQ=255) | GCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGC‑CTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGA > NC_002947/564239‑564533 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |