| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 389,850 | G→C | 7.1% | A42P (GCC→CCC) | soxA → | sarcosine oxidase subunit alpha |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 389,850 | 0 | G | C | 7.1% | 181.0 / 4.5 | 56 | A42P (GCC→CCC) | soxA | sarcosine oxidase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base G (25/27); new base C (2/2); total (27/29) | |||||||||||
| 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.89e-01 | |||||||||||
CCGTCAAAGGCCAAGGGGAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCC > NC_002947/389702‑389991 | ccGTCAAAGGCCAAGGGGAAAAAGTAAGAGCCAGACCTATCGCCTCGCCAGCGGCGGACGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAACACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTgg < 5:276264/149‑1 (MQ=255) cAAGGGGAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTc < 7:217734/149‑1 (MQ=255) aGGGGAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGa < 6:181534/149‑1 (MQ=255) aGGGGAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGa < 8:1824/149‑1 (MQ=255) gggAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACa < 6:136318/149‑1 (MQ=255) gggAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACa < 1:96206/149‑1 (MQ=255) aaaaaGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACa < 4:314981/146‑1 (MQ=255) aaaaaGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACa > 3:314981/1‑146 (MQ=255) aaaGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTc < 3:171358/149‑1 (MQ=255) aaaGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCAGGCCAACGGCGTCGACATTGTc < 6:148644/149‑1 (MQ=255) aTGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCg > 6:296966/1‑149 (MQ=255) gAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCg > 4:98920/1‑147 (MQ=255) ccAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCt < 8:2844/149‑1 (MQ=255) aCCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCaa < 5:64128/149‑1 (MQ=255) tATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTc < 3:104052/144‑1 (MQ=255) tATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTc > 4:104052/1‑144 (MQ=255) aGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGCCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCAcg > 4:287955/1‑149 (MQ=255) aGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGCCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCAcg > 4:283099/1‑149 (MQ=255) cggcggCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCAcgcg < 2:68268/149‑1 (MQ=255) aTCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATcgc > 5:80519/1‑149 (MQ=255) aaGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGGCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCAACGCCGACGGtcgcg > 8:137799/1‑144 (MQ=255) ccTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCAcgcg < 3:132119/122‑1 (MQ=255) ccTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCAcgcg > 4:132119/1‑122 (MQ=255) cTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGGCGGTGACAGGCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACTTTGTCGGCCGTAGCTTGAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGc < 4:149053/148‑1 (MQ=255) cTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGcc < 3:98920/149‑1 (MQ=255) taCTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGGCCATCATCGCCGCCGGTACCGAAGAGCCGa < 1:21115/148‑1 (MQ=255) cTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGaa < 7:313742/139‑1 (MQ=255) cTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGaa > 8:313742/1‑139 (MQ=255) cTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGa > 2:21115/1‑146 (MQ=255) cTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACg > 8:243542/1‑149 (MQ=255) ttCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGAGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTAGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGc > 1:326022/1‑149 (MQ=255) tCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGc > 1:139342/1‑148 (MQ=255) cTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATcc < 3:179431/149‑1 (MQ=255) ttCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGccgccg < 4:312337/124‑1 (MQ=255) ttCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGccgccg > 3:312337/1‑124 (MQ=255) ttCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAAc < 3:317120/141‑1 (MQ=255) ttCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAAc > 4:317120/1‑141 (MQ=255) cTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCa < 8:242817/138‑1 (MQ=255) cTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCa > 7:242817/1‑138 (MQ=255) ccAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGc > 3:201552/1‑148 (MQ=255) aGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGc > 1:136943/1‑132 (MQ=255) aGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGCCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGaa < 3:287955/149‑1 (MQ=255) aGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGCCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGaa < 3:283099/149‑1 (MQ=255) gggTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAg > 8:116129/1‑149 (MQ=255) gCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATcgccgc < 4:231083/95‑1 (MQ=255) ccGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAg < 7:243542/149‑1 (MQ=255) aCAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTAc < 7:324345/94‑1 (MQ=255) aCAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTAc > 8:324345/1‑94 (MQ=255) aCAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATcccc > 2:304795/1‑149 (MQ=255) aGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCaa > 4:225180/1‑149 (MQ=255) ccTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAg > 4:149978/1‑126 (MQ=255) ccTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAg < 3:149978/126‑1 (MQ=255) ccgccgCTTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTgcgc < 7:169374/149‑1 (MQ=255) cgccgcGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTgcgcg > 6:45116/1‑149 (MQ=255) cgcgTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCa > 6:349226/1‑149 (MQ=255) gTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCAccc < 2:139342/149‑1 (MQ=255) ttGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCAccc > 2:20467/1‑148 (MQ=255) | CCGTCAAAGGCCAAGGGGAAAAAGTATGAGCCAGACCTATCGCCTCGCCAGCGGCGGCCGTATCGACCGCAGCAAGGTCCTGAACTTCACCTTCAACGGCAAGACCTACCAGGGTTATGCCGGTGACAGCCTGGCCGCCGCGTTGCTGGCCAACGGCGTCGACATTGTCGGCCGTAGCTTCAAGTACTCGCGCCCACGCGG‑CATCATCGCCGCCGGTACCGAAGAGCCGAACGCCATCCTGCAGATCGGCTCCAGCGAAGCTACCCAGATCCCCAACGTGCGCGCCACCC > NC_002947/389702‑389991 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |