Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 222,988 | T→G | 32.9% | intergenic (+43/‑137) | lpcA → / → prfH | D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 222,988 | 0 | T | G | 32.9% | 72.1 / 48.5 | 52 | intergenic (+43/‑137) | lpcA/prfH | D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor |
Reads supporting (aligned to +/- strand): ref base T (21/13); new base G (6/10); total (27/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.36e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
GATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAG > minE/222921‑223057 | gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt < 1:1793638/70‑1 (MQ=255) attgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt < 1:786139/71‑1 (MQ=255) attgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt < 1:3412539/71‑1 (MQ=255) aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct < 1:1474638/71‑1 (MQ=255) aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct > 1:3409786/1‑70 (MQ=255) aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt < 1:2622324/60‑1 (MQ=255) aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt < 1:778494/60‑1 (MQ=255) gagaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt > 1:1312299/1‑60 (MQ=255) aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg < 1:1277879/65‑1 (MQ=255) tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcc > 1:836347/1‑69 (MQ=255) ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat > 1:38340/1‑65 (MQ=255) ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga < 1:2551613/69‑1 (MQ=255) ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg > 1:2776481/1‑64 (MQ=255) agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt > 1:3113530/1‑71 (MQ=255) gtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTg < 1:1454332/71‑1 (MQ=255) tCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt < 1:3585294/44‑1 (MQ=255) tccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAAcgc > 1:1188511/3‑71 (MQ=14) gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAtgcg > 1:67823/1‑71 (MQ=255) gTAGGCATGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg > 1:1449724/1‑56 (MQ=21) aGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt > 1:3040764/1‑61 (MQ=255) gCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATCa < 1:3162524/71‑1 (MQ=9) ccGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATc < 1:1109263/71‑1 (MQ=14) ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGc < 1:2257632/71‑1 (MQ=2) gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGcagatgcctgatgc < 1:3630928/49‑1 (MQ=11) aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgat < 1:3661025/46‑1 (MQ=255) aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc > 1:368178/1‑71 (MQ=255) aggcGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa > 1:3539000/4‑71 (MQ=11) aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt < 1:1566096/67‑1 (MQ=11) aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt < 1:1592322/67‑1 (MQ=11) aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt < 1:330049/67‑1 (MQ=11) aggcGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTa > 1:2336579/4‑71 (MQ=11) gACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa > 1:1374841/1‑70 (MQ=255) aCGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCTTGCCTACa < 1:1527663/71‑1 (MQ=255) cGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaa < 1:2412485/71‑1 (MQ=14) cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa > 1:1412261/1‑70 (MQ=255) gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa > 1:770164/1‑71 (MQ=255) gTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCag < 1:203916/62‑2 (MQ=14) tttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt < 1:2860960/71‑1 (MQ=255) ttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc > 1:2275068/1‑71 (MQ=255) tcaCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTa > 1:1025131/3‑56 (MQ=14) aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttcc > 1:2556266/1‑68 (MQ=11) aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttct < 1:3662595/71‑4 (MQ=11) gccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTa > 1:323003/1‑71 (MQ=255) ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc < 1:3931188/71‑1 (MQ=255) tgcATCCGGCA‑CTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa < 1:442082/57‑1 (MQ=14) gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc > 1:2146681/1‑71 (MQ=255) gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc > 1:1371132/1‑71 (MQ=255) aTCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGa < 1:483843/71‑1 (MQ=255) tCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAt > 1:747356/1‑60 (MQ=11) tCCGGCAC‑TGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTAt < 1:430686/46‑1 (MQ=2) ccGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt > 1:457429/1‑59 (MQ=255) gCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAtcc > 1:2705439/1‑46 (MQ=255) cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg > 1:2067074/1‑69 (MQ=255) cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg > 1:1868962/1‑69 (MQ=255) cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa < 1:2612615/71‑1 (MQ=255) aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGt < 1:381638/69‑1 (MQ=255) aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAg > 1:1276661/1‑71 (MQ=255) | GATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAG > minE/222921‑223057 |
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: ‑ |