Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 15,386 | C→T | 20.5% | intergenic (+88/+59) | dnaJ → / ← mokC | chaperone Hsp40, co‑chaperone with DnaK/regulatory protein for HokC, overlaps CDS of hokC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 15,386 | 0 | C | T | 20.5% | 104.6 / 22.4 | 54 | intergenic (+88/+59) | dnaJ/mokC | chaperone Hsp40, co‑chaperone with DnaK/regulatory protein for HokC, overlaps CDS of hokC |
Reads supporting (aligned to +/- strand): ref base C (22/19); new base T (8/3); total (31/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01 |
GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTC > minE/15318‑15454 | gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACt > 1:2462978/1‑70 (MQ=255) gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTcc > 1:3841102/1‑70 (MQ=255) agtgtgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTcc > 1:825505‑M1/66‑70 (MQ=255) tgtgctatctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCAt > 1:564670‑M1/64‑71 (MQ=255) tgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATaa > 1:3359500‑M1/62‑71 (MQ=255) tgctgtctgcgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATaa > 1:3503225‑M1/62‑71 (MQ=255) aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATaa < 1:1472933/71‑1 (MQ=255) gCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgc > 1:3355855/1‑71 (MQ=255) cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc < 1:1190796/59‑1 (MQ=255) cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAgcgc > 1:3076065/1‑71 (MQ=255) ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCt < 1:2698701/69‑1 (MQ=255) tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTc < 1:1750825/56‑1 (MQ=255) tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTa < 1:3209705/68‑1 (MQ=255) tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt > 1:3299374/1‑71 (MQ=255) cgtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTaa > 1:2377006‑M1/53‑69 (MQ=255) gtcaatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTAACtt > 1:388121‑M1/52‑71 (MQ=255) aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt > 1:3013502/1‑70 (MQ=255) caatcgcagcactggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAgcgc < 1:2074114‑M1/14‑1 (MQ=255) aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg > 1:3109025/1‑71 (MQ=255) gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt < 1:2301506/71‑1 (MQ=255) cactggtgctgtaagaccatacataaaaaagccgggtataATTCCCATAAGCGCTAACTTAAGGGTTGTgg > 1:3335450‑M1/41‑71 (MQ=255) ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt < 1:584890/71‑1 (MQ=255) tggtgctgtaagaccatacataaaaaagccggggataATTCCCATAAGCGCTa > 1:3298036‑M1/38‑53 (MQ=255) gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCAc > 1:71128/1‑37 (MQ=255) gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa > 1:1099710/1‑71 (MQ=255) gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc > 1:1724997/1‑71 (MQ=255) atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATaa > 1:1312009/1‑41 (MQ=255) atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTa > 1:3551033/1‑50 (MQ=255) tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg < 1:2917456/70‑1 (MQ=255) tGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg < 1:3897397/70‑1 (MQ=255) cGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa > 1:1283318/1‑48 (MQ=255) cGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata > 1:3952050/1‑71 (MQ=255) gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc > 1:1519043/1‑62 (MQ=255) gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat < 1:243991/71‑1 (MQ=255) ccTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt < 1:1035708/71‑1 (MQ=255) cTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg < 1:736040/50‑1 (MQ=255) ataaaaaagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa < 1:234280‑M1/53‑1 (MQ=255) gTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc > 1:3339663/1‑71 (MQ=255) tAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt < 1:340530/43‑1 (MQ=255) tAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa < 1:3302265/49‑1 (MQ=255) aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgt > 1:3893103/1‑71 (MQ=255) aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg > 1:1679547/1‑38 (MQ=255) aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtg > 1:809605/1‑71 (MQ=255) aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtg > 1:2771355/1‑71 (MQ=255) gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt < 1:1213525/46‑1 (MQ=255) aagccggggataATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgcc < 1:1899622‑M1/59‑1 (MQ=255) tGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg < 1:305575/71‑1 (MQ=255) tgCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg < 1:1060814/69‑1 (MQ=25) ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccgt > 1:1943939/1‑71 (MQ=255) gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccgtg < 1:2857892/71‑1 (MQ=255) gCGGGGAGTAATCCCATAAGCGCTAAC‑TAAGGGTTGTGGTATTACGACTGATATGATTTAAcgtgccgat > 1:3106755/1‑68 (MQ=25) ggggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt > 1:3929615/1‑42 (MQ=255) gggATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat < 1:1580968/52‑1 (MQ=255) gATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGctctc < 1:3335497/71‑1 (MQ=255) tCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtctt < 1:2917546/71‑1 (MQ=255) cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCtcttc > 1:1601470/1‑71 (MQ=255) cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa > 1:2305853/1‑35 (MQ=255) cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAcgtgccg < 1:2956801/59‑1 (MQ=255) | GTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGTGCTCTCTTC > minE/15318‑15454 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |