Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,981,093 | C→A | 8.6% | intergenic (+32/‑268) | rpsA → / → ihfB | 30S ribosomal protein S1/integration host factor subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,981,093 | 0 | C | A | 8.6% | 175.4 / 7.0 | 58 | intergenic (+32/‑268) | rpsA/ihfB | 30S ribosomal protein S1/integration host factor subunit beta |
Reads supporting (aligned to +/- strand): ref base C (29/24); new base A (2/3); total (31/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01 |
GAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC > NC_002947/1980947‑1981235 | gAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCgg < 3:125188/149‑1 (MQ=255) aaaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCggg > 2:352971/1‑149 (MQ=255) gacgacGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGCAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTCTTTGATCGGTAAAAAGGGCGACCCTCGGGTCg > 4:166583/1‑149 (MQ=255) acgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGc > 6:286375/1‑149 (MQ=255) acgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCC‑t < 7:130335/149‑1 (MQ=255) agagCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑tttttttt > 4:279329/1‑149 (MQ=255) cgcgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTc < 8:33885/149‑1 (MQ=255) cgcgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTc < 4:304747/149‑1 (MQ=255) ccAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAggt > 5:206140/1‑149 (MQ=255) ccTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTgcg > 8:120157/1‑149 (MQ=255) gCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCa < 1:352971/149‑1 (MQ=255) gAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgctg < 4:120302/149‑1 (MQ=255) gctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAgg > 5:192096/1‑120 (MQ=255) gctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAgg < 6:192096/120‑1 (MQ=255) ccGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc < 3:166583/147‑1 (MQ=255) cGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAg < 7:120157/149‑1 (MQ=255) accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑tttttt > 8:276878/1‑95 (MQ=255) accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑tttttt < 7:276878/95‑1 (MQ=255) ccaccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgct > 1:279407/1‑132 (MQ=255) ccaccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgct < 2:279407/132‑1 (MQ=255) ccaccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCtt > 6:14598/1‑147 (MQ=255) caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAg > 6:330653/1‑149 (MQ=255) caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAg > 6:283001/1‑149 (MQ=255) gctgcgctgcTGCGCGAAGCAATGACCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc < 7:341119/111‑1 (MQ=255) gctgcgctgcTGCGCGAAGCAATGACCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc > 8:341119/1‑111 (MQ=255) tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa > 2:139594/1‑149 (MQ=255) tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCg > 6:175394/1‑147 (MQ=255) cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGc < 3:312616/149‑1 (MQ=255) cgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGGGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGc < 3:279329/149‑1 (MQ=255) gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGcc > 5:132722/1‑149 (MQ=255) tgctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCg < 5:286375/149‑1 (MQ=255) gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGt > 7:221921/1‑149 (MQ=255) gcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg < 5:14598/149‑1 (MQ=255) cAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa > 2:230802/1‑132 (MQ=255) cAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa < 1:230802/132‑1 (MQ=255) aTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt < 6:206140/149‑1 (MQ=255) gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTc < 1:139594/148‑1 (MQ=255) aGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcg > 3:325879/1‑149 (MQ=255) aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgcg > 3:126453/1‑147 (MQ=255) aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCTGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAg > 3:254340/1‑149 (MQ=255) aCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGc < 5:175394/149‑1 (MQ=255) tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa < 8:181324/145‑1 (MQ=255) tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa > 7:181324/1‑145 (MQ=255) tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTGGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg > 1:88555/1‑149 (MQ=255) aGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGgc < 2:88555/149‑1 (MQ=255) tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCTGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt < 3:349603/125‑1 (MQ=255) tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt > 4:349603/1‑125 (MQ=255) gATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc < 8:221921/148‑1 (MQ=255) aaaaGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAg > 8:306731/1‑149 (MQ=255) aaaaGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAg < 5:330653/149‑1 (MQ=255) aaaGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa > 2:165264/1‑134 (MQ=255) aaaGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc > 2:208864/1‑149 (MQ=255) aaaGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCCCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCa < 1:165264/134‑1 (MQ=255) aGGGCGACCCGAGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTAt > 6:162711/1‑141 (MQ=255) aGGGCGACCCGAGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTAt < 5:162711/141‑1 (MQ=255) gggCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGGACTGCTGGTACACGCCAGCTTCAGGATGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCGCTGGGCTATCTAAGCgg > 8:115748/1‑148 (MQ=255) ggCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg > 3:143659/1‑147 (MQ=255) aCCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAc > 3:18713/1‑148 (MQ=255) | GAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGAC > NC_002947/1980947‑1981235 |
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: ‑ |