Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,765,933 | C→T | 20.0% | R345H (CGC→CAC) | PP_2418 ← | cobalamin ABC transporter substrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,765,933 | 0 | C | T | 20.0% | 40.3 / 4.8 | 20 | R345H (CGC→CAC) | PP_2418 | cobalamin ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (9/7); new base T (2/2); total (11/9) | |||||||||||
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.75e-01 |
CCGGCCAGCAGGGCGGCCAACAGAGCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGAGCTTTTTCTGGTAGTCGTCGGCGGGGAAA > NC_002947/2765808‑2766061 | ccGGCCAGCAGGGCGGCCAACAGAGCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGc > 1:507215/1‑136 (MQ=255) agcagGGCGGCCAACAGAGCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGTGGATGCCGGCCTGCAt > 1:285039/1‑136 (MQ=255) agcagGGCGGCCAACAGAGCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGTGGATGCCGGCCTGCAt > 1:212682/1‑136 (MQ=255) agCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGa < 2:507215/136‑1 (MQ=255) ccAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAatg < 2:280074/136‑1 (MQ=255) ccAGGACAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAatg < 2:211914/136‑1 (MQ=255) gggCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAatgatg < 2:390748/136‑1 (MQ=255) tcatcaCGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCgg > 2:409246/1‑136 (MQ=255) tcatcaCGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCATCACCGCGTCCATGATGCgg > 1:268454/1‑136 (MQ=255) catcaCGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCggg < 2:44040/136‑1 (MQ=255) ttACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTgc > 2:533722/1‑136 (MQ=255) tACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTgcg > 1:568831/1‑136 (MQ=255) tGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGaa < 2:534267/136‑1 (MQ=255) gAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGGTACCCGGTCACTGCGCCGGAATTCg > 2:542239/1‑136 (MQ=255) acggCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGa < 1:385472/134‑1 (MQ=255) ggACAGCGTCTGCACGCCGCGGAACAGGTGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGAGCTTTTTc < 2:285039/136‑1 (MQ=255) ggACAGCGTCTGCACGCCGCGGAACAGGTGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGAGCTTTTTc < 2:212682/136‑1 (MQ=255) ggAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTTCCGGGTCACTGCGCAGGAATTCGAGCTTTTTCTGGTAGTCGTCGGCGGGGaa > 2:291276/1‑136 (MQ=255) ggAACAGGCGGATGCCGGCCTGCATCGCGTCGGCGTCGATGATAATGATGCTGTCGTGCTTCACCTCGTCAATGATGCGGGTTACCGGGTGAAAGCTCACGAATTCAAGCGTTATATGGTAATCGTcggcgttgaa > 2:120315/1‑131 (MQ=255) gAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGAGCTTTTTCTGGTAGTCGTCGGCGGGGaaa > 2:61520/1‑136 (MQ=255) | CCGGCCAGCAGGGCGGCCAACAGAGCGGCCAGGGCAACCGCGATGCAGGGTAACAGGCGGATCATCACGGGGCTGCTTTACCGCTGGCGAAGGCCGTGGACAGCGTCTGCACGCCGCGGAACAGGCGGATGCCGGCCTGCATGGCGTCGGCGTCGAGGATAATGATGCGGTCGTGCTTCACCGCGTCCATGTTGCGGGTTACCGGGTCACTGCGCAGGAATTCGAGCTTTTTCTGGTAGTCGTCGGCGGGGAAA > NC_002947/2765808‑2766061 |
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: ‑ |