Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,650,019 | A→T | 9.1% | intergenic (+1496/+195) | PP_3214 → / ← PP_3216 | hypothetical protein/cysteine dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,650,019 | 0 | A | T | 9.1% | 145.1 / 5.5 | 55 | intergenic (+1496/+195) | PP_3214/PP_3216 | hypothetical protein/cysteine dioxygenase |
Reads supporting (aligned to +/- strand): ref base A (24/26); new base T (2/3); total (26/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 = 6.44e-01 |
GGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGG > NC_002947/3649880‑3650165 | ggCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGa < 5:175172/148‑1 (MQ=255) aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg < 1:286024/148‑3 (MQ=255) aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg > 2:286024/1‑146 (MQ=255) tCACGAAACAACTCGCCGCCCGCAGCGCACCAGCCGGCAAGGCCGCTATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCccg < 6:210983/149‑1 (MQ=255) gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg > 3:298449/1‑149 (MQ=255) gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg > 2:28997/1‑149 (MQ=255) aaaCAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCggg > 7:225120/1‑149 (MQ=255) aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc > 3:179349/1‑149 (MQ=255) aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc > 3:32290/1‑149 (MQ=255) cTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgc > 4:232681/1‑149 (MQ=255) tCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgcg > 3:260309/1‑149 (MQ=255) cgccgcCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg > 1:279982/1‑149 (MQ=255) gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAGACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGCGTGCGGGGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg > 1:235913/1‑149 (MQ=255) ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGc < 1:28997/149‑1 (MQ=255) ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGc > 3:89724/1‑149 (MQ=255) aGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGAttt < 4:32290/149‑1 (MQ=255) gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAATGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTTCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTc < 6:86345/149‑1 (MQ=255) cgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTcc < 8:191895/149‑1 (MQ=255) gcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGc > 2:36064/1‑149 (MQ=255) gCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCa < 4:26848/149‑1 (MQ=255) ggCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGtt > 2:183131/1‑149 (MQ=255) gCAAGGCCGCCATCAAGCAAGGCA‑ACATTGAAGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGAGGTGTCCCGCCGCGGAACGCGGGCTAAGATTTCCAGCTCCAGttt > 4:30051/1‑149 (MQ=255) cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc < 4:89724/149‑1 (MQ=255) cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGCCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc < 5:278516/149‑1 (MQ=255) gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg > 6:145931/1‑91 (MQ=255) gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg < 5:145931/91‑1 (MQ=255) gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACa < 3:232681/149‑1 (MQ=255) cgcgccATCAAGCAAGGCT‑ACATTGACTTAGCCCCAGGCCAGTAAACGCTGAGCGGCGTCTGCGGCGA‑GCCCCTCACCGGCATCGTAGAGTGTGAGTGCGGTGTCCCGCGGGGGAACGCGGGCGAAGATTTCCAGCTGAAGTTTCgg > 8:216217/3‑146 (MQ=255) cgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAgg > 1:54475/1‑149 (MQ=255) gccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAggg < 6:86639/149‑1 (MQ=255) aTCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTa > 8:252365/1‑149 (MQ=255) caagGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTgg < 1:183131/149‑1 (MQ=255) ggCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGGGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCcg > 6:21361/1‑149 (MQ=255) aCATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa < 1:36064/149‑1 (MQ=255) aCGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCg > 7:54203/1‑148 (MQ=255) gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg < 3:74816/148‑1 (MQ=255) ccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGcc > 2:312897/1‑149 (MQ=255) aaGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTg < 8:54203/149‑1 (MQ=255) aGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCACCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTgg < 5:309010/149‑1 (MQ=255) gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc < 7:252365/149‑1 (MQ=255) gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc < 2:279982/149‑1 (MQ=255) gTAAACGCTGCGCGGCGACTGCGGCGACGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGa > 1:14704/1‑148 (MQ=255) tGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc > 4:187597/1‑149 (MQ=255) cgcgGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc < 1:164535/147‑1 (MQ=255) cgACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGt < 2:54475/149‑1 (MQ=255) ttcggcgA‑GCCCCTCACCGTCATCGTAGAAGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGt < 3:30051/147‑1 (MQ=255) ggcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGa < 3:187597/149‑1 (MQ=255) gcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGAt < 8:225120/149‑1 (MQ=255) gA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGACGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAATGGTCTTCTTCGCGTGAGTCGAt > 2:304227/1‑147 (MQ=255) gCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg < 2:235913/149‑1 (MQ=255) ccTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGcc < 1:312897/149‑1 (MQ=255) acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa > 6:193638/2‑135 (MQ=255) acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa < 5:193638/134‑1 (MQ=255) gTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcct > 7:101060/1‑149 (MQ=255) gTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGGCTTCTTCGCGTGCGTCGATCAGGGCCAGCTCCTCCTTgg < 1:304227/149‑1 (MQ=255) | GGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGG > NC_002947/3649880‑3650165 |
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: ‑ |