Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,002,722 | C→A | 7.4% | P153P (CCC→CCA) | PP_4409 → | site‑specific recombinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,002,722 | 0 | C | A | 7.4% | 172.4 / 2.4 | 54 | P153P (CCC→CCA) | PP_4409 | site‑specific recombinase |
Reads supporting (aligned to +/- strand): ref base C (22/28); new base A (2/2); total (24/30) | |||||||||||
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.54e-01 |
GCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACGGGCCAGCTACAGTAGTCA > NC_002947/5002580‑5002871 | ccGAAGCTAAACCGCCCCCAATGGTAAGTACGGCATCCGCAAAAAAATATAACGCTCAACTCCAGCCTTTCCTCACGAGACTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAaacaa < 1:1199660/149‑1 (MQ=255) gAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAAct > 2:559629/1‑150 (MQ=255) gAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAAct > 1:1375217/1‑150 (MQ=255) tAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACtcttct < 1:200578/150‑1 (MQ=255) aCCGTTCACCATGGTACGCACGGCATCGTAAAAAAAATATCACGTTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTAt < 2:1375217/150‑1 (MQ=255) aCCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTAt < 2:461061/150‑1 (MQ=255) acGTACCCAAAGATAAGCACGGCATACTAAAAAAAATATCACTCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATc < 2:514917/149‑1 (MQ=255) gTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGa < 2:90842/150‑1 (MQ=255) ccAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGccc < 2:992680/126‑1 (MQ=255) ccAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGccc > 1:992680/1‑126 (MQ=255) ccAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTa < 1:385196/142‑1 (MQ=255) ccAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTa > 2:385196/1‑142 (MQ=255) ccAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGc > 1:119509/1‑150 (MQ=255) cAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGcc < 1:218634/150‑1 (MQ=255) cAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAATCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGcc > 2:525574/1‑150 (MQ=255) aaTGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCt < 2:1170253/150‑1 (MQ=255) ggTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCtgat < 2:582810/150‑1 (MQ=255) aaGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGaa > 2:1289864/1‑150 (MQ=255) tCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCct < 1:559629/150‑1 (MQ=255) aaaaaaTATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTAcc < 1:48702/150‑1 (MQ=255) aaTATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCa < 2:1205243/150‑1 (MQ=255) atatCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAg < 1:1289864/150‑1 (MQ=255) tatCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAgg < 2:119509/150‑1 (MQ=255) aCGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGaaa < 1:566012/114‑1 (MQ=255) aCGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGaaa > 2:566012/1‑114 (MQ=255) tCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGaaa < 1:1250329/92‑1 (MQ=255) tCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGaaa > 2:1250329/1‑92 (MQ=255) tCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgaga < 1:145037/150‑1 (MQ=255) aaGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgagaga > 2:280780/1‑150 (MQ=255) cAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGCGAGAGAGaa > 2:865375/1‑150 (MQ=255) ccTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGaaaaa < 2:1169604/150‑1 (MQ=255) tttCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCTTGATGAAAATCACCCCTACTACCATCAGGGTTCgaga > 1:1330962/1‑138 (MQ=255) tttCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCTTGATGAAAATCACCCCTACTACCATCAGGGTTCgaga < 2:1330962/138‑1 (MQ=255) tttCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGt < 2:1237624/150‑1 (MQ=255) ttCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTa < 1:865375/150‑1 (MQ=255) cTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTagca > 2:1178877/1‑150 (MQ=255) gagCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgaga > 1:97638/1‑129 (MQ=255) gagCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgaga < 2:97638/129‑1 (MQ=255) ccTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGca > 1:1140970/1‑150 (MQ=255) ccTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCAGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGca > 2:790869/1‑150 (MQ=255) ccAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCAGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGAt < 1:790869/150‑1 (MQ=255) aaaGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgag > 1:1060567/1‑100 (MQ=255) aaaGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCgag < 2:1060567/100‑1 (MQ=255) aaGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGa > 1:664734/1‑150 (MQ=255) tactGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCAc < 2:930777/71‑1 (MQ=255) tactGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCAc > 1:930777/1‑71 (MQ=255) aCGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCaaa < 2:1140970/150‑1 (MQ=255) ggCTGGTAATTGAAAAGCCAGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGc > 1:394503/1‑105 (MQ=255) ggCTGGTAATTGAAAAGCCAGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGc < 2:394503/105‑1 (MQ=255) gCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACg > 1:433316/1‑150 (MQ=255) tGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACggg > 2:266492/1‑150 (MQ=255) ggTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACGGGc < 2:664734/150‑1 (MQ=255) aTTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACgg < 2:720019/144‑1 (MQ=255) aTTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACgg > 1:720019/1‑144 (MQ=255) aaaaGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACGGGCCAGCTACa > 1:1044885/1‑150 (MQ=255) aaaaGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACGGGCCAGCTACa > 2:351213/1‑150 (MQ=255) ccGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAagc < 1:965331/94‑1 (MQ=255) ccGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAagc > 2:965331/1‑94 (MQ=255) cGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGAACGCCGTCTCGATATGTAATCAAACACGGGCCAGCTACAGTAGTCa > 2:825304/1‑150 (MQ=255) | GCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCCCGAAACAACTCTTCTTATCGAAAGCCTGATGAAAATCACCCCTACTACCATCAGGGTTCGAGAGAGAGAAAAAGTAGCAGCCAGGCAACAAGCAGCACCGATGACAGCTGCGACTGGAAATAATCAAACACGGGCCAGCTACAGTAGTCA > NC_002947/5002580‑5002871 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |