| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,467,882 | Δ1 bp | 20.8% | intergenic (+41/+87) | ubiG → / ← yfaL | bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,467,882 | 0 | C | . | 20.8% | 135.7 / 30.5 | 53 | intergenic (+41/+87) | ubiG/yfaL | bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin |
| Reads supporting (aligned to +/- strand): ref base C (20/22); new base . (7/4); total (27/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.01e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01 | |||||||||||
TGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATC > minE/1467819‑1467942 | tgctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCt > 1:3742706/1‑71 (MQ=255)tgctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCt > 1:1126061/1‑71 (MQ=255) gctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCt > 1:1741625/1‑70 (MQ=255) ctgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTa > 1:3617120/1‑71 (MQ=255) tgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTa < 1:3406434/70‑1 (MQ=255) tgcACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTAt > 1:1207719/1‑71 (MQ=255) ctgcacGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGcc > 1:2378207/4‑61 (MQ=255) cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTAt < 1:1608692/69‑1 (MQ=255) cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATc < 1:2488151/70‑1 (MQ=255) cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATc < 1:1482722/70‑1 (MQ=255) cacaCGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATc < 1:2724076/70‑1 (MQ=255) acacGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg > 1:3554321/1‑71 (MQ=255) cAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa > 1:2402202/1‑71 (MQ=255) gAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑a < 1:2607683/71‑1 (MQ=255) aTAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑aaa < 1:2899899/71‑1 (MQ=255) tAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAt < 1:3314348/71‑1 (MQ=255) tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCcgcg < 1:1597194/48‑1 (MQ=255) tGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTc < 1:1857861/71‑1 (MQ=255) gTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTAt < 1:1606085/50‑1 (MQ=255) tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTc > 1:3945823/1‑44 (MQ=255) tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCg > 1:1061835/1‑69 (MQ=255) gATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg > 1:1395669/1‑71 (MQ=255) aTGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg < 1:3460347/70‑1 (MQ=255) ttGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCg > 1:21774/1‑70 (MQ=255) gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa < 1:959386/71‑1 (MQ=255) gCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATa < 1:2080018/71‑1 (MQ=255) cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTa < 1:2018205/36‑1 (MQ=255) cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGc < 1:2404202/71‑1 (MQ=255) ggttGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGc > 1:3763323/4‑71 (MQ=2) cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCTGATAAGGCGtt < 1:105602/71‑1 (MQ=255) aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGttt < 1:3945874/71‑1 (MQ=255) aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGtt > 1:785002/1‑70 (MQ=255) aTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa > 1:552824/1‑71 (MQ=255) tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAC‑‑‑‑GCACTAC‑‑CCGTAGGTCGGATAAGGCGt > 1:155206/1‑70 (MQ=11) tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTa < 1:3159968/70‑1 (MQ=255) gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg > 1:1053630/1‑71 (MQ=255) ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTAcg < 1:58221/70‑1 (MQ=255) gATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt > 1:1678345/1‑55 (MQ=255) gATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGAt > 1:2093085/1‑55 (MQ=255) gATGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAc > 1:3139993/1‑71 (MQ=0) aTGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:3890194/1‑71 (MQ=255) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg < 1:3152114/71‑1 (MQ=0) aTGCGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcg < 1:2865534/71‑1 (MQ=0) cGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATc < 1:3065057/70‑1 (MQ=255) cGACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg < 1:1147494/71‑1 (MQ=0) gACGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTAcgccg > 1:110259/1‑70 (MQ=0) aCGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:1521629/1‑71 (MQ=0) aCGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:651183/71‑1 (MQ=0) cGCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa > 1:725712/1‑70 (MQ=0) gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc > 1:2562409/1‑71 (MQ=0) gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc > 1:3740559/1‑71 (MQ=0) gCTTG‑CGCGTCTTATCAGGCCTAC‑AAATTCCCG‑‑CACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATc > 1:2531117/1‑71 (MQ=0) cTTGCCGCGTCTTATCTGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCg < 1:3155041/36‑1 (MQ=255) gCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAg > 1:3487816/1‑36 (MQ=255) | TGCTGCACACGCAGAATAAGTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAC‑AAAT‑‑‑‑G‑‑CTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATC > minE/1467819‑1467942 |
| 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: ‑ |