Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,666,747 | A→T | 13.9% | intergenic (+21/‑45) | acnA‑II → / → prpF | aconitate hydratase/aconitate isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,666,747 | 0 | A | T | 13.9% | 76.5 / 6.2 | 36 | intergenic (+21/‑45) | acnA‑II/prpF | aconitate hydratase/aconitate isomerase |
Reads supporting (aligned to +/- strand): ref base A (12/19); new base T (3/2); total (15/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01 |
CAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGC > NC_002947/2666603‑2666889 | cAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCgg < 3:160808/149‑1 (MQ=255) aaTGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGC‑GCGCCAAGCGGc > 2:134892/1‑149 (MQ=255) gAGTGCCTGGAAGTACCGGTGACCTGCCGGCGGGACCCCGCCGAGGACGTGGCGATCAACGAAGCGGGCGGGGGGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAgg < 8:290223/149‑1 (MQ=255) aGTGCCTGGAAGTACTGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCCCGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCTGGGG‑CCGCTTTGCGGcccccga < 1:98491/149‑4 (MQ=255) tGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTa < 5:267442/149‑1 (MQ=255) cTGGAAGTACCGGGGACCTGCCGTCTGGACACCGCCGAGTAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACGGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTAttt < 1:134892/149‑1 (MQ=255) gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCa < 5:249015/138‑1 (MQ=255) ggTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTGAGCGGGCTAt > 3:26173/1‑149 (MQ=255) gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGa > 1:254923/1‑144 (MQ=255) gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGa < 2:254923/144‑1 (MQ=255) acacCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGCCCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCCGGGGATGTGCCCGGCCGAGCAGGATGCCCAAGaaaa > 6:256284/1‑149 (MQ=255) cgccgAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGATCATGAATCGGGGGTCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGt < 1:278884/149‑1 (MQ=255) gtgtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCCCAAAGCGGCCCCAGGGTATTTCACCGGCCGATAAGGATGTCCAACAAAATGTCACAATTAccc > 7:112110/1‑149 (MQ=255) gtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGa < 7:160246/130‑1 (MQ=255) gtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTAccc > 6:125673/1‑147 (MQ=255) cGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaaa > 2:53598/1‑148 (MQ=255) gcgggcggGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATcccc < 5:125673/149‑1 (MQ=255) ggggTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGGATTTCACCGGACTCTCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCAcc < 6:368914/149‑1 (MQ=255) cTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTGCCCGGTCTATAAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGcggcgg > 7:91588/1‑149 (MQ=255) cAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCa > 3:65684/1‑114 (MQ=255) cAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCa < 4:65684/114‑1 (MQ=255) ccTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaa > 5:315363/1‑108 (MQ=255) ccTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaa < 6:315363/108‑1 (MQ=255) tCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTAACCGGACGATCAGGCTGCCCAAGAAAATGTCCCATGTACCCCAGCACAAAATCCCCGCCCCCTACATCCGCGGCGGCCCCAGCCAAGGcc > 1:308013/1‑145 (MQ=255) gAAGCGACCGCCTGAGCATGACTCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCgtgt < 8:248003/147‑1 (MQ=255) gAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGtt < 1:333084/148‑1 (MQ=255) cGCCTGAGCATGAATCGGCGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt < 4:26173/149‑1 (MQ=255) ccTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCgcgg < 1:149590/117‑1 (MQ=255) ccTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTAACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTg > 3:313415/1‑148 (MQ=255) gAATCGGGGG‑ACGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCtcc < 5:203277/149‑3 (MQ=255) aaTCGGGGG‑CCGCAAAGCGGACCCAGTGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGTCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGa < 5:298528/143‑1 (MQ=255) cctgggg‑CCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt > 7:131355/4‑135 (MQ=255) cctgggg‑CCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt < 8:131355/132‑1 (MQ=255) cctgggg‑CCGCTTTGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGTACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGccc > 7:337685/4‑148 (MQ=255) cctgggg‑CCGCTTTGCGG‑CCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTAGATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGccc > 1:293151/4‑147 (MQ=255) ccGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGTCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCTGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgc < 8:112110/149‑1 (MQ=255) | CAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTG‑CCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGG‑CCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGC > NC_002947/2666603‑2666889 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |