Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,015,961 | A→T | 10.5% | 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 | 10.5% | 168.6 / 10.0 | 57 | intergenic (‑176/+58) | PP_5679/PP_4421 | hypothetical protein/aminotransferase |
Reads supporting (aligned to +/- strand): ref base A (28/23); new base T (4/2); total (32/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01 |
AATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGAG > NC_002947/5015819‑5016110 | aaTAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGAggg > 1:835214/1‑150 (MQ=255) gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑aa > 1:687107/1‑150 (MQ=255) gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑aa > 1:1346662/1‑150 (MQ=255) cctcACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑aaaa < 1:560076/147‑1 (MQ=255) ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTTTAGGCggg < 1:511791/143‑1 (MQ=255) ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTTTAGGCggg > 2:511791/1‑143 (MQ=255) gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAgg > 1:302369/1‑150 (MQ=255) tcCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAAAAAAAAggga > 2:816685/1‑149 (MQ=255) tcCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTa > 1:70759/1‑150 (MQ=255) ccGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAAt > 1:217828/1‑150 (MQ=255) ggCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAAtt > 1:1029679/1‑150 (MQ=255) gggggggCTTTTAATTTCCGGGCCTTGGAAAGTATCTGACCCGAAAGTACGTCCTTTGTAACTGATGAAAGGGTTTTGGGTTTTATTTGGAAAAAACTAAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAAtt < 1:816685/146‑1 (MQ=255) gggCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCt < 2:70759/150‑1 (MQ=255) gCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTgg < 2:302369/150‑1 (MQ=255) cgatctcAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGTGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTg < 1:1193439/144‑1 (MQ=255) cTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt < 1:1039083/150‑1 (MQ=255) ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt < 2:217828/150‑1 (MQ=255) cAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACATCGTCCTGGCTTTGATTGTGATAGCACGAC‑‑AAAAAAAAAATCCCGCCTAAAGGCGGGCTT‑ACAACGAAGAAGGTAATTAAGCTGGCCCCCTTGCCAAAGGTGATGGCTc > 2:242749/1‑150 (MQ=255) aGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAAAGAAGAGGGTAATTCAGCTGGCATCAATGCCCACGTTGaagatcg > 2:1193439/1‑144 (MQ=255) tCTGACCCGAAAGTACGTCCTTGTAAACTGATGAAAGCGTCTTGGCTTTGATTGTAATAGAACAAA‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGAc < 2:803845/150‑1 (MQ=255) ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAg > 2:738559/1‑150 (MQ=255) aaaGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTc > 2:899409/1‑145 (MQ=255) aaaGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTc < 1:899409/145‑1 (MQ=255) aattACGTCCTTGTTAACTAATGACTCCGTCTTGGCTTTGATTTTGATAGCACTA‑‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcat < 1:242749/147‑1 (MQ=255) ttGTGAACTGATGACAGCGTCTTGGCTTTGATAGTAATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGc > 1:1182203/1‑150 (MQ=255) tgAACTGATGACAGCGTCTTGGCTTTGATTGTCTTCGGATTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa < 2:1182203/150‑1 (MQ=255) agAACTAATAACAGCGTCTTGGCTTTGATTGAGATAAAAAAAA‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa < 2:835214/149‑1 (MQ=255) cTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTTTAGGCGGGc > 2:1101586/1‑68 (MQ=39) cTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTTTAGGCGGGc < 1:1101586/68‑1 (MQ=39) tgatgaCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGcaca > 1:510067/1‑150 (MQ=255) gcgtactGGCTGTGATGGAAATAGAACAAA‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtg < 2:181899/144‑1 (MQ=255) gCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtg > 2:723313/1‑150 (MQ=255) gCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtg > 2:484741/1‑150 (MQ=255) tCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTTCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGc > 2:515971/1‑150 (MQ=255) tCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGtt > 2:431500/1‑99 (MQ=255) tCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGtt < 1:431500/99‑1 (MQ=255) ttGTGATAGCACTAA‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTgatgat < 2:1346662/150‑1 (MQ=255) tAGCACTAC‑AAAAAAAAAAAGCCCGCCCAAAGGCGGGCTT‑AAAACGAAGAGGGGAATTCACCTGGCAACATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgatgagatcgga > 1:1230909/1‑143 (MQ=255) aCTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCAc > 2:377283/1‑150 (MQ=255) aaC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTc < 2:510067/149‑1 (MQ=255) aC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGGGTGCGATGCTGATGATGCGATCCACTTcc > 2:943729/1‑150 (MQ=255) aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGGGGGTAATTCAGCTGGCATTATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATTCACTTCCTg > 1:2294/1‑150 (MQ=255) aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtga > 2:1321119/1‑131 (MQ=255) aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtga < 1:1321119/131‑1 (MQ=255) aaaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGtt < 1:484741/150‑1 (MQ=255) aaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt > 2:622939/1‑72 (MQ=255) aaaaaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt < 1:622939/72‑1 (MQ=255) aaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCaa > 2:628753/1‑92 (MQ=255) aaaaaGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCaa < 1:628753/92‑1 (MQ=255) aaaGCCCGCCTAAAGGCGGGCTTAAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTc < 1:377283/150‑1 (MQ=255) gCCCGCCTTTAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGa > 1:476816/1‑150 (MQ=255) gCCCGCCTTTAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGa > 2:1248613/1‑150 (MQ=255) ccGCCTGAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACa > 2:206082/1‑150 (MQ=255) ccGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGc > 2:1287274/1‑91 (MQ=255) ccGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGc < 1:1287274/91‑1 (MQ=255) ccGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACa > 1:1386397/1‑150 (MQ=255) cGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACaa < 1:943729/150‑1 (MQ=255) aaGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAgaggag < 2:687107/150‑1 (MQ=255) | AATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAAAAAGCCCGCCTAAAGGCGGGCTT‑AAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGAG > NC_002947/5015819‑5016110 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |