Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,364,193 | G→A | 14.4% | P561P (CCG→CCA) | metG → | methionyl‑tRNA synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,364,193 | 0 | G | A | 14.4% | 69.6 / 2.3 | 41 | P561P (CCG→CCA) | metG | methionyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base G (25/10); new base A (2/4); total (27/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.57e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGAAGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCTGCGCG > minE/1364130‑1364257 | tGAAGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCt > 1:926444/1‑66 (MQ=255) tGAAGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCt > 1:1614080/1‑66 (MQ=255) aGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAg < 1:985142/68‑1 (MQ=255) aGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAg < 1:3009906/68‑1 (MQ=255) ggtggAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTTGCAGATGATcc > 1:3327713/1‑69 (MQ=255) ggtggAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGCTcc > 1:1344619/1‑69 (MQ=255) gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa > 1:2330480/1‑68 (MQ=255) gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa > 1:373435/1‑68 (MQ=255) gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa > 1:1363116/1‑68 (MQ=255) ccTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa > 1:1938370/1‑67 (MQ=255) ccTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAg > 1:1022159/1‑68 (MQ=255) tctAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGaa < 1:1430837/69‑1 (MQ=255) tctAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGaa < 1:3256921/69‑1 (MQ=255) tAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAgatg > 1:1837322/1‑53 (MQ=255) tAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAgatg > 1:2878524/1‑53 (MQ=255) tAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAgatg > 1:237984/1‑53 (MQ=255) caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa < 1:1680202/66‑1 (MQ=255) caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa < 1:3374652/66‑1 (MQ=255) caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa < 1:490429/66‑1 (MQ=255) caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa < 1:2959922/66‑1 (MQ=255) aaaGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAgatg > 1:2099852/1‑52 (MQ=255) agaagTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCAc < 1:1165001/69‑1 (MQ=255) gaagTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCAcc > 1:2405138/1‑69 (MQ=255) aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAg > 1:697978/1‑37 (MQ=255) aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac > 1:2329490/1‑69 (MQ=255) aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac < 1:2169708/69‑1 (MQ=255) aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac > 1:140623/1‑69 (MQ=255) tgTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa < 1:3291308/68‑1 (MQ=255) ggTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa < 1:2081416/69‑1 (MQ=255) ggTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa < 1:2077338/69‑1 (MQ=255) gTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa < 1:1059361/68‑1 (MQ=255) aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt > 1:3035919/1‑69 (MQ=255) aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt > 1:2714475/1‑69 (MQ=255) aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt > 1:2082150/1‑69 (MQ=255) aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt > 1:1093275/1‑69 (MQ=255) aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt > 1:1371364/1‑69 (MQ=255) ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc > 1:1229400/1‑69 (MQ=255) ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc > 1:2550720/1‑69 (MQ=255) ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc > 1:849407/1‑69 (MQ=255) gCCCACTGGCTGACGACCCGATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTg > 1:379189/1‑65 (MQ=255) gCCCACTGGCTGACGACCCGATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGcgcg > 1:3395374/1‑69 (MQ=255) | TGAAGCACTGGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCTGCGCG > minE/1364130‑1364257 |
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: ‑ |