Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,872,825 | 0 | C | . | 66.7% | 70.4 / 108.1 | 60 | intergenic (‑5/+189) | ftsB/ygbE | cell division protein FtsB/conserved inner membrane protein YgbE |
Reads supporting (aligned to +/- strand): ref base C (14/6); new base . (16/24); total (30/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.39e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACA > NC_000913/2872776‑2872874 | tAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGc > 1:8620962/1‑50 (MQ=255) gCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGccc > 1:193258/1‑50 (MQ=255) gCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGccc > 1:6146901/1‑50 (MQ=255) ccagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccc > 1:3722773/1‑50 (MQ=255) cagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccct < 1:2693553/50‑2 (MQ=255) cagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccct < 1:221311/50‑2 (MQ=255) cagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccct < 1:7069264/50‑2 (MQ=255) cagaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGccccc > 1:4914423/1‑50 (MQ=255) agaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccctg < 1:3943402/50‑3 (MQ=255) agaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccctg < 1:709191/50‑3 (MQ=255) agaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCt > 1:8375909/1‑50 (MQ=255) agaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCt < 1:4157920/50‑1 (MQ=255) gaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccctg > 1:6060060/1‑47 (MQ=255) gaccagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGcccctg > 1:7536804/1‑47 (MQ=255) cagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGaaa < 1:6260669/50‑1 (MQ=255) cagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGaaaa < 1:2921999/50‑1 (MQ=255) cagaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGaaa > 1:3078028/1‑49 (MQ=39) agaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGaaaaa > 1:1030663/1‑50 (MQ=255) agaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGaaaaa > 1:713296/1‑50 (MQ=255) gaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGaaaaa < 1:6330563/50‑1 (MQ=255) gaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAAc < 1:1299230/50‑1 (MQ=255) gaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAAc < 1:2039442/50‑1 (MQ=255) gaATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAAc < 1:4079673/50‑1 (MQ=255) aTAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACgg > 1:4564194/1‑50 (MQ=255) aTAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACgg < 1:7379881/50‑1 (MQ=255) tAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACgg > 1:5028245/1‑49 (MQ=39) tAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGc < 1:3204993/50‑1 (MQ=255) tAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGc < 1:4210802/50‑1 (MQ=255) aGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACgg > 1:7705599/1‑49 (MQ=255) gCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGca > 1:7191242/1‑50 (MQ=255) gCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGca > 1:2574552/1‑50 (MQ=255) gCCAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGca < 1:1321747/49‑1 (MQ=39) ccAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGcat > 1:2769609/1‑49 (MQ=39) ccAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGcat > 1:5798048/1‑49 (MQ=39) cAGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGcatca < 1:2576974/50‑1 (MQ=255) aGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGcatca > 1:3330947/1‑50 (MQ=255) aGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGcatca > 1:5468293/1‑50 (MQ=255) aGCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGcatcat < 1:398142/50‑1 (MQ=255) gCAACAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGcatcatc < 1:2556806/50‑1 (MQ=255) aCAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCatcatcatc < 1:2932293/50‑1 (MQ=255) aCAGCAGCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCatcatcatc > 1:6084498/1‑50 (MQ=255) cagcagCGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATcc < 1:7435472/50‑1 (MQ=255) agcagcGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATcc > 1:4184508/1‑50 (MQ=255) gcagcGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCa > 1:8086548/1‑50 (MQ=255) cagcGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCAt < 1:5317196/50‑1 (MQ=255) cagcGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCAt > 1:7778068/1‑50 (MQ=255) agcGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATg < 1:4571952/50‑1 (MQ=255) gcGTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGc < 1:1417782/50‑1 (MQ=255) cGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATg > 1:5507353/1‑49 (MQ=255) gTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCAt < 1:8469900/50‑1 (MQ=255) gTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCAt > 1:3060482/1‑50 (MQ=255) gTTAGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCAt < 1:2234983/50‑1 (MQ=255) ttAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCAt < 1:3798131/50‑1 (MQ=255) tAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATc > 1:7593497/1‑50 (MQ=255) aGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATcc < 1:2032291/49‑1 (MQ=39) aGTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATcc < 1:2026336/49‑1 (MQ=39) gTTTACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGa < 1:2135764/50‑1 (MQ=255) tttACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCg > 1:7228637/1‑49 (MQ=255) ttACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCg > 1:7502697/1‑48 (MQ=255) tACCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAga > 1:8128405/1‑50 (MQ=255) aCCCATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAga < 1:4984627/49‑1 (MQ=39) cccATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAga < 1:4164294/49‑1 (MQ=255) cccATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAgac > 1:6240904/1‑49 (MQ=39) ccATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAgac > 1:6792414/1‑49 (MQ=255) ccATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAgacga > 1:7343959/1‑50 (MQ=255) cATCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAgacgac > 1:1212968/1‑50 (MQ=255) aTCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGAct < 1:1853794/50‑1 (MQ=255) aTCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGAct < 1:132223/50‑1 (MQ=255) tCCTG‑CCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGActc > 1:6043065/1‑50 (MQ=255) ccTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGActc > 1:1091847/1‑50 (MQ=255) cTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGActc < 1:3177882/49‑1 (MQ=255) gCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTa > 1:1356269/1‑49 (MQ=255) cccccTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACa > 1:3020352/1‑50 (MQ=255) | TAGCCAGACCAGAATAGCCAGCAACAGCAGCGTTAGTTTACCCATCCTGCCCCCTGAAAAACGGCATCATCATCCCATGCATCCGAAGACGACTCTACA > NC_000913/2872776‑2872874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |