Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 743,463 | A→G | 13.3% | intergenic (+139/+135) | PP_0635 → / ← PP_0636 | group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 743,463 | 0 | A | G | 13.3% | 113.5 / 8.8 | 46 | intergenic (+139/+135) | PP_0635/PP_0636 | group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (21/18); new base G (4/2); total (25/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01 |
TACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGC > NC_002947/743334‑743598 | tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGc < 2:532620/150‑1 (MQ=255) acaccaCCGTGCGTACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGcc > 2:45758/6‑150 (MQ=255) gAACCGTACGCACGGTGGTGTGAGAGGACGGCTGGTGTGAATCCGCCTCCTACTCGATTGCGC‑AATTCATCACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCggg < 1:216067/150‑1 (MQ=37) gAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCggg > 1:427496/1‑150 (MQ=255) acggacggTGGTGTGAGAGGGCGGCGGGTGGGAATCCGCCTCCTACTCGATTGGGCGAAGTCA‑GACCTGCCGTGTGTGCCCGCTATGGGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTTCCGGGGAAacta < 2:295356/150‑1 (MQ=38) aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg < 2:427496/150‑1 (MQ=255) ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGt > 2:536368/1‑150 (MQ=255) agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg > 1:403226/1‑150 (MQ=255) gACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGa < 2:403226/150‑1 (MQ=255) aCGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg > 2:610526/1‑150 (MQ=255) ggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCa > 2:178861/1‑150 (MQ=255) cggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga < 2:26671/150‑1 (MQ=255) aaTCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg < 1:45758/150‑1 (MQ=255) gCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGa > 2:83037/1‑82 (MQ=255) gCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGa < 1:83037/82‑1 (MQ=255) ctcctACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCg > 2:283279/1‑150 (MQ=255) cctACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGCCACACAGCGTCGTGACTAATCATCGAAACTTTTTGGTCGTGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt < 1:536368/150‑1 (MQ=255) tACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa > 1:780791/1‑150 (MQ=255) tCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCTTTACATAGCGGac < 2:373836/74‑1 (MQ=255) tCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGac > 1:373836/1‑74 (MQ=255) tCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtc < 1:239545/84‑1 (MQ=255) tCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtc > 2:239545/1‑84 (MQ=255) tCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGGAACTACTTGGTCGCGGCGGGTGGAGCAGAGAATGTAGGGGAGCGTTGACAg > 1:65765/1‑150 (MQ=255) gATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTg > 2:215063/1‑150 (MQ=255) gcgAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGtt > 1:306473/1‑150 (MQ=255) gAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg < 2:343668/126‑1 (MQ=255) gAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg > 1:343668/1‑126 (MQ=255) aaTCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa < 1:610526/150‑1 (MQ=255) a‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACGTTCGGTTGCGGTGCCGAGGTAACGCACAGCGTCGTGAGTTGAAGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta < 1:178861/150‑1 (MQ=255) a‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg > 2:286202/1‑79 (MQ=255) a‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg < 1:286202/79‑1 (MQ=255) a‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta < 1:172204/150‑1 (MQ=255) tGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtgt < 1:771838/99‑1 (MQ=255) tGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtgt > 2:771838/1‑99 (MQ=255) gacgctgtgtgTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGt < 1:135341/118‑1 (MQ=17) gacgctgtgtgTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGt > 2:135341/6‑123 (MQ=17) gacgctgtgtgTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTg > 2:576036/6‑150 (MQ=21) tccGCTATGTAACGCAACCGAAGGTTGCTTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAgaga < 2:146018/97‑1 (MQ=25) tccGCTATGTAACGCAACCGAAGGTTGCTTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAgaga > 1:146018/2‑98 (MQ=21) tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCgg > 1:253043/2‑143 (MQ=25) tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCgg < 2:253043/142‑1 (MQ=25) tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCt > 2:50618/2‑150 (MQ=37) tgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAgg > 2:753044/1‑149 (MQ=255) aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGa > 1:397198/1‑150 (MQ=255) gTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATcc > 1:589257/1‑150 (MQ=255) ttGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGa > 2:144979/1‑84 (MQ=255) ttGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGa < 1:144979/84‑1 (MQ=255) ttGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCt < 2:306473/150‑1 (MQ=255) aTAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc > 2:480579/1‑143 (MQ=255) aTAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc < 1:480579/143‑1 (MQ=255) | TACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCA‑GACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGC > NC_002947/743334‑743598 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |