Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,839,707:1 | +G | 100% | coding (73/1443 nt) | ccoN‑II → | cbb3‑type cytochrome c oxidase subunit 1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,839,707 | 1 | . | G | 82.6% | 126.5 / 23.0 | 46 | coding (73/1443 nt) | ccoN‑II | cbb3‑type cytochrome c oxidase subunit 1 |
Reads supporting (aligned to +/- strand): ref base . (4/4); new base G (15/23); total (19/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01 | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCTTCGGGGGCTGT > NC_002947/4839561‑4839850 | gttAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAAAATCAAAACCCTTGAACCTTAGAATTAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCTACAGTACGCATTCAAGAGGGTGGGCT‑gg > 1:163679/3‑149 (MQ=255) gTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑ggg < 8:159108/149‑1 (MQ=255) gTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCCACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑ggg > 8:86169/1‑149 (MQ=255) ccTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑gggggga > 3:116765/1‑147 (MQ=255) ccTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAACACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑gggggga > 1:225056/1‑147 (MQ=255) ttGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑ggggggat > 1:248572/1‑145 (MQ=255) gTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTgg > 7:320621/1‑148 (MQ=255) aTCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCAt < 5:182192/148‑1 (MQ=255) gcgcCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGc < 7:86169/148‑1 (MQ=255) ccGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCt < 5:92564/149‑1 (MQ=255) tttttCGACCAACAAGAACATCAAAACCGTCGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCAt < 4:352915/149‑1 (MQ=255) tttCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTc > 6:276906/1‑145 (MQ=255) tttCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTc < 5:276906/145‑1 (MQ=255) gACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTc > 4:251285/1‑149 (MQ=255) tCAAATCCGTTGAAACTTAGAATGATCACAGCTATCAGTCCGTCTGGTTATAACTATAAGGTCGTCGGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGGCTCGTCGTCTTCTTCGCCTCGCAGCTGGTGTGGc < 4:116765/149‑1 (MQ=255) tCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGc > 1:156115/1‑149 (MQ=255) tCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGTTCCTTGGCATGGGCCTCGGCGTCt < 7:23548/124‑1 (MQ=255) tCAAAACCGTGGAACCTTAGAATGAACACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGc < 2:225056/149‑1 (MQ=255) aaaaCCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATAGCCTCGCAGCTGGTGTGGCCg > 3:238807/1‑149 (MQ=255) aaaCCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACCATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGc < 5:117361/149‑1 (MQ=255) gTGGAACCTTAGTATGATCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCTGCATTTCGACATCATGACGGTGGTCTGGGGGGATCGTTTGCATGGGCCTCGGCGTGTTCATCGCCTCGCAGCTGGGGTGGCCGCAGCTg < 2:163679/149‑1 (MQ=255) ccTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTg > 8:37819/1‑149 (MQ=255) ccTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCCTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTg < 2:279324/149‑1 (MQ=255) ttAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGa > 5:317761/1‑149 (MQ=255) aGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGa < 7:88498/147‑1 (MQ=255) aaTGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCtt < 7:37819/149‑1 (MQ=255) aGCACAGCTAGCAGGCCGGCTGCTTCTAACGATCCGGTCGGCCGCCCGTGCGCGATCATGACGGGGGGCTGGGGGGTTCGTTGGCATGGGCCTCGGCGTCTTCTGCGCCTCGCTGCTGGGGGGGCCGCAGCTGAACCTGGACCTTCCCt < 3:192483/149‑1 (MQ=255) aaTCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGc > 7:261060/1‑149 (MQ=255) aGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCg < 8:261060/149‑1 (MQ=255) gTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCgg > 7:92191/1‑149 (MQ=255) aCTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTc < 7:321988/81‑1 (MQ=255) aCTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTc > 8:321988/1‑81 (MQ=255) tataAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGc < 8:92191/149‑1 (MQ=255) ataAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGc > 3:110163/1‑145 (MQ=255) ataAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGc < 4:110163/145‑1 (MQ=255) taAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATgg > 4:111376/1‑57 (MQ=21) taAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATgg < 3:111376/57‑1 (MQ=21) tcgtcCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCaa > 8:6810/1‑149 (MQ=255) gtcCGCCAGTTCGCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAAc < 6:317761/148‑1 (MQ=255) gtcCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCT‑GACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTg > 2:316277/1‑149 (MQ=255) gCCAGTTCGCCAACATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGt > 1:121206/1‑148 (MQ=255) ccATTTCGCTATAATGACTGTGGTGTGGGGGGATCATTGGAATGGGCCTAGGCGTCTTCATCGCCTAGCTGATGGTGTGGACGGAGCTGAAACTGGATCTTACCTGGACCAGCTTCGGCCGACTGGGCGCTCTGCACTCCTACCTGGTg < 4:238807/149‑1 (MQ=255) gCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTc < 2:248572/148‑1 (MQ=255) gCCATCATGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTc < 2:156115/148‑1 (MQ=255) tcatGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGc < 7:79483/124‑1 (MQ=255) tcatGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGc > 8:79483/1‑124 (MQ=255) tcatGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTg < 2:267759/135‑1 (MQ=255) tcatGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTg > 1:267759/1‑135 (MQ=255) tcatGACGGTGGTCTGGGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCtt > 8:55961/1‑149 (MQ=255) atGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCTTCgg > 3:156283/1‑149 (MQ=255) gtggtCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCT‑GACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCTTCGGGGGCtgt < 1:316277/149‑1 (MQ=255) tggtCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCTTCGGGGGCtgt < 2:121206/149‑1 (MQ=255) | GGTAACCTTGCTGGTTATCAATACTGGCGCCGATTTTTCGACCAACAAGAACATCAAAACCGTGGAACCTTAGAATGAGCACAGCAATCAGTCCGACTGCTTATAACTATAAGGTCGTCCGCCAGTTCGCCATCATGACGGTGGTCT‑GGGGGATCCTTGGCATGGGCCTCGGCGTCTTCATCGCCTCGCAGCTGGTGTGGCCGCAGCTGAACCTGGACCTTCCCTGGACCAGCTTCGGCCGCCTGCGCCCGCTGCACACCAACCTGGTGATCTTCGCCTTCGGGGGCTGT > NC_002947/4839561‑4839850 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |