Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,045,826 | A→G | 8.1% | intergenic (+588/+207) | yqiK → / ← rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,045,826 | 0 | A | G | 8.1% | 163.7 / 6.5 | 62 | intergenic (+588/+207) | yqiK/rfaE | conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase |
Reads supporting (aligned to +/- strand): ref base A (28/29); new base G (3/2); total (31/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01 |
GCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTT > minE/2045758‑2045894 | gcgcGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaag < 1:913967/70‑1 (MQ=255) cGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTg < 1:1179650/71‑1 (MQ=255) cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat < 1:834630/71‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga < 1:3946819/52‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa < 1:1718862/71‑1 (MQ=255) tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaa < 1:1971785/71‑1 (MQ=255) caacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga < 1:3901906/51‑1 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg > 1:317614/1‑71 (MQ=255) aacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAg > 1:2129772/1‑71 (MQ=255) aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGtt < 1:4136953/71‑1 (MQ=255) aTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGt < 1:3753043/70‑1 (MQ=255) tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa > 1:1746001/1‑71 (MQ=255) tCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTa > 1:4079186/1‑71 (MQ=255) cGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat > 1:3696942/1‑71 (MQ=255) ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat < 1:722908/70‑1 (MQ=255) ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat > 1:1486467/1‑70 (MQ=255) ggggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat > 1:3568808/1‑70 (MQ=255) gggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatgatg < 1:3509603/71‑1 (MQ=255) ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTAtgatga < 1:3685976/71‑1 (MQ=255) agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat < 1:2458235/49‑1 (MQ=255) agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATaa > 1:1678597/1‑51 (MQ=255) agcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATaa > 1:50707/1‑51 (MQ=255) gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat < 1:2291861/48‑1 (MQ=255) caAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc < 1:975157/71‑1 (MQ=255) aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg > 1:1111609/1‑71 (MQ=255) aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAgg < 1:2603469/54‑1 (MQ=255) aGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt > 1:2384862/1‑71 (MQ=255) gggCTGAATCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTAGTATAAGGAAAGGGt < 1:3714754/54‑1 (MQ=18) gggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtca > 1:3981406/1‑71 (MQ=255) acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGtc < 1:1423385/69‑1 (MQ=0) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat > 1:1375826/1‑71 (MQ=255) ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat > 1:3860975/1‑71 (MQ=255) gCTGAAACGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAgg < 1:574752/45‑1 (MQ=11) gCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatga > 1:4056918/1‑57 (MQ=255) cTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc > 1:2487077/1‑70 (MQ=255) cTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTCATGATGAAGCTCGtcatca < 1:3411888/71‑1 (MQ=255) gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGta > 1:307769/1‑37 (MQ=255) gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt < 1:1798331/64‑1 (MQ=255) gAAACGGGAAAGCCCCTACCGAAGAAG‑GGGCCTTGTATAAgta > 1:77338/1‑41 (MQ=18) aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGa < 1:4084898/42‑1 (MQ=255) aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatg < 1:1446976/52‑1 (MQ=255) aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatg < 1:4055808/52‑1 (MQ=255) aaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACt < 1:83768/71‑1 (MQ=255) aCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGaaa > 1:2312814/1‑42 (MQ=14) ggAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTCATCATCGTGgtg > 1:911323/1‑70 (MQ=2) ggAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgtt > 1:552725/1‑71 (MQ=255) aaaGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGgtga < 1:3016519/71‑3 (MQ=2) aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg > 1:183808/1‑71 (MQ=255) aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg > 1:3871001/1‑71 (MQ=255) aaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgt > 1:3866165/1‑70 (MQ=255) aGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGt > 1:1868971/1‑39 (MQ=14) gCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc > 1:3552264/1‑51 (MQ=255) gCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg > 1:2698134/1‑67 (MQ=255) cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt > 1:3495394/1‑71 (MQ=255) cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt < 1:124610/71‑1 (MQ=255) ccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt < 1:355828/71‑1 (MQ=255) cTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGtc < 1:2982174/35‑2 (MQ=19) cTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatc < 1:957625/53‑1 (MQ=255) tCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat < 1:2025853/36‑1 (MQ=255) cccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg > 1:2567221/1‑63 (MQ=255) cccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt > 1:3510575/1‑71 (MQ=255) ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGaa > 1:2553070/1‑40 (MQ=255) ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat < 1:711453/49‑1 (MQ=255) ccGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt < 1:238083/70‑1 (MQ=255) cGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcatca < 1:2044391/50‑1 (MQ=255) gaagaaG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt < 1:2295558/71‑1 (MQ=255) aagaag‑GGGCCTTGTATAAGGGAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt > 1:3893925/1‑69 (MQ=255) aagaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt > 1:3216719/1‑70 (MQ=255) agaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATa > 1:2576634/1‑49 (MQ=255) agaag‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt < 1:1974021/69‑1 (MQ=255) | GCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTT > minE/2045758‑2045894 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |