| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,396,071 | G→T | 18.7% | G259G (GGG→GGT) | setB → | lactose/glucose efflux system |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,396,071 | 0 | G | T | 18.7% | 81.6 / 6.7 | 32 | G259G (GGG→GGT) | setB | lactose/glucose efflux system |
| Reads supporting (aligned to +/- strand): ref base G (21/5); new base T (2/4); total (23/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.85e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01 | |||||||||||
CATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTT > minE/1396008‑1396141 | cATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGTGAccgcc > 1:606880/1‑70 (MQ=255) cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt < 1:3109331/71‑2 (MQ=255) cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt < 1:3437316/71‑2 (MQ=255) cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt < 1:1336594/71‑1 (MQ=255) cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt < 1:2095777/71‑1 (MQ=255) tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgggg < 1:3233829/54‑1 (MQ=255) tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTgg > 1:1163031/1‑70 (MQ=255) tattatCCACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgccgc > 1:1807884/1‑60 (MQ=255) tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa < 1:2062447/70‑1 (MQ=255) tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa < 1:1376619/70‑1 (MQ=255) tCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAGATCCcg < 1:1542191/71‑1 (MQ=255) atgAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTAtcgctgccg > 1:2349334/3‑54 (MQ=255) aCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccg > 1:116538/1‑48 (MQ=255) cGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCTcgacga < 1:506966/71‑1 (MQ=255) gAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCTGGAAATTCCgatgat > 1:1785542/1‑71 (MQ=255) aCTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGt > 1:2811283/1‑71 (MQ=255) gCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGa > 1:458160/1‑71 (MQ=255) cATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTg > 1:2621546/1‑69 (MQ=255) tCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGA‑‑‑CCGCCGGGCTGGAAATcc > 1:3070076/1‑53 (MQ=255) cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAc > 1:2547902/1‑35 (MQ=255) cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGc > 1:2723217/1‑71 (MQ=255) agaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATc > 1:929507/1‑47 (MQ=255) gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGatat > 1:3039409/1‑70 (MQ=255) gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt > 1:2691213/1‑71 (MQ=255) gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt > 1:2215856/1‑71 (MQ=255) tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc > 1:2007992/1‑70 (MQ=255) tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc > 1:954128/1‑70 (MQ=255) cGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAAc > 1:1135815/1‑71 (MQ=255) tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa > 1:3357020/1‑55 (MQ=255) tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa > 1:557733/1‑55 (MQ=255) atgatgGGGACCGCCGCCGGGCTGGAAATCCcgac < 1:74870/35‑1 (MQ=255) atgatgGGGACCGCCGCCGGGCTGGAAATACCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTggg > 1:1859397/1‑71 (MQ=255) tgatgGGGGCCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTggg > 1:2012507/1‑70 (MQ=255) ggggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGc > 1:2506609/1‑71 (MQ=255) ggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGt < 1:215099/71‑1 (MQ=255) ggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGt < 1:674213/71‑1 (MQ=255) gACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGt < 1:3444370/58‑1 (MQ=255) gACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGtt < 1:3110403/71‑1 (MQ=255) | CATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTT > minE/1396008‑1396141 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |