Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,896,377 | C→T | 8.4% | intergenic (+10/‑198) | 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,377 | 0 | C | T | 8.4% | 130.4 / 2.4 | 48 | intergenic (+10/‑198) | PP_4304/sbp‑I | VIC family cation transporter/sulfate ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (21/22); new base T (2/2); total (23/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01 |
TACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGT > NC_002947/4896235‑4896524 | tACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GG < 2:316774‑M1/148‑1 (MQ=255) cgccgAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTC > 1:165568‑M1/1‑149 (MQ=255) aGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGC < 2:165568‑M1/149‑1 (MQ=255) cTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGG < 2:188011‑M1/149‑1 (MQ=255) cTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCCACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGG < 8:201690‑M1/149‑1 (MQ=255) tGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCG < 6:290280‑M1/145‑1 (MQ=255) tGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCG > 5:290280‑M1/1‑145 (MQ=255) ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCC‑GGCCTCTTCGCGGGTAAAcccgc > 2:319028‑M1/1‑144 (MQ=255) ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGc > 7:186547/1‑149 (MQ=255) gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGGTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTcc > 8:253164/1‑148 (MQ=255) caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCTCTGTGGGAATGTGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTAGCGGGTAATTACCCTACAACAGGCGCTGGCCGTgcgc > 5:228533/1‑149 (MQ=255) caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc > 7:108197/1‑149 (MQ=255) aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg < 8:108197/149‑1 (MQ=255) gATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCC‑GGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccg > 8:570‑M1/1‑118 (MQ=255) gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c < 8:264064/148‑1 (MQ=255) ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGG > 6:92311‑M1/1‑107 (MQ=255) ccGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGG < 5:92311‑M1/107‑1 (MQ=255) gACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGAGCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgta < 5:306220‑M1/149‑40 (MQ=255) aCCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCC‑GGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtag < 7:570‑M1/149‑41 (MQ=255) ccTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAgg < 4:145464/149‑1 (MQ=255) cacaAGAAGTCCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGGGCCTGCGCC‑GGCCTCTTCGCGGGTAAAcccgctccttcgggggacctcgcacggtccgcgcctgtaggagcgg < 6:258367‑M1/149‑47 (MQ=255) cacaAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGCCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg < 7:253164/149‑1 (MQ=255) acaAGAAGACCCACGAGCAGGGGGCGGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcg > 1:168079‑M1/1‑103 (MQ=255) acaAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcgga < 2:168079‑M1/149‑48 (MQ=255) agaagaCCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt > 8:186547‑M2/1‑146 (MQ=255) aagaCCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCC‑GGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttac < 1:319028‑M1/149‑53 (MQ=255) agaCCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAAc < 8:152291‑M1/97‑2 (MQ=255) agaCCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAAc > 7:152291‑M1/1‑96 (MQ=255) ccACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgc < 5:1442‑M2/148‑10 (MQ=255) cACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGGCCAATTCCCGGGTAAAACCCGACAAACCGGCCCGGACCGTGCCCGGACC‑‑‑CCCGTAGGAGCGGGTTTAccccc > 2:223560/1‑146 (MQ=255) cGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa > 8:96086/1‑109 (MQ=255) cGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa < 7:96086/109‑1 (MQ=255) cGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg < 8:174384/136‑1 (MQ=255) cGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg > 7:174384/1‑136 (MQ=255) cGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaag > 6:1442‑M2/1‑136 (MQ=255) gggggC‑GGCGTTCTGTTCGCGCTGGGGCAAGGGGTTGTTCCCCAAAGCGTGATAGCGCCTGCGCG‑GGCCTCTTGGCGGGTAAATCCGCTCCTACCGGCGGGGACGGTGCGCGTTCC‑‑‑CCCTTTGGAGCGGGTTTACCCGCGAAGAGGCCg < 2:308491/149‑1 (MQ=37) ttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCACTTAGCGGGTAAATCCGCTCCTAAAGGCGCGGACCGGGCGCGGTGC‑‑‑CCCGTAGGAGCGGGTTTAACCGCGAAGAGACCGGCACAgg > 1:269850/1‑146 (MQ=255) ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt < 7:35444‑M2/121‑4 (MQ=255) ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt > 8:35444‑M2/1‑118 (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 < 5:305389‑M2/93‑5 (MQ=255) ccAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta > 6:305389‑M2/1‑89 (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) | TACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGC‑GGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC‑GGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGC‑ATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGT > NC_002947/4896235‑4896524 |
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: ‑ |