Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,908,345 | A→G | 17.4% | intergenic (‑276/‑6) | PP_4317 ← / → PP_4318 | hypothetical protein/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,908,345 | 0 | A | G | 17.4% | 42.2 / 4.2 | 23 | intergenic (‑276/‑6) | PP_4317/PP_4318 | hypothetical protein/transposase |
Reads supporting (aligned to +/- strand): ref base A (10/9); new base G (2/2); total (12/11) | |||||||||||
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.29e-01 |
GTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCCGGAGTGGGTCGATCAGGTTTTCGAGGAGCACCGCCAACGG > NC_002947/4908197‑4908425 | gTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGCAGAAGCGGCCAAAGGGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAG > 2:749088‑M1/1‑150 (MQ=255) gTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAG > 1:600042‑M1/1‑150 (MQ=255) ttttGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAAt > 1:735278/1‑150 (MQ=32) tttGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACtt > 1:1412535/1‑150 (MQ=32) tGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGc < 2:1304009/150‑1 (MQ=32) gACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGcc < 2:600042/150‑1 (MQ=32) gACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGcc > 1:652432/1‑150 (MQ=32) aCAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGc > 1:882017/1‑139 (MQ=33) aCAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGc < 2:882017/139‑1 (MQ=33) aaGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGCCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCAAGCGTAATGGCCAAACTTGCCTTGGAg < 1:193218/150‑1 (MQ=17) tCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCa < 2:652432/150‑1 (MQ=11) ggTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccattacgctggcgg < 2:144961‑M1/137‑43 (MQ=255) ggTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccattacgctggcgg > 1:144961‑M1/1‑95 (MQ=255) gAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCACTCCGGCGcaatggcctgctccaaggcaagtttggccattacgctggcgggtgcttttttctcgaatcgcgc > 2:1221758‑M1/1‑86 (MQ=255) gCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGtaatggcctgctccaaggcaagtttggccattacgctggcgg > 1:783532‑M1/1‑78 (MQ=255) gCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGtaatggcctgctccaaggcaagtttggccattacgctggcgg < 2:783532‑M1/120‑43 (MQ=255) tCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCACTCCGGCGcaatggcctgctccaaggcaagtttggccattacgctggcgggtgcttttttctcgaatcgcgctagtacatcagcc < 1:1221758‑M1/150‑78 (MQ=255) aGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccatt < 1:192003‑M1/100‑33 (MQ=255) aGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccatt > 2:192003‑M1/1‑68 (MQ=255) gTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccattacgctggcggg < 1:1302139‑M1/110‑44 (MQ=255) gTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGcaatggcctgctccaaggcaagtttggccattacactggcggg > 2:1302139‑M1/1‑67 (MQ=255) tgGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCACTCCGGCGcaatggcctgctccaaggcaagtttggccattacgctggcgggtgctttt < 1:1103639‑M1/98‑51 (MQ=255) tgGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCACTCCGGCGcaatggcctgctccaaggcaagtttggccattacgctggcgggtgctttt > 2:1103639‑M1/1‑48 (MQ=255) | GTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACC‑‑CGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGGCCATTGCGCCGGAGTGGGTCGATCAGGTTTTCGAGGAGCACCGCCAACGG > NC_002947/4908197‑4908425 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |