Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,896,445 | G→A | 11.9% | intergenic (+78/‑130) | PP_4304 → / → sbp‑I | VIC family cation transporter/sulfate ABC transporter substrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,896,445 | 0 | G | A | 11.9% | 104.7 / 4.7 | 42 | intergenic (+78/‑130) | PP_4304/sbp‑I | VIC family cation transporter/sulfate ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base G (20/17); new base A (2/3); total (22/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTT > NC_002947/4896300‑4896588 | agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt > 8:186547‑M2/1‑146 (MQ=255) ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgc < 5:1442‑M2/148‑10 (MQ=255) cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGGCCAATTCCCGGGTAAAACCCGACAAACCGGCCCGGACCGTGCCCGGACC‑‑‑CCCGTAGGAGCGGGTTTAccccc > 2:223560/1‑146 (MQ=255) cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg > 7:174384/1‑136 (MQ=255) cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg < 8:174384/136‑1 (MQ=255) cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaag > 6:1442‑M2/1‑136 (MQ=255) gcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAaga > 4:176819‑M1/68‑149 (MQ=255) ctataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAc > 7:169467‑M1/64‑136 (MQ=255) ctataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAc < 8:169467‑M1/73‑1 (MQ=255) gggggCGGCGTTCTGTTCGCGCTGGGGCAAGGGGTTGTTCCCCAAAGCGTGATAGCGCCTGCGCG‑GGCCTCTTGGCGGGTAAATCCGCTCCTACCGGCGGGGACGGTGCGCGTTCC‑‑‑CCCTTTGGAGCGGGTTTACCCGCGAAGAGGCCg < 2:308491/149‑1 (MQ=37) taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa > 8:31902‑M1/61‑140 (MQ=255) taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa < 7:31902‑M1/80‑1 (MQ=255) taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCC‑GCCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCg > 2:99861‑M1/61‑147 (MQ=255) ttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCACTTAGCGGGTAAATCCGCTCCTAAAGGCGCGGACCGGGCGCGGTGC‑‑‑CCCGTAGGAGCGGGTTTAACCGCGAAGAGACCGGCACAgg > 1:269850/1‑146 (MQ=255) ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt > 8:35444‑M2/1‑118 (MQ=255) ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt < 7:35444‑M2/121‑4 (MQ=255) gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACTGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc > 6:251500‑M2/1‑116 (MQ=255) gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc < 5:251500‑M2/137‑22 (MQ=255) aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCACATAGCCGGTAAATACGCTACAAAACGCCCGGACCCTGCCCGGTCC‑‑‑CCCCTACGAGCGGGTTTAACCCCCAAGAGGCCGGCACAAGCAAAAGC‑AAAACCAAAACaat > 6:297037/1‑147 (MQ=255) tcgcAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCtt < 5:297037/145‑1 (MQ=255) ccAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta > 6:305389‑M2/1‑89 (MQ=255) ccAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta < 5:305389‑M2/93‑5 (MQ=255) cAAAGCGTGATAGCGCCTGCGCC‑GGCCTAATAGCGGGTAAAAACCCACCAACAGGGGCGGGCCCGGGGCGGTGC‑‑‑CCCGGATGGGCGGGGTTTACAGCCAAAAGGCCCGCACAAAGAGAAGCAAAAACCAAAAAATtttttttttata > 8:14087/1‑143 (MQ=255) cAAAGCGGTATAGCGCCTGAGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGCCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGTAGCGGGTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTc < 7:14087/143‑1 (MQ=38) aaaGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAAcc > 7:164824/1‑149 (MQ=255) cgcC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGt < 1:223560/148‑1 (MQ=255) GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAGCCCGTAATAACCCACGCGTTATAGt < 5:267789‑M1/149‑1 (MQ=255) gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtctt < 1:319028‑M2/149‑98 (MQ=255) ATCCGCTCCTACAGGCGCTGACCGTTCGCGGTCC‑‑‑CCCGTAGTAGCGGGTTTACCCGCGAAGAGGCCGTCACAGGCAGAAGC‑ATAACCAAACCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc > 6:228533‑M1/1‑149 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt < 2:168079‑M2/149‑103 (MQ=255) ccGCTCCTACAGGCGCGGACCGTGCGAGGTCC‑‑‑CCCGAAGGAGCGGGtttacccgcgaagaggccggcgcaggccctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggacttcttgtg < 6:258367‑M2/149‑104 (MQ=255) cGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgcccgccccctgctcgtgggtcttcttgt > 1:168079‑M2/1‑45 (MQ=255) cTACAGGCGCGGACCGTGCGCGGTCG‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta < 4:185116/121‑1 (MQ=255) cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggt < 7:570‑M2/149‑110 (MQ=255) cTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta > 3:185116/1‑121 (MQ=255) tACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagctcgccgccccctgctcgtgggtcttcttgtggcaggtc < 5:306220‑M2/149‑111 (MQ=255) aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca > 7:120710/1‑142 (MQ=255) aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAaca < 8:120710/142‑1 (MQ=255) cgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatc > 8:570‑M2/1‑31 (MQ=255) cgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc < 2:269850/149‑1 (MQ=255) cGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGa > 7:97208/1‑149 (MQ=255) gAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAActct > 5:263261/1‑149 (MQ=255) aGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt < 6:263261/149‑1 (MQ=255) gCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctccccacgcatgg > 2:319028‑M2/1‑5 (MQ=255) gtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtct > 7:152291‑M2/1‑1 (MQ=255) gtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtct < 8:152291‑M2/97‑97 (MQ=255) | AGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTT > NC_002947/4896300‑4896588 |
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: ‑ |