Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,180,012 | C→G | 9.0% | intergenic (‑240/‑43) | PP_4557 ← / → PP_4558 | hypothetical protein/GNAT family acetyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,180,012 | 0 | C | G | 9.0% | 151.3 / 4.9 | 53 | intergenic (‑240/‑43) | PP_4557/PP_4558 | hypothetical protein/GNAT family acetyltransferase |
Reads supporting (aligned to +/- strand): ref base C (22/26); new base G (2/3); total (24/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.54e-01 |
GTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTA > NC_002947/5179868‑5180159 | gTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCctct < 5:132169/149‑1 (MQ=255) tCCCGCGTCGTACTGCCAGGGGGGATGAACAAACTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTCGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga < 5:241774/149‑3 (MQ=255) cGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGtt > 1:2490/1‑148 (MQ=255) cTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTACGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCa < 1:253740/149‑1 (MQ=255) cGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAgg < 4:12300/148‑1 (MQ=255) tACTGCCTGGGGGTATGAACAAAATTCTACTCCCGGCAAGCAATAGTTTAACCCCGTCGTCGGATATTCACCCTTGCCGGTTTTGCAAAAATTTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTa < 1:93952/149‑1 (MQ=255) aCTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACACTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCTCAAATCATCAGGCCACAGCCTATGGCCTCGTTTCAGCGTAg > 6:101701/1‑149 (MQ=255) ccAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt < 4:111575/135‑1 (MQ=255) ccAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt > 3:111575/1‑135 (MQ=255) aGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCa > 8:165366/1‑149 (MQ=255) aaCAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTg < 8:136594/148‑1 (MQ=255) cAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGAGACAAAGTTGCCGCAAATCATCAGGCCAGACCCTCTTGCCACCTTTCAGGGTAGGGC‑TCACCAATTCTGaa > 3:115887/1‑148 (MQ=255) aaaaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGCCCAGAGCCTCTGGCCACGTTTCAGGGAAGGGC‑TCACCATTTCTGAAc > 7:44065/1‑148 (MQ=255) ttCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAg < 8:168111/121‑1 (MQ=255) ttCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAg > 7:168111/1‑121 (MQ=255) tACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga < 1:214605/112‑3 (MQ=255) tACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga > 2:214605/1‑110 (MQ=255) cTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCAGCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGGCGCAAATGATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTa < 8:139848/146‑1 (MQ=255) cTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTa < 2:225092/146‑1 (MQ=255) ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGtt < 4:174562/148‑1 (MQ=255) ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTc > 6:295845/1‑149 (MQ=255) ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTc > 2:247822/1‑149 (MQ=255) gCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGAc > 7:235198/1‑149 (MQ=255) ccAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGAc < 8:235198/148‑1 (MQ=255) aaTAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTc < 2:37339/149‑1 (MQ=255) aaTAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTc > 7:43947/1‑149 (MQ=255) gTTGAACCGCTTAGTGGGATTTTCACCCTTGGCGGTTTTGCAAGAATCTGCGACTAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTTGTGTTCGACGCTTTCGGa < 4:115887/148‑1 (MQ=255) ttGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGaca < 5:295845/149‑1 (MQ=255) ccccGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGtt > 2:268123/1‑129 (MQ=255) cccGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAAc > 1:95061/1‑149 (MQ=255) cccGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAAc < 8:43947/149‑1 (MQ=255) cGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGc > 1:319940/1‑149 (MQ=255) gTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGtt < 1:136908/149‑1 (MQ=255) tCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGttt > 1:222235/1‑149 (MQ=255) cccTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACt > 5:9118/1‑149 (MQ=255) ccTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCTCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACtt < 2:95061/149‑1 (MQ=255) ccGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTAc < 2:222235/148‑1 (MQ=255) cGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCg > 5:87034/1‑149 (MQ=255) gTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGc > 4:222616/1‑148 (MQ=255) ttGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAg < 3:222616/148‑1 (MQ=255) ttGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAg < 2:319940/148‑1 (MQ=255) acaaTCTTCGACAAAGTGGCCGCAAATCATGTGGCGAGAGCCTCTGGCATCGTTTCAGGTTAGGGC‑TCAGCATTTCTGAACTTTTGTTCTGTTCGACGCTTTCGGACAGAGATCGCCGTTTGCTGGAGCACTTGTACCGCCAGCACAgg < 8:44065/149‑2 (MQ=255) caaTCTGCGACAAATTGTCCGCAAATCATCAGGCCAGAGCCTCTTGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGATCAATTGTACTGTTCGACGCTTTCGGACACAGAACTCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCt < 5:101701/149‑1 (MQ=255) aCAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgc < 6:87034/149‑1 (MQ=255) tcatcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGc > 5:147429/1‑149 (MQ=255) aaacgaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGc > 2:173190/6‑149 (MQ=255) aGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc > 1:293494/1‑149 (MQ=255) ggCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCGTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc > 1:220618/1‑149 (MQ=255) ggCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACt > 3:318995/1‑149 (MQ=255) cAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgtg < 2:293494/149‑1 (MQ=255) agagCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg < 4:318995/149‑1 (MQ=255) agagCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg < 2:2490/149‑1 (MQ=255) gCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTa < 1:247822/149‑1 (MQ=255) | GTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGC‑TCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTA > NC_002947/5179868‑5180159 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |