Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,691,815 | C→G | 13.8% | intergenic (+40/‑131) | PP_1487 → / → PP_1488 | hypothetical protein/methyl‑accepting chemotaxis transducer |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,691,815 | 0 | C | G | 13.8% | 129.9 / 8.0 | 48 | intergenic (+40/‑131) | PP_1487/PP_1488 | hypothetical protein/methyl‑accepting chemotaxis transducer |
Reads supporting (aligned to +/- strand): ref base C (21/20); new base G (2/5); total (23/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTG > NC_002947/1691667‑1691963 | ggAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCg > 2:940963/1‑150 (MQ=255) gAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑a > 1:553136/1‑150 (MQ=255) ccTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGTCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTTTGAGATCG‑AGCGCCGCACGGGCGGCGCt < 1:431719/150‑1 (MQ=255) gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTcgatctca < 1:533186/150‑8 (MQ=38) cGAGCTGGCCCGCCAGCTGCGGCTGGACCCTGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTcaca < 2:553136/150‑1 (MQ=255) cGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTCGATCTcaca < 2:614149/150‑1 (MQ=38) gAGCTGGCCCGCCCGCTGAGGCTGGACCCGGGCTTGGGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACCGACTCATCAGCGACTGTCAGATCG‑ATCACCGCCCGTGCGCCTCTCCATCTCACAg < 2:945335/150‑1 (MQ=255) gAGCTGGCCCGCCAGCTGCGGCTGGAACCGGGCTTGCGCCAGGAGCTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAg > 2:400053/1‑150 (MQ=255) aGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAgg < 1:190023/150‑1 (MQ=255) gCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCTGGGCTTGCCACGGCCGATTGAGCGCCTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTCGATCTCACAGgc < 2:161148/150‑1 (MQ=37) gCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaa > 1:424658/1‑150 (MQ=255) ccGCCAGCTGCGGCTGGACCCGGGCTTGCACCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaaaa > 1:254645/1‑150 (MQ=255) tGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCCGGTGCGGTTGGGCGTGGCGTTTTGGCGGTGCTTGCCACGGCCTATTGAGCGACTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAAt < 2:750375/150‑1 (MQ=38) ccGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACAt > 1:845482/1‑150 (MQ=255) cTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGAC‑ACTACGCTCTTCCGATCTGGTAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATa < 2:1118905/150‑1 (MQ=255) ccAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTaa > 1:1197859/1‑150 (MQ=255) ggAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGAc < 2:424658/150‑1 (MQ=255) ggAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGAc > 1:697515/1‑150 (MQ=255) aaCTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACct > 1:1144032/1‑150 (MQ=255) cTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACctct < 2:845482/150‑1 (MQ=255) aaGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCt < 2:1144032/150‑1 (MQ=255) ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCCAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCtata > 2:954337/1‑150 (MQ=255) ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCtata > 1:124396/1‑150 (MQ=255) gTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTc < 1:1064275/150‑1 (MQ=255) gTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCGGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAAGACCGTAAGGCCTCTGGGCTATACTTCGGCGATTTTTCCCGCTc < 1:400053/150‑1 (MQ=255) gTGGCGGTGCGTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCCATCCCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGtt < 2:254645/150‑1 (MQ=255) gCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgtg > 2:885936/1‑150 (MQ=255) gCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCt > 2:1072558/1‑78 (MQ=255) gCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCt < 1:1072558/78‑1 (MQ=255) gCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGCCTCACAGGCGCAAAATTCCTGACGACAATcaca < 1:754242/81‑1 (MQ=255) gCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGCCTCACAGGCGCAAAATTCCTGACGACAATcaca > 2:754242/1‑81 (MQ=255) cTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTga < 2:194042/68‑1 (MQ=255) cTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTga > 1:194042/1‑68 (MQ=255) tATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt > 2:495856/1‑150 (MQ=255) aTTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtg > 1:341315/1‑150 (MQ=255) tGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgtg < 2:1197859/150‑1 (MQ=255) gcCTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCTAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTgaagaa > 2:224142/1‑150 (MQ=37) ccTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCTGCTCGTTGTGAATTCCTGa > 2:1357614/1‑132 (MQ=37) ccTGTGAGACGG‑AGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGa < 1:1357614/132‑1 (MQ=37) cTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgt < 1:1150226/121‑1 (MQ=255) cTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgt > 2:1150226/1‑121 (MQ=255) tgAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTCAGCGGCCAGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa > 1:247146/1‑150 (MQ=255) tgAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa > 1:1048042/1‑150 (MQ=255) gagaTCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGAc < 1:954337/150‑1 (MQ=255) agaTCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTTAATTCCTGAGGGCTGAATGTGAAGAACTGGAcc < 2:854753/150‑1 (MQ=255) aTCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCt > 1:1184797/1‑129 (MQ=255) aTCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCt < 2:1184797/129‑1 (MQ=255) tCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTg < 1:885936/150‑1 (MQ=255) | GGAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCG‑AGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTG > NC_002947/1691667‑1691963 |
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: ‑ |