| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,040,713 | C→T | 35.5% | D40D (GAC→GAT) | yqiH → | predicted periplasmic pilin chaperone |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,040,713 | 0 | C | T | 35.5% | 62.0 / 55.6 | 62 | D40D (GAC→GAT) | yqiH | predicted periplasmic pilin chaperone |
| Reads supporting (aligned to +/- strand): ref base C (19/21); new base T (14/8); total (33/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.90e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGA > minE/2040647‑2040777 | tATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATaa < 1:3331950/69‑1 (MQ=255) aTGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACa > 1:1475685/1‑67 (MQ=255) agcagcATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCg > 1:880186/1‑70 (MQ=255) ttGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGAc < 1:2616074/56‑1 (MQ=255) ttGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGAc < 1:1925165/56‑1 (MQ=255) ttGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCAtt > 1:77271/1‑70 (MQ=255) ttGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCAtt > 1:1785474/1‑70 (MQ=255) gCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa < 1:1210164/70‑1 (MQ=255) gCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa < 1:852780/70‑1 (MQ=255) ccTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCa < 1:11608/65‑1 (MQ=255) ccTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCa < 1:936889/65‑1 (MQ=255) tGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACaa > 1:964716/1‑53 (MQ=255) gggTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACt < 1:1210492/69‑1 (MQ=255) gCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGc < 1:2403360/51‑1 (MQ=255) ccATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAAc > 1:338026/1‑63 (MQ=255) atcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATc > 1:584852/1‑70 (MQ=255) atcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa < 1:3027920/71‑1 (MQ=255) atcatcCCGGATCGTACCCGTATCATAATCATTGAGTCAGACAAAGTCGAAGCATTAAAACTTACTAATCa < 1:3739917/71‑1 (MQ=255) atcatcCCGGATCGTACCCGTACCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTAGTAATCa < 1:2310369/71‑1 (MQ=255) tcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCg < 1:3644124/48‑1 (MQ=255) tcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa > 1:1061463/1‑70 (MQ=255) tcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa > 1:981772/1‑70 (MQ=255) tcatcCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa > 1:2451896/1‑70 (MQ=255) tcatcCCGGATCGTACCCGTATCATAATGAATGAGTAAGATAAACGCGAAGCATTAAAACTTACTAATCa > 1:192430/1‑70 (MQ=255) catcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa > 1:3583942/1‑69 (MQ=255) catcCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaaa > 1:2649355/1‑57 (MQ=255) atcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATaaa > 1:351446/1‑42 (MQ=255) atcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCg > 1:1684885/1‑46 (MQ=255) atcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAg > 1:427509/1‑71 (MQ=255) atcCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAg > 1:3770674/1‑71 (MQ=255) cccGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaaaa > 1:2631456/1‑55 (MQ=255) ccGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAg > 1:3792444/1‑46 (MQ=255) ccGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTa < 1:3034996/70‑1 (MQ=255) ggATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaa < 1:2824514/71‑1 (MQ=255) ggATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaa < 1:1653441/71‑1 (MQ=255) gATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaaa < 1:1898103/71‑1 (MQ=255) gATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaa < 1:3007946/70‑1 (MQ=255) aTCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaa < 1:3502298/69‑1 (MQ=255) tCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa < 1:3318924/47‑1 (MQ=255) tCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa < 1:3578785/47‑1 (MQ=255) tCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa < 1:3809984/47‑1 (MQ=255) gTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCt < 1:3316470/70‑1 (MQ=255) gTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaaaa > 1:546415/1‑47 (MQ=255) ccGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAg < 1:1233845/57‑1 (MQ=255) cGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTaa > 1:3565404/1‑40 (MQ=255) cGTATCATAATCAATGAGTCAGACAAAGGCGAAGCATTa > 1:3013809/1‑39 (MQ=255) gTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTAt > 1:3849047/1‑71 (MQ=255) tATCATAATGAATGAGTCAGACAAAGGCGAAGCAtt < 1:1574015/36‑1 (MQ=255) tATCATAATGAATGAGTCAGACAAAGGCGAAGCAtt < 1:2678665/36‑1 (MQ=255) tATCATAATGAATGAGTCAGACAAAGGCGAAGCAtt < 1:2766415/36‑1 (MQ=255) aTCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCt > 1:334887/1‑71 (MQ=255) tCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAAc > 1:2761965/1‑59 (MQ=255) aTAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAg < 1:2124153/50‑1 (MQ=255) atgaatgaGTCAGACAAAGGCGAAGCATTAAAACTTa > 1:3337552/1‑37 (MQ=255) atgaatgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCa < 1:3072353/44‑1 (MQ=255) tgaatgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGccc > 1:945614/1‑69 (MQ=255) tgaatgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCAGGccc > 1:2754158/1‑69 (MQ=255) gaatgaGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaa > 1:2671942/1‑50 (MQ=255) gaatgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCaaa > 1:2613447/1‑71 (MQ=255) atgaGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAGCCTGCCCTATCTGGc > 1:1478172/1‑64 (MQ=255) tgaGTCAGACAAACGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACAt < 1:2221155/71‑1 (MQ=255) tCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTa < 1:378080/53‑1 (MQ=255) aGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGa > 1:3336423/1‑68 (MQ=255) aGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATc > 1:3429538/1‑53 (MQ=255) | TATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGACAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGA > minE/2040647‑2040777 |
| 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: ‑ |