| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,896,445 | G→A | 9.8% | 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 | 9.8% | 140.3 / 4.3 | 51 | 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 (23/23); new base A (2/3); total (25/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01 | |||||||||||
ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCG > NC_002947/4896297‑4896593 | acaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg > 7:170114/1‑149 (MQ=255) caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGt < 1:289679/149‑1 (MQ=255) caAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACATAGCGGGTAAATACGCAACTAACTGCGCGGACCCTGCCCGGTAC‑‑‑CCCGGAGGAGCGGGt > 8:271925/1‑149 (MQ=255) ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCcgc > 1:217585/1‑148 (MQ=255) gCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAgg < 2:16953/149‑1 (MQ=255) gggggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTCTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCTGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCg < 8:170114/149‑1 (MQ=255) gcgTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcagg > 2:217585‑M2/1‑145 (MQ=255) ggTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCAttttt < 5:247375/148‑1 (MQ=255) gcgcCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAaccc < 7:271925/149‑1 (MQ=255) ggCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGCAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAg > 6:259537/1‑149 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCAAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGtt < 3:228886/141‑1 (MQ=14) tctTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGATTTACCCGCAAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGtt > 4:228886/1‑141 (MQ=14) ttCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa > 8:101508/1‑129 (MQ=255) ttCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa < 7:101508/129‑1 (MQ=255) tCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTgg > 3:133609/1‑149 (MQ=255) cgcgGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgg < 4:81633‑M2/149‑71 (MQ=255) ggTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATc > 8:148013/1‑149 (MQ=255) gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtctt < 5:75825‑M2/149‑76 (MQ=255) aaTCCGCTCCTACAGGCCCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACTCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGc < 1:58502/148‑1 (MQ=255) aTCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttg < 6:196309‑M2/149‑80 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt < 3:133404‑M2/149‑81 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt > 4:133404‑M2/1‑69 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt < 1:270981‑M2/149‑81 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc > 1:214867/1‑148 (MQ=255) tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt < 6:91682‑M2/149‑81 (MQ=255) tCCGCGCCTGTAGGAGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc > 5:37674/1‑148 (MQ=255) gCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCaa > 8:6033/1‑148 (MQ=255) tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATaaa > 2:231366/1‑112 (MQ=255) tCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcag < 6:196117‑M2/149‑86 (MQ=255) aCAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAACCCCCCCCGTTATAGACGCTGGCAAACTGCCCCAACACCCtt > 6:63463/1‑149 (MQ=255) cAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgg > 2:270981‑M2/1‑59 (MQ=255) tacgggggACCGCGCACGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTc > 2:317440/7‑148 (MQ=255) aGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcggg < 6:187721‑M2/149‑92 (MQ=255) ggagcgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa < 6:37674/146‑1 (MQ=255) ggCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa < 7:6033/149‑1 (MQ=255) ggCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCaa < 7:277318/149‑1 (MQ=255) ggCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCa > 8:277318/1‑148 (MQ=255) gggggACCGCGCACGGTCC‑‑‑GCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTGTAGTCGCTGGCAATCTGCCCCAACACCCTTCAAg > 1:21553/4‑148 (MQ=39) cgGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGATTTACCCGCGAAGAGGCCGGCGcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgtggcaggtcgggcaatc > 5:91682‑M2/1‑53 (MQ=255) gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAac < 7:223347/134‑1 (MQ=255) gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAac > 8:223347/1‑134 (MQ=255) gACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAc > 2:60469/1‑147 (MQ=255) ccGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAgg < 7:148013/149‑1 (MQ=255) cgGTCC‑‑‑CCCGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAATCATTACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAAGAAGGACGCAAc < 6:24337/149‑1 (MQ=39) gTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGt > 8:120043/1‑149 (MQ=255) c‑‑‑cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGa > 8:53416/1‑148 (MQ=255) cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAc > 6:234106/1‑143 (MQ=255) cccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAAc < 5:234106/143‑1 (MQ=255) ccGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGaaaa < 2:214867/149‑1 (MQ=255) gTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGa > 1:174805/1‑134 (MQ=255) gTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGa < 2:174805/134‑1 (MQ=255) gAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAActct < 1:60469/149‑1 (MQ=255) ggTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCAcc > 3:41151/1‑149 (MQ=255) gTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCg > 6:142873/1‑149 (MQ=255) | ACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCCCAACACCCTTCAAGAACAAGGACGCAACGTGAAAAAACTCTTCACCG > NC_002947/4896297‑4896593 |
| 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: ‑ |