| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,896,445 | G→A | 21.0% | 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 | 21.0% | 78.0 / 12.5 | 38 | 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 (18/12); new base A (5/3); total (23/15) | |||||||||||
| 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.78e-01 | |||||||||||
ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACC > NC_002947/4896297‑4896592 | acaAGAAGACCCACTAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGTGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 3:241914/149‑1 (MQ=255)acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 4:146774/149‑1 (MQ=255)acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 3:81573/149‑1 (MQ=255)acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 3:335173/149‑1 (MQ=255)acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 3:323177/149‑1 (MQ=255)acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCggg < 3:167470/149‑1 (MQ=255) agaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacc > 5:182562‑M2/1‑143 (MQ=255) aCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGT‑‑‑CCGCGCCTGTAGGAGCGGAtttacccg > 8:355362‑M2/1‑141 (MQ=255) ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgc < 1:390910‑M2/148‑10 (MQ=255) gAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaaga > 2:390910‑M2/1‑135 (MQ=255) gttctCGCTGTGGGAATGGATTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGt > 4:12284‑M2/1‑116 (MQ=255) gttctCGCTGTGGGAATGGATTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGt < 3:12284‑M2/117‑2 (MQ=255) ttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCAACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGgcagcagca > 6:87948/1‑149 (MQ=255) aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagag > 6:73308‑M2/1‑102 (MQ=255) aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagag < 5:73308‑M2/117‑16 (MQ=255) ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttaccc < 1:199569‑M2/103‑8 (MQ=255) ttgttTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttaccc > 2:199569‑M2/1‑96 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTACGCG‑‑‑CCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGc > 8:139344/1‑149 (MQ=14) ctTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt > 2:137148/1‑149 (MQ=255) gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctattacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg > 7:370868‑M2/1‑54 (MQ=255) gggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcgcaggcgctattacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg < 8:370868‑M2/146‑93 (MQ=255) gggTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGT‑‑‑CCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGCAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAAt > 7:187851/1‑149 (MQ=14) ggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCCCCCGTTATAGTCGCTGGCAATc > 1:96790/1‑149 (MQ=255) ggTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGT‑‑‑CCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc > 6:112545/1‑149 (MQ=17) gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt > 7:289335/1‑149 (MQ=255) gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt > 7:222210/1‑149 (MQ=255) gTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGT‑‑‑CCGCGCCTGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGCAGGCGGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCt > 1:29141/1‑149 (MQ=14) tCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAATCCCACCCGTTATAGTCGCTGGCAATCTGccc > 2:236610/1‑148 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCAGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaccccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt < 7:124982‑M2/149‑103 (MQ=255) gCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtgg < 8:360562‑M2/148‑105 (MQ=255) gCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtgg < 2:391828‑M2/148‑105 (MQ=255) gCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCGCCaa > 7:165642/1‑148 (MQ=255) gCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa > 3:105833/1‑148 (MQ=255) cTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc > 6:317931/1‑129 (MQ=255) cTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTc < 5:317931/129‑1 (MQ=255) cTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACAcc < 1:137148/148‑1 (MQ=255) aCAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcg < 7:329085‑M2/149‑112 (MQ=255) ggCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa < 4:105833/149‑1 (MQ=255) aCCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGTGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAg < 8:289335/149‑1 (MQ=255) ggT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACg > 2:82768/1‑149 (MQ=255) ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt > 8:102272/1‑102 (MQ=255) ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCt < 7:102272/102‑1 (MQ=255) gAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatcatgctgcagctgctccccacgcat > 1:391828‑M2/1‑7 (MQ=255) aGCGGATTTACCCGCGAAGAGGCCGGCGCAGGCGGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCtt < 2:29141/149‑1 (MQ=255) ggTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCAcc < 8:165642/149‑1 (MQ=255) | ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGT‑‑‑CC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACC > NC_002947/4896297‑4896592 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |