Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,348,163 | C→A | 24.3% | E173* (GAG→TAG) ‡ | PP_3820 ← | group II intron‑encoding maturase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,348,163 | 0 | C | A | 24.3% | 63.8 / 18.1 | 37 | E173* (GAG→TAG) ‡ | PP_3820 | group II intron‑encoding maturase |
Reads supporting (aligned to +/- strand): ref base C (14/14); new base A (4/5); total (18/19) | |||||||||||
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.33e-01 |
GCGGCGAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGC > NC_002947/4348014‑4348256 | gcggcgAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCAtatc < 2:1068248/150‑1 (MQ=17) gcgAATAAGCCTGAGTACGCACTTGTCTTTGATCCGCCGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAt < 1:1007563/150‑1 (MQ=17) gcgAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAt > 2:857818/1‑150 (MQ=32) aaTAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGc > 1:588120/1‑150 (MQ=32) aTAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCt > 2:328013/1‑150 (MQ=32) aGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGa < 2:1349783/150‑1 (MQ=32) aGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtcgt < 1:1011718/150‑4 (MQ=255) gCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtcgtg > 1:524090/1‑146 (MQ=255) cTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtcgtgct < 2:524090/150‑7 (MQ=255) cTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtc < 2:151640/145‑2 (MQ=12) cTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtc > 1:151640/1‑144 (MQ=12) gAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGt < 2:771157/150‑1 (MQ=32) gATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGc < 1:1281406/130‑1 (MQ=33) gATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGc > 2:1281406/1‑130 (MQ=33) ccGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtc < 1:42620/122‑2 (MQ=12) ccGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtc > 2:42549/1‑121 (MQ=255) ccGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtc < 1:42549/122‑2 (MQ=255) ccGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCACGCTGTCACATGAGCACGAGCCATCGAGATGGCtc > 2:42620/1‑121 (MQ=255) ccATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGc > 2:913286/1‑120 (MQ=33) ccATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGc < 1:913286/120‑1 (MQ=33) gATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTgg < 1:922200/150‑1 (MQ=32) cccGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtg > 2:267663/1‑150 (MQ=32) cGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGtgtg > 1:950447/1‑150 (MQ=32) gATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGa < 2:1397090/150‑1 (MQ=32) aTCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtgt < 2:973370/147‑1 (MQ=32) aTCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtgt > 1:973370/1‑147 (MQ=32) cGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAgct > 2:1020087/1‑150 (MQ=32) aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGccc < 1:627563/93‑1 (MQ=34) aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGccc > 2:627563/1‑93 (MQ=34) aaaaCTTCTCCGGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAgctgct < 1:306498/150‑1 (MQ=17) ttCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCtt > 1:1104686/1‑70 (MQ=37) ttCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCtt < 2:1104686/70‑1 (MQ=37) ctcCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCtgctgctgccg < 2:950447/150‑3 (MQ=17) cAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtgc > 2:722433/1‑150 (MQ=32) gATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGcc < 1:1017125/79‑1 (MQ=35) gATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGcc > 2:1017125/1‑79 (MQ=35) gAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGt > 1:934933/1‑109 (MQ=32) gAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGt < 2:934933/109‑1 (MQ=32) | GCGGCGAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGC > NC_002947/4348014‑4348256 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |