Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 624,589 | 0 | T | C | 42.9% | 37.7 / 96.5 | 70 | intergenic (‑79/‑296) | fepB/entC | ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC |
Reads supporting (aligned to +/- strand): ref base T (31/9); new base C (16/14); total (47/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.20e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCG > NC_000913/624541‑624633 | gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt > 1:1409326/1‑50 (MQ=255) gggctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCctt > 1:933435‑M2/35‑50 (MQ=255) gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg < 1:1904402/50‑1 (MQ=255) gctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTgg > 1:355722‑M2/33‑50 (MQ=255) gctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTgg > 1:3148003‑M2/33‑50 (MQ=255) cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgg > 1:2238034/1‑50 (MQ=255) cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgg < 1:172929/50‑1 (MQ=255) gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg > 1:3545119/1‑50 (MQ=255) ctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg > 1:99981‑M2/32‑50 (MQ=255) ctgcaagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTggg < 1:2706045‑M2/19‑1 (MQ=255) tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgggg < 1:2400097/50‑1 (MQ=255) caagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga < 1:3004847‑M2/21‑1 (MQ=255) aTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg > 1:2534671/1‑47 (MQ=255) aagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga > 1:2955058‑M2/28‑48 (MQ=255) aagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga > 1:3210264‑M2/28‑48 (MQ=255) aagttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGga > 1:2969948‑M2/28‑48 (MQ=255) gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag > 1:429897/1‑50 (MQ=255) gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag > 1:623956/1‑50 (MQ=255) agttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGagag < 1:1944157‑M2/24‑1 (MQ=255) agttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGagag < 1:1686071‑M2/24‑1 (MQ=255) agttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGagag < 1:1382644‑M2/24‑1 (MQ=255) aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg > 1:290592/1‑50 (MQ=255) ttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAggg > 1:2481175‑M2/25‑50 (MQ=255) ttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAggg > 1:2746657‑M2/25‑50 (MQ=255) ttgcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAggg > 1:93735‑M2/25‑50 (MQ=255) cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAggg > 1:1621947/1‑50 (MQ=255) aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:172039/1‑50 (MQ=255) aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:478499/1‑50 (MQ=255) gcagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGt < 1:161675‑M2/27‑1 (MQ=255) cagaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGt > 1:2383063‑M2/22‑48 (MQ=255) aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:1550934/1‑48 (MQ=255) agaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg > 1:3044611‑M2/21‑50 (MQ=255) gaccgttacgagattgacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg > 1:3436646‑M2/20‑50 (MQ=255) gaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg > 1:996739‑M2/20‑50 (MQ=255) gaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg > 1:1988357‑M2/20‑50 (MQ=255) gaccgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg > 1:448073‑M2/20‑50 (MQ=255) cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg > 1:1686414/1‑50 (MQ=255) cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg < 1:1321720/49‑1 (MQ=255) cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg < 1:4204429/49‑1 (MQ=255) accgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg < 1:1729967‑M2/32‑1 (MQ=255) gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt > 1:207075/1‑50 (MQ=255) cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGt > 1:106873‑M2/17‑49 (MQ=255) cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTg < 1:1570815‑M2/34‑1 (MQ=255) cgttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTg < 1:201179‑M2/34‑1 (MQ=255) gttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGa < 1:3290864‑M2/35‑1 (MQ=255) gttacgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGa < 1:4021692‑M2/35‑1 (MQ=255) cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGCGGGATAGGGTgc > 1:1778069/1‑49 (MQ=18) cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt > 1:945468/1‑48 (MQ=255) cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt > 1:3712634/1‑48 (MQ=255) aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGat > 1:40674/1‑49 (MQ=255) aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAg < 1:4035954/50‑1 (MQ=255) aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg > 1:78157/1‑50 (MQ=255) aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg < 1:1261574/50‑1 (MQ=255) aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg < 1:3374940/50‑1 (MQ=255) acgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg < 1:3013144‑M2/37‑1 (MQ=255) aGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg > 1:2605628/1‑50 (MQ=255) aGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg > 1:2782963/1‑50 (MQ=255) cgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgggg < 1:2782246‑M2/39‑1 (MQ=255) cgagatggacgGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggg < 1:1649660‑M2/38‑1 (MQ=255) gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa > 1:3707241/1‑50 (MQ=255) acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa > 1:3989751/1‑48 (MQ=255) acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAca < 1:3340481/50‑1 (MQ=255) caATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcag > 1:2251372/1‑50 (MQ=255) aaTCCGTCCCCTCGCCCCTTTGGGGCGAGGGTTAGGGTGAGGGGAAcagc > 1:4042617/1‑50 (MQ=255) aaTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc > 1:752196/1‑50 (MQ=255) cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag > 1:128049/1‑49 (MQ=255) gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCa > 1:3302913/1‑50 (MQ=255) gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCa > 1:991631/1‑50 (MQ=255) tCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg > 1:294912/1‑50 (MQ=255) tCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg > 1:3771232/1‑50 (MQ=255) cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGc > 1:1520512/1‑50 (MQ=255) ccccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg > 1:359831/1‑49 (MQ=255) | GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCG > NC_000913/624541‑624633 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |