Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,338,538 | A→T | 13.3% | intergenic (+114/+28) | PP_3812 → / ← PP_5632 | lipase/allo‑isoleucine metabolism protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,338,538 | 0 | A | T | 13.3% | 71.9 / 5.7 | 30 | intergenic (+114/+28) | PP_3812/PP_5632 | lipase/allo‑isoleucine metabolism protein |
Reads supporting (aligned to +/- strand): ref base A (11/15); new base T (2/2); total (13/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01 |
GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTA > NC_002947/4338395‑4338676 | ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCcgc < 2:125610/150‑1 (MQ=255) gCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCcgcg > 1:506638/1‑150 (MQ=255) aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGaa > 2:347788/1‑150 (MQ=255) aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGaaa < 2:506638/150‑1 (MQ=255) cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTTTACCCGCGAAAg > 2:427700/1‑150 (MQ=255) cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAg < 1:347788/150‑1 (MQ=255) aCTTAAAGGGCATGTGCGGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTTAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc < 2:415454/150‑1 (MQ=255) ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg < 1:427700/150‑1 (MQ=255) cGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaa < 1:334565/150‑1 (MQ=255) cAATAGTCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAGCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc < 1:7407/150‑1 (MQ=255) aaTAGACCCGCCTGGGCGGGTCTTTTGCGTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATGCTGCTCAGGGCTCAAACCCTTcc > 1:395589/1‑150 (MQ=255) gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGCATCTGCC‑TCGAT‑CTGCTGTTTGCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc < 1:391510/150‑1 (MQ=255) gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc > 2:350593/1‑150 (MQ=255) cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCTTCGCCAcc > 1:452935/1‑150 (MQ=255) ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCa < 2:69452/106‑1 (MQ=255) ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCa > 1:69452/1‑106 (MQ=255) tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg < 1:350593/150‑1 (MQ=255) ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCCCAAACCCTTCCACAATGAGCACCTCAGCCTTGGCCACCCCTGGGCGCTg > 1:380162/1‑150 (MQ=255) cAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt > 1:355448/1‑105 (MQ=255) cAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt < 2:355448/105‑1 (MQ=255) ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGc < 1:11977/150‑1 (MQ=255) aaTGAATCTGCCAGCGGTGCGGCAGGT‑‑CACGTCATTCGGCGTGTAAACCTGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACa < 1:494196/150‑1 (MQ=255) gCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAGTCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc > 1:392496/1‑150 (MQ=255) ccAGCGGTGCGGCAGGC‑‑CCGGCCCTTTCGCGCGTAAACCCGCGCGAGGAGCAGGTAAGCCTGCTCAGGGCTCAACCCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg < 2:395589/150‑1 (MQ=255) aGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCcaca > 2:15786/1‑83 (MQ=255) aGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCcaca < 1:15786/83‑1 (MQ=255) agat‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACg < 1:266778/147‑1 (MQ=255) tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt > 2:15665/1‑66 (MQ=37) tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt < 1:15665/66‑1 (MQ=37) tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa > 2:496872/1‑150 (MQ=255) | GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTA > NC_002947/4338395‑4338676 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |