Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,145,454 | A→G | 17.7% | intergenic (‑115/+26) | hemO ← / ← PP_1006 | heme oxygenase/heme receptor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,145,454 | 0 | A | G | 17.7% | 67.6 / 8.1 | 34 | intergenic (‑115/+26) | hemO/PP_1006 | heme oxygenase/heme receptor |
Reads supporting (aligned to +/- strand): ref base A (13/15); new base G (3/3); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01 |
GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTGC > NC_002947/1145307‑1145603 | gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCaaa < 1:304462/150‑1 (MQ=255) gcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAaga < 2:643803/150‑1 (MQ=255) ggggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGCCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAgg < 2:613256/149‑1 (MQ=255) aGGGCGGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAg > 2:673922/1‑150 (MQ=255) gggggAGTCCTTGGAAAGAGGGGGTCTGGATACAAGACGAAACAGCGGGAGGGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCTTCTTCGCGGGCGTGCCCTCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATcc < 2:149581/150‑1 (MQ=255) gggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCACCGATTAGGCCGGTACAGACGTAACGCCTCAGATATcc < 1:839752/150‑1 (MQ=255) gAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAgg < 2:720717/117‑1 (MQ=38) gAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAgg > 1:720717/1‑117 (MQ=38) ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑At > 1:585150/1‑150 (MQ=255) tGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑Att < 2:851829/150‑1 (MQ=255) gagagggggCGTGTAAAAACGAAACGAAACAGCGGTAACGAGACAGTAAAAAAAGCCGGTATGTACGGTGCCAGTACCGGCCTCTTCGCGTGCGTGCCTGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGa < 2:799264/144‑1 (MQ=255) aGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTgcgc < 1:673922/150‑1 (MQ=255) gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCAt > 2:721122/1‑150 (MQ=255) aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGc < 1:721122/150‑1 (MQ=255) gCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg > 2:144966/1‑150 (MQ=255) gCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt < 1:144966/150‑1 (MQ=255) gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGAc < 1:74621/80‑1 (MQ=255) gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGAc > 2:74621/1‑80 (MQ=255) gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTATTCGAGGGAGTGCCCGCAAAAAGGCCGTTACAGACGTAAATGCCCAAATATCCCAAACCAG‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt > 2:736539/1‑150 (MQ=255) gaaaGTCAAAAGAGCCGGTATTTGCTTTTCCGGTACCGGCCTCGTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGttt > 2:692735/4‑149 (MQ=255) gaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg < 2:577545/150‑1 (MQ=255) tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac < 2:52459/84‑1 (MQ=255) tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac > 1:52459/1‑84 (MQ=255) tctgTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGATAACAGTCGCCTCAccata > 2:819560/4‑147 (MQ=37) cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTgg > 1:572314/1‑150 (MQ=255) aGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGCAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGc < 2:572314/150‑1 (MQ=255) gCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa > 1:551277/1‑150 (MQ=255) gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc < 1:46624/105‑1 (MQ=255) gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc > 2:46624/1‑105 (MQ=255) ctctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCACACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc > 2:797911/1‑150 (MQ=255) tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt > 2:353389/1‑150 (MQ=255) ttCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAATCGACGAC‑AATTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc < 1:692735/150‑1 (MQ=255) tgcgcgGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc < 2:378370/148‑1 (MQ=255) gCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGAGTCACAgccgc < 2:38161/112‑1 (MQ=255) gCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGAGTCACAgccgc > 1:38161/1‑112 (MQ=255) aaaGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTGc > 1:238578/1‑150 (MQ=255) | GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCAC‑ATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCACTGC > NC_002947/1145307‑1145603 |
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: ‑ |