Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,252,857 | T→A | 13.6% | intergenic (+70/‑49) | leuD → / → PP_1987 | 3‑isopropylmalate dehydratase small subunit/UbiE/COQ5 family methlytransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,252,857 | 0 | T | A | 13.6% | 116.7 / 10.2 | 44 | intergenic (+70/‑49) | leuD/PP_1987 | 3‑isopropylmalate dehydratase small subunit/UbiE/COQ5 family methlytransferase |
Reads supporting (aligned to +/- strand): ref base T (21/17); new base A (2/4); total (23/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01 |
GGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCA > NC_002947/2252715‑2252996 | ggACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTAAACCcgcg < 2:117665/150‑1 (MQ=255) aaCAGCGCCATGAGCAAGGCCTGCTAAGCCATGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTAAACCcgcgc < 2:341646/149‑2 (MQ=255) aCAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTAAACCcgcgcc < 1:326861/150‑3 (MQ=255) aCAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGCGTAAACCcgcgcc < 1:239157/150‑3 (MQ=255) gACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAACTGTAC‑CGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGc > 1:338968/1‑150 (MQ=255) aCGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGcc > 1:183142/1‑150 (MQ=255) aTCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGcac > 2:464345/1‑150 (MQ=255) tCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCGGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGcaca > 1:109953/1‑150 (MQ=255) ttCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGcaacaa > 2:390900/1‑150 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTAc < 1:242646/115‑1 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTAc > 2:242646/1‑115 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGc > 1:161029/1‑127 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGc < 2:161029/127‑1 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaa > 1:84709/1‑150 (MQ=255) cAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTAGCTGTATTCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaa > 1:591601/1‑150 (MQ=255) aGCACCGCGCCAGCCAGCCTTGGCTGTTCAGCGATGGCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAg < 2:338968/150‑1 (MQ=255) aGCACCGCGCCAGCCAGCCTGGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAg < 2:118962/150‑1 (MQ=255) gcgcCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAACGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACga < 2:520310/150‑1 (MQ=255) aaccatccatGGCAGGTCCGTGAAGCCTGAAATTCGTAGGTGGCAATGCGGGCTGCCTCGCTGCAACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACAt < 2:109953/148‑1 (MQ=37) tGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGAccagcacccag < 2:218003/150‑1 (MQ=255) gTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATc > 1:10186/1‑150 (MQ=255) aTGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAt < 2:523374/150‑1 (MQ=255) ccTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAtgtg < 2:10186/150‑1 (MQ=255) ccTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAtgtg < 1:390900/150‑1 (MQ=255) cTGAAGTTCGTAGGTGACAGTGCGGCCCTGCTCGCTCCTAGAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTgg < 1:608969/150‑1 (MQ=255) gTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGcc > 1:518004/1‑150 (MQ=255) tGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCTGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTc > 1:193620/1‑150 (MQ=255) tGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCg < 1:248416/70‑1 (MQ=255) tGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCg > 2:248416/1‑70 (MQ=255) ctcctcGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCa > 1:214521/1‑150 (MQ=255) ctcctcGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCa > 2:187326/1‑150 (MQ=255) tcGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaa < 2:330902/80‑1 (MQ=255) tcGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaa > 1:330902/1‑80 (MQ=255) cTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCt > 2:541300/1‑150 (MQ=255) tACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAAGTAGTCCAACGCCAAATCTGCGCGGAGACCAACACATACCTc > 2:101392/1‑150 (MQ=37) aCAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGtt < 1:399748/124‑1 (MQ=255) aCAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGtt > 2:399748/1‑124 (MQ=255) cAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACgaga > 1:490683/1‑81 (MQ=255) cAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACgaga < 2:490683/81‑1 (MQ=255) gAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTcagca < 2:84709/150‑1 (MQ=255) tGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGc < 2:200459/142‑1 (MQ=255) tGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGc > 1:200459/1‑142 (MQ=255) tcttcgcgGGTAAACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAgg > 2:142922/5‑150 (MQ=255) ttcgcgGGTAAACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCa > 1:287780/3‑150 (MQ=255) | GGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGT‑TACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCA > NC_002947/2252715‑2252996 |
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: ‑ |