Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,155,592 | C→T | 18.7% | intergenic (+103/+83) | dacB → / ← obgE | D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,155,592 | 0 | C | T | 18.7% | 177.0 / 29.7 | 81 | intergenic (+103/+83) | dacB/obgE | D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair |
Reads supporting (aligned to +/- strand): ref base C (45/20); new base T (8/7); total (54/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.64e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01 |
CAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC > minE/2155523‑2155659 | cAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg > 1:4382492/1‑71 (MQ=255) attgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg > 1:474442/1‑70 (MQ=255) tgttGACAACGTACGCATTGTTCATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAACAGTGCa > 1:534654/3‑70 (MQ=11) tgttGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc > 1:339141/3‑69 (MQ=255) tgctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCa > 1:4908578/4‑70 (MQ=11) gctGACAAATTGCGCTTTATGT‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt > 1:2276917/3‑71 (MQ=11) ctGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAt > 1:4521451/2‑71 (MQ=255) aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt > 1:297689/3‑71 (MQ=11) gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt > 1:1553547/1‑70 (MQ=255) aCAAATTGCGCTTTATGT‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa > 1:2350845/1‑71 (MQ=11) aCAAATTGCGCTTTATGT‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa > 1:4820455/1‑71 (MQ=11) aCAAATTGCGCTTTATGT‑TGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa > 1:4762957/1‑71 (MQ=11) aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc > 1:1072952/1‑70 (MQ=39) aacaGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaa > 1:4747307/4‑71 (MQ=11) aaaGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaat < 1:422614/71‑2 (MQ=255) aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaatac < 1:4900152/71‑4 (MQ=9) tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:2710053/1‑66 (MQ=255) tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata < 1:4553164/68‑3 (MQ=11) tGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCaata < 1:4563233/68‑3 (MQ=11) cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa < 1:970542/55‑1 (MQ=255) ttGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt > 1:1546328/1‑71 (MQ=255) tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat > 1:5692738/1‑64 (MQ=40) tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat > 1:276101/1‑64 (MQ=40) aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa < 1:2334484/54‑1 (MQ=255) tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTg < 1:4561491/71‑1 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc > 1:444215/1‑56 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCa > 1:3612334/1‑57 (MQ=255) tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc‑t < 1:891784/60‑1 (MQ=255) tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCt < 1:5120457/71‑1 (MQ=255) cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt > 1:253713/1‑70 (MQ=11) cggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg < 1:2459078/71‑1 (MQ=14) cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTg < 1:695466/59‑1 (MQ=255) cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGa > 1:2469852/1‑68 (MQ=255) cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGAt > 1:2812655/1‑69 (MQ=255) cgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa < 1:1772080/71‑1 (MQ=14) tAAACGCCTTATCCGGCCTACAAAATCGTGCAAAt < 1:614982/35‑1 (MQ=40) tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc < 1:3992907/37‑1 (MQ=40) tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc < 1:5697768/47‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTg < 1:1200361/57‑1 (MQ=255) aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc < 1:79329/36‑1 (MQ=40) aaaCGCCTTATCCGGCCTACAAAATCGTGCAAATTc < 1:3829410/36‑1 (MQ=40) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:1040965/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:909019/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:726113/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:1087824/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:1094565/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:2181420/1‑37 (MQ=37) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTa > 1:3049852/1‑57 (MQ=255) aaCGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCg > 1:4316326/1‑71 (MQ=11) aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAACATATTGCAATTC‑TCTTg < 1:3097118/54‑1 (MQ=255) aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa > 1:3463578/1‑36 (MQ=33) aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTg < 1:5603832/54‑1 (MQ=255) aCGCCTTATCCGGCCTACAAAATCCTGCAAATTCAACATATTGCAATTC‑TCTTg < 1:3922463/54‑1 (MQ=255) cGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCgca < 1:775285/71‑3 (MQ=11) cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAgcg > 1:1525361/1‑70 (MQ=255) aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc‑t > 1:4672749/1‑43 (MQ=255) aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:2995949/1‑70 (MQ=255) aTCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt > 1:2438123/1‑71 (MQ=11) tCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc‑t > 1:1891833/1‑42 (MQ=255) tCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc < 1:218615/71‑1 (MQ=11) ccGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGt > 1:4635591/1‑46 (MQ=255) gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc < 1:2681351/71‑1 (MQ=11) gCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc < 1:4763606/71‑1 (MQ=11) cTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAg > 1:5233610/1‑58 (MQ=255) cTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAg > 1:5235795/1‑58 (MQ=255) tACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat > 1:5363373/1‑68 (MQ=0) aCAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt > 1:381283/1‑71 (MQ=255) cataaTCTTGAAAATTCAATATATTGCAGTTA‑TCTGGTAGGCCAGGTAAGCGTAGCGCATc > 1:2272322/4‑61 (MQ=2) cataaTCTTGAAAATTCAATATATTGCAGTTA‑TCTGGTAGGCCAGGTAAGCGTAGCGCATc > 1:3137942/4‑61 (MQ=2) cataaTCTTGAAAATTCAATATATTGCAGTTA‑TCTGGTAGGCCAGGTAAGCGTAGCGCATc > 1:4678403/4‑61 (MQ=2) cAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt > 1:623795/1‑71 (MQ=255) cAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat < 1:1641788/69‑4 (MQ=0) aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt > 1:1547663/1‑70 (MQ=255) aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt > 1:1090193/1‑70 (MQ=255) aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg > 1:3142566/1‑71 (MQ=255) aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg > 1:3810780/1‑71 (MQ=255) aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGt > 1:3145127/1‑35 (MQ=255) aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:2970338/1‑57 (MQ=255) aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc > 1:4833084/1‑71 (MQ=255) aTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGc > 1:724251/1‑38 (MQ=255) tCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt < 1:1816328/65‑1 (MQ=255) | CAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC > minE/2155523‑2155659 |
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: ‑ |