| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,254,204 | 0 | G | T | 23.3% | 151.6 / 29.1 | 61 | A945D (GCT→GAT) | rpoC | RNA polymerase, beta prime subunit |
| Reads supporting (aligned to +/- strand): ref base G (26/20); new base T (3/11); total (30/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.20e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACAACCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCAC > minE/2254137‑2254267 | cacaACCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGCTGCTGAATCCAGCATcc < 1:2200912/70‑1 (MQ=38)cacaACCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCg > 1:999026/1‑71 (MQ=255)cacaACCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCg > 1:2747043/1‑71 (MQ=255) aaCCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCAc > 1:2219384/1‑71 (MQ=255) aaCCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCAc > 1:1540913/1‑71 (MQ=255) aaCCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCAc > 1:1208028/1‑71 (MQ=255) ccGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCAc > 1:2109304/1‑69 (MQ=255) aCTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGAt > 1:2300634/1‑71 (MQ=255) aCTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGAt > 1:2056226/1‑71 (MQ=255) ttCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGc > 1:2943021/1‑71 (MQ=255) cACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCg > 1:2850624/1‑71 (MQ=255) cGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGca > 1:2771707/1‑71 (MQ=255) ttGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGcac > 1:80095/1‑70 (MQ=255) ttGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGcac > 1:2238788/1‑70 (MQ=255) tGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCaccacc > 1:2508866/1‑70 (MQ=255) aGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCg < 1:2672828/69‑1 (MQ=255) aGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGAt < 1:1915963/71‑1 (MQ=255) gCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGAtg > 1:1436999/1‑71 (MQ=255) cTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGAtgt > 1:941644/1‑71 (MQ=255) gCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCAc < 1:2587085/44‑1 (MQ=255) tACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTa > 1:2275633/1‑71 (MQ=255) ccTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTa < 1:734014/69‑1 (MQ=255) ccTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACg < 1:2453784/71‑1 (MQ=255) ccTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACg < 1:1260025/71‑1 (MQ=255) ccTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACg < 1:2058822/71‑1 (MQ=255) cTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACg < 1:1594986/70‑1 (MQ=255) cTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGc < 1:3109899/71‑1 (MQ=255) tttGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCa > 1:1160520/1‑71 (MQ=255) tttCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGt > 1:1682246/1‑69 (MQ=255) tttCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGt > 1:2693242/1‑69 (MQ=255) tttCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGt > 1:1931192/1‑69 (MQ=255) ttCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAg > 1:785859/1‑71 (MQ=255) tCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAg < 1:2596198/70‑1 (MQ=255) aCTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGcaccac < 1:2340154/43‑1 (MQ=255) cTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCtgt < 1:225775/71‑1 (MQ=255) cTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCtgt < 1:3139459/71‑1 (MQ=255) cTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCtgt < 1:3031897/71‑1 (MQ=255) gATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGc > 1:781164/1‑55 (MQ=255) tGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTAc > 1:229230/1‑51 (MQ=255) tGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTAc > 1:2058769/1‑51 (MQ=255) tGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTAccc < 1:377223/71‑1 (MQ=255) tGCTGAATCCAGCATCCGCACTAGATGCCGCACCACCGATGTGGAAc < 1:1302868/47‑1 (MQ=25) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATTTGGAAc < 1:408230/47‑1 (MQ=25) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:237802/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:123219/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:179109/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:3045452/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:1774096/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:654594/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:3139495/47‑1 (MQ=37) tGCTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAAc < 1:379821/47‑1 (MQ=37) gCTGGATTCAGCACCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCg > 1:18846/1‑71 (MQ=255) cTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAACGt > 1:1059291/1‑47 (MQ=38) cTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAACGt > 1:759218/1‑47 (MQ=38) cTGAATCCAGCATCCGCACGAGATGCCGCACCACCGATGTGGAACGt > 1:2138428/1‑47 (MQ=38) tGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTAc < 1:187701/48‑1 (MQ=255) tGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTAc < 1:1264824/48‑1 (MQ=255) aTTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCa < 1:1700865/71‑1 (MQ=255) aTTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCa < 1:664892/71‑1 (MQ=255) ttCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCAc < 1:1032687/71‑1 (MQ=255) ttCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCAc > 1:230149/1‑71 (MQ=255) agcagcCGCACGAGATGCCGCACCACCGATGTGGAACGTa > 1:1560574/1‑40 (MQ=255) agcagcCGCACGAGATGCCGCACCACCGATGTGGAACGTa > 1:144474/1‑40 (MQ=255) | CACAACCGACTTCACGTTGCTGAGCTTGATGCTACCTTTGTTTTTCACTTGGATGCTGGATTCAGCAGCCGCACGAGATGCCGCACCACCGATGTGGAACGTACGCATGGTCAGCTGTGTACCCGGTTCAC > minE/2254137‑2254267 |
| 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: ‑ |