Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,975,345 | G→T | 20.0% | intergenic (‑109/+285) | fliC ← / ← PP_4379 | flagellin/3‑oxoacyl‑ACP synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,975,345 | 0 | G | T | 20.0% | 49.1 / 2.5 | 20 | intergenic (‑109/+285) | fliC/PP_4379 | flagellin/3‑oxoacyl‑ACP synthase |
Reads supporting (aligned to +/- strand): ref base G (8/8); new base T (2/2); total (10/10) | |||||||||||
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.75e-01 |
GTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTG > NC_002947/4975213‑4975477 | gTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGtt > 1:114824/1‑136 (MQ=255) ggTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCaaaaaaa > 1:519472/1‑136 (MQ=255) aaaGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAAt < 1:120646/136‑1 (MQ=255) tGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGGTTGATGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTgg > 1:71777/1‑135 (MQ=255) tGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTga > 1:545520/1‑136 (MQ=255) ggATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAAt < 1:178243/136‑1 (MQ=255) gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTGGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATg < 2:1041/108‑1 (MQ=255) gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTa > 1:211171/1‑136 (MQ=255) gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTa > 1:473862/1‑136 (MQ=255) gCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATg > 1:1041/1‑108 (MQ=255) ccTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGttt > 1:311625/1‑136 (MQ=255) aagatAACCTTAGTCATGGTTTTCGTCGGGGGCGCGGGTTGCTTGTGGGTTATTTCAAATAAAATATTCTGGCATGGGGCCTGGCCTGAGAAATCAAGGGGTTTATGTCGGTTTATTGGCGTCATAAACTTTTTGa < 2:405706/135‑1 (MQ=255) tGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTTGGCTGAATGCCt < 2:114824/136‑1 (MQ=255) aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCa > 2:511836/1‑103 (MQ=255) aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCa < 1:511836/103‑1 (MQ=255) aTCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATcc > 1:138878/1‑136 (MQ=255) gggTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCgg < 2:472492/136‑1 (MQ=255) gCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAg < 1:447023/136‑1 (MQ=255) gAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTg < 2:211171/136‑1 (MQ=255) gAGTGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTg < 2:473862/136‑1 (MQ=255) | GTTGGTGTTAACGGTCAAAGCCATGACGAATTCCTCGTTGTATGGATACTGCGGCTTCCGGCCCTGGCGACAGCCGGGTGTGGCCTAGAGAACCTTCGTAATGGTTATCGTCGGGTGCGCGGGTTGCTTGAGGGTTTTTTCAAAAAAAAAATACTGGCATGGGGCCAGCCCTGAGAAATCAAGGGGTTAAGGGCGGTTTATTGGCGTCAAAAACTGTTTGAACTGGGCTGAATGCCTCATCCCGTCGGTGGAGCGGGTGAGGCTG > NC_002947/4975213‑4975477 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |