Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,554,502 | G→C | 11.4% | L25L (CTG→CTC) | PP_3139 → | group 1 family glycosyl transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,554,502 | 0 | G | C | 11.4% | 126.4 / 7.5 | 44 | L25L (CTG→CTC) | PP_3139 | group 1 family glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base G (18/21); new base C (2/3); 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 = 9.11e-01 |
GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGT > NC_002947/3554354‑3554648 | gCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTTGCGCTgg > 1:464524/1‑150 (MQ=255) aCAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGCCGCGAAATACTGGCGCTCCAGCGCCagtat < 1:157772/150‑5 (MQ=255) cTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGa > 1:452481/1‑150 (MQ=255) tgaagtgaagGTTGCCAGGGGAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAGAGATCTCGTCAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg < 1:61488/150‑1 (MQ=255) tgaagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg < 2:570496/150‑1 (MQ=255) gTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGta > 2:211332/1‑150 (MQ=255) aGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGccc > 1:700132/1‑150 (MQ=255) ttGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCagta < 2:317155/119‑4 (MQ=255) ttGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCagta > 1:317155/1‑116 (MQ=255) atcaAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAAGCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGACGTGCA‑ACGTATAGCCCTGCGTGGGTGGGAc < 1:473933/150‑1 (MQ=255) caAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGATGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGc > 1:61993/1‑150 (MQ=255) aaGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGcc < 1:757596/150‑1 (MQ=255) gCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg > 2:714630/1‑115 (MQ=255) gCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg < 1:714630/115‑1 (MQ=255) ggTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAg < 1:371455/104‑1 (MQ=255) ggTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAg > 2:371455/1‑104 (MQ=255) gTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACg < 2:25756/119‑1 (MQ=255) gTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACg > 1:25756/1‑119 (MQ=255) aTTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAg < 2:464524/150‑1 (MQ=255) aTTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGac > 2:832194/1‑149 (MQ=255) cTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACggg < 2:845960/150‑1 (MQ=255) cTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCaa > 1:254222/1‑150 (MQ=255) ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAg < 1:296475/129‑1 (MQ=255) ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGGGGGACGCCGAAGTGCAGGATGCCGAg > 2:296475/1‑129 (MQ=255) atcaaacaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGGTGAGGTGCA‑ACGTATCGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAcc > 2:269215/1‑150 (MQ=255) tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAGATCCTCGCGCTGGAACGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCAAAAGTGCAGGATGCGGAGGACGCAAGCGAACGGGGCAAGAccc > 2:122398/1‑150 (MQ=255) tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg < 1:481531/74‑1 (MQ=255) tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg > 2:481531/1‑74 (MQ=255) aatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGc < 1:678100/150‑1 (MQ=255) tcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGGATTTCAGGATGACGAGGACGCAAGCGAACGGGGCAAGACCCGCTa > 2:393912/1‑150 (MQ=255) aCCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAt < 2:700132/150‑1 (MQ=255) cccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAGGCGATGGGGGCAAGACCCGCTACGTAtt < 1:542885/150‑1 (MQ=255) ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCaa < 1:517512/134‑1 (MQ=255) ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCaa > 2:517512/1‑134 (MQ=255) aaaGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCGCGGGGGATTGTTTCAGGCGTGTGGCGGT‑CCTGGGTGGTGGGCGGTAGTGAAGGCTGAAGAGCACGAAATCGAACGGGGCAAGACCCGCTACGTATTGCaa < 1:832194/150‑1 (MQ=255) aGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTggg < 2:254222/150‑1 (MQ=255) caGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGCCGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTAGGTATTCCAATGTggagtcca > 2:370746/1‑144 (MQ=255) gTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAgg < 1:211332/150‑1 (MQ=255) cATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgtt < 1:103532/150‑1 (MQ=255) cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGa > 2:394666/3‑72 (MQ=255) cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGa < 1:394666/70‑1 (MQ=255) ggCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAg > 1:5527/1‑76 (MQ=255) ggCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAg < 2:5527/76‑1 (MQ=255) cgcTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAgg < 1:370746/150‑1 (MQ=255) gcTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGt > 2:760863/1‑150 (MQ=255) | GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCA‑ACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGT > NC_002947/3554354‑3554648 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |