Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,962,626:1 | +T | 9.3% | coding (780/1152 nt) | PP_3492 → | acyl‑CoA dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,962,626 | 1 | . | T | 9.3% | 153.3 / 11.6 | 43 | coding (780/1152 nt) | PP_3492 | acyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base . (19/20); new base T (2/2); total (21/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01 |
TGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCAGTGCC > NC_002947/3962477‑3962773 | ggcgcGGGCTCCTGCGACAACTGCCGCATTCCTGCCGCCAATATCCCCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGaa < 2:768252/149‑1 (MQ=255) gcgcGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGCCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑g < 1:1044557/150‑1 (MQ=255) gcgGTCACCCTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑gcg < 2:598564/150‑1 (MQ=255) cACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcg < 2:985846/150‑2 (MQ=255) aCAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGtgatgc > 2:855503/1‑145 (MQ=255) cGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGt > 2:64862/1‑137 (MQ=255) cGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGt < 1:64862/137‑1 (MQ=255) ccgccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCTCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGt > 2:235918/1‑150 (MQ=255) cgccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGtt < 2:215020/150‑1 (MQ=255) gccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGttt < 1:279193/150‑1 (MQ=255) gccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGttt < 1:537370/150‑1 (MQ=255) cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTgg > 2:891585/1‑150 (MQ=255) cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcgg < 2:946114/121‑3 (MQ=255) cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcgg > 1:946114/1‑119 (MQ=255) tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc < 2:279050/150‑1 (MQ=255) gAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGGCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAAc < 1:446304/150‑1 (MQ=255) ggCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAacac < 2:542919/150‑1 (MQ=255) aGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCaga < 1:1036795/150‑1 (MQ=255) ggCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc > 1:503196/1‑125 (MQ=255) ggCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc < 2:503196/125‑1 (MQ=255) cATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCaa < 1:939547/150‑1 (MQ=255) gTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGc < 1:746948/150‑1 (MQ=255) cAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGa > 2:736571/1‑150 (MQ=255) aGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACAtgcaggt > 2:867589/1‑150 (MQ=255) cGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGcaggtgcaggtg < 2:133414/150‑1 (MQ=255) tATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGaa > 1:1127293/1‑150 (MQ=255) tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc > 1:869097/1‑90 (MQ=255) tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc < 2:869097/90‑1 (MQ=255) tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATgcaggtgcaggt > 2:1085632/1‑145 (MQ=255) tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATgcaggtgcaggt < 1:1085632/145‑1 (MQ=255) aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAAcgccgc < 2:626708/149‑1 (MQ=255) aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAAcgccgc > 1:626708/1‑149 (MQ=255) aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCa > 2:1161316/1‑150 (MQ=255) gccgccCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGt < 1:1161316/150‑1 (MQ=255) cccAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGc < 1:736571/150‑1 (MQ=255) ccTTGGGCATCGCCCGCGCCGCATTCGCA‑GCGGGCCTGGTATAATCGCGTTACTGCAG‑ACAGTTTGGCAAGCCGATCAACCACGACCATAGCATGGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCACACGTGGGCTGATcc > 1:931335/1‑150 (MQ=255) cTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc > 2:875462/1‑71 (MQ=255) cTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc < 1:875462/71‑1 (MQ=255) ttGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTg > 1:641361/1‑150 (MQ=255) cagcgcCGC‑TTCGAATGCGGCGCTGGTCTACTCGCGTGACCGCATCACCG‑TTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGACATACAGGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGccgcccg > 1:874282/3‑150 (MQ=255) cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCt < 2:402813/124‑1 (MQ=255) cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCt > 1:402813/1‑124 (MQ=255) cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCa > 2:441436/1‑150 (MQ=255) gAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACAtgcaggtgcag > 2:868207/1‑99 (MQ=255) gAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACAtgcaggtgcag < 1:868207/99‑1 (MQ=255) gAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCAGTGcc > 1:851888/1‑150 (MQ=255) | TGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCAGTGCC > NC_002947/3962477‑3962773 |
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: ‑ |