Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,015,961 | A→T | 25.0% | intergenic (‑176/+58) | PP_5679 ← / ← PP_4421 | hypothetical protein/aminotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,015,961 | 0 | A | T | 25.0% | 31.6 / 9.4 | 20 | intergenic (‑176/+58) | PP_5679/PP_4421 | hypothetical protein/aminotransferase |
Reads supporting (aligned to +/- strand): ref base A (5/10); new base T (2/3); total (7/13) | |||||||||||
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.67e-01 |
AGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGT > NC_002947/5015823‑5016099 | aGCCAGCACTGGTCTCCAGCGTTGGTTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTTTAGGCGGGCttt < 1:566940/150‑2 (MQ=255) ccGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGTTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTaa < 1:579091/150‑1 (MQ=255) ccGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGAAATAAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTa < 1:395650/150‑1 (MQ=255) ttttACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGca > 2:193057/1‑150 (MQ=255) cTTCCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACCGATGACAGCTTCTTGGCTTTGATTTTGATAGCACTAA‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt < 1:79173/150‑1 (MQ=255) ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGa > 2:382721/1‑150 (MQ=255) gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg < 1:479213/150‑1 (MQ=255) cacGAAAGTATGTACTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGAACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGTGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCa < 2:2297/148‑1 (MQ=255) cGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg < 1:193057/150‑1 (MQ=255) ccTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAACGCCCGCCTCAAGGATGGTTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGAc < 1:382721/150‑1 (MQ=255) ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGc > 1:390529/1‑75 (MQ=255) ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGc < 2:390529/75‑1 (MQ=255) gaCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCgg < 1:59308/61‑1 (MQ=255) gaCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCgg > 2:59308/1‑61 (MQ=255) aTAGCACTAC‑AAAAAAAAAAAGCCCCCCTAAAGGCGGGCTTAAAAAGAAGAGGGAAATTCACCTGGAATCATTGCCAACGTTTATGGCTCAGACCTTCAGCTCATCAAAGACCCTCAGCACCGCTTCGGGTGCGATGCTGATGTTGCGAt > 2:80660/1‑150 (MQ=255) aaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGt < 1:80660/150‑1 (MQ=255) gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGa > 2:174122/1‑90 (MQ=255) gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGa < 1:174122/90‑1 (MQ=255) gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa > 1:31294/1‑96 (MQ=255) gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa < 2:31294/96‑1 (MQ=255) | AGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGT > NC_002947/5015823‑5016099 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |